Geometry/macros/ECALMaterialBudgetCompare.C

0001 #include "TFile.h"
0002 #include "TTree.h"
0003 #include "TText.h"
0004
0005 void ECALMaterialBudgetCompare()
0006 {
0007
0008  gROOT ->Reset();
0009  const char*  rfilename = "old_matbdg_ECAL.root";
0010  const char*  sfilename = "new_matbdg_ECAL.root";
0011
0012  int rcolor = 2;
0013  int scolor = 4;
0014
0015  delete gROOT->GetListOfFiles()->FindObject(rfilename);
0016  delete gROOT->GetListOfFiles()->FindObject(sfilename);
0017
0018  TText* te = new TText();
0019  te->SetTextSize(0.1);
0020  TFile * rfile = new TFile(rfilename);
0021  TFile * sfile = new TFile(sfilename);
0022
0023  Char_t histo[200];
0024
0025  gStyle->SetOptStat("n");
0026
0027
0028 // Global Eta and Phi plots
0029
0030  if (1) {
0031    TCanvas * Ecal = new TCanvas("Ecal","Ecal",800,1000);
0032    Ecal->Divide(1,2);
0033
0034    TH1* eta_;
0035    rfile->GetObject("10;1",eta_);
0036    eta_;
0037    eta_->SetLineColor(rcolor);
0038
0039    TH1* neweta_;
0040    sfile->GetObject("10;1",neweta_);
0041    neweta_;
0042    neweta_->SetLineColor(scolor);
0043
0044    TH1* phi_;
0045    rfile->GetObject("20;1",phi_);
0046    phi_;
0047    phi_->SetLineColor(rcolor);
0048
0049    TH1* newphi_;
0050    sfile->GetObject("20;1",newphi_);
0051    newphi_;
0052    newphi_->SetLineColor(scolor);
0053
0054    Ecal->cd(1);
0055    eta_->Draw();
0056    neweta_->Draw("same");
0057    Ecal->cd(2);
0058    phi_->Draw();
0059    newphi_->Draw("same");
0060    Ecal->Print("Global_MB_compare.eps");
0061  }
0062
0063  // ECAL barrel plots
0064
0065
0066   if (1) {
0067     TCanvas * Ecal = new TCanvas("Ecal","Ecal",800,1000);
0068     Ecal->Divide(1,3);
0069
0070
0071     TH1* eta_;
0072     rfile->GetObject("1001;1",eta_);
0073     eta_;
0074     eta_->SetLineColor(rcolor);
0075
0076     TH1* neweta_;
0077     sfile->GetObject("1001;1",neweta_);
0078     neweta_;
0079     neweta_->SetLineColor(scolor);
0080
0081     TH1* phi_;
0082     rfile->GetObject("1002;1",phi_);
0083     phi_;
0084     phi_->SetLineColor(rcolor);
0085
0086     TH1* newphi_;
0087     sfile->GetObject("1002;1",newphi_);
0088     newphi_;
0089     newphi_->SetLineColor(scolor);
0090
0091     TH1* smphi_;
0092     rfile->GetObject("1003;1",smphi_);
0093     smphi_;
0094     smphi_->SetLineColor(rcolor);
0095
0096     TH1* newsmphi_;
0097     sfile->GetObject("1003;1",newsmphi_);
0098     newsmphi_;
0099     newsmphi_->SetLineColor(scolor);
0100
0101     Ecal->cd(1);
0102     eta_->Draw();
0103     neweta_->Draw("same");
0104     Ecal->cd(2);
0105     phi_->Draw();
0106     newphi_->Draw("same");
0107     Ecal->cd(3);
0108     smphi_->Draw();
0109     newsmphi_->Draw("same");
0110     Ecal->Print("ECAL_Barrel_MB_compare.eps");
0111
0112   }
0113
0114   if (1) {
0115     TCanvas * Ecal = new TCanvas("Ecal","Ecal",800,1000);
0116     Ecal->Divide(2,2);
0117
0118
0119     TH1* phimod1_;
0120     rfile->GetObject("1004;1",phimod1_);
0121     phimod1_;
0122     phimod1_->SetLineColor(rcolor);
0123
0124     TH1* newphimod1_;
0125     sfile->GetObject("1004;1",newphimod1_);
0126     newphimod1_;
0127     newphimod1_->SetLineColor(scolor);
0128
0129     TH1* phimod2_;
0130     rfile->GetObject("1005;1",phimod2_);
0131     phimod2_;
0132     phimod2_->SetLineColor(rcolor);
0133
0134     TH1* newphimod2_;
0135     sfile->GetObject("1005;1",newphimod2_);
0136     newphimod2_;
0137     newphimod2_->SetLineColor(scolor);
0138
0139     TH1* phimod3_;
0140     rfile->GetObject("1006;1",phimod3_);
0141     phimod3_;
0142     phimod3_->SetLineColor(rcolor);
0143
0144     TH1* newphimod3_;
0145     sfile->GetObject("1006;1",newphimod3_);
0146     newphimod3_;
0147     newphimod3_->SetLineColor(scolor);
0148
0149     TH1* phimod4_;
0150     rfile->GetObject("1007;1",phimod4_);
0151     phimod4_;
0152     phimod4_->SetLineColor(rcolor);
0153
0154     TH1* newphimod4_;
0155     sfile->GetObject("1007;1",newphimod4_);
0156     newphimod4_;
0157     newphimod4_->SetLineColor(scolor);
0158
0159
0160     Ecal->cd(1);
0161     phimod1_->Draw();
0162     newphimod1_->Draw("same");
0163     Ecal->cd(2);
0164     phimod2_->Draw();
0165     newphimod2_->Draw("same");
0166     Ecal->cd(3);
0167     phimod3_->Draw();
0168     newphimod3_->Draw("same");
0169     Ecal->cd(4);
0170     phimod4_->Draw();
0171     newphimod4_->Draw("same");
0172     Ecal->Print("ECAL_Barrel_modules_MB_compare.eps");
0173
0174   }
0175
0176  // ECAL preshower plots
0177
0178  if (1) {
0179    TCanvas * Ecal = new TCanvas("Ecal","Ecal",800,1000);
0180    Ecal->Divide(2,2);
0181
0182
0183    TH1* etazp_;
0184    rfile->GetObject("1011;1",etazp_);
0185    etazp_;
0186    etazp_->SetLineColor(rcolor);
0187
0188    TH1* newetazp_;
0189    sfile->GetObject("1011;1",newetazp_);
0190    newetazp_;
0191    newetazp_->SetLineColor(scolor);
0192
0193    TH1* phizp_;
0194    rfile->GetObject("1012;1",phizp_);
0195    phizp_;
0196    phizp_->SetLineColor(rcolor);
0197
0198    TH1* newphizp_;
0199    sfile->GetObject("1012;1",newphizp_);
0200    newphizp_;
0201    newphizp_->SetLineColor(scolor);
0202
0203    TH1* etazm_;
0204    rfile->GetObject("1013;1",etazm_);
0205    etazm_;
0206    etazm_->SetLineColor(rcolor);
0207
0208    TH1* newetazm_;
0209    sfile->GetObject("1013;1",newetazm_);
0210    newetazm_;
0211    newetazm_->SetLineColor(scolor);
0212
0213    TH1* phizm_;
0214    rfile->GetObject("1014;1",phizm_);
0215    phizm_;
0216    phizm_->SetLineColor(rcolor);
0217
0218    TH1* newphizm_;
0219    sfile->GetObject("1014;1",newphizm_);
0220    newphizm_;
0221    newphizm_->SetLineColor(scolor);
0222
0223
0224    Ecal->cd(1);
0225    etazp_->Draw();
0226    newetazp_->Draw("same");
0227    Ecal->cd(2);
0228    phizp_->Draw();
0229    newphizp_->Draw("same");
0230    Ecal->cd(3);
0231    etazm_->Draw();
0232    newetazm_->Draw("same");
0233    Ecal->cd(4);
0234    phizm_->Draw();
0235    newphizm_->Draw("same");
0236    Ecal->Print("ECAL_Preshower_MB_compare.eps");
0237
0238  }
0239
0240 }
0241