{
TH1F* hjt20 = J20_dist_pt.Clone();
TH1F* hjt20trig = J20_dist_L1_pt.Clone();
TH1F* hjt20_int = J20_dist_pt.Clone();
TH1F* hjt20trig_int = J20_dist_L1_pt.Clone();
for (i=1;i<=60;i++) {
Float_t sum = 0;
for (ii=60;ii>=i;ii--) sum += hjt20trig->GetBinContent(ii);
hjt20trig_int->SetBinContent(i,sum);
}
Float_t sum = 0;
for (i=1;i<=60;i++) {
Float_t sum = 0;
for (ii=60;ii>=i;ii--) sum += hjt20->GetBinContent(ii);
hjt20_int->SetBinContent(i,sum);
}
hjt20trig_int->Sumw2(); hjt20_int->Sumw2();
TH1F* heffint = (TH1F*) hjt20_int->Clone("heffint");
heffint->Divide(hjt20,hjt20trig,1,1,"B");
heffint->SetTitle("Trigger efficiency for JT20");
heffint->GetXaxis()->SetTitle("pT of matched jet, GeV");
heffint->Draw();
}
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>