#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::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 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 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 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 void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Attributes
std::unique_ptr< TrackDetectorAssociator >	theAssociator

std::unique_ptr< TrackAssociatorParameters >	theAssociatorParameters

double	theDR_Max

double	theDR_Veto

bool	theExcludeMuonVeto

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	theFieldToken

edm::EDGetTokenT< reco::CaloJetCollection >	theJetCollectionToken

bool	thePrintTimeReport

std::string	thePropagatorName

std::unique_ptr< MuonServiceProxy >	theService
	the event setup proxy, it takes care the services update More...

double	theThreshold

Detailed Description

Definition at line 37 of file JetExtractor.h.

Constructor & Destructor Documentation

◆ JetExtractor() [1/2]

JetExtractor::JetExtractor ( )

Definition at line 27 of file JetExtractor.cc.

27 {}

◆ JetExtractor() [2/2]

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

Definition at line 29 of file JetExtractor.cc.

     : theJetCollectionToken(iC.consumes<CaloJetCollection>(par.getParameter<edm::InputTag>("JetCollectionLabel"))),
       thePropagatorName(par.getParameter<std::string>("PropagatorName")),
       theFieldToken(iC.esConsumes()),
       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 = std::make_unique<MuonServiceProxy>(serviceParameters, edm::ConsumesCollector(iC));
  
   //  theAssociatorParameters = new TrackAssociatorParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"), iC_);
   theAssociatorParameters = std::make_unique<TrackAssociatorParameters>();
   theAssociatorParameters->loadParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"), iC);
   theAssociator = std::make_unique<TrackDetectorAssociator>();
 }

References edm::ParameterSet::getParameter(), theAssociator, theAssociatorParameters, and theService.

◆ ~JetExtractor()

JetExtractor::~JetExtractor ( )

override

Definition at line 49 of file JetExtractor.cc.

49 {}

Member Function Documentation

◆ deposit()

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 57 of file JetExtractor.cc.

                                                                                                           {
   theService->update(eventSetup);
   theAssociator->setPropagator(&*(theService->propagator(thePropagatorName)));
  
   typedef IsoDeposit::Veto Veto;
   IsoDeposit::Direction muonDir(muon.eta(), muon.phi());
  
   IsoDeposit depJet(muonDir);
  
   auto const& bField = eventSetup.getData(theFieldToken);
  
   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;
 }

References Calorimetry_cff::bField, TrackDetMatchInfo::crossedTowerIds, TrackDetMatchInfo::crossedTowers, reco::deltaR(), edm::EventSetup::getData(), reco::TransientTrack::initialFreeState(), theAssociator, theAssociatorParameters, theDR_Max, theDR_Veto, theExcludeMuonVeto, theFieldToken, theJetCollectionToken, thePropagatorName, theService, theThreshold, and TrackDetMatchInfo::trkGlobPosAtHcal.

◆ fillVetos()

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 51 of file JetExtractor.cc.

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

Member Data Documentation

◆ theAssociator

std::unique_ptr<TrackDetectorAssociator> muonisolation::JetExtractor::theAssociator

private

Definition at line 68 of file JetExtractor.h.

Referenced by deposit(), and JetExtractor().

◆ theAssociatorParameters

std::unique_ptr<TrackAssociatorParameters> muonisolation::JetExtractor::theAssociatorParameters

private

Definition at line 67 of file JetExtractor.h.

Referenced by deposit(), and JetExtractor().

◆ theDR_Max

double muonisolation::JetExtractor::theDR_Max

private

Definition at line 59 of file JetExtractor.h.

Referenced by deposit().

◆ theDR_Veto

double muonisolation::JetExtractor::theDR_Veto

private

Definition at line 58 of file JetExtractor.h.

Referenced by deposit().

◆ theExcludeMuonVeto

bool muonisolation::JetExtractor::theExcludeMuonVeto

private

Definition at line 62 of file JetExtractor.h.

Referenced by deposit().

◆ theFieldToken

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> muonisolation::JetExtractor::theFieldToken

private

Definition at line 54 of file JetExtractor.h.

Referenced by deposit().

◆ theJetCollectionToken

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

private

Definition at line 50 of file JetExtractor.h.

Referenced by deposit().

◆ thePrintTimeReport

bool muonisolation::JetExtractor::thePrintTimeReport

private

Definition at line 70 of file JetExtractor.h.

◆ thePropagatorName

std::string muonisolation::JetExtractor::thePropagatorName

private

Definition at line 52 of file JetExtractor.h.

Referenced by deposit().

◆ theService

std::unique_ptr<MuonServiceProxy> muonisolation::JetExtractor::theService

private

the event setup proxy, it takes care the services update

Definition at line 65 of file JetExtractor.h.

Referenced by deposit(), and JetExtractor().

◆ theThreshold

double muonisolation::JetExtractor::theThreshold

private

Definition at line 57 of file JetExtractor.h.

Referenced by deposit().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ JetExtractor() [1/2]

◆ JetExtractor() [2/2]

◆ ~JetExtractor()

Member Function Documentation

◆ deposit()

◆ fillVetos()

Member Data Documentation

◆ theAssociator

◆ theAssociatorParameters

◆ theDR_Max

◆ theDR_Veto

◆ theExcludeMuonVeto

◆ theFieldToken

◆ theJetCollectionToken

◆ thePrintTimeReport

◆ thePropagatorName

◆ theService

◆ theThreshold