macros/leptoquarks_analysis/cone_energies.C

0001 /**
0002  * Plot energies in cones around jet axis
0003  *
0004  * Written by nils.feege@stonybrook.edu
0005  */
0006
0007 int cone_energies()
0008 {
0009   gStyle->SetOptStat(0);
0010
0011   unsigned col1 = kOrange+7;
0012   unsigned col2 = kBlue+2;
0013
0014   /* open inout files and merge trees */
0015   TChain chain("candidates");
0016   chain.Add("data_3pions/p250_e20_0events_file1093_LeptoAna_r05.root");
0017   chain.Add("data_3pions/p250_e20_0events_file1096_LeptoAna_r05.root");
0018   chain.Add("data_3pions/p250_e20_0events_file1101_LeptoAna_r05.root");
0019   chain.Add("data_3pions/p250_e20_0events_file1115_LeptoAna_r05.root");
0020   chain.Add("data_3pions/p250_e20_0events_file1122_LeptoAna_r05.root");
0021   chain.Add("data_3pions/p250_e20_0events_file1127_LeptoAna_r05.root");
0022   chain.Add("data_3pions/p250_e20_0events_file1131_LeptoAna_r05.root");
0023   chain.Add("data_3pions/p250_e20_0events_file1164_LeptoAna_r05.root");
0024
0025   /* Create temporary canvas */
0026   TCanvas* ctemp = new TCanvas();
0027
0028   ctemp->cd();
0029
0030   TH1F* h_uds = new TH1F( "prof_uds", "", 5, 0, 0.5);
0031   h_uds->GetXaxis()->SetNdivisions(505);
0032   TH1F* h_tau = new TH1F( "prof_tau", "", 5, 0, 0.5);
0033
0034   for ( int r = 1; r <= 5; r++ )
0035     {
0036       /* particle selection */
0037       TCut select_true_uds( Form( "jetshape_econe_r0%d * ( evtgen_is_uds==1 && abs(evtgen_uds_eta)<0.8 )", r ) );
0038       TCut select_true_tau( Form( "jetshape_econe_r0%d * ( evtgen_is_tau==1 && evtgen_tau_decay_prong==3 && abs(evtgen_tau_eta)<0.8 )", r ) );
0039
0040       chain.Draw( Form("0.9*(%f) >>+ prof_uds", r*0.1), select_true_uds );
0041       chain.Draw( Form("0.9*(%f) >>+ prof_tau", r*0.1), select_true_tau );
0042     }
0043
0044   h_uds->Scale(1./h_uds->Integral());
0045   h_uds->GetXaxis()->SetTitle( "#DeltaR" );
0046   h_uds->GetYaxis()->SetRangeUser(0, 0.5);
0047   h_uds->SetLineColor(col1);
0048   h_uds->SetFillColor(col1);
0049   h_uds->SetFillStyle(1001);
0050
0051   h_tau->Scale(1./h_tau->Integral());
0052   h_tau->SetLineColor(col2);
0053   h_tau->SetFillColor(col2);
0054   h_tau->SetFillStyle(0);
0055
0056   /* create Canvas and draw histograms */
0057   TCanvas *c1 = new TCanvas();
0058
0059   h_uds->DrawClone("");
0060   h_tau->DrawClone("sames");
0061
0062   gPad->RedrawAxis();
0063
0064   c1->Print("plots/cone_energies.eps");
0065   c1->Print("plots/cone_energies.png");
0066
0067   return 0;
0068 }