Skip to content

Commit

Permalink
Implement btag cut + new commands
Browse files Browse the repository at this point in the history
  • Loading branch information
Norraphat committed Dec 19, 2013
1 parent b06a806 commit a08ac86
Show file tree
Hide file tree
Showing 12 changed files with 2,274 additions and 182 deletions.
79 changes: 55 additions & 24 deletions NtupleAnalyzer/bin/AnaMonoJet.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -26,17 +26,20 @@ using namespace std;

int main(int argc, char ** argv)
{
if(argc < 6){
if(argc < 7){
cerr << "Program need more than this parameter " << endl;
cerr << "Example: Analysis sampleName(ex:wjets or met or ...) maxEvents isMC(0 or 1) cuts(jet,met,jetmet) folder(./results/...)" << endl;
cerr << "Example: Analysis sampleName(ex:wjets or met or ...) maxEvents isMC(0 or 1) cuts(jet,met,jetmet,btagjet,btagmet,btagjetmet) jetThreshold folder(./results/...)" << endl;
return 1;
}

//cout<<"argc = "<<argc<<" | argv[0] = "<<argv[0]<<", argv[4] = "<<argv[4]<<endl;

string anaout = Constants::outputDir;
anaout += "results/";
anaout += argv[5];
anaout += argv[6];

int secJetCut;
sscanf(argv[5], "%d", &secJetCut);

int varycutIndex = 0;
string varycut = "";
Expand All @@ -45,8 +48,10 @@ int main(int argc, char ** argv)
varycutIndex += 1;
if(varycut.find("jet") != std::string::npos || varycut.find("Jet") != std::string::npos || varycut.find("JET") != std::string::npos)
varycutIndex += 2;
if(varycut.find("btag") != std::string::npos || varycut.find("BTag") != std::string::npos || varycut.find("BTAG") != std::string::npos)
varycutIndex += 4;
if(varycutIndex<=0){
cerr << "Check the defined cuts. It should be [jet,met,jetmet]" << endl;
cerr << "Check the defined cuts. It should be [jet,met,jetmet,btagjet,btagmet,btagjetmet]" << endl;
return 1;
}

Expand All @@ -56,21 +61,43 @@ int main(int argc, char ** argv)
int nev;
if ( argc >= 3 ) sscanf ( argv[2], "%d", &nev );
else nev = 100000000;

int wjet_zjet = 0;
std::string dataset = "";
dataset += argv[1];
if(dataset.find("zjets") != std::string::npos ||dataset.find("w4jets") != std::string::npos||dataset.find("wjets") != std::string::npos){
wjet_zjet += 1;
cout<<"Including Gamma Cut"<<endl;
}

std::string logFileName = anaout + "/";
logFileName += argv[1];
logFileName += "_";
logFileName += argv[4];
logFileName += "_";
logFileName += argv[5];
logFileName += ".log";

Manager manager(logFileName);

string histFile = anaout + "/";
histFile += argv[1];
histFile += "_";
histFile += argv[4];
histFile += "_";
histFile += argv[5];

cout << "Starting to Run ........" << endl;
cout << "Vary cuts by ";
cout << " - Vary cuts by ";
if(varycutIndex==1) cout<<"MET";
if(varycutIndex==2) cout<<"Jet";
if(varycutIndex==3) cout<<"JetMET";
if(varycutIndex==4) cout<<"BTag";
if(varycutIndex==5) cout<<"BTag + MET";
if(varycutIndex==6) cout<<"BTag + Jet";
if(varycutIndex==7) cout<<"BTag + JetMET";
cout<<endl;
cout << " - Count only jets which pt > "<<secJetCut<<" GeV";
cout<<endl;

//---------------------------------------CUTS & Fill Hist--------------------------------------------------------------
Expand All @@ -84,10 +111,12 @@ int main(int argc, char ** argv)

CutHLT CHLT0(0); // Primary vertex
CutHLT CHLT1(1); // HLT
CutAbnormalEvents CAbnormalEvents(abrun, abevt);
CutAbnormalEvents CAbnormalEvents(abrun, abevt);
CutGamma CGamma(5.0);
manager.Add(&CHLT1);
manager.Add(&CHLT0);
manager.Add(&CAbnormalEvents);
if(wjet_zjet==1) manager.Add(&CGamma);
hDataMcMatching DataMcMatching0(histFile+"_AnaMonoJet_0.root");
manager.Add(&DataMcMatching0);

Expand Down Expand Up @@ -132,8 +161,10 @@ int main(int argc, char ** argv)

CutJet1 CJet1(110 , 2.4, 0.02, 0.98);
manager.Add(&CJet1);
//CutJet2 CJet20(30,4.7);
//manager.Add(&CJet20);
CutJet1BTag CJet1BTag(0.679);
if(varycutIndex==4 || varycutIndex==5 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet1BTag);
}
hDataMcMatching DataMcMatching3(histFile+"_AnaMonoJet_3.root");
manager.Add(&DataMcMatching3);
hHLTEff HLTEff1(histFile+"_HLTEff1.root");
Expand Down Expand Up @@ -174,10 +205,10 @@ int main(int argc, char ** argv)

CutMet CMet2(250);
CutJet1 CJet2(250, 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet2);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CMet2);
manager.Add(&CJet2);
}
Expand All @@ -187,10 +218,10 @@ int main(int argc, char ** argv)

