{ TFile *file = new TFile("test.root"); gROOT->SetStyle("Plain"); gStyle->SetOptStat(0000); gStyle->SetPalette(1,0); gStyle->SetPaperSize(8,8); gStyle->SetPaperSize(20,26); /* //--------------------------------------------------------------------------------------------- // add relevant histograms //--------------------------------------------------------------------------------------------- TH2D *f1_acc_e = (TH2D*)file1->Get("fh_acc_e"); TH2D *f2_acc_e = (TH2D*)file2->Get("fh_acc_e"); TH2D *f3_acc_e = (TH2D*)file3->Get("fh_acc_e"); TH2D *f4_acc_e = (TH2D*)file4->Get("fh_acc_e"); TH2D *f5_acc_e = (TH2D*)file5->Get("fh_acc_e"); TH2D *f6_acc_e = (TH2D*)file6->Get("fh_acc_e"); TH2D *f7_acc_e = (TH2D*)file7->Get("fh_acc_e"); TH2D *f8_acc_e = (TH2D*)file8->Get("fh_acc_e"); TH2D *f9_acc_e = (TH2D*)file9->Get("fh_acc_e"); TH2D *f10_acc_e = (TH2D*)file10->Get("fh_acc_e"); TH2D *f11_acc_e = (TH2D*)file11->Get("fh_acc_e"); TH2D *f12_acc_e = (TH2D*)file12->Get("fh_acc_e"); fh_acc_e = new TH2D("fh_acc_e","geometrical acceptance e",50,0,5,40,0,2); fh_acc_e->Add(f1_acc_e,f2_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f3_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f4_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f5_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f6_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f7_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f8_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f9_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f10_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f11_acc_e,1,1); fh_acc_e->Add(fh_acc_e,f12_acc_e,1,1); fh_acc_e->Add(fh_acc_e,fh_acc_e,0.0833,0); TH2D *f1_acc_pi = (TH2D*)file1->Get("fh_acc_pi"); TH2D *f2_acc_pi = (TH2D*)file2->Get("fh_acc_pi"); TH2D *f3_acc_pi = (TH2D*)file3->Get("fh_acc_pi"); TH2D *f4_acc_pi = (TH2D*)file4->Get("fh_acc_pi"); TH2D *f5_acc_pi = (TH2D*)file5->Get("fh_acc_pi"); TH2D *f6_acc_pi = (TH2D*)file6->Get("fh_acc_pi"); TH2D *f7_acc_pi = (TH2D*)file7->Get("fh_acc_pi"); TH2D *f8_acc_pi = (TH2D*)file8->Get("fh_acc_pi"); TH2D *f9_acc_pi = (TH2D*)file9->Get("fh_acc_pi"); TH2D *f10_acc_pi = (TH2D*)file10->Get("fh_acc_pi"); TH2D *f11_acc_pi = (TH2D*)file11->Get("fh_acc_pi"); TH2D *f12_acc_pi = (TH2D*)file12->Get("fh_acc_pi"); fh_acc_pi = new TH2D("fh_acc_pi","geometrical acceptance pi",50,0,5,40,0,2); fh_acc_pi->Add(f1_acc_pi,f2_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f3_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f4_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f5_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f6_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f7_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f8_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f9_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f10_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f11_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,f12_acc_pi,1,1); fh_acc_pi->Add(fh_acc_pi,fh_acc_pi,0.0833,0); TH2D *f1_meas_e = (TH2D*)file1->Get("fh_meas_e"); TH2D *f2_meas_e = (TH2D*)file2->Get("fh_meas_e"); TH2D *f3_meas_e = (TH2D*)file3->Get("fh_meas_e"); TH2D *f4_meas_e = (TH2D*)file4->Get("fh_meas_e"); TH2D *f5_meas_e = (TH2D*)file5->Get("fh_meas_e"); TH2D *f6_meas_e = (TH2D*)file6->Get("fh_meas_e"); TH2D *f7_meas_e = (TH2D*)file7->Get("fh_meas_e"); TH2D *f8_meas_e = (TH2D*)file8->Get("fh_meas_e"); TH2D *f9_meas_e = (TH2D*)file9->Get("fh_meas_e"); TH2D *f10_meas_e = (TH2D*)file10->Get("fh_meas_e"); TH2D *f11_meas_e = (TH2D*)file11->Get("fh_meas_e"); TH2D *f12_meas_e = (TH2D*)file12->Get("fh_meas_e"); fh_meas_e = new TH2D("fh_meas_e","measured e",50,0,5,40,0,2); fh_meas_e->Add(f1_meas_e,f2_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f3_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f4_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f5_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f6_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f7_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f8_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f9_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f10_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f11_meas_e,1,1); fh_meas_e->Add(fh_meas_e,f12_meas_e,1,1); TH2D *f1_meas_pi = (TH2D*)file1->Get("fh_meas_pi"); TH2D *f2_meas_pi = (TH2D*)file2->Get("fh_meas_pi"); TH2D *f3_meas_pi = (TH2D*)file3->Get("fh_meas_pi"); TH2D *f4_meas_pi = (TH2D*)file4->Get("fh_meas_pi"); TH2D *f5_meas_pi = (TH2D*)file5->Get("fh_meas_pi"); TH2D *f6_meas_pi = (TH2D*)file6->Get("fh_meas_pi"); TH2D *f7_meas_pi = (TH2D*)file7->Get("fh_meas_pi"); TH2D *f8_meas_pi = (TH2D*)file8->Get("fh_meas_pi"); TH2D *f9_meas_pi = (TH2D*)file9->Get("fh_meas_pi"); TH2D *f10_meas_pi = (TH2D*)file10->Get("fh_meas_pi"); TH2D *f11_meas_pi = (TH2D*)file11->Get("fh_meas_pi"); TH2D *f12_meas_pi = (TH2D*)file12->Get("fh_meas_pi"); fh_meas_pi = new TH2D("fh_meas_pi","measured pi",50,0,5,40,0,2); fh_meas_pi->Add(f1_meas_pi,f2_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f3_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f4_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f5_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f6_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f7_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f8_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f9_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f10_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f11_meas_pi,1,1); fh_meas_pi->Add(fh_meas_pi,f12_meas_pi,1,1); TH2D *f1_ring_nsts_e = (TH2D*)file1->Get("fh_ring_nsts_e"); TH2D *f2_ring_nsts_e = (TH2D*)file2->Get("fh_ring_nsts_e"); TH2D *f3_ring_nsts_e = (TH2D*)file3->Get("fh_ring_nsts_e"); TH2D *f4_ring_nsts_e = (TH2D*)file4->Get("fh_ring_nsts_e"); TH2D *f5_ring_nsts_e = (TH2D*)file5->Get("fh_ring_nsts_e"); TH2D *f6_ring_nsts_e = (TH2D*)file6->Get("fh_ring_nsts_e"); TH2D *f7_ring_nsts_e = (TH2D*)file7->Get("fh_ring_nsts_e"); TH2D *f8_ring_nsts_e = (TH2D*)file8->Get("fh_ring_nsts_e"); TH2D *f9_ring_nsts_e = (TH2D*)file9->Get("fh_ring_nsts_e"); TH2D *f10_ring_nsts_e = (TH2D*)file10->Get("fh_ring_nsts_e"); TH2D *f11_ring_nsts_e = (TH2D*)file11->Get("fh_ring_nsts_e"); TH2D *f12_ring_nsts_e = (TH2D*)file12->Get("fh_ring_nsts_e"); fh_ring_nsts_e = new TH2D("fh_ring_nsts_e","found prim ring e",50,0,5,40,0,2); fh_ring_nsts_e->Add(f1_ring_nsts_e,f2_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f3_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f4_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f5_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f6_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f7_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f8_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f9_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f10_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f11_ring_nsts_e,1,1); fh_ring_nsts_e->Add(fh_ring_nsts_e,f12_ring_nsts_e,1,1); TH2D *f1_ring_nsts_pi = (TH2D*)file1->Get("fh_ring_nsts_pi"); TH2D *f2_ring_nsts_pi = (TH2D*)file2->Get("fh_ring_nsts_pi"); TH2D *f3_ring_nsts_pi = (TH2D*)file3->Get("fh_ring_nsts_pi"); TH2D *f4_ring_nsts_pi = (TH2D*)file4->Get("fh_ring_nsts_pi"); TH2D *f5_ring_nsts_pi = (TH2D*)file5->Get("fh_ring_nsts_pi"); TH2D *f6_ring_nsts_pi = (TH2D*)file6->Get("fh_ring_nsts_pi"); TH2D *f7_ring_nsts_pi = (TH2D*)file7->Get("fh_ring_nsts_pi"); TH2D *f8_ring_nsts_pi = (TH2D*)file8->Get("fh_ring_nsts_pi"); TH2D *f9_ring_nsts_pi = (TH2D*)file9->Get("fh_ring_nsts_pi"); TH2D *f10_ring_nsts_pi = (TH2D*)file10->Get("fh_ring_nsts_pi"); TH2D *f11_ring_nsts_pi = (TH2D*)file11->Get("fh_ring_nsts_pi"); TH2D *f12_ring_nsts_pi = (TH2D*)file12->Get("fh_ring_nsts_pi"); fh_ring_nsts_pi = new TH2D("fh_ring_nsts_pi","found prim ring pi",50,0,5,40,0,2); fh_ring_nsts_pi->Add(f1_ring_nsts_pi,f2_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f3_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f4_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f5_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f6_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f7_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f8_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f9_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f10_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f11_ring_nsts_pi,1,1); fh_ring_nsts_pi->Add(fh_ring_nsts_pi,f12_ring_nsts_pi,1,1); TH2D *f1_ringmatchtrue_e = (TH2D*)file1->Get("fh_ringmatchtrue_e"); TH2D *f2_ringmatchtrue_e = (TH2D*)file2->Get("fh_ringmatchtrue_e"); TH2D *f3_ringmatchtrue_e = (TH2D*)file3->Get("fh_ringmatchtrue_e"); TH2D *f4_ringmatchtrue_e = (TH2D*)file4->Get("fh_ringmatchtrue_e"); TH2D *f5_ringmatchtrue_e = (TH2D*)file5->Get("fh_ringmatchtrue_e"); TH2D *f6_ringmatchtrue_e = (TH2D*)file6->Get("fh_ringmatchtrue_e"); TH2D *f7_ringmatchtrue_e = (TH2D*)file7->Get("fh_ringmatchtrue_e"); TH2D *f8_ringmatchtrue_e = (TH2D*)file8->Get("fh_ringmatchtrue_e"); TH2D *f9_ringmatchtrue_e = (TH2D*)file9->Get("fh_ringmatchtrue_e"); TH2D *f10_ringmatchtrue_e = (TH2D*)file10->Get("fh_ringmatchtrue_e"); TH2D *f11_ringmatchtrue_e = (TH2D*)file11->Get("fh_ringmatchtrue_e"); TH2D *f12_ringmatchtrue_e = (TH2D*)file12->Get("fh_ringmatchtrue_e"); fh_ringmatchtrue_e = new TH2D("fh_ringmatchtrue_e","found prim ring e",50,0,5,40,0,2); fh_ringmatchtrue_e->Add(f1_ringmatchtrue_e,f2_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f3_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f4_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f5_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f6_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f7_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f8_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f9_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f10_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f11_ringmatchtrue_e,1,1); fh_ringmatchtrue_e->Add(fh_ringmatchtrue_e,f12_ringmatchtrue_e,1,1); TH2D *f1_ringmatchtrue_pi = (TH2D*)file1->Get("fh_ringmatchtrue_pi"); TH2D *f2_ringmatchtrue_pi = (TH2D*)file2->Get("fh_ringmatchtrue_pi"); TH2D *f3_ringmatchtrue_pi = (TH2D*)file3->Get("fh_ringmatchtrue_pi"); TH2D *f4_ringmatchtrue_pi = (TH2D*)file4->Get("fh_ringmatchtrue_pi"); TH2D *f5_ringmatchtrue_pi = (TH2D*)file5->Get("fh_ringmatchtrue_pi"); TH2D *f6_ringmatchtrue_pi = (TH2D*)file6->Get("fh_ringmatchtrue_pi"); TH2D *f7_ringmatchtrue_pi = (TH2D*)file7->Get("fh_ringmatchtrue_pi"); TH2D *f8_ringmatchtrue_pi = (TH2D*)file8->Get("fh_ringmatchtrue_pi"); TH2D *f9_ringmatchtrue_pi = (TH2D*)file9->Get("fh_ringmatchtrue_pi"); TH2D *f10_ringmatchtrue_pi = (TH2D*)file10->Get("fh_ringmatchtrue_pi"); TH2D *f11_ringmatchtrue_pi = (TH2D*)file11->Get("fh_ringmatchtrue_pi"); TH2D *f12_ringmatchtrue_pi = (TH2D*)file12->Get("fh_ringmatchtrue_pi"); fh_ringmatchtrue_pi = new TH2D("fh_ringmatchtrue_pi","found prim ring pi",50,0,5,40,0,2); fh_ringmatchtrue_pi->Add(f1_ringmatchtrue_pi,f2_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f3_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f4_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f5_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f6_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f7_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f8_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f9_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f10_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f11_ringmatchtrue_pi,1,1); fh_ringmatchtrue_pi->Add(fh_ringmatchtrue_pi,f12_ringmatchtrue_pi,1,1); TH2D *f1_ringmatchfalse_track_e = (TH2D*)file1->Get("fh_ringmatchfalse_track_e"); TH2D *f2_ringmatchfalse_track_e = (TH2D*)file2->Get("fh_ringmatchfalse_track_e"); TH2D *f3_ringmatchfalse_track_e = (TH2D*)file3->Get("fh_ringmatchfalse_track_e"); TH2D *f4_ringmatchfalse_track_e = (TH2D*)file4->Get("fh_ringmatchfalse_track_e"); TH2D *f5_ringmatchfalse_track_e = (TH2D*)file5->Get("fh_ringmatchfalse_track_e"); TH2D *f6_ringmatchfalse_track_e = (TH2D*)file6->Get("fh_ringmatchfalse_track_e"); TH2D *f7_ringmatchfalse_track_e = (TH2D*)file7->Get("fh_ringmatchfalse_track_e"); TH2D *f8_ringmatchfalse_track_e = (TH2D*)file8->Get("fh_ringmatchfalse_track_e"); TH2D *f9_ringmatchfalse_track_e = (TH2D*)file9->Get("fh_ringmatchfalse_track_e"); TH2D *f10_ringmatchfalse_track_e = (TH2D*)file10->Get("fh_ringmatchfalse_track_e"); TH2D *f11_ringmatchfalse_track_e = (TH2D*)file11->Get("fh_ringmatchfalse_track_e"); TH2D *f12_ringmatchfalse_track_e = (TH2D*)file12->Get("fh_ringmatchfalse_track_e"); fh_ringmatchfalse_track_e = new TH2D("fh_ringmatchfalse_track_e","found prim ring e",50,0,5,40,0,2); fh_ringmatchfalse_track_e->Add(f1_ringmatchfalse_track_e,f2_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f3_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f4_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f5_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f6_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f7_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f8_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f9_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f10_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f11_ringmatchfalse_track_e,1,1); fh_ringmatchfalse_track_e->Add(fh_ringmatchfalse_track_e,f12_ringmatchfalse_track_e,1,1); TH2D *f1_ringmatchfalse_track_pi = (TH2D*)file1->Get("fh_ringmatchfalse_track_pi"); TH2D *f2_ringmatchfalse_track_pi = (TH2D*)file2->Get("fh_ringmatchfalse_track_pi"); TH2D *f3_ringmatchfalse_track_pi = (TH2D*)file3->Get("fh_ringmatchfalse_track_pi"); TH2D *f4_ringmatchfalse_track_pi = (TH2D*)file4->Get("fh_ringmatchfalse_track_pi"); TH2D *f5_ringmatchfalse_track_pi = (TH2D*)file5->Get("fh_ringmatchfalse_track_pi"); TH2D *f6_ringmatchfalse_track_pi = (TH2D*)file6->Get("fh_ringmatchfalse_track_pi"); TH2D *f7_ringmatchfalse_track_pi = (TH2D*)file7->Get("fh_ringmatchfalse_track_pi"); TH2D *f8_ringmatchfalse_track_pi = (TH2D*)file8->Get("fh_ringmatchfalse_track_pi"); TH2D *f9_ringmatchfalse_track_pi = (TH2D*)file9->Get("fh_ringmatchfalse_track_pi"); TH2D *f10_ringmatchfalse_track_pi = (TH2D*)file10->Get("fh_ringmatchfalse_track_pi"); TH2D *f11_ringmatchfalse_track_pi = (TH2D*)file11->Get("fh_ringmatchfalse_track_pi"); TH2D *f12_ringmatchfalse_track_pi = (TH2D*)file12->Get("fh_ringmatchfalse_track_pi"); fh_ringmatchfalse_track_pi = new TH2D("fh_ringmatchfalse_track_pi","found prim ring pi",50,0,5,40,0,2); fh_ringmatchfalse_track_pi->Add(f1_ringmatchfalse_track_pi,f2_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f3_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f4_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f5_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f6_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f7_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f8_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f9_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f10_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f11_ringmatchfalse_track_pi,1,1); fh_ringmatchfalse_track_pi->Add(fh_ringmatchfalse_track_pi,f12_ringmatchfalse_track_pi,1,1); TH2D *f1_ringmatchtrue_mv_e = (TH2D*)file1->Get("fh_ringmatchtrue_mv_e"); TH2D *f2_ringmatchtrue_mv_e = (TH2D*)file2->Get("fh_ringmatchtrue_mv_e"); TH2D *f3_ringmatchtrue_mv_e = (TH2D*)file3->Get("fh_ringmatchtrue_mv_e"); TH2D *f4_ringmatchtrue_mv_e = (TH2D*)file4->Get("fh_ringmatchtrue_mv_e"); TH2D *f5_ringmatchtrue_mv_e = (TH2D*)file5->Get("fh_ringmatchtrue_mv_e"); TH2D *f6_ringmatchtrue_mv_e = (TH2D*)file6->Get("fh_ringmatchtrue_mv_e"); TH2D *f7_ringmatchtrue_mv_e = (TH2D*)file7->Get("fh_ringmatchtrue_mv_e"); TH2D *f8_ringmatchtrue_mv_e = (TH2D*)file8->Get("fh_ringmatchtrue_mv_e"); TH2D *f9_ringmatchtrue_mv_e = (TH2D*)file9->Get("fh_ringmatchtrue_mv_e"); TH2D *f10_ringmatchtrue_mv_e = (TH2D*)file10->Get("fh_ringmatchtrue_mv_e"); TH2D *f11_ringmatchtrue_mv_e = (TH2D*)file11->Get("fh_ringmatchtrue_mv_e"); TH2D *f12_ringmatchtrue_mv_e = (TH2D*)file12->Get("fh_ringmatchtrue_mv_e"); fh_ringmatchtrue_mv_e = new TH2D("fh_ringmatchtrue_mv_e","found prim ring e",50,0,5,40,0,2); fh_ringmatchtrue_mv_e->Add(f1_ringmatchtrue_mv_e,f2_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f3_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f4_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f5_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f6_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f7_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f8_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f9_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f10_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f11_ringmatchtrue_mv_e,1,1); fh_ringmatchtrue_mv_e->Add(fh_ringmatchtrue_mv_e,f12_ringmatchtrue_mv_e,1,1); TH2D *f1_ringmatchtrue_mv_pi = (TH2D*)file1->Get("fh_ringmatchtrue_mv_pi"); TH2D *f2_ringmatchtrue_mv_pi = (TH2D*)file2->Get("fh_ringmatchtrue_mv_pi"); TH2D *f3_ringmatchtrue_mv_pi = (TH2D*)file3->Get("fh_ringmatchtrue_mv_pi"); TH2D *f4_ringmatchtrue_mv_pi = (TH2D*)file4->Get("fh_ringmatchtrue_mv_pi"); TH2D *f5_ringmatchtrue_mv_pi = (TH2D*)file5->Get("fh_ringmatchtrue_mv_pi"); TH2D *f6_ringmatchtrue_mv_pi = (TH2D*)file6->Get("fh_ringmatchtrue_mv_pi"); TH2D *f7_ringmatchtrue_mv_pi = (TH2D*)file7->Get("fh_ringmatchtrue_mv_pi"); TH2D *f8_ringmatchtrue_mv_pi = (TH2D*)file8->Get("fh_ringmatchtrue_mv_pi"); TH2D *f9_ringmatchtrue_mv_pi = (TH2D*)file9->Get("fh_ringmatchtrue_mv_pi"); TH2D *f10_ringmatchtrue_mv_pi = (TH2D*)file10->Get("fh_ringmatchtrue_mv_pi"); TH2D *f11_ringmatchtrue_mv_pi = (TH2D*)file11->Get("fh_ringmatchtrue_mv_pi"); TH2D *f12_ringmatchtrue_mv_pi = (TH2D*)file12->Get("fh_ringmatchtrue_mv_pi"); fh_ringmatchtrue_mv_pi = new TH2D("fh_ringmatchtrue_mv_pi","found prim ring pi",50,0,5,40,0,2); fh_ringmatchtrue_mv_pi->Add(f1_ringmatchtrue_mv_pi,f2_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f3_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f4_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f5_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f6_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f7_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f8_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f9_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f10_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f11_ringmatchtrue_mv_pi,1,1); fh_ringmatchtrue_mv_pi->Add(fh_ringmatchtrue_mv_pi,f12_ringmatchtrue_mv_pi,1,1); TH2D *f1_ringmatchfalse_track_mv_e = (TH2D*)file1->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f2_ringmatchfalse_track_mv_e = (TH2D*)file2->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f3_ringmatchfalse_track_mv_e = (TH2D*)file3->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f4_ringmatchfalse_track_mv_e = (TH2D*)file4->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f5_ringmatchfalse_track_mv_e = (TH2D*)file5->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f6_ringmatchfalse_track_mv_e = (TH2D*)file6->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f7_ringmatchfalse_track_mv_e = (TH2D*)file7->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f8_ringmatchfalse_track_mv_e = (TH2D*)file8->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f9_ringmatchfalse_track_mv_e = (TH2D*)file9->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f10_ringmatchfalse_track_mv_e = (TH2D*)file10->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f11_ringmatchfalse_track_mv_e = (TH2D*)file11->Get("fh_ringmatchfalse_track_mv_e"); TH2D *f12_ringmatchfalse_track_mv_e = (TH2D*)file12->Get("fh_ringmatchfalse_track_mv_e"); fh_ringmatchfalse_track_mv_e = new TH2D("fh_ringmatchfalse_track_mv_e","found prim ring e",50,0,5,40,0,2); fh_ringmatchfalse_track_mv_e->Add(f1_ringmatchfalse_track_mv_e,f2_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f3_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f4_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f5_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f6_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f7_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f8_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f9_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f10_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f11_ringmatchfalse_track_mv_e,1,1); fh_ringmatchfalse_track_mv_e->Add(fh_ringmatchfalse_track_mv_e,f12_ringmatchfalse_track_mv_e,1,1); TH2D *f1_ringmatchfalse_track_mv_fake_e = (TH2D*)file1->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f2_ringmatchfalse_track_mv_fake_e = (TH2D*)file2->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f3_ringmatchfalse_track_mv_fake_e = (TH2D*)file3->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f4_ringmatchfalse_track_mv_fake_e = (TH2D*)file4->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f5_ringmatchfalse_track_mv_fake_e = (TH2D*)file5->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f6_ringmatchfalse_track_mv_fake_e = (TH2D*)file6->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f7_ringmatchfalse_track_mv_fake_e = (TH2D*)file7->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f8_ringmatchfalse_track_mv_fake_e = (TH2D*)file8->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f9_ringmatchfalse_track_mv_fake_e = (TH2D*)file9->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f10_ringmatchfalse_track_mv_fake_e = (TH2D*)file10->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f11_ringmatchfalse_track_mv_fake_e = (TH2D*)file11->Get("fh_ringmatchfalse_track_mv_fake_e"); TH2D *f12_ringmatchfalse_track_mv_fake_e = (TH2D*)file12->Get("fh_ringmatchfalse_track_mv_fake_e"); fh_ringmatchfalse_track_mv_fake_e = new TH2D("fh_ringmatchfalse_track_mv_fake_e","found prim ring e, fakes",50,0,5,40,0,2); fh_ringmatchfalse_track_mv_fake_e->Add(f1_ringmatchfalse_track_mv_fake_e,f2_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f3_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f4_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f5_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f6_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f7_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f8_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f9_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f10_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f11_ringmatchfalse_track_mv_fake_e,1,1); fh_ringmatchfalse_track_mv_fake_e->Add(fh_ringmatchfalse_track_mv_fake_e,f12_ringmatchfalse_track_mv_fake_e,1,1); TH2D *f1_ringmatchfalse_track_mv_pi = (TH2D*)file1->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f2_ringmatchfalse_track_mv_pi = (TH2D*)file2->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f3_ringmatchfalse_track_mv_pi = (TH2D*)file3->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f4_ringmatchfalse_track_mv_pi = (TH2D*)file4->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f5_ringmatchfalse_track_mv_pi = (TH2D*)file5->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f6_ringmatchfalse_track_mv_pi = (TH2D*)file6->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f7_ringmatchfalse_track_mv_pi = (TH2D*)file7->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f8_ringmatchfalse_track_mv_pi = (TH2D*)file8->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f9_ringmatchfalse_track_mv_pi = (TH2D*)file9->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f10_ringmatchfalse_track_mv_pi = (TH2D*)file10->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f11_ringmatchfalse_track_mv_pi = (TH2D*)file11->Get("fh_ringmatchfalse_track_mv_pi"); TH2D *f12_ringmatchfalse_track_mv_pi = (TH2D*)file12->Get("fh_ringmatchfalse_track_mv_pi"); fh_ringmatchfalse_track_mv_pi = new TH2D("fh_ringmatchfalse_track_mv_pi","found prim ring pi",50,0,5,40,0,2); fh_ringmatchfalse_track_mv_pi->Add(f1_ringmatchfalse_track_mv_pi,f2_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f3_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f4_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f5_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f6_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f7_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f8_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f9_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f10_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f11_ringmatchfalse_track_mv_pi,1,1); fh_ringmatchfalse_track_mv_pi->Add(fh_ringmatchfalse_track_mv_pi,f12_ringmatchfalse_track_mv_pi,1,1); TH2D *f1_ring_mv_e = (TH2D*)file1->Get("fh_ring_mv_e"); TH2D *f2_ring_mv_e = (TH2D*)file2->Get("fh_ring_mv_e"); TH2D *f3_ring_mv_e = (TH2D*)file3->Get("fh_ring_mv_e"); TH2D *f4_ring_mv_e = (TH2D*)file4->Get("fh_ring_mv_e"); TH2D *f5_ring_mv_e = (TH2D*)file5->Get("fh_ring_mv_e"); TH2D *f6_ring_mv_e = (TH2D*)file6->Get("fh_ring_mv_e"); TH2D *f7_ring_mv_e = (TH2D*)file7->Get("fh_ring_mv_e"); TH2D *f8_ring_mv_e = (TH2D*)file8->Get("fh_ring_mv_e"); TH2D *f9_ring_mv_e = (TH2D*)file9->Get("fh_ring_mv_e"); TH2D *f10_ring_mv_e = (TH2D*)file10->Get("fh_ring_mv_e"); TH2D *f11_ring_mv_e = (TH2D*)file11->Get("fh_ring_mv_e"); TH2D *f12_ring_mv_e = (TH2D*)file12->Get("fh_ring_mv_e"); fh_ring_mv_e = new TH2D("fh_ring_mv_e","found prim ring e",50,0,5,40,0,2); fh_ring_mv_e->Add(f1_ring_mv_e,f2_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f3_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f4_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f5_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f6_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f7_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f8_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f9_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f10_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f11_ring_mv_e,1,1); fh_ring_mv_e->Add(fh_ring_mv_e,f12_ring_mv_e,1,1); TH2D *f1_ring_mv_pi = (TH2D*)file1->Get("fh_ring_mv_pi"); TH2D *f2_ring_mv_pi = (TH2D*)file2->Get("fh_ring_mv_pi"); TH2D *f3_ring_mv_pi = (TH2D*)file3->Get("fh_ring_mv_pi"); TH2D *f4_ring_mv_pi = (TH2D*)file4->Get("fh_ring_mv_pi"); TH2D *f5_ring_mv_pi = (TH2D*)file5->Get("fh_ring_mv_pi"); TH2D *f6_ring_mv_pi = (TH2D*)file6->Get("fh_ring_mv_pi"); TH2D *f7_ring_mv_pi = (TH2D*)file7->Get("fh_ring_mv_pi"); TH2D *f8_ring_mv_pi = (TH2D*)file8->Get("fh_ring_mv_pi"); TH2D *f9_ring_mv_pi = (TH2D*)file9->Get("fh_ring_mv_pi"); TH2D *f10_ring_mv_pi = (TH2D*)file10->Get("fh_ring_mv_pi"); TH2D *f11_ring_mv_pi = (TH2D*)file11->Get("fh_ring_mv_pi"); TH2D *f12_ring_mv_pi = (TH2D*)file12->Get("fh_ring_mv_pi"); fh_ring_mv_pi = new TH2D("fh_ring_mv_pi","found prim ring pi",50,0,5,40,0,2); fh_ring_mv_pi->Add(f1_ring_mv_pi,f2_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f3_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f4_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f5_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f6_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f7_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f8_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f9_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f10_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f11_ring_mv_pi,1,1); fh_ring_mv_pi->Add(fh_ring_mv_pi,f12_ring_mv_pi,1,1); */ //------------------------------------------------------------------------------------------------ cout << " electrons: true ring-track, wrong " << fh_ringmatchtrue_mv_e->GetEntries()<< " " << fh_ringmatchfalse_track_mv_e->GetEntries() << " " << fh_ringmatchfalse_track_mv_fake_e->GetEntries() << endl; cout << " pions: true ring-track, wrong " << fh_ringmatchtrue_mv_pi->GetEntries()<< " " << fh_ringmatchfalse_track_mv_pi->GetEntries() << endl; //-------------------------------------------------------------------------------------------------- TCanvas *c1 = new TCanvas("c1","electrons",1200,800); c1->Divide(3,2); TF1 f1_e("f1_e","sqrt(1*1-0.000511*0.000511*sinh(x)*sinh(x))/cosh(x)",0.,5.); TF1 f05_e("f05_e","sqrt(0.5*0.5-0.000511*0.000511*sinh(x)*sinh(x))/cosh(x)",0.,5.); TF1 f2_e("f2_e","sqrt(2*2-0.000511*0.000511*sinh(x)*sinh(x))/cosh(x)",0.,5.); TF1 f12_e("f12_e","sqrt(12*12-0.000511*0.000511*sinh(x)*sinh(x))/cosh(x)",0.,5.); c1->cd(1); fh_acc_e->SetMaximum(1); fh_acc_e->GetXaxis()->SetTitle("Ylab"); fh_acc_e->GetYaxis()->SetTitle("Pt[GeV/c]"); fh_acc_e->Draw("colz"); // f1_e->Draw("same"); f05_e->Draw("same"); // f2_e->Draw("same"); f12_e->Draw("same"); c1->cd(2); ring_eff_nsts_e = new TH2D("ring_nsts_eff_e","ring_nsts_eff_e",50,0,5,40,0,2); ring_eff_nsts_e->GetXaxis()->SetTitle("Ylab"); ring_eff_nsts_e->GetYaxis()->SetTitle("Pt[GeV/c]"); ring_eff_nsts_e->Divide(fh_ring_nsts_e,fh_meas_e,1,1); ring_eff_nsts_e->SetMaximum(1); ring_eff_nsts_e->Draw("colz"); // f1_e->Draw("same"); f05_e->Draw("same"); // f2_e->Draw("same"); f12_e->Draw("same"); c1->cd(3); ring_match_nsts_e = new TH2D("ring_match_nsts_e","ring_match_nsts_e",50,0,5,40,0,2); ring_match_nsts_e->GetXaxis()->SetTitle("Ylab"); ring_match_nsts_e->GetYaxis()->SetTitle("Pt[GeV/c]"); ring_match_nsts_e->Divide(fh_ringmatchtrue_e,fh_ring_nsts_e,1,1); ring_match_nsts_e->SetMaximum(1); ring_match_nsts_e->Draw("colz"); // f1_e->Draw("same"); f05_e->Draw("same"); // f2_e->Draw("same"); f12_e->Draw("same"); // c1->cd(4); // ring_match_false_e = new TH2D("ring_match_false_e","ring_match_false_e",50,0,5,40,0,2); // ring_match_false_e->Divide(fh_ringmatchfalse_track_e,fh_ring_nsts_e,1,1); // ring_match_false_e->SetMaximum(1); // ring_match_false_e->Draw("colz"); c1->cd(4); ring_match_mv_e = new TH2D("ring_match_mv_e","ring_match_mv_e",50,0,5,40,0,2); ring_match_mv_e->GetXaxis()->SetTitle("Ylab"); ring_match_mv_e->GetYaxis()->SetTitle("Pt[GeV/c]"); // ring_match_mv_e->Divide(fh_ringmatchtrue_mv_e,fh_ring_nsts_e,1,1); ring_match_mv_e->Divide(fh_ringmatchtrue_mv_e,fh_ring_mv_e,1,1); ring_match_mv_e->SetMaximum(1); ring_match_mv_e->Draw("colz"); c1->cd(5); ring_match_mv_false_e = new TH2D("ring_match_mv_false_e","ring_match_mv_false_e",50,0,5,40,0,2); ring_match_mv_false_e->GetXaxis()->SetTitle("Ylab"); ring_match_mv_false_e->GetYaxis()->SetTitle("Pt[GeV/c]"); // ring_match_mv_false_e->Divide(fh_ringmatchfalse_track_mv_e,fh_ring_nsts_e,1,1); ring_match_mv_false_e->Divide(fh_ringmatchfalse_track_mv_e,fh_ring_mv_e,1,1); ring_match_mv_false_e->SetMaximum(1); ring_match_mv_false_e->Draw("colz"); c1->cd(6); ring_match_mv_false_fake_e = new TH2D("ring_match_mv_false_fake_e","ring_match_mv_false_fake_e",50,0,5,40,0,2); ring_match_mv_false_fake_e->GetXaxis()->SetTitle("Ylab"); ring_match_mv_false_fake_e->GetYaxis()->SetTitle("Pt[GeV/c]"); // ring_match_mv_false_e->Divide(fh_ringmatchfalse_track_mv_e,fh_ring_nsts_e,1,1); ring_match_mv_false_fake_e->Divide(fh_ringmatchfalse_track_mv_fake_e,fh_ring_mv_e,1,1); ring_match_mv_false_fake_e->SetMaximum(1); ring_match_mv_false_fake_e->Draw("colz"); /* c1->cd(4); acc_all_e = new TH2D("acc_all_e","acc_all_e",50,0,5,40,0,2); acc_all_e->Divide(fh_spot_e,fh_gen_e,1,1); acc_all_e->SetMaximum(1); acc_all_e->Draw("colz"); c1->cd(5); ring_eff_e = new TH2D("ring_eff_e","ring_eff_e",50,0,5,40,0,2); ring_eff_e->Divide(fh_ring_e,fh_spot_e,1,1); ring_eff_e->SetMaximum(1); ring_eff_e->Draw("colz"); c1->cd(6); spot_match_e = new TH2D("spot_match_e","spot_match_e",50,0,5,40,0,2); spot_match_e->Divide(fh_spotmatchtrue_e,fh_spot_e,1,1); spot_match_e->SetMaximum(1); spot_match_e->Draw("colz"); */ //------------------------------------------------------------------------------------------------ TCanvas *c2 = new TCanvas("c2","pions",1200,800); c2->Divide(3,2); c2->cd(1); fh_acc_pi->SetMaximum(1); fh_acc_pi->GetXaxis()->SetTitle("Ylab"); fh_acc_pi->GetYaxis()->SetTitle("Pt[GeV/c]"); fh_acc_pi->Draw("colz"); c2->cd(2); ring_eff_nsts_pi = new TH2D("ring_nsts_eff_pi","ring_nsts_eff_pi",50,0,5,40,0,2); ring_eff_nsts_pi->GetXaxis()->SetTitle("Ylab"); ring_eff_nsts_pi->GetYaxis()->SetTitle("Pt[GeV/c]"); ring_eff_nsts_pi->Divide(fh_ring_nsts_pi,fh_meas_pi,1,1); ring_eff_nsts_pi->SetMaximum(1); ring_eff_nsts_pi->Draw("colz"); c2->cd(3); ring_match_nsts_pi = new TH2D("ring_match_nsts_pi","ring_match_nsts_pi",50,0,5,40,0,2); ring_match_nsts_pi->GetXaxis()->SetTitle("Ylab"); ring_match_nsts_pi->GetYaxis()->SetTitle("Pt[GeV/c]"); ring_match_nsts_pi->Divide(fh_ringmatchtrue_pi,fh_ring_nsts_pi,1,1); ring_match_nsts_pi->SetMaximum(1); ring_match_nsts_pi->Draw("colz"); c2->cd(4); ring_match_false_pi = new TH2D("ring_match_false_pi","ring_match_false_pi",50,0,5,40,0,2); ring_match_false_pi->GetXaxis()->SetTitle("Ylab"); ring_match_false_pi->GetYaxis()->SetTitle("Pt[GeV/c]"); ring_match_false_pi->Divide(fh_ringmatchfalse_track_pi,fh_ring_nsts_pi,1,1); ring_match_false_pi->SetMaximum(1); ring_match_false_pi->Draw("colz"); c2->cd(5); ring_match_mv_pi = new TH2D("ring_match_mv_pi","ring_match_mv_pi",50,0,5,40,0,2); ring_match_mv_pi->GetXaxis()->SetTitle("Ylab"); ring_match_mv_pi->GetYaxis()->SetTitle("Pt[GeV/c]"); // ring_match_mv_pi->Divide(fh_ringmatchtrue_mv_pi,fh_ring_nsts_pi,1,1); ring_match_mv_pi->Divide(fh_ringmatchtrue_mv_pi,fh_ring_mv_pi,1,1); ring_match_mv_pi->SetMaximum(1); ring_match_mv_pi->Draw("colz"); c2->cd(6); ring_match_mv_false_pi = new TH2D("ring_match_mv_false_pi","ring_match_mv_false_pi",50,0,5,40,0,2); ring_match_mv_false_pi->GetXaxis()->SetTitle("Ylab"); ring_match_mv_false_pi->GetYaxis()->SetTitle("Pt[GeV/c]"); // ring_match_mv_false_pi->Divide(fh_ringmatchfalse_track_mv_pi,fh_ring_nsts_pi,1,1); ring_match_mv_false_pi->Divide(fh_ringmatchfalse_track_mv_pi,fh_ring_mv_pi,1,1); ring_match_mv_false_pi->SetMaximum(1); ring_match_mv_false_pi->Draw("colz"); /* c2->cd(4); acc_all_pi = new TH2D("acc_all_pi","acc_all_pi",50,0,5,40,0,2); acc_all_pi->Divide(fh_spot_pi,fh_gen_pi,1,1); acc_all_pi->SetMaximum(1); acc_all_pi->Draw("colz"); c2->cd(5); ring_eff_pi = new TH2D("ring_eff_pi","ring_eff_pi",50,0,5,40,0,2); ring_eff_pi->Divide(fh_ring_pi,fh_spot_pi,1,1); ring_eff_pi->SetMaximum(1); ring_eff_pi->Draw("colz"); c2->cd(6); spot_match_pi = new TH2D("spot_match_pi","spot_match_pi",50,0,5,40,0,2); spot_match_pi->Divide(fh_spotmatchtrue_pi,fh_spot_pi,1,1); spot_match_pi->SetMaximum(1); spot_match_pi->Draw("colz"); */ //--------------------------------------------------------------------------------------------------------- TCanvas *c3 = new TCanvas("c3","not reconstructed rings/ fakes",800,800); c3->Divide(2,2); // Detector l1a = new TLine(-158,75.0,158,75.0); l1a->SetLineWidth(3); l2a = new TLine(-158,250.0,158,250.0); l2a->SetLineWidth(3); l3a = new TLine(-158,75.0,-158,250.0); l3a->SetLineWidth(3); l4a = new TLine(158,75.0,158,250.0); l4a->SetLineWidth(3); l1b = new TLine(-158,-75.0,158,-75.0); l1b->SetLineWidth(3); l2b = new TLine(-158,-250.0,158,-250.0); l2b->SetLineWidth(3); l3b = new TLine(-158,-75.0,-158,-250.0); l3b->SetLineWidth(3); l4b = new TLine(158,-75.0,158,-250.0); l4b->SetLineWidth(3); c3->cd(1); fh_fake_ring->Draw("colz"); l1a->Draw(); l2a->Draw(); l3a->Draw(); l4a->Draw(); l1b->Draw(); l2b->Draw(); l3b->Draw(); l4b->Draw(); c3->cd(2); fh_spot_prim_notrec_e->Draw("colz"); l1a->Draw(); l2a->Draw(); l3a->Draw(); l4a->Draw(); l1b->Draw(); l2b->Draw(); l3b->Draw(); l4b->Draw(); c3->cd(3); fh_spot_prim_notrec_pi->Draw("colz"); l1a->Draw(); l2a->Draw(); l3a->Draw(); l4a->Draw(); l1b->Draw(); l2b->Draw(); l3b->Draw(); l4b->Draw(); //--------------------------------------------------------------------------------------------------------- TCanvas *c4 = new TCanvas("c4","distances ring - track: electrons",800,800); c4->Divide(2,2); cutline = new TLine(0.8,0.,0.8,1000.); fh_disttrue_e->Project3D("xz"); fh_disttrue_e->Project3D("yz"); fh_disttrue_e->Project3D("z"); fh_distfalse_e->Project3D("z"); fh_distfalse_fake_e->Project3D("z"); fh_disttrue_mv_e->Project3D("z"); fh_distfalse_mv_e->Project3D("z"); fh_distfalse_mv_fake_e->Project3D("z"); c4->cd(1); fh_disttrue_e_xz->Draw("colz"); c4->cd(2); fh_disttrue_e_yz->Draw("colz"); c4->cd(3); fh_disttrue_e_z->Draw(); fh_distfalse_e_z->Draw("same"); fh_distfalse_fake_e_z->Draw("same"); // cutline->Draw(); c4->cd(4); fh_disttrue_mv_e_z->Draw(); Double_t true_e = fh_disttrue_mv_e_z->GetEntries(); fh_distfalse_mv_e_z->Draw("same"); Double_t wrong_e = fh_distfalse_mv_e_z->GetEntries(); fh_distfalse_mv_fake_e_z->Draw("same"); Double_t wrong_fake_e = fh_distfalse_mv_fake_e_z->GetEntries(); cout << true_e << " " << wrong_e << " " << wrong_fake_e << endl; //------------------------------------------------------- TCanvas *c4a = new TCanvas("c4a","distances ring - track: pions",800,800); c4a->Divide(2,2); fh_disttrue_pi->Project3D("xz"); fh_disttrue_pi->Project3D("yz"); fh_disttrue_pi->Project3D("z"); fh_distfalse_pi->Project3D("z"); fh_disttrue_mv_pi->Project3D("z"); fh_distfalse_mv_pi->Project3D("z"); c4a->cd(1); fh_disttrue_pi_xz->Draw("colz"); c4a->cd(2); fh_disttrue_pi_yz->Draw("colz"); c4a->cd(3); fh_disttrue_pi_z->Draw(); fh_distfalse_pi_z->Draw("same"); // cutline->Draw(); c4a->cd(4); fh_disttrue_mv_pi_z->Draw(); fh_distfalse_mv_pi_z->Draw("same"); //--------------------------------------------------------------------------------------------------------- TCanvas *c5 = new TCanvas("c5","radius",1200,800); c5->Divide(3,2); fh_radius->Project3D("xz"); fh_radius->Project3D("yz"); fh_radius->Project3D("z"); c5->cd(1); fh_radius_xz->GetXaxis()->SetTitle("Radius [sm]"); fh_radius_xz->GetYaxis()->SetTitle("P[GeV/c]"); fh_radius_xz->Draw("colz"); c5->cd(2); fh_radius_yz->GetXaxis()->SetTitle("Radius [sm]"); fh_radius_yz->GetYaxis()->SetTitle("P[GeV/c]"); fh_radius_yz->Draw("colz"); c5->cd(3); fh_radius_z->GetXaxis()->SetTitle("Radius [sm]"); fh_radius_z->GetYaxis()->SetTitle("P[GeV/c]"); fh_radius_z->Draw(); c5->cd(4); fh_fake_p_radius->GetXaxis()->SetTitle("Radius [sm]"); fh_fake_p_radius->GetYaxis()->SetTitle("P[GeV/c]"); fh_fake_p_radius->Draw("colz"); c5->cd(5); fh_r_points_true->GetXaxis()->SetTitle("Number of Hits"); fh_r_points_true->GetYaxis()->SetTitle("Radius [sm]"); fh_r_points_true->Draw("colz"); c5->cd(6); fh_r_points_fake->GetXaxis()->SetTitle("Number of Hits"); fh_r_points_fake->GetYaxis()->SetTitle("Radius [sm]"); fh_r_points_fake->Draw("colz"); }