예제 #1
0
void makePlots() {
   // turn on/off batch mode
   gROOT->SetBatch(kTRUE);
   // set ROOT style
   style();
   //********************************************
   // root files
   //********************************************
   // *** MC ***
   string MC_noise = "OffsetAnalysis_histos_NoiseOnly_STARTUP3X_V8I.root";
   string MC_MB900 = "OffsetAnalysis_histos_MinBias900GeV_STARTUP3X_V8I.root";
   string MC_MB2360 = "OffsetAnalysis_histos_MinBias2360GeV_STARTUP3X_V8D.root";
   // *** data ***
   string data_ZB_900 = "OffsetAnalysis_histos_ZeroBias_R124022_900GeV.root";
   string data_Coll_900 = "OffsetAnalysis_histos_Collisions_R124022_900GeV.root";
   string data_Coll_2360 = "OffsetAnalysis_histos_Collisions_R124120_2360GeV.root";
   //********************************************
   // make plots
   //********************************************
   overlay_plots(MC_noise, data_ZB_900, 0.365, "p_AvgEinC5", "NoiseOnly MC, STARTUP3X_V8I","ZeroBias, Run 124022", "p_AvgEinC5_NoiseMC_vs_ZB_R124022_900GeV.png");
   overlay_plots(MC_noise, data_ZB_900, 0.365, "p_AvgpTinC5", "NoiseOnly MC, STARTUP3X_V8I","ZeroBias, Run 124022", "p_AvgpTinC5_NoiseMC_vs_ZB_R124022_900GeV.png");
   stack_in_cone_contribs(MC_noise, 0.365, "NoiseOnly MC, STARTUP3X_V8I", "Noise_MC.png");
   stack_in_cone_contribs(data_ZB_900, 0.365, "ZeroBias, Run 124022, 900 GeV", "ZB_R124022_900GeV.png");
   
   overlay_plots(MC_MB900, data_Coll_900, 2.2, "p_AvgEinC5", "MinBias900GeV MC, STARTUP3X_V8I","Collisions, Run 124022", "p_AvgEinC5_MinBias900GeVMC_vs_Coll_R124022_900GeV.png");
   overlay_plots(MC_MB900, data_Coll_900, 0.41, "p_AvgpTinC5", "MinBias900GeV MC, STARTUP3X_V8I","Collisions, Run 124022", "p_AvgpTinC5_MinBias900GeVMC_vs_Coll_R124022_900GeV.png");
   stack_in_cone_contribs(MC_MB900, 2.2, "MinBias900GeV MC, STARTUP3X_V8I", "MinBias900GeV_MC.png");
   stack_in_cone_contribs(data_Coll_900, 2.2, "Collisions, Run 124022, 900 GeV", "Coll_R124022_900GeV.png");
   
   overlay_plots(MC_MB2360, data_Coll_2360, 4.2, "p_AvgEinC5", "MinBias2360GeV MC, STARTUP3X_V8D","Collisions, Run 124120", "p_AvgEinC5_MinBias2360GeVMC_vs_Coll_R124120_2360GeV.png");
   overlay_plots(MC_MB2360, data_Coll_2360, 0.5, "p_AvgpTinC5", "MinBias2360GeV MC, STARTUP3X_V8D","Collisions, Run 124120", "p_AvgpTinC5_MinBias2360GeVMC_vs_Coll_R124120_2360GeV.png");
   stack_in_cone_contribs(MC_MB2360, 4.2, "MinBias2360GeV MC, STARTUP3X_V8D", "MinBias2360GeV_MC.png");
   stack_in_cone_contribs(data_Coll_2360, 4.2, "Collisions, Run 124120, 2360 GeV", "Coll_R124120_2360GeV.png");
}
예제 #2
0
void makePlots() {
   // turn on/off batch mode
   gROOT->SetBatch(kTRUE);
   // set ROOT style
   style();
   //********************************************
   // root files
   //********************************************
   // *** MC ***
   //string MC = "output_MC900GeV.root";
   string MC = "output_MC2360GeV.root";
   // *** data ***
   //string data = "output_DATA900GeV.root";
   string data = "output_DATA2360GeV.root";
   //********************************************
   // make plots
   //********************************************

   /*
   overlay_plots(MC, data, "g_calometPtMean_vs_ieta", -41, 41, -0.1, 0.7, "Calotower ieta", "CaloMET Mean [GeV]", "g_calometPtMean_vs_ieta_900.png");
   overlay_plots(MC, data, "g_calometPtRMS_vs_ieta", -41, 41, 0, 0.7, "Calotower ieta", "CaloMET RMS [GeV]", "g_calometPtRMS_vs_ieta_900.png");

   overlay_plots(MC, data, "g_caloSumetMean_vs_ieta", -41, 41, -0.1, 0.7,"Calotower ieta", "CaloSumET [GeV]", "g_caloSumetMean_vs_ieta_900.png");
   overlay_plots(MC, data, "g_caloSumetRMS_vs_ieta", -41, 41, 0, 0.7,"Calotower ieta", "CaloSumET RMS [GeV]", "g_caloSumetRMS_vs_ieta_900.png");

   overlay_plots(MC, data, "g_calometPxMean_vs_ieta", -41, 41, -0.1, 0.1, "Calotower ieta", "CaloMETPx Mean [GeV]", "g_calometPxMean_vs_ieta_900.png");
   overlay_plots(MC, data, "g_calometPxRMS_vs_ieta", -41, 41, 0, 0.7, "Calotower ieta", "CaloMETPx RMS [GeV]", "g_calometPxRMS_vs_ieta_900.png");

   overlay_plots(MC, data, "g_calometPyMean_vs_ieta", -41, 41, -0.1, 0.1, "Calotower ieta", "CaloMETPy Mean [GeV]", "g_calometPyMean_vs_ieta_900.png");
   overlay_plots(MC, data, "g_calometPyRMS_vs_ieta", -41, 41, 0., 0.7, "Calotower ieta", "CaloMETPy RMS [GeV]", "g_calometPyRMS_vs_ieta_900.png");
   */

   
   overlay_plots(MC, data, "g_calometPtMean_vs_ieta", -41, 41, -0.1, 0.7, "Calotower ieta", "CaloMET Mean [GeV]", "g_calometPtMean_vs_ieta_2360.png");
   overlay_plots(MC, data, "g_calometPtRMS_vs_ieta", -41, 41, 0, 0.7, "Calotower ieta", "CaloMET RMS [GeV]", "g_calometPtRMS_vs_ieta_2360.png");

   overlay_plots(MC, data, "g_caloSumetMean_vs_ieta", -41, 41, -0.1, 0.7,"Calotower ieta", "CaloSumET [GeV]", "g_caloSumetMean_vs_ieta_2360.png");
   overlay_plots(MC, data, "g_caloSumetRMS_vs_ieta", -41, 41, 0, 0.7,"Calotower ieta", "CaloSumET RMS [GeV]", "g_caloSumetRMS_vs_ieta_2360.png");

   overlay_plots(MC, data, "g_calometPxMean_vs_ieta", -41, 41, -0.1, 0.1, "Calotower ieta", "CaloMETPx Mean [GeV]", "g_calometPxMean_vs_ieta_2360.png");
   overlay_plots(MC, data, "g_calometPxRMS_vs_ieta", -41, 41, 0, 0.7, "Calotower ieta", "CaloMETPx RMS [GeV]", "g_calometPxRMS_vs_ieta_2360.png");

   overlay_plots(MC, data, "g_calometPyMean_vs_ieta", -41, 41, -0.1, 0.1, "Calotower ieta", "CaloMETPy Mean [GeV]", "g_calometPyMean_vs_ieta_2360.png");
   overlay_plots(MC, data, "g_calometPyRMS_vs_ieta", -41, 41, 0., 0.7, "Calotower ieta", "CaloMETPy RMS [GeV]", "g_calometPyRMS_vs_ieta_2360.png");
   

}
예제 #3
0
void makePlots() {
   // turn on/off batch mode
   gROOT->SetBatch(kTRUE);
   // set ROOT style
   style();
   //********************************************
   // root files
   //********************************************
   // *** MC ***
   string MC = "output_MC.root";
   // *** data ***
   string data = "output_DATA.root";
   //********************************************
   // make plots
   //********************************************

   overlay_plots(MC, data, "h_caloSumet", 0, 100, "CaloSumET [GeV]", "h_caloSumet.png",1);
   overlay_plots(MC, data, "h_calometPt", 0, 100, "CaloMET [GeV]", "h_calometPt.png",1);
   overlay_plots(MC, data, "h_calometPx", -50, 50, "CaloMETX [GeV]", "h_calometPx.png",1);
   overlay_plots(MC, data, "h_calometPy", -50, 50, "CaloMETY [GeV]", "h_calometPy.png",1);
   overlay_plots(MC, data, "h_calometPhi", -4, 4, "CaloMETPhi (rad.)", "h_calometPhi.png",1);
   overlay_plots(MC, data, "h_calometSig", 0, 10, "CaloMET/sqrt(SumET)", "h_calometSig.png",1);

   overlay_plots(MC, data, "h_caloSumetEB", 0, 100, "CaloSumET - EB [GeV]", "h_caloSumetEB.png",1);
   overlay_plots(MC, data, "h_caloSumetEE", 0, 100, "CaloSumET - EE [GeV]", "h_caloSumetEE.png",1);
   overlay_plots(MC, data, "h_caloSumetHB", 0, 100, "CaloSumET - HB [GeV]", "h_caloSumetHB.png",1);
   overlay_plots(MC, data, "h_caloSumetHE", 0, 100, "CaloSumET - HE [GeV]", "h_caloSumetHE.png",1);
   overlay_plots(MC, data, "h_caloSumetHadHF", 0, 100, "CaloSumET - had HF [GeV]", "h_caloSumetHadHF.png",1);
   overlay_plots(MC, data, "h_caloSumetEmHF", 0, 100, "CaloSumET - em HF [GeV]", "h_caloSumetEmHF.png",1);

}
예제 #4
0
void makePlots() {
   // turn on/off batch mode
   gROOT->SetBatch(kTRUE);
   // set ROOT style
   style();
   //********************************************
   // root files
   //********************************************
   // *** MC ***
   //string file0 = "$AXIGLUONDATA/axigluons_enujj/5fb-1_Summer11MC_AxigluonW_enujj_13092011/output_cutTable_axigluons_enujj/analysisClass_axigluons_enujj_plots.root";
   string file0 = "$AXIGLUONDATA/axigluons_enujj/5fb-1_Summer11MC_AxigluonW_enujj_Mjj500cut_13092011_2/output_cutTable_axigluons_enujj/analysisClass_axigluons_enujj_plots.root";
   // *** signal ***
   string file1 = file0;
   string file2 = file0;
   string file3 = file0;
   string file4 = file0;
//    string file1 = "$AXIGLUONDATA/axigluons_lnujj/AxigluonW_GenStudies_13092011/AG150WGen.root";
//    string file2 = "$AXIGLUONDATA/axigluons_lnujj/AxigluonW_GenStudies_13092011/AG500WGen.root";
//    string file3 = "$AXIGLUONDATA/axigluons_lnujj/AxigluonW_GenStudies_13092011/AG1000WGen.root";
//    string file4 = "$AXIGLUONDATA/axigluons_lnujj/AxigluonW_GenStudies_13092011/AG1500WGen.root";

   //********************************************
   // make plots
   //********************************************

   TFile file0_(file0.c_str());
   TFile file1_(file1.c_str());

   
   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Ele1_Pt" , 
		 "histo1D__AxigluonW_M150__Ele1_Pt", "histo1D__AxigluonW_M500__Ele1_Pt", "histo1D__AxigluonW_M1000__Ele1_Pt", "histo1D__AxigluonW_M1500__Ele1_Pt" ,
		 0, 600, "Ele1 Pt (GeV)", "Ele1_Pt.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Ele1_Eta" , 
		 "histo1D__AxigluonW_M150__Ele1_Eta", "histo1D__AxigluonW_M500__Ele1_Eta", "histo1D__AxigluonW_M1000__Ele1_Eta", "histo1D__AxigluonW_M1500__Ele1_Eta" ,
		 -3, 3, "Ele1 \\eta (GeV)", "Ele1_Eta.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet1_Pt" , 
		 "histo1D__AxigluonW_M150__Jet1_Pt", "histo1D__AxigluonW_M500__Jet1_Pt", "histo1D__AxigluonW_M1000__Jet1_Pt", "histo1D__AxigluonW_M1500__Jet1_Pt" ,
		 0, 1000, "Jet1 Pt (GeV)", "Jet1_Pt.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet2_Pt" , 
		 "histo1D__AxigluonW_M150__Jet2_Pt", "histo1D__AxigluonW_M500__Jet2_Pt", "histo1D__AxigluonW_M1000__Jet2_Pt", "histo1D__AxigluonW_M1500__Jet2_Pt" ,
		 0, 1000, "Jet2 Pt (GeV)", "Jet2_Pt.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet1_Eta" , 
		 "histo1D__AxigluonW_M150__Jet1_Eta", "histo1D__AxigluonW_M500__Jet1_Eta", "histo1D__AxigluonW_M1000__Jet1_Eta", "histo1D__AxigluonW_M1500__Jet1_Eta" ,
		 -3, 3, "Jet1 \\eta (GeV)", "Jet1_Eta.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet2_Eta" , 
		 "histo1D__AxigluonW_M150__Jet2_Eta", "histo1D__AxigluonW_M500__Jet2_Eta", "histo1D__AxigluonW_M1000__Jet2_Eta", "histo1D__AxigluonW_M1500__Jet2_Eta" ,
		 -3, 3, "Jet2 \\eta (GeV)", "Jet2_Eta.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__mDEta_Jet1Jet2" , 
		 "histo1D__AxigluonW_M150__mDEta_Jet1Jet2", "histo1D__AxigluonW_M500__mDEta_Jet1Jet2", "histo1D__AxigluonW_M1000__mDEta_Jet1Jet2", "histo1D__AxigluonW_M1500__mDEta_Jet1Jet2" ,
		 0, 5, "|\\Delta\\eta (j1,j2)|", "mDEta_Jet1Jet2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__mDPhi_Jet1Jet2" , 
		 "histo1D__AxigluonW_M150__mDPhi_Jet1Jet2", "histo1D__AxigluonW_M500__mDPhi_Jet1Jet2", "histo1D__AxigluonW_M1000__mDPhi_Jet1Jet2", "histo1D__AxigluonW_M1500__mDPhi_Jet1Jet2" ,
		 0, 3.1416, "|\\Delta\\phi (j1,j2)|", "mDPhi_Jet1Jet2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__mDphi_BosonJet1" , 
		 "histo1D__AxigluonW_M150__mDphi_BosonJet1", "histo1D__AxigluonW_M500__mDphi_BosonJet1", "histo1D__AxigluonW_M1000__mDphi_BosonJet1", "histo1D__AxigluonW_M1500__mDphi_BosonJet1" ,
		 0, 3.1416, "|\\Delta\\phi (boson,j1)|", "mDphi_BosonJet1.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__mDphi_BosonJet2" , 
		 "histo1D__AxigluonW_M150__mDphi_BosonJet2", "histo1D__AxigluonW_M500__mDphi_BosonJet2", "histo1D__AxigluonW_M1000__mDphi_BosonJet2", "histo1D__AxigluonW_M1500__mDphi_BosonJet2" ,
		 0, 3.1416, "|\\Delta\\phi (boson,j2)|", "mDphi_BosonJet2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet1_Pt_over_Mj1j2" , 
		 "histo1D__AxigluonW_M150__Jet1_Pt_over_Mj1j2", "histo1D__AxigluonW_M500__Jet1_Pt_over_Mj1j2", "histo1D__AxigluonW_M1000__Jet1_Pt_over_Mj1j2", "histo1D__AxigluonW_M1500__Jet1_Pt_over_Mj1j2" ,
		 0, 10, "Jet1 Pt / M_j1j2", "Jet1_Pt_over_Mj1j2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet2_Pt_over_Mj1j2" , 
		 "histo1D__AxigluonW_M150__Jet2_Pt_over_Mj1j2", "histo1D__AxigluonW_M500__Jet2_Pt_over_Mj1j2", "histo1D__AxigluonW_M1000__Jet2_Pt_over_Mj1j2", "histo1D__AxigluonW_M1500__Jet2_Pt_over_Mj1j2" ,
		 0, 10, "Jet2 Pt / M_j1j2", "Jet2_Pt_over_Mj1j2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__M_j1j2" , 
		 "histo1D__AxigluonW_M150__M_j1j2", "histo1D__AxigluonW_M500__M_j1j2", "histo1D__AxigluonW_M1000__M_j1j2", "histo1D__AxigluonW_M1500__M_j1j2" ,
		 0, 2000, "M_j1j2 (GeV)", "M_j1j2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__MT_Ele1MET" , 
		 "histo1D__AxigluonW_M150__MT_Ele1MET", "histo1D__AxigluonW_M500__MT_Ele1MET", "histo1D__AxigluonW_M1000__MT_Ele1MET", "histo1D__AxigluonW_M1500__MT_Ele1MET" ,
		 0, 500, "M_{T}(e,MET) (GeV)", "MT_Ele1MET.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__MET_Pt" , 
		 "histo1D__AxigluonW_M150__MET_Pt", "histo1D__AxigluonW_M500__MET_Pt", "histo1D__AxigluonW_M1000__MET_Pt", "histo1D__AxigluonW_M1500__MET_Pt" ,
		 0, 600, "PFMET (GeV)", "MET_Pt.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Pt_j1j2" , 
		 "histo1D__AxigluonW_M150__Pt_j1j2", "histo1D__AxigluonW_M500__Pt_j1j2", "histo1D__AxigluonW_M1000__Pt_j1j2", "histo1D__AxigluonW_M1500__Pt_j1j2" ,
		 0, 600, "Pt(j1,j2) (GeV)", "Pt_j1j2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet1_E_over_Mj1j2" , 
		 "histo1D__AxigluonW_M150__Jet1_E_over_Mj1j2", "histo1D__AxigluonW_M500__Jet1_E_over_Mj1j2", "histo1D__AxigluonW_M1000__Jet1_E_over_Mj1j2", "histo1D__AxigluonW_M1500__Jet1_E_over_Mj1j2" ,
		 0, 10, "Jet1 Energy / M_j1j2", "Jet1_E_over_Mj1j2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet1_E_over_Jet1_Pt" , 
		 "histo1D__AxigluonW_M150__Jet1_E_over_Jet1_Pt", "histo1D__AxigluonW_M500__Jet1_E_over_Jet1_Pt", "histo1D__AxigluonW_M1000__Jet1_E_over_Jet1_Pt", "histo1D__AxigluonW_M1500__Jet1_E_over_Jet1_Pt" ,
		 0, 10, "Jet1 Energy / Jet1 Pt", "Jet1_E_over_Jet1_Pt.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet2_E_over_Mj1j2" , 
		 "histo1D__AxigluonW_M150__Jet2_E_over_Mj1j2", "histo1D__AxigluonW_M500__Jet2_E_over_Mj1j2", "histo1D__AxigluonW_M1000__Jet2_E_over_Mj1j2", "histo1D__AxigluonW_M1500__Jet2_E_over_Mj1j2" ,
		 0, 10, "Jet2 Energy / M_j1j2", "Jet2_E_over_Mj1j2.png",0, 1, -99, -99,  "");

   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__Jet2_E_over_Jet2_Pt" , 
		 "histo1D__AxigluonW_M150__Jet2_E_over_Jet2_Pt", "histo1D__AxigluonW_M500__Jet2_E_over_Jet2_Pt", "histo1D__AxigluonW_M1000__Jet2_E_over_Jet2_Pt", "histo1D__AxigluonW_M1500__Jet2_E_over_Jet2_Pt" ,
		 0, 10, "Jet2 Energy / Jet2 Pt", "Jet2_E_over_Jet2_Pt.png",0, 1, -99, -99,  "");


  
   

   /* 
   overlay_plots(file0, file1, file2, file3, file4, "histo1D__ALLBKG__M_j1j2" , 
		 "h1_Mass_AG", "h1_Mass_AG", "h1_Mass_AG", "h1_Mass_AG" ,
		 0, 2000, "Axigluon Mass (GeV)", "M_j1j2.png",0, 1, -99, -99,  "");

   */



}