macros/diskinematics_analysis/evaluate_cluster_track_matching.C

0001 int
0002 evaluate_cluster_track_matching()
0003 {
0004   TString fname( "eicana_disreco_test.root" );
0005
0006   TFile *fin = new TFile( fname , "OPEN" );
0007
0008   TTree *event = (TTree*)fin->Get("event");
0009
0010   event->Print();
0011
0012   TCanvas *ctemp = new TCanvas();
0013   event->Draw("1","em_cluster_e>0");
0014   cout << "EM candidates total: " << htemp->GetEntries() << endl;
0015   event->Draw("1","em_cluster_e>1 && em_cluster_prob > 0.01");
0016   cout << "EM candidates E_cluster > 1 GeV && Prob_EM > 0.01: " << htemp->GetEntries() << endl;
0017   event->Draw("1","em_cluster_e>1 && em_cluster_prob > 0.01 && em_track_cluster_dr < 0.01");
0018   cout << "EM candidates E_cluster > 1 GeV && Prob_EM > 0.01 && dr_cluster_track < 0.01: " << htemp->GetEntries() << endl;
0019
0020   TCanvas *c0 = new TCanvas();
0021   event->Draw("em_track_cluster_dr",
0022           "em_cluster_e > 1");
0023
0024   TCanvas *c1 = new TCanvas();
0025   event->Draw("1",//(-1*log( tan( em_cluster_theta / 2. ) ) )",// - em_evtgen_eta",
0026           "em_cluster_e>1 && em_cluster_prob > 0.01");
0027
0028   TCanvas *c2 = new TCanvas();
0029   event->Draw("(-1*log( tan( em_track_theta / 2. ) ) ) - em_evtgen_eta : em_evtgen_eta",
0030           "em_cluster_e>1 && em_cluster_prob > 0.01 && em_track_cluster_dr < 0.01");
0031
0032 //.x Fun4All_EICAnalysis_DISReco.C
0033 //TFile *_file0 = TFile::Open()
0034 //event->Scan("-1*log(tan(em_cluster_theta/2.)):em_cluster_e:-1*log(tan(em_track_theta/2.)):em_track_ptotal","em_cluster_e>2")
0035 //event->Scan("-1*log(tan(em_cluster_theta/2.)):em_cluster_e:-1*log(tan(em_track_theta/2.)):em_track_ptotal","em_cluster_e>2&&em_evtgen_pid==-11")
0036 //event->Scan("-1*log(tan(em_cluster_theta/2.)):em_cluster_e:-1*log(tan(em_track_theta/2.)):em_track_ptotal","em_cluster_e>2&&em_evtgen_pid==11")
0037 //event->Draw("-1*log(tan(em_cluster_theta/2.)):-1*log(tan(em_track_theta/2.))","em_cluster_e>2&&em_evtgen_pid==11")
0038 //event->Draw("-1*log(tan(em_cluster_theta/2.)):em_evtgen_theta","em_cluster_e>2&&em_evtgen_pid==11")
0039 //event->Draw("-1*log(tan(em_cluster_theta/2.)):em_evtgen_eta","em_cluster_e>2&&em_evtgen_pid==11")
0040 //event->Draw("-1*log(tan(em_cluster_theta/2.))-em_evtgen_eta","em_cluster_e>2&&em_evtgen_pid==11")
0041 //event->Draw("-1*log(tan(em_track_theta/2.))-em_evtgen_eta","em_cluster_e>2&&em_evtgen_pid==11")
0042 //event->Draw("-1*log(tan(em_track_theta/2.))-em_evtgen_eta>>h1(100,0,1)","em_cluster_e>2&&em_evtgen_pid==11")
0043 //event->Draw("-1*log(tan(em_track_theta/2.))-em_evtgen_eta>>h1(100,-0.2,0.2)","em_cluster_e>2&&em_evtgen_pid==11")
0044 //event->Draw("-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))","em_cluster_e>0")
0045 //event->Draw("-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))","em_cluster_e>0&& abs(-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))) < 0.1")
0046 //event->Draw("-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))","em_cluster_e>0&& abs(-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))) < 0.1")
0047 //event->Draw("-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))","em_cluster_e>0&& abs(-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))) < 0.01")
0048 //event->Draw("(-1*log(tan(em_track_theta/2.)))-em_evtgen_eta","em_cluster_e>0&& abs(-1*log(tan(em_cluster_theta/2.))-(-1*log(tan(em_track_theta/2.)))) < 0.01")
0049
0050 //.q
0051
0052   return 0;
0053 }