Skip to content

Dev/11 0 x #5

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.

Sign up for GitHub

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

Open
selvaggi wants to merge 4 commits into
base: dev/11_0_X
Choose a base branch
from
Open

Dev/11 0 x #5

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
41e25d1
Remove buggy BulkGrav samples.
hqucms Nov 22, 2017
54d9c74
Update README.md
hqucms Nov 22, 2017
a548cdb
Merge remote-tracking branch 'upstream/master' into dev/11_0_X
selvaggi Jan 3, 2020
ac1de32
added generic targetmass to fatjet ntuples
selvaggi Jan 10, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 8 additions & 4 deletions Ntupler/src/FatJetInfoFiller.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -91,10 +91,14 @@ void FatJetInfoFiller::book() {
data.add<float>("fj_genOverReco_mass_null", 0); // default to 0 if not gen-matched
data.add<float>("fj_genjet_sdmass", 0);
data.add<float>("fj_genjet_sdmass_sqrt", 0);
data.add<float>("fj_genjet_targetmass", 0);
data.add<float>("fj_genOverReco_sdmass", 1); // default to 1 if not gen-matched
data.add<float>("fj_genOverReco_sdmass_null", 0); // default to 0 if not gen-matched
// ----------------------------------

// --- target mass for energy regression ---
data.add<float>("fj_targetgenmass", 0);

// fatjet kinematics
data.add<float>("fj_pt", 0);
data.add<float>("fj_eta", 0);
Expand Down Expand Up @@ -142,7 +146,6 @@ void FatJetInfoFiller::book() {
data.add<float>("fj_relptdiff", 0);
data.add<float>("fj_sdn2", 0);


//double-b
data.add<float>("fj_doubleb", 0);

Expand Down Expand Up @@ -246,7 +249,8 @@ bool FatJetInfoFiller::fill(const pat::Jet& jet, size_t jetidx, const JetHelper&
// gen-matched particle (top/W/etc.)
data.fill<float>("fj_gen_pt", fjlabel.second ? fjlabel.second->pt() : -999);
data.fill<float>("fj_gen_eta", fjlabel.second ? fjlabel.second->eta() : -999);
data.fill<float>("fj_gen_mass", (fjlabel.first < FatJetMatching::QCD_all && fjlabel.second) ? fjlabel.second->mass() : 0);
float gen_mass = (fjlabel.first < FatJetMatching::QCD_all && fjlabel.second) ? fjlabel.second->mass() : 0.;
data.fill<float>("fj_gen_mass", gen_mass);
data.fill<float>("fj_gen_deltaR", fjlabel.second ? reco::deltaR(jet, fjlabel.second->p4()) : 999);

// ----------------------------------
Expand Down Expand Up @@ -323,13 +327,13 @@ bool FatJetInfoFiller::fill(const pat::Jet& jet, size_t jetidx, const JetHelper&
if (sdgenjet){
// jet here points to the uncorrected jet
auto pos = [](double x){ return x<0 ? 0 : x; };
data.fill<float>("fj_genjet_sdmass", pos(sdgenjet->mass()));
data.fill<float>("fj_genjet_sdmass", pos(sdgenjet->mass()));
data.fill<float>("fj_genjet_targetmass",(gen_mass > 6.) ? gen_mass : pos(sdgenjet->mass()) );
data.fill<float>("fj_genjet_sdmass_sqrt", std::sqrt(pos(sdgenjet->mass())));
data.fill<float>("fj_genOverReco_sdmass", catchInfs(pos(sdgenjet->mass()) / pos(msd_uncorr), 1));
data.fill<float>("fj_genOverReco_sdmass_null", catchInfs(pos(sdgenjet->mass()) / pos(msd_uncorr), 0));
}


// --------
// double-b

Expand Down