macros/diskinematics_analysis/makePlot_reco_check.C

0001 int
0002 makePlot_reco_check()
0003 {
0004   gStyle->SetOptStat(1);
0005
0006   TFile *fin = new TFile("../../data/Sample_DISReco_ep.root", "OPEN");
0007   //TFile *fin = new TFile("../../data/DISReco_Pythia6_DIS_20x250_1k.root", "OPEN");
0008
0009   TTree *t_reco = (TTree*)fin->Get("event_cluster");
0010
0011   TCanvas *c1 = new TCanvas();
0012   t_reco->Draw( "em_cluster_e:em_cluster_eta", "em_evtgen_pid==11", "colz" );
0013   htemp->GetXaxis()->SetTitle("#eta_{cluster}");
0014   htemp->GetYaxis()->SetTitle("E_{cluster} (GeV)");
0015   c1->Print("plots/reco_check_c1.eps");
0016
0017   TCanvas *c2 = new TCanvas();
0018   t_reco->Draw( "em_evtgen_ptotal:em_evtgen_eta", "em_evtgen_pid==11", "colz" );
0019   htemp->GetXaxis()->SetTitle("#eta_{truth}");
0020   htemp->GetYaxis()->SetTitle("|p_{truth}| (GeV)");
0021   c2->Print("plots/reco_check_c2.eps");
0022
0023   TCanvas *c3 = new TCanvas();
0024   t_reco->Draw( "em_cluster_eta:em_evtgen_eta", "em_evtgen_pid==11", "colz" );
0025   htemp->GetXaxis()->SetTitle("#eta_{truth}");
0026   htemp->GetYaxis()->SetTitle("#eta_{cluster}");
0027   c3->Print("plots/reco_check_c3.eps");
0028
0029   TCanvas *c4 = new TCanvas();
0030   t_reco->Draw( "em_cluster_e:em_evtgen_ptotal", "em_evtgen_pid==11", "colz" );
0031   htemp->GetXaxis()->SetTitle("|p_{truth}| (GeV)");
0032   htemp->GetYaxis()->SetTitle("E_{cluster} (GeV)");
0033   c4->Print("plots/reco_check_c4.eps");
0034
0035   TCanvas *c5 = new TCanvas();
0036   t_reco->Draw( "em_cluster_phi:em_evtgen_phi", "em_evtgen_pid==11", "colz" );
0037   htemp->GetXaxis()->SetTitle("#phi_{truth}");
0038   htemp->GetYaxis()->SetTitle("#phi_{cluster}");
0039   c5->Print("plots/reco_check_c5.eps");
0040
0041   return 1;
0042 }