#include <CaloExtractorByAssociator.h>

Inheritance diagram for muonisolation::CaloExtractorByAssociator:

Public Member Functions
	CaloExtractorByAssociator ()
	constructors More...

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

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
	no-op: by design of this extractor the deposits are pulled out all at a time More...

virtual std::vector < reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
	return deposits for 3 calorimeter subdetectors (ecal, hcal, ho) – in this order More...

virtual void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks)
	allows to set extra vetoes (in addition to the muon) – no-op at this point More...

virtual	~CaloExtractorByAssociator ()
	destructor More...

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::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	~IsoDepositExtractor ()
	Destructor. More...

Private Member Functions
double	noiseEcal (const CaloTower &tower) const
	Determine noise for HCAL and ECAL (take some defaults for the time being) More...

double	noiseHcal (const CaloTower &tower) const

double	noiseHOcal (const CaloTower &tower) const

double	noiseRecHit (const DetId &detId) const

Private Attributes
TrackDetectorAssociator *	theAssociator

TrackAssociatorParameters *	theAssociatorParameters
	associator, its' parameters and the propagator More...

bool	theCenterConeOnCalIntersection

std::vector< std::string >	theDepositInstanceLabels
	multiple deposits: labels – expect 3 labels beginning with "e", "h", "ho" More...

std::string	theDepositLabel
	Label of deposit – suggest to set to "" (all info is in collection name anyways) More...

double	theDR_Max
	max cone size in which towers are considered More...

double	theDR_Veto_E
	cone sizes inside which the Et (towers) are not counted More...

double	theDR_Veto_H

double	theDR_Veto_HO

double	theNoise_EB

double	theNoise_EE

double	theNoise_HB

double	theNoise_HE

double	theNoise_HO

double	theNoiseTow_EB

double	theNoiseTow_EE

bool	thePrintTimeReport
	flag to turn on/off printing of a time report More...

std::string	thePropagatorName
	propagator name to feed into the track associator More...

MuonServiceProxy *	theService
	the event setup proxy, it takes care the services update More...

double	theThreshold_E

double	theThreshold_H

double	theThreshold_HO

bool	theUseRecHitsFlag
	use towers or rec hits More...

std::vector< DetId >	theVetoCollection
	Vector of calo Ids to veto – not used. More...

Detailed Description

Definition at line 34 of file CaloExtractorByAssociator.h.

Constructor & Destructor Documentation

muonisolation::CaloExtractorByAssociator::CaloExtractorByAssociator ( )

inline

constructors

Definition at line 39 of file CaloExtractorByAssociator.h.

39 {};

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

Definition at line 35 of file CaloExtractorByAssociator.cc.

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

                                                                                                         :
   theUseRecHitsFlag(par.getParameter<bool>("UseRecHitsFlag")),
   theDepositLabel(par.getUntrackedParameter<string>("DepositLabel")),
   theDepositInstanceLabels(par.getParameter<std::vector<std::string> >("DepositInstanceLabels")),
   thePropagatorName(par.getParameter<std::string>("PropagatorName")),
   theThreshold_E(par.getParameter<double>("Threshold_E")),
   theThreshold_H(par.getParameter<double>("Threshold_H")),
   theThreshold_HO(par.getParameter<double>("Threshold_HO")),
   theDR_Veto_E(par.getParameter<double>("DR_Veto_E")),
   theDR_Veto_H(par.getParameter<double>("DR_Veto_H")),
   theDR_Veto_HO(par.getParameter<double>("DR_Veto_HO")),
   theCenterConeOnCalIntersection(par.getParameter<bool>("CenterConeOnCalIntersection")),
   theDR_Max(par.getParameter<double>("DR_Max")),
   theNoise_EB(par.getParameter<double>("Noise_EB")),
   theNoise_EE(par.getParameter<double>("Noise_EE")),
   theNoise_HB(par.getParameter<double>("Noise_HB")),
   theNoise_HE(par.getParameter<double>("Noise_HE")),
   theNoise_HO(par.getParameter<double>("Noise_HO")),
   theNoiseTow_EB(par.getParameter<double>("NoiseTow_EB")),
   theNoiseTow_EE(par.getParameter<double>("NoiseTow_EE")),
   theService(0),
   theAssociator(0),
   thePrintTimeReport(par.getUntrackedParameter<bool>("PrintTimeReport"))
 {
   ParameterSet serviceParameters = par.getParameter<ParameterSet>("ServiceParameters");
   theService = new MuonServiceProxy(serviceParameters);
 
   theAssociatorParameters = new TrackAssociatorParameters(par.getParameter<edm::ParameterSet>("TrackAssociatorParameters"));
   theAssociator = new TrackDetectorAssociator();
 }

