macros/leptoquarks_analysis/event_counts.C

0001 int event_counts()
0002 {
0003   gStyle->SetOptStat(0);
0004
0005   unsigned col1 = kOrange+7;
0006   unsigned col2 = kBlue+2;
0007
0008   /* open inout files and merge trees */
0009   TChain chain_event("ntp_event");
0010   chain_event.Add("data_3pions/p250_e20_0events_file1093_LeptoAna_TruthEvent.root");
0011   chain_event.Add("data_3pions/p250_e20_0events_file1096_LeptoAna_TruthEvent.root");
0012   chain_event.Add("data_3pions/p250_e20_0events_file1101_LeptoAna_TruthEvent.root");
0013   chain_event.Add("data_3pions/p250_e20_0events_file1115_LeptoAna_TruthEvent.root");
0014   chain_event.Add("data_3pions/p250_e20_0events_file1122_LeptoAna_TruthEvent.root");
0015   chain_event.Add("data_3pions/p250_e20_0events_file1127_LeptoAna_TruthEvent.root");
0016   chain_event.Add("data_3pions/p250_e20_0events_file1131_LeptoAna_TruthEvent.root");
0017   chain_event.Add("data_3pions/p250_e20_0events_file1164_LeptoAna_TruthEvent.root");
0018
0019   TChain chain("candidates");
0020   chain.Add("data_3pions/p250_e20_0events_file1093_LeptoAna_r05.root");
0021   chain.Add("data_3pions/p250_e20_0events_file1096_LeptoAna_r05.root");
0022   chain.Add("data_3pions/p250_e20_0events_file1101_LeptoAna_r05.root");
0023   chain.Add("data_3pions/p250_e20_0events_file1115_LeptoAna_r05.root");
0024   chain.Add("data_3pions/p250_e20_0events_file1122_LeptoAna_r05.root");
0025   chain.Add("data_3pions/p250_e20_0events_file1127_LeptoAna_r05.root");
0026   chain.Add("data_3pions/p250_e20_0events_file1131_LeptoAna_r05.root");
0027   chain.Add("data_3pions/p250_e20_0events_file1164_LeptoAna_r05.root");
0028
0029   /* particle selection */
0030   cout << "TOTAL events: " << chain_event.GetEntries() << endl;
0031   cout << "  && tau within abs(eta)<1: " << chain_event.GetEntries("abs(tau_eta) < 1.0") << endl;
0032
0033   cout << "TOTAL tau found: " << chain.GetEntries("evtgen_is_tau == 1") << endl;
0034   cout << "  && tau within abs(jet_eta)<1 && jet_pT > 5 GeV: " << chain.GetEntries("evtgen_is_tau == 1 && abs(jet_eta) < 1.0 && jet_ptrans > 5") << endl;
0035   cout << "  && Delte_R < 0.1: " << chain.GetEntries("evtgen_is_tau == 1 && abs(jet_eta) < 1.0 && jet_ptrans>5 && sqrt( (evtgen_tau_eta-jet_eta)*(evtgen_tau_eta-jet_eta) + (evtgen_tau_phi-jet_phi)*(evtgen_tau_phi-jet_phi) ) < 0.1") << endl;
0036
0037   cout << endl;
0038   cout << endl;
0039   cout << endl;
0040
0041   cout << "OTHER JETS found within acceptance: " << chain.GetEntries("evtgen_is_tau == 0 && abs(jet_eta) < 1.0 && jet_ptrans>5") << endl;
0042
0043
0044
0045   /* Plot tau angle reconstruction quality */
0046   TCanvas *c1 = new TCanvas();
0047   chain.Draw("sqrt( (evtgen_tau_eta-jet_eta)*(evtgen_tau_eta-jet_eta) + (evtgen_tau_phi-jet_phi)*(evtgen_tau_phi-jet_phi) ) >> h1(50,0,0.5)", "evtgen_is_tau == 1 && abs(jet_eta) < 1.0 && jet_ptrans > 5" );
0048   h1->GetXaxis()->SetTitle("#DeltaR (#tau_{jet} - #tau_{true})");
0049   c1->Print("plots/event_counts_deltaR.eps");
0050   c1->Print("plots/event_counts_deltaR.png");
0051
0052
0053   return 0;
0054 }