#include <TtSemiLepSignalSelMVAComputer.h>

Inheritance diagram for TtSemiLepSignalSelMVAComputer:

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

	~TtSemiLepSignalSelMVAComputer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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

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

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

Private Attributes
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

edm::ESGetToken< PhysicsTools::Calibration::MVAComputerContainer, TtSemiLepSignalSelMVARcd >	mvaToken_

edm::EDPutTokenT< double >	putToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 21 of file TtSemiLepSignalSelMVAComputer.h.

Constructor & Destructor Documentation

◆ TtSemiLepSignalSelMVAComputer()

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

explicit

Definition at line 13 of file TtSemiLepSignalSelMVAComputer.cc.

     : mvaToken_(esConsumes()),
       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"))),
       putToken_(produces("DiscSel")) {}

◆ ~TtSemiLepSignalSelMVAComputer()

TtSemiLepSignalSelMVAComputer::~TtSemiLepSignalSelMVAComputer ( )

override

Definition at line 21 of file TtSemiLepSignalSelMVAComputer.cc.

21 {}

Member Function Documentation

◆ DeltaPhi()

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

private

Definition at line 110 of file TtSemiLepSignalSelMVAComputer.cc.

References HLT_2022v12_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;
 }

◆ DeltaR()

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

private

Definition at line 117 of file TtSemiLepSignalSelMVAComputer.cc.

References DeltaPhi(), HLT_2022v12_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;
 }

◆ produce()

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

overrideprivate

Definition at line 23 of file TtSemiLepSignalSelMVAComputer.cc.

References DeltaR(), HLT_2022v12_cff::dRmin, pwdgSkimBPark_cfi::electrons, pat::Flags::Overlap::Electrons, electronsToken_, edm::Event::emplace(), evaluateTtSemiLepSignalSel(), edm::Event::getByToken(), edm::Ref< C, T, F >::isNull(), edm::HandleBase::isValid(), metsig::jet, singleTopDQM_cfi::jets, jetsToken_, METsToken_, PDWG_BPHSkim_cff::muons, muonsToken_, mvaComputer, mvaToken_, B2GMonitoring_cff::nelectrons, B2GMonitoring_cff::nmuons, putToken_, corrVsCorr::selection, singleTopDQM_cfi::setup, pat::Flags::test(), and PhysicsTools::MVAComputerCache::update().

                                                                                      {
   mvaComputer.update(&setup.getData(mvaToken_), "ttSemiLepSignalSelMVA");
 
   //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);
   }
 
   evt.emplace(putToken_, discrim);
 }