CaloExtractorByAssociator::~CaloExtractorByAssociator ( )

virtual

destructor

Definition at line 66 of file CaloExtractorByAssociator.cc.

References gather_cfg::cout, TimingReport::current(), TimingReport::dump(), theAssociator, theAssociatorParameters, thePrintTimeReport, and theService.

                                                      {
   if (thePrintTimeReport) TimingReport::current()->dump(std::cout);
   if (theAssociatorParameters) delete theAssociatorParameters;
   if (theService) delete theService;
   if (theAssociator) delete theAssociator;
 }

Member Function Documentation

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

virtual

no-op: by design of this extractor the deposits are pulled out all at a time

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 81 of file CaloExtractorByAssociator.cc.

References reco::TrackBase::eta(), and reco::TrackBase::phi().

 {
   IsoDeposit::Direction muonDir(muon.eta(), muon.phi());
   IsoDeposit dep(muonDir );
 
 //   LogWarning("CaloExtractorByAssociator")
 //     <<"single deposit is not an option here\n"
 //     <<"use ::deposits --> extract all and reweight as necessary";
 
   return dep;
 
 }

std::vector< IsoDeposit > CaloExtractorByAssociator::deposits	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

virtual

return deposits for 3 calorimeter subdetectors (ecal, hcal, ho) – in this order

Make separate deposits: for ECAL, HCAL, HO.

check configuration consistency could've been made at construction stage (fix later?)

this should be (eventually) set to the eta-phi of the crossing point of a straight line tangent to a muon at IP and the calorimeter

each deposit type veto is at the point of intersect with that detector

do things based on rec-hits here too much copy-pasting now (refactor later?)

first check if the hit is inside the veto cone by dR-alone

and now pitch those in the crossed list

first check if the hit is inside the veto cone by dR-alone

and now pitch those in the crossed list

first check if the hit is inside the veto cone by dR-alone

and now pitch those in the crossed list

use calo towers

first check if the tower is inside the veto cone by dR-alone

and now pitch those in the crossed list

add the Et of the tower to deposits if it's not a vetoed; put into muonEnergy otherwise

Reimplemented from reco::isodeposit::IsoDepositExtractor.

Definition at line 96 of file CaloExtractorByAssociator.cc.

