muonisolation::JetExtractor Class Reference

#include <JetExtractor.h>

Inheritance diagram for muonisolation::JetExtractor:

Public Member Functions
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const override
void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks) override
	JetExtractor ()
	JetExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)
	~JetExtractor () override
Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
	make single IsoDeposit based on a CandidateBaseRef as input More...
virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const
virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
virtual void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...
virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Attributes
TrackDetectorAssociator *	theAssociator
TrackAssociatorParameters *	theAssociatorParameters
double	theDR_Max
double	theDR_Veto
bool	theExcludeMuonVeto
edm::EDGetTokenT< reco::CaloJetCollection >	theJetCollectionToken
bool	thePrintTimeReport
std::string	thePropagatorName
MuonServiceProxy *	theService
	the event setup proxy, it takes care the services update More...
double	theThreshold

Detailed Description

Definition at line 34 of file JetExtractor.h.

Constructor & Destructor Documentation

muonisolation::JetExtractor::JetExtractor ( )

inline

Definition at line 36 of file JetExtractor.h.

References deposit(), ev, fillVetos(), HLT_2018_cff::track, PDWG_EXOHSCP_cff::tracks, and ~JetExtractor().

{};

JetExtractor::JetExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &&	iC
	)

Definition at line 32 of file JetExtractor.cc.

References edm::ParameterSet::getParameter(), TrackAssociatorParameters::loadParameters(), MuonServiceProxy_cff::MuonServiceProxy, theAssociator, theAssociatorParameters, theService, and PDWG_EXOHSCP_cff::TrackAssociatorParameters.

    : theJetCollectionToken(iC.consumes<CaloJetCollection>(par.getParameter<edm::InputTag>("JetCollectionLabel"))),
      thePropagatorName(par.getParameter<std::string>("PropagatorName")),
      theThreshold(par.getParameter<double>("Threshold")),
      theDR_Veto(par.getParameter<double>("DR_Veto")),
      theDR_Max(par.getParameter<double>("DR_Max")),
      theExcludeMuonVeto(par.getParameter<bool>("ExcludeMuonVeto")),
      theService(nullptr),
      theAssociator(nullptr),
      thePrintTimeReport(par.getUntrackedParameter<bool>("PrintTimeReport")) {
  ParameterSet serviceParameters = par.getParameter<ParameterSet>("ServiceParameters");
  theService = new MuonServiceProxy(serviceParameters);

  //  theAssociatorParameters = new TrackAssociatorParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"), iC_);
  theAssociatorParameters = new TrackAssociatorParameters();
  theAssociatorParameters->loadParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"), iC);
  theAssociator = new TrackDetectorAssociator();
}

JetExtractor::~JetExtractor ( )

override

Definition at line 51 of file JetExtractor.cc.

References theAssociator, theAssociatorParameters, and theService.

Referenced by JetExtractor().

                            {
  if (theAssociatorParameters)
    delete theAssociatorParameters;
  if (theService)
    delete theService;
  if (theAssociator)
    delete theAssociator;
}

Member Function Documentation

IsoDeposit JetExtractor::deposit ( const edm::Event &  ev, const edm::EventSetup &  evSetup, const reco::Track &  track  ) const

overridevirtual

make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 66 of file JetExtractor.cc.

References TrackDetectorAssociator::associate(), Calorimetry_cff::bField, TrackDetMatchInfo::crossedTowerIds, TrackDetMatchInfo::crossedTowers, reco::deltaR(), reco::TrackBase::eta(), edm::EventSetup::get(), reco::TransientTrack::initialFreeState(), reco::TrackBase::phi(), TrackDetectorAssociator::setPropagator(), theAssociator, theAssociatorParameters, theDR_Max, theDR_Veto, theExcludeMuonVeto, theJetCollectionToken, thePropagatorName, theService, theThreshold, and TrackDetMatchInfo::trkGlobPosAtHcal.

