Return to jeturnon_tau_L2.C CVS log

Up to [Nicadd] / ttbar / p20_taujets_note / tau_jettrig

Annotation of ttbar/p20_taujets_note/tau_jettrig/jeturnon_tau_L2.C, revision 1.1.1.1

1.1       uid12904    1: void jeturnon_tau_L2()
                      2: {
                      3:   gROOT->SetStyle("Plain");
                      4:   gStyle->SetOptStat(0);
                      5:   
                      6:   TFile* jettrigeff = new TFile("jettrigeff.root");
                      7:   TFile* jettrigeff_tau = new TFile("jettrigeff_tau.root");
                      8: 
                      9:   jettrigeff->cd();
                     10:   TH1F* MUJ2_2J20LM3DR3_1_dist_L1_pt = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_1_dist_L1_pt")->Clone();
                     11:   TH1F* MUJ2_2J20LM3DR3_1_dist_L2_pt = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_1_dist_L2_pt")->Clone();
                     12:   
                     13:   jettrigeff_tau->cd();
                     14:   TH1F* MUJ2_2J20LM3DR3_1_dist_L1_pt_tau = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_1_dist_L1_pt")->Clone();
                     15:   TH1F* MUJ2_2J20LM3DR3_1_dist_L2_pt_tau = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_1_dist_L2_pt")->Clone();
                     16:   
                     17:  
                     18:   Int_t n = MUJ2_2J20LM3DR3_1_dist_L1_pt->GetNbinsX();
                     19:   Double_t xmin = MUJ2_2J20LM3DR3_1_dist_L1_pt->GetXaxis()->GetXmin();
                     20:   Double_t xmax = MUJ2_2J20LM3DR3_1_dist_L1_pt->GetXaxis()->GetXmax();
                     21:   
                     22:   
                     23:   TH1F* h_eff = new TH1F("eff","L2 turnon (JET(0,5)8)",n,xmin,xmax);
                     24:   TH1F* h_eff_tau = new TH1F("eff_tau","L2 turnon (JET(0,5)8)",n,xmin,xmax);
                     25: 
                     26:   h_eff->Divide(MUJ2_2J20LM3DR3_1_dist_L2_pt,MUJ2_2J20LM3DR3_1_dist_L1_pt,1,1,"B");
                     27:   h_eff->SetMinimum(0.0);
                     28:   h_eff->SetMaximum(1.05);
                     29:   h_eff->SetMarkerStyle(20);
                     30:   h_eff->SetMarkerColor(2);
                     31:   h_eff->GetXaxis()->SetTitle("pT of matched jet, GeV");
                     32:   h_eff->Draw("pe0");
                     33: 
                     34:   h_eff_tau->Divide(MUJ2_2J20LM3DR3_1_dist_L2_pt_tau,MUJ2_2J20LM3DR3_1_dist_L1_pt_tau,1,1,"B");
                     35:   h_eff_tau->SetMinimum(0.0);
                     36:   h_eff_tau->SetMaximum(1.05);
                     37:   h_eff_tau->GetXaxis()->SetTitle("pT of matched jet, GeV");
                     38:   h_eff_tau->Draw("pe0same"); //jettrigeff->Close();
                     39:   //jettrigeff_tau->Close();
                     40: }