macros/diskinematics_analysis/compare_eictree_evtgenana.C

0001 int
0002 compare_eictree_evtgenana()
0003 {
0004   /* INPUT: EICTree created with eic-smear */
0005   TString fname_eictree("/gpfs/mnt/gpfs04/sphenix/user/nfeege/data/copy_eic_data_PYTHIA_ep/TREES/pythia.ep.10x250.1Mevents.RadCor=0.root");
0006
0007   /* INPUT: tree generated in Fun4All with DISKinematicsReco (truth) */
0008   TString fname_evtgenana("/gpfs/mnt/gpfs04/sphenix/user/nfeege/sphenixsw/devel/analysis/EICAnalysis/data/EventGenAna_pythiaeRHIC_DIS_10x250_1M.root");
0009
0010   /* get Fun4All tree from input file */
0011   TFile *file_evtgenana = new TFile (fname_evtgenana, "OPEN");
0012   TTree *tree_evtgenana = (TTree*)file_evtgenana->Get("event_truth");
0013
0014   /* add eic-smear tree as friend */
0015   tree_evtgenana->AddFriend("EICTree", fname_eictree);
0016
0017   /* plots correlations */
0018   TString drawmode("colz");
0019   //TString drawmode("");
0020
0021   TCanvas *cy = new TCanvas();
0022   tree_evtgenana->Draw("evtgen_y:EICTree.trueY","",drawmode);
0023
0024   TCanvas *cx = new TCanvas();
0025   tree_evtgenana->Draw("evtgen_x:EICTree.trueX","",drawmode);
0026
0027   TCanvas *cq2 = new TCanvas();
0028   tree_evtgenana->Draw("evtgen_Q2:EICTree.trueQ2","",drawmode);
0029
0030   TCanvas *cs = new TCanvas();
0031   tree_evtgenana->Draw("evtgen_s:EICTree.sHat","",drawmode);
0032
0033   TCanvas *cw = new TCanvas();
0034   tree_evtgenana->Draw("evtgen_W*evtgen_W:EICTree.trueW2","",drawmode);
0035
0036 }