void jeturnon_tau()
{
gROOT->SetStyle("Plain");
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
Int_t icol[4] = {1,2,4,3};
Char_t ch_name[] = "t000";
const Int_t nf = 1;
TFile* jettrigeff = new TFile("jettrigeff.root");
TFile* jettrigeff_tau = new TFile("jettrigeff_tau.root");
jettrigeff->cd();
TH1F* J20_dist_pt = (TH1F*)gROOT->FindObject("J20_dist_pt")->Clone();
TH1F* J20_dist_L1_pt = (TH1F*)gROOT->FindObject("J20_dist_L1_pt")->Clone();
jettrigeff_tau->cd();
TH1F* J20_dist_pt_tau = (TH1F*)gROOT->FindObject("J20_dist_pt")->Clone();
TH1F* J20_dist_L1_pt_tau = (TH1F*)gROOT->FindObject("J20_dist_L1_pt")->Clone();
Int_t n = J20_dist_pt->GetNbinsX();
Double_t xmin = J20_dist_pt->GetXaxis()->GetXmin();
Double_t xmax = J20_dist_pt->GetXaxis()->GetXmax();
//cout<<xmin<<endl;
TH1F* h_eff = new TH1F("L1 eff","eff",n,xmin,xmax);
TH1F* h_eff_tau = new TH1F("L1 eff for tau-matched ","eff_tau",n,xmin,xmax);
h_eff->Divide(J20_dist_pt,J20_dist_L1_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(J20_dist_pt_tau,J20_dist_L1_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();
}
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>