CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
TtSemiLepJetCombWMassDeltaTopMass Class Reference

#include <TtSemiLepJetCombWMassDeltaTopMass.h>

Inheritance diagram for TtSemiLepJetCombWMassDeltaTopMass:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 TtSemiLepJetCombWMassDeltaTopMass (const edm::ParameterSet &)
 
 ~TtSemiLepJetCombWMassDeltaTopMass () override
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
SerialTaskQueueglobalLuminosityBlocksQueue ()
 
SerialTaskQueueglobalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducer () override
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
 ~ProducerBase () noexcept(false) override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void beginJob () override
 
void endJob () override
 
bool isValid (const int &idx, const edm::Handle< std::vector< pat::Jet > > &jets)
 
void produce (edm::Event &evt, const edm::EventSetup &setup) override
 

Private Attributes

std::string bTagAlgorithm_
 
edm::EDGetTokenT< std::vector< pat::Jet > > jetsToken_
 
edm::EDGetTokenT< edm::View< reco::RecoCandidate > > lepsToken_
 
double maxBDiscLightJets_
 
int maxNJets_
 
edm::EDGetTokenT< std::vector< pat::MET > > metsToken_
 
double minBDiscBJets_
 
int neutrinoSolutionType_
 
bool useBTagging_
 
double wMass_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::ProducerBase
ProducesCollector producesCollector ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 13 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Constructor & Destructor Documentation

TtSemiLepJetCombWMassDeltaTopMass::TtSemiLepJetCombWMassDeltaTopMass ( const edm::ParameterSet cfg)
explicit

Definition at line 7 of file TtSemiLepJetCombWMassDeltaTopMass.cc.

References Exception, and maxNJets_.

8  : jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
10  metsToken_(consumes<std::vector<pat::MET> >(cfg.getParameter<edm::InputTag>("mets"))),
11  maxNJets_(cfg.getParameter<int>("maxNJets")),
12  wMass_(cfg.getParameter<double>("wMass")),
13  useBTagging_(cfg.getParameter<bool>("useBTagging")),
14  bTagAlgorithm_(cfg.getParameter<std::string>("bTagAlgorithm")),
15  minBDiscBJets_(cfg.getParameter<double>("minBDiscBJets")),
16  maxBDiscLightJets_(cfg.getParameter<double>("maxBDiscLightJets")),
17  neutrinoSolutionType_(cfg.getParameter<int>("neutrinoSolutionType")) {
18  if (maxNJets_ < 4 && maxNJets_ != -1)
19  throw cms::Exception("WrongConfig") << "Parameter maxNJets can not be set to " << maxNJets_ << ". \n"
20  << "It has to be larger than 4 or can be set to -1 to take all jets.";
21 
22  produces<std::vector<std::vector<int> > >();
23  produces<int>("NumberOfConsideredJets");
24 }
edm::EDGetTokenT< std::vector< pat::Jet > > jetsToken_
T getParameter(std::string const &) const
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
edm::EDGetTokenT< edm::View< reco::RecoCandidate > > lepsToken_
edm::EDGetTokenT< std::vector< pat::MET > > metsToken_
TtSemiLepJetCombWMassDeltaTopMass::~TtSemiLepJetCombWMassDeltaTopMass ( )
override

Definition at line 26 of file TtSemiLepJetCombWMassDeltaTopMass.cc.

26 {}

Member Function Documentation

void TtSemiLepJetCombWMassDeltaTopMass::beginJob ( void  )
inlineoverrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 19 of file TtSemiLepJetCombWMassDeltaTopMass.h.

References produce(), and singleTopDQM_cfi::setup.

19 {};
void TtSemiLepJetCombWMassDeltaTopMass::endJob ( void  )
inlineoverrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 21 of file TtSemiLepJetCombWMassDeltaTopMass.h.

21 {};
bool TtSemiLepJetCombWMassDeltaTopMass::isValid ( const int &  idx,
const edm::Handle< std::vector< pat::Jet > > &  jets 
)
inlineprivate

Definition at line 23 of file TtSemiLepJetCombWMassDeltaTopMass.h.

References singleTopDQM_cfi::jets.

Referenced by ntupleDataFormat._Object::_checkIsValid(), produce(), and core.AutoHandle.AutoHandle::ReallyLoad().

23  {
24  return (0 <= idx && idx < (int)jets->size());
25  };
void TtSemiLepJetCombWMassDeltaTopMass::produce ( edm::Event evt,
const edm::EventSetup setup 
)
overrideprivate

Definition at line 28 of file TtSemiLepJetCombWMassDeltaTopMass.cc.

