Return to jeturnon_tau.C CVS log

Up to [Nicadd] / ttbar / p20_taujets_note / tau_jettrig

Annotation of ttbar/p20_taujets_note/tau_jettrig/jeturnon_tau.C, revision 1.1

1.1     ! uid12904    1: void jeturnon_tau()
        !             2: {
        !             3:   gROOT->SetStyle("Plain");
        !             4:   gStyle->SetOptStat(0);
        !             5:   gStyle->SetOptTitle(0);
        !             6: 
        !             7:   Int_t icol[4] = {1,2,4,3};
        !             8: 
        !             9:   Char_t ch_name[] = "t000";
        !            10:   
        !            11:   const Int_t nf = 1;
        !            12:   
        !            13:   TFile* jettrigeff = new TFile("jettrigeff.root");
        !            14:   TFile* jettrigeff_tau = new TFile("jettrigeff_tau.root");
        !            15: 
        !            16:   jettrigeff->cd();
        !            17:   TH1F* J20_dist_pt = (TH1F*)gROOT->FindObject("J20_dist_pt")->Clone();
        !            18:   TH1F* J20_dist_L1_pt = (TH1F*)gROOT->FindObject("J20_dist_L1_pt")->Clone();
        !            19:   
        !            20:   jettrigeff_tau->cd();
        !            21:   TH1F* J20_dist_pt_tau = (TH1F*)gROOT->FindObject("J20_dist_pt")->Clone();
        !            22:   TH1F* J20_dist_L1_pt_tau = (TH1F*)gROOT->FindObject("J20_dist_L1_pt")->Clone();
        !            23:   
        !            24:  
        !            25:   Int_t n = J20_dist_pt->GetNbinsX();
        !            26:   Double_t xmin = J20_dist_pt->GetXaxis()->GetXmin();
        !            27:   Double_t xmax = J20_dist_pt->GetXaxis()->GetXmax();
        !            28:   //cout<<xmin<<endl;
        !            29:   TH1F* h_eff = new TH1F("L1 eff","eff",n,xmin,xmax);
        !            30:   TH1F* h_eff_tau = new TH1F("L1 eff for tau-matched ","eff_tau",n,xmin,xmax);
        !            31: 
        !            32:   h_eff->Divide(J20_dist_pt,J20_dist_L1_pt,1,1,"B");
        !            33:   h_eff->SetMinimum(0.0);
        !            34:   h_eff->SetMaximum(1.05);
        !            35:   h_eff->SetMarkerStyle(20);
        !            36:   h_eff->SetMarkerColor(2);
        !            37:   h_eff->GetXaxis()->SetTitle("pT of matched jet, GeV");
        !            38:   h_eff->Draw("pe0");
        !            39: 
        !            40:   h_eff_tau->Divide(J20_dist_pt_tau,J20_dist_L1_pt_tau,1,1,"B");
        !            41:   h_eff_tau->SetMinimum(0.0);
        !            42:   h_eff_tau->SetMaximum(1.05);
        !            43:   h_eff_tau->GetXaxis()->SetTitle("pT of matched jet, GeV");
        !            44:   h_eff_tau->Draw("pe0same"); //jettrigeff->Close();
        !            45:   //jettrigeff_tau->Close();
        !            46: }