TtSemiLepSignalSelMVAComputer Class Reference

#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 ( 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 ( )

override

Definition at line 21 of file TtSemiLepSignalSelMVAComputer.cc.

{}

Member Function Documentation

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

private

Definition at line 110 of file TtSemiLepSignalSelMVAComputer.cc.

References 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 117 of file TtSemiLepSignalSelMVAComputer.cc.

References DeltaPhi().

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::produce ( edm::Event &  evt, const edm::EventSetup &  setup  )

overrideprivate

Definition at line 23 of file TtSemiLepSignalSelMVAComputer.cc.

References edm::View< T >::begin(), DeltaR(), pat::Flags::Overlap::Electrons, HI_PhotonSkim_cff::electrons, electronsToken_, edm::Event::emplace(), edm::View< T >::end(), evaluateTtSemiLepSignalSel(), edm::Event::getByToken(), edm::EventSetup::getData(), edm::Ref< C, T, F >::isNull(), edm::HandleBase::isValid(), metsig::jet, fwrapper::jets, jetsToken_, METsToken_, HLT_FULL_cff::muon, patZpeak::muons, muonsToken_, mvaComputer, mvaToken_, putToken_, corrVsCorr::selection, 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);
}

Member Data Documentation

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

private

Definition at line 36 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().

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

private

Definition at line 34 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().

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

private

Definition at line 35 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().

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

private

Definition at line 33 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().

PhysicsTools::MVAComputerCache TtSemiLepSignalSelMVAComputer::mvaComputer

private

Definition at line 39 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().

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

private

Definition at line 32 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().

edm::EDPutTokenT<double> TtSemiLepSignalSelMVAComputer::putToken_

private

Definition at line 37 of file TtSemiLepSignalSelMVAComputer.h.

Referenced by produce().