References reco::IsoDeposit::addCandEnergy(), reco::IsoDeposit::addDeposit(), TrackDetectorAssociator::associate(), ecalTB2006H4_GenSimDigiReco_cfg::bField, compare_using_db::compare, TrackDetMatchInfo::crossedEcalIds, TrackDetMatchInfo::crossedHcalIds, TrackDetMatchInfo::crossedHOIds, TrackDetMatchInfo::crossedTowerIds, reco::deltaR(), CaloRecHit::detid(), reco::IsoDeposit::Veto::dR, TrackDetMatchInfo::ecalRecHits, CaloTower::emEnergy(), CaloTower::emEt(), CaloRecHit::energy(), relval_parameters_module::energy, PV3DBase< T, PVType, FrameType >::eta(), reco::TrackBase::eta(), reco::LeafCandidate::eta(), edm::EventSetup::get(), CaloTower::hadEnergy(), CaloTower::hadEt(), TrackDetMatchInfo::hcalRecHits, TrackDetMatchInfo::hoRecHits, CaloTower::id(), reco::TransientTrack::initialFreeState(), LogDebug, noiseEcal(), noiseHcal(), noiseHOcal(), noiseRecHit(), CaloTower::outerEnergy(), CaloTower::outerEt(), PV3DBase< T, PVType, FrameType >::phi(), reco::TrackBase::phi(), reco::LeafCandidate::phi(), DetId::rawId(), TrackDetectorAssociator::setPropagator(), reco::IsoDeposit::setVeto(), AlCaHLTBitMon_QueryRunRegistry::string, theAssociator, theAssociatorParameters, theCenterConeOnCalIntersection, theDepositInstanceLabels, theDR_Max, theDR_Veto_E, theDR_Veto_H, theDR_Veto_HO, thePropagatorName, theService, theThreshold_E, theThreshold_H, theThreshold_HO, theUseRecHitsFlag, TrackDetMatchInfo::towers, TrackDetMatchInfo::trkGlobPosAtEcal, TrackDetMatchInfo::trkGlobPosAtHcal, TrackDetMatchInfo::trkGlobPosAtHO, reco::IsoDeposit::veto(), and reco::IsoDeposit::Veto::vetoDir.

 {
   theService->update(eventSetup);
   theAssociator->setPropagator(&*(theService->propagator(thePropagatorName)));
 
   if (theDepositInstanceLabels.size() != 3){
     LogError("MuonIsolation")<<"Configuration is inconsistent: Need 3 deposit instance labels";
   }
   if (! theDepositInstanceLabels[0].compare(0,1, std::string("e")) == 0
       || ! theDepositInstanceLabels[1].compare(0,1, std::string("h")) == 0
       || ! theDepositInstanceLabels[2].compare(0,2, std::string("ho")) == 0){
     LogWarning("MuonIsolation")<<"Deposit instance labels do not look like  (e*, h*, ho*):"
                    <<"proceed at your own risk. The extractor interprets lab0=from ecal; lab1=from hcal; lab2=from ho";
   }
 
   typedef IsoDeposit::Veto Veto;
   IsoDeposit::Direction muonDir(muon.eta(), muon.phi());
 
   IsoDeposit depEcal(muonDir);
   IsoDeposit depHcal(muonDir);
   IsoDeposit depHOcal(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);
 
   depEcal.setVeto(Veto(reco::isodeposit::Direction(mInfo.trkGlobPosAtEcal.eta(), mInfo.trkGlobPosAtEcal.phi()),
                theDR_Veto_E));
   depHcal.setVeto(Veto(reco::isodeposit::Direction(mInfo.trkGlobPosAtHcal.eta(), mInfo.trkGlobPosAtHcal.phi()),
                theDR_Veto_H));
   depHOcal.setVeto(Veto(reco::isodeposit::Direction(mInfo.trkGlobPosAtHO.eta(), mInfo.trkGlobPosAtHO.phi()),
             theDR_Veto_HO));
 
   if (theCenterConeOnCalIntersection){
     reco::isodeposit::Direction dirTmp = depEcal.veto().vetoDir;
     double dRtmp = depEcal.veto().dR;
     depEcal = IsoDeposit(dirTmp); depEcal.setVeto(Veto(dirTmp, dRtmp));
 
     dirTmp = depHcal.veto().vetoDir;
     dRtmp = depHcal.veto().dR;
     depHcal = IsoDeposit(dirTmp); depHcal.setVeto(Veto(dirTmp, dRtmp));
 
     dirTmp = depHOcal.veto().vetoDir;
     dRtmp = depHOcal.veto().dR;
     depHOcal = IsoDeposit(dirTmp); depHOcal.setVeto(Veto(dirTmp, dRtmp));
   }
 
   if (theUseRecHitsFlag){
     edm::ESHandle<CaloGeometry> caloGeom;
     eventSetup.get<CaloGeometryRecord>().get(caloGeom);
 
     //Ecal
     std::vector<const EcalRecHit*>::const_iterator eHitCI = mInfo.ecalRecHits.begin();
     for (; eHitCI != mInfo.ecalRecHits.end(); ++eHitCI){
       const EcalRecHit* eHitCPtr = *eHitCI;
       GlobalPoint eHitPos = caloGeom->getPosition(eHitCPtr->detid());
       double deltar0 = reco::deltaR(muon, eHitPos);
       double cosTheta = 1./cosh(eHitPos.eta());
       double energy = eHitCPtr->energy();
       double et = energy*cosTheta;
       if (deltar0 > theDR_Max
       || ! (et > theThreshold_E && energy > 3*noiseRecHit(eHitCPtr->detid()))) continue;
 
       bool vetoHit = false;
       double deltar = reco::deltaR(mInfo.trkGlobPosAtEcal, eHitPos);
       if (deltar < theDR_Veto_E ){
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Veto ECAL hit: Calo deltaR= " << deltar;
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Calo eta phi ethcal: " << eHitPos.eta() << " " << eHitPos.phi() << " " << et;
     vetoHit = true;
       }
       if (! vetoHit){
     for (unsigned int iH = 0; iH< mInfo.crossedEcalIds.size() && ! vetoHit; ++iH){
       if (mInfo.crossedEcalIds[iH].rawId() == eHitCPtr->detid().rawId()) vetoHit = true;
     }
       }
 
       if (vetoHit ){
     depEcal.addCandEnergy(et);
       } else {
     depEcal.addDeposit(reco::isodeposit::Direction(eHitPos.eta(), eHitPos.phi()), et);
       }
     }
 
     //Hcal
     std::vector<const HBHERecHit*>::const_iterator hHitCI = mInfo.hcalRecHits.begin();
     for (; hHitCI != mInfo.hcalRecHits.end(); ++hHitCI){
       const HBHERecHit* hHitCPtr = *hHitCI;
       GlobalPoint hHitPos = caloGeom->getPosition(hHitCPtr->detid());
       double deltar0 = reco::deltaR(muon, hHitPos);
       double cosTheta = 1./cosh(hHitPos.eta());
       double energy = hHitCPtr->energy();
       double et = energy*cosTheta;
       if (deltar0 > theDR_Max
       || ! (et > theThreshold_H && energy > 3*noiseRecHit(hHitCPtr->detid()))) continue;
 
       bool vetoHit = false;
       double deltar = reco::deltaR(mInfo.trkGlobPosAtHcal, hHitPos);
       if (deltar < theDR_Veto_H ){
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Veto HBHE hit: Calo deltaR= " << deltar;
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Calo eta phi ethcal: " << hHitPos.eta() << " " << hHitPos.phi() << " " << et;
     vetoHit = true;
       }
       if (! vetoHit){
     for (unsigned int iH = 0; iH< mInfo.crossedHcalIds.size() && ! vetoHit; ++iH){
       if (mInfo.crossedHcalIds[iH].rawId() == hHitCPtr->detid().rawId()) vetoHit = true;
     }
       }
 
       if (vetoHit ){
     depHcal.addCandEnergy(et);
       } else {
     depHcal.addDeposit(reco::isodeposit::Direction(hHitPos.eta(), hHitPos.phi()), et);
       }
     }
 
     //HOcal
     std::vector<const HORecHit*>::const_iterator hoHitCI = mInfo.hoRecHits.begin();
     for (; hoHitCI != mInfo.hoRecHits.end(); ++hoHitCI){
       const HORecHit* hoHitCPtr = *hoHitCI;
       GlobalPoint hoHitPos = caloGeom->getPosition(hoHitCPtr->detid());
       double deltar0 = reco::deltaR(muon, hoHitPos);
       double cosTheta = 1./cosh(hoHitPos.eta());
       double energy = hoHitCPtr->energy();
       double et = energy*cosTheta;
       if (deltar0 > theDR_Max
       || ! (et > theThreshold_HO && energy > 3*noiseRecHit(hoHitCPtr->detid()))) continue;
 
       bool vetoHit = false;
       double deltar = reco::deltaR(mInfo.trkGlobPosAtHO, hoHitPos);
       if (deltar < theDR_Veto_HO ){
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Veto HO hit: Calo deltaR= " << deltar;
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Calo eta phi ethcal: " << hoHitPos.eta() << " " << hoHitPos.phi() << " " << et;
     vetoHit = true;
       }
       if (! vetoHit){
     for (unsigned int iH = 0; iH< mInfo.crossedHOIds.size() && ! vetoHit; ++iH){
       if (mInfo.crossedHOIds[iH].rawId() == hoHitCPtr->detid().rawId()) vetoHit = true;
     }
       }
 
       if (vetoHit ){
     depHOcal.addCandEnergy(et);
       } else {
     depHOcal.addDeposit(reco::isodeposit::Direction(hoHitPos.eta(), hoHitPos.phi()), et);
       }
     }
 
 
   } else {
     std::vector<const CaloTower*>::const_iterator calCI = mInfo.towers.begin();
     for (; calCI != mInfo.towers.end(); ++calCI){
       const CaloTower* calCPtr = *calCI;
       double deltar0 = reco::deltaR(muon,*calCPtr);
       if (deltar0>theDR_Max) continue;
 
       //even more copy-pasting .. need to refactor
       double etecal = calCPtr->emEt();
       double eecal = calCPtr->emEnergy();
       bool doEcal = etecal>theThreshold_E && eecal>3*noiseEcal(*calCPtr);
       double ethcal = calCPtr->hadEt();
       double ehcal = calCPtr->hadEnergy();
       bool doHcal = ethcal>theThreshold_H && ehcal>3*noiseHcal(*calCPtr);
       double ethocal = calCPtr->outerEt();
       double ehocal = calCPtr->outerEnergy();
       bool doHOcal = ethocal>theThreshold_HO && ehocal>3*noiseHOcal(*calCPtr);
       if ((!doEcal) && (!doHcal) && (!doHcal)) continue;
 
       bool vetoTowerEcal = false;
       double deltarEcal = reco::deltaR(mInfo.trkGlobPosAtEcal, *calCPtr);
       if (deltarEcal < theDR_Veto_E ){
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Veto ecal tower: Calo deltaR= " << deltarEcal;
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Calo eta phi ethcal: " << calCPtr->eta() << " " << calCPtr->phi() << " " << ethcal;
     vetoTowerEcal = true;
       }
       bool vetoTowerHcal = false;
       double deltarHcal = reco::deltaR(mInfo.trkGlobPosAtHcal, *calCPtr);
       if (deltarHcal < theDR_Veto_H ){
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Veto hcal tower: Calo deltaR= " << deltarHcal;
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Calo eta phi ethcal: " << calCPtr->eta() << " " << calCPtr->phi() << " " << ethcal;
     vetoTowerHcal = true;
       }
       bool vetoTowerHOCal = false;
       double deltarHOcal = reco::deltaR(mInfo.trkGlobPosAtHO, *calCPtr);
       if (deltarHOcal < theDR_Veto_HO ){
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Veto HO tower: Calo deltaR= " << deltarHOcal;
     LogDebug("RecoMuon|CaloExtractorByAssociator")
       << " >>> Calo eta phi ethcal: " << calCPtr->eta() << " " << calCPtr->phi() << " " << ethcal;
     vetoTowerHOCal = true;
       }
 
       if (! (vetoTowerHOCal && vetoTowerHcal &&  vetoTowerEcal )){
     for (unsigned int iH = 0; iH< mInfo.crossedTowerIds.size(); ++iH){
       if (mInfo.crossedTowerIds[iH].rawId() == calCPtr->id().rawId()){
         vetoTowerEcal = true;
         vetoTowerHcal = true;
         vetoTowerHOCal = true;
         break;
       }
     }
       }
 
       reco::isodeposit::Direction towerDir(calCPtr->eta(), calCPtr->phi());
       if (doEcal){
     if (vetoTowerEcal) depEcal.addCandEnergy(etecal);
     else depEcal.addDeposit(towerDir, etecal);
       }
       if (doHcal){
     if (vetoTowerHcal) depHcal.addCandEnergy(ethcal);
     else depHcal.addDeposit(towerDir, ethcal);
       }
       if (doHOcal){
     if (vetoTowerHOCal) depHOcal.addCandEnergy(ethocal);
     else depHOcal.addDeposit(towerDir, ethocal);
       }
     }
   }
 
   std::vector<IsoDeposit> resultDeps;
   resultDeps.push_back(depEcal);
   resultDeps.push_back(depHcal);
   resultDeps.push_back(depHOcal);
 
   return resultDeps;
 
 }

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

