-
Notifications
You must be signed in to change notification settings - Fork 20
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
New scripts for DY stitching #88
Conversation
if (DEBUG) cout << "NBs = " << nbs << endl; | ||
|
||
float ht = theBigTree.lheHt; | ||
|
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
should we add here a protection like this?
if (ht >= 10000)
{
cout << --> ev: " << iEvent << " ht is " << ht << " lowered to 9999" << endl;
ht = 9999.;
}
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Done!
test/computDYsplitting_HT.cpp
Outdated
|
||
TFile* fOut = new TFile (outputFileName.c_str(), "recreate"); | ||
TH3F* h_nJets_nBs = new TH3F ("h_nJets_nBs" , ";nJets;nBs;ht", 5, xbins, 5, ybins, 9, zbins); | ||
TH3I* hINT_nJets_nBs = new TH3I ("hINT_nJets_nBs", ";nJets;nBs;ht", 5, xbins, 5, ybins, 9, zbins); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
can you change the names of the TH3 histograms here? something like h_nJets_nBs_HT
and hINT_nJets_nBs_HT
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Done!
Ciao @dzuolo, thanks! It looks ok to me, but did you try to test the result of the stitching? |
@dzuolo @camendola the file |
Created two new scripts to stitch DY samples, including the HT binned ones.
I'm quite confident on the changes made to test/computDYsplitting_HT.cpp but to the ones made to scripts/stitchWeightsDY_HT.py, I'm not sure if I understood correctly the meaning of fractionXS, please review these changes.