CutMet CMet3(300);
CutJet1 CJet3(300, 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet3);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet3);
}
hDataMcMatching DataMcMatching10(histFile+"_AnaMonoJet_10.root");
Expand All @@ -201,10 +232,10 @@ int main(int argc, char ** argv)

CutMet CMet4(350);
CutJet1 CJet4(350 , 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet4);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet4);
}
hDataMcMatching DataMcMatching11(histFile+"_AnaMonoJet_11.root");
Expand All @@ -213,10 +244,10 @@ int main(int argc, char ** argv)

CutMet CMet5(400);
CutJet1 CJet5(400, 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet5);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet5);
}
hDataMcMatching DataMcMatching12(histFile+"_AnaMonoJet_12.root");
Expand All @@ -225,10 +256,10 @@ int main(int argc, char ** argv)

CutMet CMet6(450);
CutJet1 CJet6(450, 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet6);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet6);
}
hDataMcMatching DataMcMatching13(histFile+"_AnaMonoJet_13.root");
Expand All @@ -237,10 +268,10 @@ int main(int argc, char ** argv)

CutMet CMet7(500);
CutJet1 CJet7(500, 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet7);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet7);
}
hDataMcMatching DataMcMatching14(histFile+"_AnaMonoJet_14.root");
Expand All @@ -249,20 +280,20 @@ int main(int argc, char ** argv)

