| 1 |
#include "../MitHtt/NtupleDefs/interface/HiggsAnaDefs.hh"
|
| 2 |
|
| 3 |
plotMuFR(char* fname="NonMCBackground/root/musel-m10.root" ) {
|
| 4 |
|
| 5 |
TChain * t = new TChain("mutuple");
|
| 6 |
t->Add("NonMCBackground/root/musel-a05.root");
|
| 7 |
t->Add("NonMCBackground/root/musel-m10.root");
|
| 8 |
t->Add("NonMCBackground/root/musel-prv4.root");
|
| 9 |
t->Add("NonMCBackground/root/musel-prv6.root");
|
| 10 |
|
| 11 |
|
| 12 |
TH1F * hdenom = new TH1F("hdenom","hdenom",50,0,200);
|
| 13 |
TH1F * hnumer = new TH1F("hnumer","hnumer",50,0,200);
|
| 14 |
|
| 15 |
char cutstr[256];
|
| 16 |
sprintf( cutstr, "(trigger&0x%x||trigger&0x%x)&&(type==0x%x)", kHLT_Mu8, kHLT_Mu15, kTracker );
|
| 17 |
// sprintf( cutstr, "(trigger&0x%x||trigger&0x%x)&&(type&0x%x)", kHLT_Mu8, kHLT_Mu15, kGlobal );
|
| 18 |
t->Draw("pt>>hdenom", cutstr , "goff");
|
| 19 |
sprintf( cutstr, "%s&&(!failsID&&!failsIso)", cutstr );
|
| 20 |
t->Draw("pt>>hnumer", cutstr, "goff");
|
| 21 |
|
| 22 |
hdenom->Sumw2();
|
| 23 |
hnumer->Sumw2();
|
| 24 |
hdenom->SetLineWidth(2);
|
| 25 |
hnumer->SetLineWidth(2);
|
| 26 |
hnumer->SetLineColor(kRed);
|
| 27 |
|
| 28 |
|
| 29 |
TCanvas * c1 = new TCanvas("c1","c1", 0, 0, 500, 400);
|
| 30 |
TH1F * c1frame = (TH1F*)(c1->DrawFrame(0,0.0001,200,1.5*hdenom->GetMaximum()));
|
| 31 |
hdenom->Draw("histsame");
|
| 32 |
hnumer->Draw("histsame");
|
| 33 |
// c1->SetLogy();
|
| 34 |
|
| 35 |
|
| 36 |
TCanvas * c2 = new TCanvas("c2","c2", 500, 0, 500, 400);
|
| 37 |
// TH1F * c2frame = (TH1F*)(c2->DrawFrame( 0, 0, 200, 1.1*hnumer->GetMaximum()));
|
| 38 |
TH1F * mufr = (TH1F*)(hnumer->Clone());
|
| 39 |
mufr->Divide(mufr, hdenom, 1., 1., "B");
|
| 40 |
mufr->SetMarkerStyle(20);
|
| 41 |
mufr->Draw("pe");
|
| 42 |
|
| 43 |
}
|