virtual

allows to set extra vetoes (in addition to the muon) – no-op at this point

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 73 of file CaloExtractorByAssociator.cc.

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

double CaloExtractorByAssociator::noiseEcal ( const CaloTower & tower ) const

private

Determine noise for HCAL and ECAL (take some defaults for the time being)

Definition at line 356 of file CaloExtractorByAssociator.cc.

References reco::LeafCandidate::eta(), eta(), theNoiseTow_EB, and theNoiseTow_EE.

Referenced by deposits().

                                                                         {
       double noise = theNoiseTow_EB;
       double eta = tower.eta();
       if (fabs(eta)>1.479) noise = theNoiseTow_EE;
       return noise;
 }

double CaloExtractorByAssociator::noiseHcal ( const CaloTower & tower ) const

private

Definition at line 363 of file CaloExtractorByAssociator.cc.

References reco::LeafCandidate::eta(), theNoise_HB, and theNoise_HE.

Referenced by deposits().

                                                                         {
   double noise = fabs(tower.eta())> 1.479 ? theNoise_HE : theNoise_HB;
   return noise;
 }

double CaloExtractorByAssociator::noiseHOcal ( const CaloTower & tower ) const

private

Definition at line 368 of file CaloExtractorByAssociator.cc.

References theNoise_HO.

