Return to jeturnon_tau_L3.C CVS log

Up to [Nicadd] / ttbar / p20_taujets_note / tau_jettrig

Initial revision

void jeturnon_tau_L3()
{
  gROOT->SetStyle("Plain");
  gStyle->SetOptStat(0);
  
  TFile* jettrigeff = new TFile("jettrigeff.root");
  TFile* jettrigeff_tau = new TFile("jettrigeff_tau.root");

  jettrigeff->cd();
  TH1F* MUJ2_2J20LM3DR3_2_dist_L2_pt = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_2_dist_L2_pt")->Clone();
  TH1F* MUJ2_2J20LM3DR3_2_dist_L2runL3_pt = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_2_dist_L2runL3_pt")->Clone();
  
  jettrigeff_tau->cd();
  TH1F* MUJ2_2J20LM3DR3_2_dist_L2_pt_tau = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_2_dist_L2_pt")->Clone();
  TH1F* MUJ2_2J20LM3DR3_2_dist_L2runL3_pt_tau = (TH1F*)gROOT->FindObject("MUJ2_2J20LM3DR3_2_dist_L2runL3_pt")->Clone();
  
 
  Int_t n = MUJ2_2J20LM3DR3_2_dist_L2_pt->GetNbinsX();
  Double_t xmin = MUJ2_2J20LM3DR3_2_dist_L2_pt->GetXaxis()->GetXmin();
  Double_t xmax = MUJ2_2J20LM3DR3_2_dist_L2_pt->GetXaxis()->GetXmax();
  //cout<<xmin<<endl;
  TH1F* h_eff = new TH1F("eff","L3 turnon (SC5JET_9_PV3 20)",n,xmin,xmax);
  TH1F* h_eff_tau = new TH1F("eff_tau","L3 turnon (SC5JET_9_PV3 20)",n,xmin,xmax);

  h_eff->Divide(MUJ2_2J20LM3DR3_2_dist_L2runL3_pt,MUJ2_2J20LM3DR3_2_dist_L2_pt,1,1,"B");
  h_eff->SetMinimum(0.0);
  h_eff->SetMaximum(1.05);
  h_eff->SetMarkerStyle(20);
  h_eff->SetMarkerColor(2);
  h_eff->GetXaxis()->SetTitle("pT of matched jet, GeV");
  h_eff->Draw("pe0");

  h_eff_tau->Divide(MUJ2_2J20LM3DR3_2_dist_L2runL3_pt_tau,MUJ2_2J20LM3DR3_2_dist_L2_pt_tau,1,1,"B");
  h_eff_tau->SetMinimum(0.0);
  h_eff_tau->SetMaximum(1.05);
  h_eff_tau->GetXaxis()->SetTitle("pT of matched jet, GeV");
  h_eff_tau->Draw("pe0same"); //jettrigeff->Close();
  //jettrigeff_tau->Close();
}