Skip to content

Fix next prod #36

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
swertz wants to merge 5 commits into
base: master
Choose a base branch
from
Open

Fix next prod #36

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
d3fe644
Tweaked neutrino solver
Apr 5, 2016
9a13ba2
Fixed object cuts, added option to remove ECAL gap
Apr 5, 2016
e5a03fa
Fixed jet matching; removed b/non-b separation
Apr 11, 2016
ef1a021
Separate HLT bit and match
Apr 11, 2016
cb437d6
Store the HLT SF in diLepton objects (pending adding MuonEG in the fr…
Apr 11, 2016
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
8 changes: 4 additions & 4 deletions interface/Defines.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2,10 +2,10 @@

// Debugging macros

// #define _TT_DEBUG_
#define TT_MTT_DEBUG (false)
#define TT_HLT_DEBUG (false)
#define TT_GEN_DEBUG (false)
//#define _TT_DEBUG_
//#define TT_MTT_DEBUG
//#define TT_HLT_DEBUG
//#define TT_GEN_DEBUG


#if TT_GEN_DEBUG
Expand Down
34 changes: 19 additions & 15 deletions interface/Indices.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -46,35 +46,39 @@ namespace TTAnalysis {
}

// Combination of Jet IDs for two jets (NOTE: NOT USED FOR NOW)
uint16_t JetJetID(const JetID::JetID& id1, const JetID::JetID& id2);
std::string JetJetIDStr(const JetID::JetID& id1, const JetID::JetID& id2);
//uint16_t JetJetID(const JetID::JetID& id1, const JetID::JetID& id2);
//std::string JetJetIDStr(const JetID::JetID& id1, const JetID::JetID& id2);

// B-tagging working points
namespace BWP {
enum BWP{ L, M, T, Count };
const std::array<BWP, Count> it = {{ L, M, T }};
const std::map<BWP, std::string> map = { {L, "L"}, {M, "M"}, {T, "T"} };
enum BWP{ A, L, M, Count };
const std::array<BWP, Count> it = {{ A, L, M }};
const std::map<BWP, std::string> map = { {A, "A"}, {L, "L"}, {M, "M"} };
}

// Combination of Jet ID and B-tagging working point (NOTE: NOT USED FOR NOW)
uint16_t JetIDBWP(const JetID::JetID& id, const BWP::BWP& wp);
std::string JetIDBWPStr(const JetID::JetID& id, const BWP::BWP& wp);
//uint16_t JetIDBWP(const JetID::JetID& id, const BWP::BWP& wp);
//std::string JetIDBWPStr(const JetID::JetID& id, const BWP::BWP& wp);

// Combination of Lepton ID + Lepton Isolation (one lepton) and B-tagging working point for one jet
uint16_t LepIDIsoJetBWP(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp);
std::string LepIDIsoJetBWPStr(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp);
// Combination of Lepton ID + Lepton Isolation (on lepton) and B-tagging working point for one jet (NOTE: NOT USED FOR NOW)
//uint16_t LepIDIsoJetBWP(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp);
//std::string LepIDIsoJetBWPStr(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp);

// Combination of Lepton ID + Lepton Isolation (two leptons) and B-tagging working point for one jet
uint16_t LepLepIDIsoJetBWP(const LepID::LepID& id1, const LepIso::LepIso& iso1, const LepID::LepID& id2, const LepIso::LepIso& iso2, const BWP::BWP& wp);
std::string LepLepIDIsoJetBWPStr(const LepID::LepID& id1, const LepIso::LepIso& iso1, const LepID::LepID& id2, const LepIso::LepIso& iso2, const BWP::BWP& wp);

// Combination of B-tagging working points for two jets
uint16_t JetJetBWP(const BWP::BWP& wp1, const BWP::BWP& wp2);
std::string JetJetBWPStr(const BWP::BWP& wp1, const BWP::BWP& wp2);

// Combination of Jet ID and B-tagging working points for two jets (NOTE: NOT USED FOR NOW)
uint16_t JetJetIDBWP(const JetID::JetID& id1, const BWP::BWP& wp1, const JetID::JetID& id2, const BWP::BWP& wp2);
std::string JetJetIDBWPStr(const JetID::JetID& id1, const BWP::BWP wp1, const JetID::JetID& id2, const BWP::BWP wp2);
//uint16_t JetJetIDBWP(const JetID::JetID& id1, const BWP::BWP& wp1, const JetID::JetID& id2, const BWP::BWP& wp2);
//std::string JetJetIDBWPStr(const JetID::JetID& id1, const BWP::BWP wp1, const JetID::JetID& id2, const BWP::BWP wp2);

// Combination of Lepton ID + Lepton Isolation (one lepton) and B-tagging working points for two jets
uint16_t LepIDIsoJetJetBWP(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp1, const BWP::BWP& wp2);
std::string LepIDIsoJetJetBWPStr(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp1, const BWP::BWP& wp2);
// Combination of Lepton ID + Lepton Isolation (one lepton) and B-tagging working points for two jets (NOTE: NOT USED FOR NOW)
//uint16_t LepIDIsoJetJetBWP(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp1, const BWP::BWP& wp2);
//std::string LepIDIsoJetJetBWPStr(const LepID::LepID& id, const LepIso::LepIso& iso, const BWP::BWP& wp1, const BWP::BWP& wp2);

// Combination of Lepton ID, Lepton Isolation, and B-tagging working points for a two-lepton-two-b-jets object
uint16_t LepLepIDIsoJetJetBWP(const LepID::LepID& id1, const LepIso::LepIso& iso1, const LepID::LepID& id2, const LepIso::LepIso& iso2, const BWP::BWP& wp1, const BWP::BWP& wp2);
Expand Down
11 changes: 4 additions & 7 deletions interface/NeutrinosSolver.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -24,18 +24,15 @@ class NeutrinosSolver {
public:
using LorentzVector = ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>>;

NeutrinosSolver(float top_mass, float w_mass):
t_mass(top_mass), w_mass(w_mass) {
NeutrinosSolver(){
// Empty
}

std::vector<std::pair<LorentzVector, LorentzVector>> getNeutrinos(const LorentzVector& lepton1_p4,
const LorentzVector& lepton2_p4,
const LorentzVector& bjet1_p4,
const LorentzVector& bjet2_p4,
const LorentzVector& met);

private:
float t_mass = 172.5;
float w_mass = 80.4;
const LorentzVector& met,
const double wplus_mass=80.4, const double wminus_mass=80.4,
const double t_mass=172.5, const double tbar_mass=172.5);
};
87 changes: 47 additions & 40 deletions interface/TTAnalyzer.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,6 +8,7 @@
#include <cp3_llbb/Framework/interface/MuonsProducer.h>
#include <cp3_llbb/Framework/interface/JetsProducer.h>
#include <cp3_llbb/Framework/interface/Analyzer.h>
#include <cp3_llbb/Framework/interface/BinnedValuesJSONParser.h>

#include <cp3_llbb/TTAnalysis/interface/Types.h>
#include <cp3_llbb/TTAnalysis/interface/Tools.h>
Expand All @@ -24,7 +25,8 @@ class TTAnalyzer: public Framework::Analyzer {
m_met_producer(config.getParameter<std::string>("metProducer")),

m_electronPtCut( config.getUntrackedParameter<double>("electronPtCut", 20) ),
m_electronEtaCut( config.getUntrackedParameter<double>("electronEtaCut", 2.5) ),
m_electronEtaCut( config.getUntrackedParameter<double>("electronEtaCut", 2.4) ),
m_electronRemoveGap( config.getUntrackedParameter<bool>("electronRemoveGap", true) ),
m_electronVetoIDName( config.getUntrackedParameter<std::string>("electronVetoIDName") ),
m_electronLooseIDName( config.getUntrackedParameter<std::string>("electronLooseIDName") ),
m_electronMediumIDName( config.getUntrackedParameter<std::string>("electronMediumIDName") ),
Expand All @@ -36,19 +38,27 @@ class TTAnalyzer: public Framework::Analyzer {
m_muonTightIsoCut( config.getUntrackedParameter<double>("muonTightIsoCut", 0.12) ),

m_jetPtCut( config.getUntrackedParameter<double>("jetPtCut", 30) ),
m_jetEtaCut( config.getUntrackedParameter<double>("jetEtaCut", 2.5) ),
m_jetEtaCut( config.getUntrackedParameter<double>("jetEtaCut", 2.4) ),
m_bJetEtaCut( config.getUntrackedParameter<double>("bJetEtaCut", 2.4) ),
m_jetPUID( config.getUntrackedParameter<double>("jetPUID", std::numeric_limits<float>::min()) ),
m_jetPUID( config.getUntrackedParameter<double>("jetPUID", std::numeric_limits<double>::min()) ),
m_jetDRleptonCut( config.getUntrackedParameter<double>("jetDRleptonCut", 0.3) ),
m_jetID( config.getUntrackedParameter<std::string>("jetID", "loose") ),
m_jetCSVv2Name( config.getUntrackedParameter<std::string>("jetCSVv2Name", "pfCombinedInclusiveSecondaryVertexV2BJetTags") ),
m_jetCSVv2L( config.getUntrackedParameter<double>("jetCSVv2L", 0.605) ),
m_jetCSVv2M( config.getUntrackedParameter<double>("jetCSVv2M", 0.89) ),
m_jetCSVv2T( config.getUntrackedParameter<double>("jetCSVv2T", 0.97) ),

m_hltDRCut( config.getUntrackedParameter<double>("hltDRCut", std::numeric_limits<float>::max()) ),
m_hltDPtCut( config.getUntrackedParameter<double>("hltDPtCut", std::numeric_limits<float>::max()) )
m_hltDRCut( config.getUntrackedParameter<double>("hltDRCut", std::numeric_limits<double>::max()) ),
m_hltDPtCut( config.getUntrackedParameter<double>("hltDPtCut", std::numeric_limits<double>::max()) )
{
if(config.exists("hltScaleFactors")){
const edm::ParameterSet& hltSFPSet = config.getUntrackedParameter<edm::ParameterSet>("hltScaleFactors");
std::vector<std::string> hltSFNames = hltSFPSet.getParameterNames();
for(const std::string& sf: hltSFNames){
BinnedValuesJSONParser parser(hltSFPSet.getUntrackedParameter<edm::FileInPath>(sf).fullPath());
m_hltSF.emplace(sf, std::move(parser.get_values()));
}
}
}

virtual void analyze(const edm::Event&, const edm::EventSetup&, const ProducersManager&, const AnalyzersManager&, const CategoryManager&) override;
Expand All @@ -64,33 +74,24 @@ class TTAnalyzer: public Framework::Analyzer {
BRANCH(diLeptons_IDIso, std::vector<std::vector<uint16_t>>);

BRANCH(selJets, std::vector<TTAnalysis::Jet>);
BRANCH(selJets_selID, std::vector<uint16_t>);
// ex.: selectedJets_..._DRCut[X][0] is the highest Pt selected jet with minDRjl>0.3 taking into account ID/Iso-X Leptons
BRANCH(selJets_selID_DRCut, std::vector<std::vector<uint16_t>>);
// ex.: selectedBJets_..._PtOrdered[X][0] is the highest Pt selected jet with minDRjl>0.3 taking into account ID/Iso/Btag-X combination
BRANCH(selBJets_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>);
BRANCH(selBJets_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);
// ex.: selJets_..._PtOrdered[X][0] is the highest Pt selected jet with minDRjl>0.3 taking into account leptonID/Iso pair & Btag-X combination
BRANCH(selJets_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>);
BRANCH(selJets_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);

BRANCH(diJets, std::vector<TTAnalysis::DiJet>);
// ex.: diJets_DRCut[X][0] is first diJet with minDRjl>0.3 taking into account ID/Iso-X Leptons
BRANCH(diJets_DRCut, std::vector<std::vector<uint16_t>>);
// ex.: diBJets_..._CSVv2Ordered[X][0] is the b-jet pair with highest CSVv2 values and with minDRjl>0.3 taking into account the leptonID/Iso/Btag-X combination
BRANCH(diBJets_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>);
BRANCH(diBJets_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);
// ex.: diJets_..._CSVv2Ordered[X][0] is the jet pair with highest CSVv2 values and with minDRjl>0.3 taking into account the leptonID/Iso pair & Btag-X combination
BRANCH(diJets_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>);
BRANCH(diJets_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);

// For all the following: indices are combinations of LeptonID/LeptonIso/(B-tagging working point)

BRANCH(diLepDiJets, std::vector<TTAnalysis::DiLepDiJet>);
BRANCH(diLepDiJets_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>); // di-leptons of combined ID/Iso with di-jets built out of jets having minDRjl>cut taking into account lepton ID/Iso corresponding to the chosen combination of the two leptons of the object
BRANCH(diLepDiJets_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);

BRANCH(diLepDiJets_DRCut, std::vector<std::vector<uint16_t>>); // di-leptons of combined ID/Iso with di-jets built out of jets having minDRjl>cut taking into account lepton ID/Iso corresponding to the loosest combination of the two leptons of the object
BRANCH(diLepDiBJets_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>);
BRANCH(diLepDiBJets_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);

BRANCH(diLepDiJetsMet, std::vector<TTAnalysis::DiLepDiJetMet>);

BRANCH(diLepDiJetsMet_DRCut, std::vector<std::vector<uint16_t>>);
BRANCH(diLepDiBJetsMet_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>);
BRANCH(diLepDiBJetsMet_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);
BRANCH(diLepDiJetsMet_DRCut_BWP_PtOrdered, std::vector<std::vector<uint16_t>>); // di-leptons of combined ID/Iso with di-jets built out of jets having minDRjl>cut taking into account lepton ID/Iso corresponding to the chosen combination of the two leptons of the object
BRANCH(diLepDiJetsMet_DRCut_BWP_CSVv2Ordered, std::vector<std::vector<uint16_t>>);

BRANCH(ttbar, std::vector<std::vector<std::vector<TTAnalysis::TTBar>>>);

Expand Down Expand Up @@ -140,28 +141,32 @@ class TTAnalyzer: public Framework::Analyzer {
BRANCH(gen_b_lepton_t_deltaR, float); // DeltaR between the b quark and the lepton coming from the top decay chain
BRANCH(gen_bbar_lepton_tbar_deltaR, float); // DeltaR between the b quark and the lepton coming from the top decay chain

// These two vectors are indexed wrt LepLepId enum
BRANCH(gen_b_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen b coming from the top decay and each selected jets. Indexed as `selectedJets_tightID_DRcut` array
BRANCH(gen_bbar_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen bbar coming from the anti-top decay chain and each selected jets. Indexed as `selectedJets_tightID_DRcut` array

// These two vectors are indexed wrt LepLepId enum
BRANCH(gen_b_beforeFSR_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen b coming from the top decay and each selected jets. Indexed as `selectedJets_tightID_DRcut` array
BRANCH(gen_bbar_beforeFSR_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen bbar coming from the anti-top decay chain and each selected jets. Indexed as `selectedJets_tightID_DRcut` array

BRANCH(gen_lepton_t_deltaR, std::vector<float>); // DeltaR between the gen lepton coming from the top decay chain and each selected lepton. Indexed as `leptons` array
BRANCH(gen_lepton_tbar_deltaR, std::vector<float>); // DeltaR between the gen lepton coming from the anti-top decay chain and each selected lepton. Indexed as `leptons` array

// These two vectors are indexed wrt LepLepId enum
BRANCH(gen_matched_b, std::vector<int8_t>); // Index inside the `selectedJets_tightID_DRcut` collection of the jet with the smallest deltaR with the gen b coming from the top decay
BRANCH(gen_matched_bbar, std::vector<int8_t>); // Index inside the `selectedJets_tightID_DRcut` collection of the jet with the smallest deltaR with the gen bbar coming from the anti-top decay

// These two vectors are indexed wrt LepLepId enum
BRANCH(gen_matched_b_beforeFSR, std::vector<int8_t>); // Index inside the `selectedJets_tightID_DRcut` collection of the jet with the smallest deltaR with the gen b coming from the top decay
BRANCH(gen_matched_bbar_beforeFSR, std::vector<int8_t>); // Index inside the `selectedJets_tightID_DRcut` collection of the jet with the smallest deltaR with the gen bbar coming from the anti-top decay
BRANCH(gen_lepton_beforeFSR_t_deltaR, std::vector<float>); // DeltaR between the gen lepton coming from the top decay chain and each selected lepton. Indexed as `leptons` array
BRANCH(gen_lepton_beforeFSR_tbar_deltaR, std::vector<float>); // DeltaR between the gen lepton coming from the anti-top decay chain and each selected lepton. Indexed as `leptons` array

BRANCH(gen_matched_lepton_t, int16_t); // Index inside the `leptons` collection of the lepton with the smallest deltaR with the gen lepton coming from the top decay chain
BRANCH(gen_matched_lepton_tbar, int16_t); // Index inside the `leptons` collection of the lepton with the smallest deltaR with the gen lepton coming from the anti-top decay chain

BRANCH(gen_matched_lepton_beforeFSR_t, int16_t); // Index inside the `leptons` collection of the lepton with the smallest deltaR with the gen lepton coming from the top decay chain
BRANCH(gen_matched_lepton_beforeFSR_tbar, int16_t); // Index inside the `leptons` collection of the lepton with the smallest deltaR with the gen lepton coming from the anti-top decay chain

// The following vectors are indexed wrt. LepLepIdJetJetBWP enum (as, for instance, selJets_DRCut_BWP_PtOrdered), giving indices pointing to the `selJets` collection
BRANCH(gen_b_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen b coming from the top decay and each selected jets.
BRANCH(gen_bbar_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen bbar coming from the anti-top decay chain and each selected jets.

BRANCH(gen_b_beforeFSR_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen b coming from the top decay and each selected jets.
BRANCH(gen_bbar_beforeFSR_deltaR, std::vector<std::vector<float>>); // DeltaR between the gen bbar coming from the anti-top decay chain and each selected jets.

// The following vectors are indexed wrt. LepLepIdJetJetBWP enum (as, for instance, selJets_DRCut_BWP_PtOrdered)
BRANCH(gen_matched_b, std::vector<int8_t>); // Index inside the `selJets` collection of the jet with the smallest deltaR with the gen b coming from the top decay
BRANCH(gen_matched_bbar, std::vector<int8_t>); // Index inside the `selJets` collection of the jet with the smallest deltaR with the gen bbar coming from the anti-top decay

BRANCH(gen_matched_b_beforeFSR, std::vector<int8_t>); // Index inside the `selJets` collection of the jet with the smallest deltaR with the gen b coming from the top decay
BRANCH(gen_matched_bbar_beforeFSR, std::vector<int8_t>); // Index inside the `selJets` collection of the jet with the smallest deltaR with the gen bbar coming from the anti-top decay

private:

// Producers name
Expand All @@ -171,6 +176,7 @@ class TTAnalyzer: public Framework::Analyzer {
const std::string m_met_producer;

const float m_electronPtCut, m_electronEtaCut;
const bool m_electronRemoveGap;
const std::string m_electronVetoIDName;
const std::string m_electronLooseIDName;
const std::string m_electronMediumIDName;
Expand All @@ -183,6 +189,7 @@ class TTAnalyzer: public Framework::Analyzer {
const float m_jetCSVv2L, m_jetCSVv2M, m_jetCSVv2T;

const float m_hltDRCut, m_hltDPtCut;
std::map<std::string, BinnedValues> m_hltSF;

std::shared_ptr<NeutrinosSolver> m_neutrinos_solver;

Expand Down
Loading