Referenced by deposits().

                                                                          {
       double noise = theNoise_HO;
       return noise;
 }

double CaloExtractorByAssociator::noiseRecHit ( const DetId & detId ) const

private

Definition at line 374 of file CaloExtractorByAssociator.cc.

References DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, DetId::subdetId(), theNoise_EB, theNoise_EE, theNoise_HB, theNoise_HE, and theNoise_HO.

Referenced by deposits().

                                                                       {
   double  noise = 100;
   DetId::Detector det = detId.det();
   if (det == DetId::Ecal){
     EcalSubdetector subDet = (EcalSubdetector)(detId.subdetId());
     if (subDet == EcalBarrel){
       noise = theNoise_EB;
     } else if (subDet == EcalEndcap){
       noise = theNoise_EE;
     }
   } else if (det == DetId::Hcal){
     HcalSubdetector subDet = (HcalSubdetector)(detId.subdetId());
     if (subDet == HcalBarrel){
       noise = theNoise_HB;
     } else if (subDet == HcalEndcap){
       noise = theNoise_HE;
     } else if (subDet == HcalOuter){
       noise = theNoise_HO;
     }
   }
   return noise;
 }

Member Data Documentation

TrackDetectorAssociator* muonisolation::CaloExtractorByAssociator::theAssociator