References bTagAlgorithm_, MEzCalculator::Calculate(), Exception, edm::Event::getByToken(), TtSemiLepEvtPartons::HadB, mps_fire::i, training_settings::idx, isValid(), singleTopDQM_cfi::jets, jetsToken_, TtSemiLepEvtPartons::LepB, TtSemiLepHitFitProducer_Electrons_cfi::leps, lepsToken_, TtSemiLepEvtPartons::LightQ, TtSemiLepEvtPartons::LightQBar, match(), maxBDiscLightJets_, SUSYDQMAnalyzer_cfi::maxNJets, maxNJets_, singleTopDQM_cfi::mets, metsToken_, minBDiscBJets_, eostools::move(), neutrinoSolutionType_, edm::Event::put(), MEzCalculator::SetLepton(), MEzCalculator::SetMET(), mathSSE::sqrt(), useBTagging_, and wMass_.

Referenced by beginJob().

28  {
29  std::unique_ptr<std::vector<std::vector<int> > > pOut(new std::vector<std::vector<int> >);
30  std::unique_ptr<int> pJetsConsidered(new int);
31 
32  std::vector<int> match;
33  for (unsigned int i = 0; i < 4; ++i)
34  match.push_back(-1);
35 
36  // get jets
38  evt.getByToken(jetsToken_, jets);
39 
40  // get leptons
42  evt.getByToken(lepsToken_, leps);
43 
44  // get MET
46  evt.getByToken(metsToken_, mets);
47 
48  // skip events without lepton candidate or less than 4 jets or no MET
49  if (leps->empty() || jets->size() < 4 || mets->empty()) {
50  pOut->push_back(match);
51  evt.put(std::move(pOut));
52  *pJetsConsidered = jets->size();
53  evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
54  return;
55  }
56 
57  unsigned maxNJets = maxNJets_;
58  if (maxNJets_ == -1 || (int)jets->size() < maxNJets_)
59  maxNJets = jets->size();
60  *pJetsConsidered = maxNJets;
61  evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
62 
63  std::vector<bool> isBJet;
64  std::vector<bool> isLJet;
65  int cntBJets = 0;
66  if (useBTagging_) {
67  for (unsigned int idx = 0; idx < maxNJets; ++idx) {
68  isBJet.push_back(((*jets)[idx].bDiscriminator(bTagAlgorithm_) > minBDiscBJets_));
69  isLJet.push_back(((*jets)[idx].bDiscriminator(bTagAlgorithm_) < maxBDiscLightJets_));
70  if ((*jets)[idx].bDiscriminator(bTagAlgorithm_) > minBDiscBJets_)
71  cntBJets++;
72  }
73  }
74 
75  // -----------------------------------------------------
76  // associate those jets that get closest to the W mass
77  // with their invariant mass to the hadronic W boson
78  // -----------------------------------------------------
79  double wDist = -1.;
80  std::vector<int> closestToWMassIndices;
81  closestToWMassIndices.push_back(-1);
82  closestToWMassIndices.push_back(-1);
83  for (unsigned idx = 0; idx < maxNJets; ++idx) {
84  if (useBTagging_ && (!isLJet[idx] || (cntBJets <= 2 && isBJet[idx])))
85  continue;
86  for (unsigned jdx = (idx + 1); jdx < maxNJets; ++jdx) {
87  if (useBTagging_ &&
88  (!isLJet[jdx] || (cntBJets <= 2 && isBJet[jdx]) || (cntBJets == 3 && isBJet[idx] && isBJet[jdx])))
89  continue;
90  reco::Particle::LorentzVector sum = (*jets)[idx].p4() + (*jets)[jdx].p4();
91  if (wDist < 0. || wDist > fabs(sum.mass() - wMass_)) {
92  wDist = fabs(sum.mass() - wMass_);
93  closestToWMassIndices.clear();
94  closestToWMassIndices.push_back(idx);
95  closestToWMassIndices.push_back(jdx);
96  }
97  }
98  }
99 
100  // -----------------------------------------------------
101  // build a leptonic W boson from the lepton and the MET
102  // -----------------------------------------------------
103  double neutrino_pz = 0.;
104  if (neutrinoSolutionType_ != -1) {
105  MEzCalculator mez;
106  mez.SetMET(*mets->begin());
107  if (dynamic_cast<const reco::Muon*>(&(leps->front())))
108  mez.SetLepton((*leps)[0], true);
109  else if (dynamic_cast<const reco::GsfElectron*>(&(leps->front())))
110  mez.SetLepton((*leps)[0], false);
111  else
112  throw cms::Exception("UnimplementedFeature")
113  << "Type of lepton given together with MET for solving neutrino kinematics is neither muon nor electron.\n";
114  neutrino_pz = mez.Calculate(neutrinoSolutionType_);
115  }
116  const math::XYZTLorentzVector neutrino(mets->begin()->px(),
117  mets->begin()->py(),
118  neutrino_pz,
119  sqrt(mets->begin()->px() * mets->begin()->px() +
120  mets->begin()->py() * mets->begin()->py() + neutrino_pz * neutrino_pz));
121  const reco::Particle::LorentzVector lepW = neutrino + leps->front().p4();
122 
123  // -----------------------------------------------------
124  // associate those jets to the hadronic and the leptonic
125  // b quark that minimize the difference between
126  // hadronic and leptonic top mass
127  // -----------------------------------------------------
128  double deltaTop = -1.;
129  int hadB = -1;
130  int lepB = -1;
131  if (isValid(closestToWMassIndices[0], jets) && isValid(closestToWMassIndices[1], jets)) {
132  const reco::Particle::LorentzVector hadW =
133  (*jets)[closestToWMassIndices[0]].p4() + (*jets)[closestToWMassIndices[1]].p4();
134  // find hadronic b candidate
135  for (unsigned idx = 0; idx < maxNJets; ++idx) {
136  if (useBTagging_ && !isBJet[idx])
137  continue;
138  // make sure it's not used up already from the hadronic W
139  if ((int)idx != closestToWMassIndices[0] && (int)idx != closestToWMassIndices[1]) {
140  reco::Particle::LorentzVector hadTop = hadW + (*jets)[idx].p4();
141  // find leptonic b candidate
142  for (unsigned jdx = 0; jdx < maxNJets; ++jdx) {
143  if (useBTagging_ && !isBJet[jdx])
144  continue;
145  // make sure it's not used up already from the hadronic branch
146  if ((int)jdx != closestToWMassIndices[0] && (int)jdx != closestToWMassIndices[1] && jdx != idx) {
147  reco::Particle::LorentzVector lepTop = lepW + (*jets)[jdx].p4();
148  if (deltaTop < 0. || deltaTop > fabs(hadTop.mass() - lepTop.mass())) {
149  deltaTop = fabs(hadTop.mass() - lepTop.mass());
150  hadB = idx;
151  lepB = jdx;
152  }
153  }
154  }
155  }
156  }
157  }
158 
159  match[TtSemiLepEvtPartons::LightQ] = closestToWMassIndices[0];
160  match[TtSemiLepEvtPartons::LightQBar] = closestToWMassIndices[1];
161  match[TtSemiLepEvtPartons::HadB] = hadB;
162  match[TtSemiLepEvtPartons::LepB] = lepB;
163 
164  pOut->push_back(match);
165  evt.put(std::move(pOut));
166 }
edm::EDGetTokenT< std::vector< pat::Jet > > jetsToken_
void SetLepton(const pat::Particle &lepton, bool isMuon=true)
Set lepton.
Definition: MEzCalculator.h:31
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:131
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
void SetMET(const pat::MET &MET)
Set MET.
Definition: MEzCalculator.h:25
double Calculate(int type=1)
member functions
T sqrt(T t)
Definition: SSEVec.h:19
bool isValid(const int &idx, const edm::Handle< std::vector< pat::Jet > > &jets)
edm::EDGetTokenT< edm::View< reco::RecoCandidate > > lepsToken_
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10
edm::EDGetTokenT< std::vector< pat::MET > > metsToken_
math::XYZTLorentzVector LorentzVector
Lorentz vector.
Definition: Particle.h:21
def move(src, dest)
Definition: eostools.py:511

