#include <TtSemiLepSignalSelMVAComputer.h>

Inheritance diagram for TtSemiLepSignalSelMVAComputer:

Public Member Functions
	TtSemiLepSignalSelMVAComputer (const edm::ParameterSet &)

	~TtSemiLepSignalSelMVAComputer () override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

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< ConsumesInfo >	consumesInfo () 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

EDConsumerBase &	operator= (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

double	DeltaPhi (const math::XYZTLorentzVector &v1, const math::XYZTLorentzVector &v2)

double	DeltaR (const math::XYZTLorentzVector &v1, const math::XYZTLorentzVector &v2)

void	endJob () override

void	produce (edm::Event &evt, const edm::EventSetup &setup) override

Private Attributes
double	DiscSel

edm::EDGetTokenT< edm::View< pat::Electron > >	electronsToken_

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

edm::EDGetTokenT< edm::View< pat::MET > >	METsToken_

edm::EDGetTokenT< edm::View< pat::Muon > >	muonsToken_

PhysicsTools::MVAComputerCache	mvaComputer

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 21 of file TtSemiLepSignalSelMVAComputer.h.

Constructor & Destructor Documentation

TtSemiLepSignalSelMVAComputer::TtSemiLepSignalSelMVAComputer ( const edm::ParameterSet & cfg )

explicit

Definition at line 13 of file TtSemiLepSignalSelMVAComputer.cc.

     : muonsToken_(consumes<edm::View<pat::Muon> >(cfg.getParameter<edm::InputTag>("muons"))),
       jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
       METsToken_(consumes<edm::View<pat::MET> >(cfg.getParameter<edm::InputTag>("mets"))),
       electronsToken_(consumes<edm::View<pat::Electron> >(cfg.getParameter<edm::InputTag>("elecs"))) {
   produces<double>("DiscSel");
 }

TtSemiLepSignalSelMVAComputer::~TtSemiLepSignalSelMVAComputer ( )

override

Definition at line 21 of file TtSemiLepSignalSelMVAComputer.cc.

21 {}

Member Function Documentation

void TtSemiLepSignalSelMVAComputer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 123 of file TtSemiLepSignalSelMVAComputer.cc.

123 {}

double TtSemiLepSignalSelMVAComputer::DeltaPhi	(	const math::XYZTLorentzVector &	v1,
		const math::XYZTLorentzVector &	v2
	)

private

Definition at line 127 of file TtSemiLepSignalSelMVAComputer.cc.

References HLT_2018_cff::dPhi, and Pi.

Referenced by DeltaR().

                                                                                                                {
   double dPhi = fabs(v1.Phi() - v2.Phi());
   if (dPhi > TMath::Pi())
     dPhi = 2 * TMath::Pi() - dPhi;
   return dPhi;
 }

double TtSemiLepSignalSelMVAComputer::DeltaR	(	const math::XYZTLorentzVector &	v1,
		const math::XYZTLorentzVector &	v2
	)

private

Definition at line 134 of file TtSemiLepSignalSelMVAComputer.cc.

References DeltaPhi(), HLT_2018_cff::dPhi, and HGC3DClusterGenMatchSelector_cfi::dR.

Referenced by produce().

                                                                                                              {
   double dPhi = DeltaPhi(v1, v2);
   double dR = TMath::Sqrt((v1.Eta() - v2.Eta()) * (v1.Eta() - v2.Eta()) + dPhi * dPhi);
   return dR;
 }

void TtSemiLepSignalSelMVAComputer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 125 of file TtSemiLepSignalSelMVAComputer.cc.

125 {}

void TtSemiLepSignalSelMVAComputer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	setup
	)

overrideprivate

Definition at line 23 of file TtSemiLepSignalSelMVAComputer.cc.

