macros/leptoquarks_analysis/tracks.C

0001 tracks()
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_track("ntp_track");
0010   chain_track.Add("data_3pions/p250_e20_0events_file1093_LeptoAna_r05.root");
0011   chain_track.Add("data_3pions/p250_e20_0events_file1096_LeptoAna_r05.root");
0012   chain_track.Add("data_3pions/p250_e20_0events_file1101_LeptoAna_r05.root");
0013   chain_track.Add("data_3pions/p250_e20_0events_file1115_LeptoAna_r05.root");
0014   chain_track.Add("data_3pions/p250_e20_0events_file1122_LeptoAna_r05.root");
0015   chain_track.Add("data_3pions/p250_e20_0events_file1127_LeptoAna_r05.root");
0016   chain_track.Add("data_3pions/p250_e20_0events_file1131_LeptoAna_r05.root");
0017   chain_track.Add("data_3pions/p250_e20_0events_file1164_LeptoAna_r05.root");
0018
0019   TCut select_tracks_tau("1 * (is_tau==1 && abs(jet_eta) < 1.0 && jet_etotal>5 && sqrt( (tau_eta-jet_eta)*(tau_eta-jet_eta) + (tau_phi-jet_phi)*(tau_phi-jet_phi) ) < 0.1)");
0020
0021   TCut select_tracks_uds("1 * (is_tau==0 && abs(jet_eta) < 1.0 && jet_etotal>5)");
0022
0023   TH2F* h2_tracks_tau = new TH2F("h2_tracks_tau", ";#Delta#phi;#Delta#eta;N/#SigmaN", 40, -0.2, 0.2, 40, -0.2, 0.2 );
0024   h2_tracks_tau->SetLineColor( col2 );
0025   h2_tracks_tau->GetXaxis()->SetNdivisions(505);
0026   h2_tracks_tau->GetYaxis()->SetNdivisions(505);
0027   h2_tracks_tau->GetZaxis()->SetNdivisions(505);
0028   h2_tracks_tau->GetZaxis()->SetTitleOffset(1.35);
0029   chain_track.Draw("( track_eta - jet_eta ) : ( track_phi - jet_phi ) >> h2_tracks_tau",select_tracks_tau);
0030   h2_tracks_tau->Scale( 1. / h2_tracks_tau->Integral() );
0031
0032   TH2F* h2_tracks_uds = new TH2F("h2_tracks_uds", ";#Delta#phi;#Delta#eta;N/#SigmaN", 40, -0.2, 0.2, 40, -0.2, 0.2 );
0033   h2_tracks_uds->SetLineColor( col1 );
0034   h2_tracks_uds->GetXaxis()->SetNdivisions(505);
0035   h2_tracks_uds->GetYaxis()->SetNdivisions(505);
0036   h2_tracks_uds->GetZaxis()->SetNdivisions(505);
0037   h2_tracks_uds->GetZaxis()->SetTitleOffset(1.35);
0038   chain_track.Draw("( track_eta - jet_eta ) : ( track_phi - jet_phi ) >> h2_tracks_uds",select_tracks_uds);
0039   h2_tracks_uds->Scale( 1. / h2_tracks_uds->Integral() );
0040
0041   /* plot histograms */
0042   TCanvas *c1 = new TCanvas();
0043   h2_tracks_tau->Draw("lego");
0044   c1->Print("plots/tracks_2d_tau.eps");
0045   c1->Print("plots/tracks_2d_tau.png");
0046
0047   TCanvas *c2 = new TCanvas();
0048   h2_tracks_uds->Draw("lego");
0049   c2->Print("plots/tracks_2d_uds.eps");
0050   c2->Print("plots/tracks_2d_uds.png");
0051 }