CutMet CMet8(550);
CutJet1 CJet8(550 , 2.4, 0.02, 0.98);
if(varycutIndex==1 || varycutIndex==3){
if(varycutIndex==1 || varycutIndex==3 || varycutIndex==5 || varycutIndex==7){
manager.Add(&CMet8);
}
if(varycutIndex==2 || varycutIndex==3){
if(varycutIndex==2 || varycutIndex==3 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet8);
}
hDataMcMatching DataMcMatching15(histFile+"_AnaMonoJet_15.root");
manager.Add(&DataMcMatching15);

//-------------------------------------------------------------------------------------------------------------------------

cout << "Running over sample " << argv[1] << endl;
cout << " - Running over sample " << argv[1] << endl;

EventData eventData(argv[1], nev, isMC);
EventData eventData(argv[1], nev, isMC, (float)secJetCut);

// Loop over events
manager.Run(eventData);
Expand Down
51 changes: 42 additions & 9 deletions NtupleAnalyzer/bin/AnaMonoJetCutEff.cc
View file
Original file line number Diff line number Diff line change
Expand Up @@ -25,23 +25,28 @@ using namespace std;

int main(int argc, char ** argv)
{
if(argc < 7){
if(argc < 8){
cerr << "Program need more than this parameter " << endl;
cerr << "Example: Analysis sampleName(ex:wjets or met or ...) maxEvents cutNumber isMC(0 or 1) cuts(jet,met,jetmet) folder(./results/...)" << endl;
cerr << "Example: Analysis sampleName(ex:wjets or met or ...) maxEvents cutNumber isMC(0 or 1) cuts(jet,met,jetmet,btagjet,btagmet,btagjetmet) jetThreshold folder(./results/...)" << endl;
return 1;
}

string anaout = Constants::outputDir;
anaout += "results/";
anaout += argv[6];
anaout += argv[7];

int secJetCut;
sscanf(argv[6], "%d", &secJetCut);

int varycutIndex = 0;
string varycut = "";
varycut += argv[5];
if(varycut.find("met") != std::string::npos || varycut.find("Met") != std::string::npos || varycut.find("MET") != std::string::npos)
varycutIndex += 1;
if(varycut.find("jet") != std::string::npos || varycut.find("Jet") != std::string::npos || varycut.find("JET") != std::string::npos)
varycutIndex += 2;
if(varycut.find("btag") != std::string::npos || varycut.find("BTag") != std::string::npos || varycut.find("BTAG") != std::string::npos)
varycutIndex += 4;
if(varycutIndex<=0){
cerr << "Check the defined cuts. It should be [jet,met,jetmet]" << endl;
return 1;
Expand All @@ -56,22 +61,44 @@ int main(int argc, char ** argv)

int nev;
if ( argc >= 5 ) sscanf ( argv[2], "%d", &nev );
else nev = 100000000;
else nev = 100000000;

int wjet_zjet = 0;
std::string dataset = "";
dataset += argv[1];
if(dataset.find("zjets") != std::string::npos ||dataset.find("w4jets") != std::string::npos||dataset.find("wjets") != std::string::npos){
wjet_zjet += 1;
cout<<"Including Gamma Cut"<<endl;
}

std::string logFileName = anaout + "/";
logFileName += argv[1];
logFileName += "_";
logFileName += argv[5];
logFileName += "_";
logFileName += argv[6];
logFileName += ".log";

Manager manager(logFileName);

string histFile = anaout + "/";
histFile += argv[1];
histFile += argv[1];
histFile += "_";
histFile += argv[5];
histFile += "_";
histFile += argv[6];

cout << "Starting to Run ........" << endl;
cout << "Vary cuts by ";
cout << " - Vary cuts by ";
if(varycutIndex==1) cout<<"MET";
if(varycutIndex==2) cout<<"Jet";
if(varycutIndex==3) cout<<"JetMET";
if(varycutIndex==4) cout<<"BTag";
if(varycutIndex==5) cout<<"BTag + MET";
if(varycutIndex==6) cout<<"BTag + Jet";
if(varycutIndex==7) cout<<"BTag + JetMET";
cout<<endl;
cout << " - Count only jets which pt > "<<secJetCut<<" GeV";
cout<<endl;


Expand All @@ -81,6 +108,7 @@ int main(int argc, char ** argv)
vector<int> abevt;

CutAbnormalEvents CAbnormalEvents(abrun , abevt);
CutGamma CGamma(5.0);
CutHLT CHLT(0);
CutHLT CHLT1(1); //HLT
CutHLT CHLT2(2); //No scraping
Expand All @@ -99,6 +127,7 @@ int main(int argc, char ** argv)
CutNoiseClean CNoiseClean(0.95 , 0.98, 1, 0.01, 0.99);
CutJet1 CJet1(110 , 2.4, 0.02, 0.98);
CutJet1 CJet2(250 , 2.4, 0.02, 0.98);
CutJet1BTag CJet1BTag(0.679);
CutNJet CNJet(3);
CutDeltaPhi3 CDeltaPhi3(2.5);
CutMet CMet(250);
Expand All @@ -107,7 +136,8 @@ int main(int argc, char ** argv)
NoPFElec CNoPFElec(10., 66.);

manager.Add(&CAbnormalEvents);
manager.Add(&CHLT);
manager.Add(&CHLT);
if(wjet_zjet==1) manager.Add(&CGamma);
/*
manager.Add(&CHLT1); //HLT
manager.Add(&CHLT2); //No scraping
Expand All @@ -130,6 +160,9 @@ int main(int argc, char ** argv)
manager.Add(&CJet1);
else
manager.Add(&CJet2);
if(varycutIndex==4 || varycutIndex==5 || varycutIndex==6 || varycutIndex==7){
manager.Add(&CJet1BTag);
}
}
if(sel!=3) manager.Add(&CNJet);
if(sel!=4) manager.Add(&CDeltaPhi3);
Expand All @@ -142,9 +175,9 @@ int main(int argc, char ** argv)

//-------------------------------------------------------------------------------------------------------------------------

cout << "Running over sample " << argv[1] << endl;
cout << " - Running over sample " << argv[1] << endl;

EventData eventData(argv[1], nev, isMC);
EventData eventData(argv[1], nev, isMC, (float)secJetCut);

// Loop over events
manager.Run(eventData);
Expand Down
Loading

0 comments on commit a08ac86

Please sign in to comment.