Member Data Documentation

std::string TtSemiLepJetCombWMassDeltaTopMass::bTagAlgorithm_
private

Definition at line 33 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

edm::EDGetTokenT<std::vector<pat::Jet> > TtSemiLepJetCombWMassDeltaTopMass::jetsToken_
private

Definition at line 25 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

edm::EDGetTokenT<edm::View<reco::RecoCandidate> > TtSemiLepJetCombWMassDeltaTopMass::lepsToken_
private

Definition at line 28 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

double TtSemiLepJetCombWMassDeltaTopMass::maxBDiscLightJets_
private

Definition at line 35 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

int TtSemiLepJetCombWMassDeltaTopMass::maxNJets_
private
edm::EDGetTokenT<std::vector<pat::MET> > TtSemiLepJetCombWMassDeltaTopMass::metsToken_
private

Definition at line 29 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

double TtSemiLepJetCombWMassDeltaTopMass::minBDiscBJets_
private

Definition at line 34 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

int TtSemiLepJetCombWMassDeltaTopMass::neutrinoSolutionType_
private

Definition at line 36 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

bool TtSemiLepJetCombWMassDeltaTopMass::useBTagging_
private

Definition at line 32 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().

double TtSemiLepJetCombWMassDeltaTopMass::wMass_
private

Definition at line 31 of file TtSemiLepJetCombWMassDeltaTopMass.h.

Referenced by produce().