macros/diskinematics_analysis/makePlot_reco_xQ2.C

0001 /* em_cluster_caloid values:
0002    1 = CEMC
0003    4 = EEMC
0004    5 = FEMC
0005 */
0006
0007 #include "plot_commons.h"
0008
0009 int
0010 makePlot_reco_xQ2()
0011 {
0012   gStyle->SetOptStat(0);
0013
0014   TFile *fin = new TFile("../../data/Sample_DISReco_ep.root", "OPEN");
0015   //TFile *fin = new TFile("../../data/DISReco_Pythia6_DIS_20x250_1k.root", "OPEN");
0016
0017   TTree *t_reco = (TTree*)fin->Get("event_cluster");
0018
0019   /* Q2: Create histogram */
0020   TH2F* hQ2_reco_truth = make_TH2D_xylog( 60, 0, 3, 60, 0, 3 );
0021   hQ2_reco_truth->SetName("hQ2_reco_truth");
0022   hQ2_reco_truth->GetXaxis()->SetTitle("Q^{2}_{truth} (GeV^{2})");
0023   hQ2_reco_truth->GetYaxis()->SetTitle("Q^{2}_{reco} (GeV^{2})");
0024
0025   /* Q2: Fill histogram */
0026   TCanvas *c1 = new TCanvas();
0027   c1->SetLogx(1);
0028   c1->SetLogy(1);
0029   c1->SetLogz(1);
0030   t_reco->Draw("em_reco_q2_e:evtgen_Q2 >> hQ2_reco_truth","em_evtgen_pid==11", "colz");
0031   c1->Print("plots/reco_xQ2_1.eps");
0032
0033   /* x: Create histogram */
0034   TH2F* hx_reco_truth = make_TH2D_xylog( 40, -4, 0, 40, -4, 0 );
0035   hx_reco_truth->SetName("hx_reco_truth");
0036   hx_reco_truth->GetXaxis()->SetTitle("x_{truth}");
0037   hx_reco_truth->GetYaxis()->SetTitle("x_{reco}");
0038
0039   /* x: Fill histogram */
0040   TCanvas *c2 = new TCanvas();
0041   c2->SetLogx(1);
0042   c2->SetLogy(1);
0043   c2->SetLogz(1);
0044   t_reco->Draw("em_reco_q2_e:evtgen_x >> hx_reco_truth","em_evtgen_pid==11", "colz");
0045   c2->Print("plots/reco_xQ2_2.eps");
0046
0047   /* W: Create histogram */
0048   TH2F* hW_reco_truth = new TH2F("hW_reco_truth","W_{truth};W_{reco}",50,0,100,50,0,100);
0049   hW_reco_truth->SetName("hW_reco_truth");
0050   hW_reco_truth->GetXaxis()->SetTitle("W_{truth}");
0051   hW_reco_truth->GetYaxis()->SetTitle("W_{reco}");
0052
0053   /* W: Fill histogram */
0054   TCanvas *c3 = new TCanvas();
0055   //c3->SetLogx(1);
0056   //c3->SetLogy(1);
0057   c3->SetLogz(1);
0058   t_reco->Draw("em_reco_w_e:evtgen_W >> hW_reco_truth","em_evtgen_pid==11", "colz");
0059   c3->Print("plots/reco_xQ2_3.eps");
0060
0061   /* y: Create histogram */
0062   TH2F* hy_reco_truth = new TH2F("hy_reco_truth","y_{truth};y_{reco}",100,0,1,100,0,1);
0063   hy_reco_truth->SetName("hy_reco_truth");
0064   hy_reco_truth->GetXaxis()->SetTitle("y_{truth}");
0065   hy_reco_truth->GetYaxis()->SetTitle("y_{reco}");
0066
0067   /* y: Fill histogram */
0068   TCanvas *c4 = new TCanvas();
0069   //c4->SetLogx(1);
0070   //c4->SetLogy(1);
0071   c3->SetLogz(1);
0072   t_reco->Draw("em_reco_y_e:evtgen_y >> hy_reco_truth","em_evtgen_pid==11", "colz");
0073   c4->Print("plots/reco_xQ2_4.eps");
0074
0075   return 1;
0076 }