Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-04-21 02:53:35

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