private

Definition at line 103 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), deposits(), and ~CaloExtractorByAssociator().

TrackAssociatorParameters* muonisolation::CaloExtractorByAssociator::theAssociatorParameters

private

associator, its' parameters and the propagator

Definition at line 102 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), deposits(), and ~CaloExtractorByAssociator().

bool muonisolation::CaloExtractorByAssociator::theCenterConeOnCalIntersection

private

centers the cone on the veto direction – makes more sense for very displaced tracks like in cosmics

Definition at line 80 of file CaloExtractorByAssociator.h.

Referenced by deposits().

std::vector<std::string> muonisolation::CaloExtractorByAssociator::theDepositInstanceLabels

private

multiple deposits: labels – expect 3 labels beginning with "e", "h", "ho"

Definition at line 63 of file CaloExtractorByAssociator.h.

Referenced by deposits().

std::string muonisolation::CaloExtractorByAssociator::theDepositLabel

private

Label of deposit – suggest to set to "" (all info is in collection name anyways)

Definition at line 60 of file CaloExtractorByAssociator.h.

double muonisolation::CaloExtractorByAssociator::theDR_Max

private

max cone size in which towers are considered

Definition at line 82 of file CaloExtractorByAssociator.h.

Referenced by deposits().

double muonisolation::CaloExtractorByAssociator::theDR_Veto_E