References edm::View< T >::begin(), DeltaR(), DiscSel, pwdgSkimBPark_cfi::electrons, pat::Flags::Overlap::Electrons, electronsToken_, edm::View< T >::end(), evaluateTtSemiLepSignalSel(), edm::eventsetup::EventSetupRecordImplementation< T >::get(), edm::Event::getByToken(), edm::Ref< C, T, F >::isNull(), edm::HandleBase::isValid(), metsig::jet, singleTopDQM_cfi::jets, jetsToken_, METsToken_, eostools::move(), PDWG_BPHSkim_cff::muons, muonsToken_, mvaComputer, B2GMonitoring_cff::nelectrons, B2GMonitoring_cff::nmuons, edm::Event::put(), corrVsCorr::selection, singleTopDQM_cfi::setup, pat::Flags::test(), and PhysicsTools::MVAComputerCache::update().

                                                                                      {
   std::unique_ptr<double> pOutDisc(new double);
 
   mvaComputer.update<TtSemiLepSignalSelMVARcd>(setup, "ttSemiLepSignalSelMVA");
 
   // read name of the last processor in the MVA calibration
   // (to be used as meta information)
   edm::ESHandle<PhysicsTools::Calibration::MVAComputerContainer> calibContainer;
   setup.get<TtSemiLepSignalSelMVARcd>().get(calibContainer);
   std::vector<PhysicsTools::Calibration::VarProcessor*> processors =
       (calibContainer->find("ttSemiLepSignalSelMVA")).getProcessors();
 
   //make your preselection! This must!! be the same one as in TraintreeSaver.cc
   edm::Handle<edm::View<pat::MET> > MET_handle;
   evt.getByToken(METsToken_, MET_handle);
   if (!MET_handle.isValid())
     return;
   const edm::View<pat::MET> MET = *MET_handle;
 
   edm::Handle<std::vector<pat::Jet> > jet_handle;
   evt.getByToken(jetsToken_, jet_handle);
   if (!jet_handle.isValid())
     return;
   const std::vector<pat::Jet> jets = *jet_handle;
   unsigned int nJets = 0;
   std::vector<pat::Jet> seljets;
   for (std::vector<pat::Jet>::const_iterator it = jets.begin(); it != jets.end(); it++) {
     if (!(pat::Flags::test(*it, pat::Flags::Overlap::Electrons)))
       continue;
     if (it->pt() > 30. && fabs(it->eta()) < 2.4) {
       seljets.push_back(*it);
       nJets++;
     }
   }
 
   edm::Handle<edm::View<pat::Muon> > muon_handle;
   evt.getByToken(muonsToken_, muon_handle);
   if (!muon_handle.isValid())
     return;
   const edm::View<pat::Muon> muons = *muon_handle;
   int nmuons = 0;
   std::vector<pat::Muon> selMuons;
   for (edm::View<pat::Muon>::const_iterator it = muons.begin(); it != muons.end(); it++) {
     reco::TrackRef gltr = it->track();       // global track
     reco::TrackRef trtr = it->innerTrack();  // tracker track
     if (it->pt() > 30 && fabs(it->eta()) < 2.1 && (it->pt() / (it->pt() + it->trackIso() + it->caloIso())) > 0.95 &&
         it->isGlobalMuon()) {
       if (gltr.isNull())
         continue;  //temporary problems with dead trackrefs
       if ((gltr->chi2() / gltr->ndof()) < 10 && trtr->numberOfValidHits() > 11) {
         double dRmin = 9999.;
         for (std::vector<pat::Jet>::const_iterator ajet = seljets.begin(); ajet != seljets.end(); ajet++) {
           math::XYZTLorentzVector jet = ajet->p4();
           math::XYZTLorentzVector muon = it->p4();
           double tmpdR = DeltaR(muon, jet);
           if (tmpdR < dRmin)
             dRmin = tmpdR;
         }
         if (dRmin > 0.3) {  //temporary problems with muon isolation
           nmuons++;
           selMuons.push_back(*it);
         }
       }
     }
   }
 
   edm::Handle<edm::View<pat::Electron> > electron_handle;
   evt.getByToken(electronsToken_, electron_handle);
   if (!electron_handle.isValid())
     return;
   const edm::View<pat::Electron> electrons = *electron_handle;
   int nelectrons = 0;
   for (edm::View<pat::Electron>::const_iterator it = electrons.begin(); it != electrons.end(); it++) {
     if (it->pt() > 30 && fabs(it->eta()) < 2.4 && (it->pt() / (it->pt() + it->trackIso() + it->caloIso())) > 0.95 &&
         it->isElectronIDAvailable("eidTight")) {
       if (it->electronID("eidTight") == 1)
         nelectrons++;
     }
   }
 
   double discrim;
   // discriminator output for events which do not pass the preselection is set to -1
   if (nmuons != 1 || nJets < 4 || nelectrons > 0)
     discrim = -1.;  //std::cout<<"nJets: "<<seljets.size()<<" numLeptons: "<<nleptons<<std::endl;}
   else {
     //check wheter a event was already selected (problem with duplicated events)
     math::XYZTLorentzVector muon = selMuons.begin()->p4();
 
     TtSemiLepSignalSel selection(seljets, muon, MET);
 
     discrim = evaluateTtSemiLepSignalSel(mvaComputer, selection);
   }
 
   *pOutDisc = discrim;
 
   evt.put(std::move(pOutDisc), "DiscSel");
 
   DiscSel = discrim;
 }