Referenced by JetExtractor().

                                                                                                          {
  theService->update(eventSetup);
  theAssociator->setPropagator(&*(theService->propagator(thePropagatorName)));

  typedef IsoDeposit::Veto Veto;
  IsoDeposit::Direction muonDir(muon.eta(), muon.phi());

  IsoDeposit depJet(muonDir);

  edm::ESHandle<MagneticField> bField;
  eventSetup.get<IdealMagneticFieldRecord>().get(bField);

  reco::TransientTrack tMuon(muon, &*bField);
  FreeTrajectoryState iFTS = tMuon.initialFreeState();
  TrackDetMatchInfo mInfo = theAssociator->associate(event, eventSetup, iFTS, *theAssociatorParameters);

  reco::isodeposit::Direction vetoDirection(mInfo.trkGlobPosAtHcal.eta(), mInfo.trkGlobPosAtHcal.phi());
  depJet.setVeto(Veto(vetoDirection, theDR_Veto));

  edm::Handle<CaloJetCollection> caloJetsH;
  event.getByToken(theJetCollectionToken, caloJetsH);

  //use calo towers
  CaloJetCollection::const_iterator jetCI = caloJetsH->begin();
  for (; jetCI != caloJetsH->end(); ++jetCI) {
    double deltar0 = reco::deltaR(muon, *jetCI);
    if (deltar0 > theDR_Max)
      continue;
    if (jetCI->et() < theThreshold)
      continue;

    //should I make a separate config option for this?
    std::vector<CaloTowerPtr> jetConstituents = jetCI->getCaloConstituents();

    std::vector<DetId>::const_iterator crossedCI = mInfo.crossedTowerIds.begin();
    std::vector<CaloTowerPtr>::const_iterator jetTowCI = jetConstituents.begin();

    double sumEtExcluded = 0;
    for (; jetTowCI != jetConstituents.end(); ++jetTowCI) {
      bool isExcluded = false;
      double deltaRLoc = reco::deltaR(vetoDirection, *jetCI);
      if (deltaRLoc < theDR_Veto) {
        isExcluded = true;
      }
      for (; !isExcluded && crossedCI != mInfo.crossedTowerIds.end(); ++crossedCI) {
        if (crossedCI->rawId() == (*jetTowCI)->id().rawId()) {
          isExcluded = true;
        }
      }
      if (isExcluded)
        sumEtExcluded += (*jetTowCI)->et();
    }
    if (theExcludeMuonVeto) {
      if (jetCI->et() - sumEtExcluded < theThreshold)
        continue;
    }

    double depositEt = jetCI->et();
    if (theExcludeMuonVeto)
      depositEt = depositEt - sumEtExcluded;

    reco::isodeposit::Direction jetDir(jetCI->eta(), jetCI->phi());
    depJet.addDeposit(jetDir, depositEt);
  }

  std::vector<const CaloTower*>::const_iterator crossedCI = mInfo.crossedTowers.begin();
  double muSumEt = 0;
  for (; crossedCI != mInfo.crossedTowers.end(); ++crossedCI) {
    muSumEt += (*crossedCI)->et();
  }
  depJet.addCandEnergy(muSumEt);

  return depJet;
}

void JetExtractor::fillVetos ( const edm::Event &  ev, const edm::EventSetup &  evSetup, const reco::TrackCollection &  tracks  )

overridevirtual

fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 60 of file JetExtractor.cc.

Referenced by JetExtractor().

                                                                                                                 {
  //   LogWarning("JetExtractor")
  //     <<"fillVetos does nothing now: IsoDeposit provides enough functionality\n"
  //     <<"to remove a deposit at/around given (eta, phi)";
}

Member Data Documentation

TrackDetectorAssociator* muonisolation::JetExtractor::theAssociator

private

Definition at line 63 of file JetExtractor.h.

Referenced by deposit(), JetExtractor(), and ~JetExtractor().

TrackAssociatorParameters* muonisolation::JetExtractor::theAssociatorParameters

private

Definition at line 62 of file JetExtractor.h.

Referenced by deposit(), JetExtractor(), and ~JetExtractor().

double muonisolation::JetExtractor::theDR_Max

private

Definition at line 54 of file JetExtractor.h.

Referenced by deposit().

double muonisolation::JetExtractor::theDR_Veto

private

Definition at line 53 of file JetExtractor.h.

Referenced by deposit().

bool muonisolation::JetExtractor::theExcludeMuonVeto

private

Definition at line 57 of file JetExtractor.h.

Referenced by deposit().

edm::EDGetTokenT<reco::CaloJetCollection> muonisolation::JetExtractor::theJetCollectionToken

private

Definition at line 47 of file JetExtractor.h.

Referenced by deposit().

bool muonisolation::JetExtractor::thePrintTimeReport

private

Definition at line 65 of file JetExtractor.h.

std::string muonisolation::JetExtractor::thePropagatorName

private

Definition at line 49 of file JetExtractor.h.

Referenced by deposit().

MuonServiceProxy* muonisolation::JetExtractor::theService

private

the event setup proxy, it takes care the services update

Definition at line 60 of file JetExtractor.h.

Referenced by deposit(), JetExtractor(), and ~JetExtractor().

double muonisolation::JetExtractor::theThreshold

private

Definition at line 52 of file JetExtractor.h.

Referenced by deposit().