private

cone sizes inside which the Et (towers) are not counted

Definition at line 75 of file CaloExtractorByAssociator.h.

Referenced by deposits().

double muonisolation::CaloExtractorByAssociator::theDR_Veto_H

private

Definition at line 76 of file CaloExtractorByAssociator.h.

Referenced by deposits().

double muonisolation::CaloExtractorByAssociator::theDR_Veto_HO

private

Definition at line 77 of file CaloExtractorByAssociator.h.

Referenced by deposits().

double muonisolation::CaloExtractorByAssociator::theNoise_EB

private

the noise "sigmas" for a hit or tower to be considered consider if Energy > 3.*sigma

Definition at line 86 of file CaloExtractorByAssociator.h.

Referenced by noiseRecHit().

double muonisolation::CaloExtractorByAssociator::theNoise_EE

private

Definition at line 87 of file CaloExtractorByAssociator.h.

Referenced by noiseRecHit().

double muonisolation::CaloExtractorByAssociator::theNoise_HB

private

Definition at line 88 of file CaloExtractorByAssociator.h.

Referenced by noiseHcal(), and noiseRecHit().

double muonisolation::CaloExtractorByAssociator::theNoise_HE

private

Definition at line 89 of file CaloExtractorByAssociator.h.

Referenced by noiseHcal(), and noiseRecHit().

double muonisolation::CaloExtractorByAssociator::theNoise_HO

private

Definition at line 90 of file CaloExtractorByAssociator.h.

Referenced by noiseHOcal(), and noiseRecHit().

double muonisolation::CaloExtractorByAssociator::theNoiseTow_EB

private

Definition at line 91 of file CaloExtractorByAssociator.h.

Referenced by noiseEcal().

double muonisolation::CaloExtractorByAssociator::theNoiseTow_EE

private

Definition at line 92 of file CaloExtractorByAssociator.h.

Referenced by noiseEcal().

bool muonisolation::CaloExtractorByAssociator::thePrintTimeReport

private

flag to turn on/off printing of a time report

Definition at line 106 of file CaloExtractorByAssociator.h.

Referenced by ~CaloExtractorByAssociator().

std::string muonisolation::CaloExtractorByAssociator::thePropagatorName

private

propagator name to feed into the track associator

Definition at line 66 of file CaloExtractorByAssociator.h.

Referenced by deposits().

MuonServiceProxy* muonisolation::CaloExtractorByAssociator::theService

private

the event setup proxy, it takes care the services update

Definition at line 98 of file CaloExtractorByAssociator.h.

Referenced by CaloExtractorByAssociator(), deposits(), and ~CaloExtractorByAssociator().

double muonisolation::CaloExtractorByAssociator::theThreshold_E

private

Cone cuts and thresholds min values of Et to be included in deposits

Definition at line 70 of file CaloExtractorByAssociator.h.

Referenced by deposits().

double muonisolation::CaloExtractorByAssociator::theThreshold_H

private

Definition at line 71 of file CaloExtractorByAssociator.h.

Referenced by deposits().

double muonisolation::CaloExtractorByAssociator::theThreshold_HO

private

Definition at line 72 of file CaloExtractorByAssociator.h.

Referenced by deposits().

bool muonisolation::CaloExtractorByAssociator::theUseRecHitsFlag

private

use towers or rec hits

Definition at line 57 of file CaloExtractorByAssociator.h.

Referenced by deposits().

std::vector<DetId> muonisolation::CaloExtractorByAssociator::theVetoCollection

private

Vector of calo Ids to veto – not used.

Definition at line 95 of file CaloExtractorByAssociator.h.

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation