CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
muonisolation::CaloExtractor Class Reference

#include <CaloExtractor.h>

Inheritance diagram for muonisolation::CaloExtractor:
reco::isodeposit::IsoDepositExtractor

Public Member Functions

 CaloExtractor ()
 
 CaloExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)
 
virtual reco::IsoDeposit deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
 
virtual void fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks)
 
virtual ~CaloExtractor ()
 
- 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 ~IsoDepositExtractor ()
 Destructor. More...
 

Static Public Member Functions

static GlobalPoint MuonAtCaloPosition (const reco::Track &muon, const double bz, const GlobalPoint &endpos, bool fixVxy=false, bool fixVz=false)
 Extrapolate muons to calorimeter-object positions. More...
 

Private Member Functions

double noiseEcal (const CaloTower &tower) const
 
double noiseHcal (const CaloTower &tower) const
 

Private Attributes

edm::EDGetTokenT
< CaloTowerCollection > 		theCaloTowerCollectionToken
 
std::string theDepositLabel
 
double theDR_Max
 
double theDR_Veto_E
 
double theDR_Veto_H
 
double theThreshold_E
 
double theThreshold_H
 
std::vector< DetId > theVetoCollection
 
double theWeight_E
 
double theWeight_H
 
bool vertexConstraintFlag_XY
 
bool vertexConstraintFlag_Z
 

Detailed Description

Definition at line 20 of file CaloExtractor.h.

Constructor & Destructor Documentation

muonisolation::CaloExtractor::CaloExtractor ( )
inline

Definition at line 24 of file CaloExtractor.h.

24 {};
CaloExtractor::CaloExtractor ( const edm::ParameterSet &  par,
edm::ConsumesCollector &&  iC 
)

Definition at line 25 of file CaloExtractor.cc.

25  :
26  theCaloTowerCollectionToken(iC.consumes<CaloTowerCollection>(par.getParameter<edm::InputTag>("CaloTowerCollectionLabel"))),
27  theDepositLabel(par.getUntrackedParameter<string>("DepositLabel")),
28  theWeight_E(par.getParameter<double>("Weight_E")),
29  theWeight_H(par.getParameter<double>("Weight_H")),
30  theThreshold_E(par.getParameter<double>("Threshold_E")),
31  theThreshold_H(par.getParameter<double>("Threshold_H")),
32  theDR_Veto_E(par.getParameter<double>("DR_Veto_E")),
33  theDR_Veto_H(par.getParameter<double>("DR_Veto_H")),
34  theDR_Max(par.getParameter<double>("DR_Max")),
35  vertexConstraintFlag_XY(par.getParameter<bool>("Vertex_Constraint_XY")),
36  vertexConstraintFlag_Z(par.getParameter<bool>("Vertex_Constraint_Z"))
37 {
38 }
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:41
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
muonisolation::CaloExtractor::theCaloTowerCollectionToken
edm::EDGetTokenT< CaloTowerCollection > theCaloTowerCollectionToken
Definition: CaloExtractor.h:37
virtual muonisolation::CaloExtractor::~CaloExtractor ( )
inlinevirtual

Definition at line 27 of file CaloExtractor.h.

27 {}

Member Function Documentation

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

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

make this abit faster

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 92 of file CaloExtractor.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::bField, reco::deltaR(), PV3DBase< T, PVType, FrameType >::eta(), reco::TrackBase::eta(), spr::find(), edm::EventSetup::get(), LogDebug, MuonAtCaloPosition(), noiseEcal(), noiseHcal(), PV3DBase< T, PVType, FrameType >::phi(), reco::TrackBase::phi(), reco::TrackBase::pt(), theCaloTowerCollectionToken, theDR_Max, theDR_Veto_E, theDR_Veto_H, theThreshold_E, theThreshold_H, theVetoCollection, theWeight_E, theWeight_H, vertexConstraintFlag_XY, and vertexConstraintFlag_Z.

93 {
94  IsoDeposit dep(muon.eta(), muon.phi() );
95  LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
96  << " >>> Muon: pt " << muon.pt()
97  << " eta " << muon.eta()
98  << " phi " << muon.phi();
99 
100  Handle<CaloTowerCollection> towers;
101  event.getByToken(theCaloTowerCollectionToken,towers);
102 
103  edm::ESHandle<CaloGeometry> caloGeom;
104  eventSetup.get<CaloGeometryRecord>().get(caloGeom);
105 
106  edm::ESHandle<MagneticField> bField;
107  eventSetup.get<IdealMagneticFieldRecord>().get(bField);
108  double bz = bField->inInverseGeV(GlobalPoint(0.,0.,0.)).z();
109 
110  CaloTowerCollection::const_iterator cal;
111  CaloTowerCollection::const_iterator calEnd(towers->end());
112  for ( cal = towers->begin(); cal != calEnd; ++cal ) {
114  double dEta = fabs(muon.eta()-cal->eta());
115  if (fabs(dEta) > theDR_Max) continue;
116 
117  double deltar0 = reco::deltaR(muon,*cal);
118  if (deltar0>theDR_Max) continue;
119 
120  double etecal = cal->emEt();
121  double eecal = cal->emEnergy();
122  bool doEcal = theWeight_E>0 && etecal>theThreshold_E && eecal>3*noiseEcal(*cal);
123  double ethcal = cal->hadEt();
124  double ehcal = cal->hadEnergy();
125  bool doHcal = theWeight_H>0 && ethcal>theThreshold_H && ehcal>3*noiseHcal(*cal);
126  if ((!doEcal) && (!doHcal)) continue;
127 
128  DetId calId = cal->id();
129  GlobalPoint endpos = caloGeom->getPosition(calId);
130  GlobalPoint muatcal = MuonAtCaloPosition(muon,bz,endpos,vertexConstraintFlag_XY, vertexConstraintFlag_Z);
131  double deltar = reco::deltaR(muatcal,endpos);
132 
133  if (deltar<theDR_Veto_H) {
134  dep.setVeto(IsoDeposit::Veto(reco::isodeposit::Direction(muatcal.eta(), muatcal.phi()), theDR_Veto_H));
135  }
136 
137  if (doEcal) {
138  if (deltar<theDR_Veto_E) {
139  double calodep = theWeight_E*etecal;
140  if (doHcal) calodep += theWeight_H*ethcal;
141  dep.addCandEnergy(calodep);
142  LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
143  << " >>> Calo deposit inside veto (with ECAL): deltar " << deltar
144  << " calodep " << calodep
145  << " ecaldep " << etecal
146  << " hcaldep " << ethcal
147  << " eta " << cal->eta()
148  << " phi " << cal->phi();
149  continue;
150  }
151  } else {
152  if (deltar<theDR_Veto_H) {
153  dep.addCandEnergy(theWeight_H*ethcal);
154  LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
155  << " >>> Calo deposit inside veto (no ECAL): deltar " << deltar
156  << " calodep " << theWeight_H*ethcal
157  << " eta " << cal->eta()
158  << " phi " << cal->phi();
159  continue;
160  }
161  }
162 
163  if (std::find(theVetoCollection.begin(), theVetoCollection.end()
164  , calId)!=theVetoCollection.end()) {
165  LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
166  << " >>> Deposits belongs to other track: deltar, etecal, ethcal= "
167  << deltar << ", " << etecal << ", " << ethcal;
168  continue;
169  }
170 
171  if (doEcal) {
172  if (deltar>theDR_Veto_E) {
173  double calodep = theWeight_E*etecal;
174  if (doHcal) calodep += theWeight_H*ethcal;
175  dep.addDeposit(reco::isodeposit::Direction(endpos.eta(), endpos.phi()),calodep);
176  LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
177  << " >>> Calo deposit (with ECAL): deltar " << deltar
178  << " calodep " << calodep
179  << " ecaldep " << etecal
180  << " hcaldep " << ethcal
181  << " eta " << cal->eta()
182  << " phi " << cal->phi();
183  }
184  } else {
185  if (deltar>theDR_Veto_H) {
186  dep.addDeposit(reco::isodeposit::Direction(endpos.eta(), endpos.phi()),theWeight_H*ethcal);
187  LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
188  << " >>> Calo deposit (no ECAL): deltar " << deltar
189  << " calodep " << theWeight_H*ethcal
190  << " eta " << cal->eta()
191  << " phi " << cal->phi();
192  }
193  }
194  }
195 
196  return dep;
197 
198 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
muonisolation::CaloExtractor::theCaloTowerCollectionToken
edm::EDGetTokenT< CaloTowerCollection > theCaloTowerCollectionToken
Definition: CaloExtractor.h:37
muonisolation::CaloExtractor::noiseEcal
double noiseEcal(const CaloTower &tower) const
Definition: CaloExtractor.cc:293
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:82
reco::deltaR
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
trackingTruthProducerFastSim_cfi.muon
tuple muon
Definition: trackingTruthProducerFastSim_cfi.py:10
ecalTB2006H4_GenSimDigiReco_cfg.bField
tuple bField
Definition: ecalTB2006H4_GenSimDigiReco_cfg.py:235
muonisolation::CaloExtractor::noiseHcal
double noiseHcal(const CaloTower &tower) const
Definition: CaloExtractor.cc:300
DetId
Definition: DetId.h:18
muonisolation::CaloExtractor::MuonAtCaloPosition
static GlobalPoint MuonAtCaloPosition(const reco::Track &muon, const double bz, const GlobalPoint &endpos, bool fixVxy=false, bool fixVz=false)
Extrapolate muons to calorimeter-object positions.
Definition: CaloExtractor.cc:200
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
muonisolation::CaloExtractor::theVetoCollection
std::vector< DetId > theVetoCollection
Definition: CaloExtractor.h:54
void CaloExtractor::fillVetos ( const edm::Event &  ev,
const edm::EventSetup &  evSetup,
const reco::TrackCollection &  tracks 
)
virtual

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

make this abit faster

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 40 of file CaloExtractor.cc.

References ecalTB2006H4_GenSimDigiReco_cfg::bField, reco::deltaR(), edm::EventSetup::get(), RPCpg::mu, MuonAtCaloPosition(), noiseEcal(), noiseHcal(), theCaloTowerCollectionToken, theDR_Max, theDR_Veto_E, theDR_Veto_H, theThreshold_E, theThreshold_H, theVetoCollection, theWeight_E, theWeight_H, vertexConstraintFlag_XY, and vertexConstraintFlag_Z.

41 {
42  theVetoCollection.clear();
43 
44  Handle<CaloTowerCollection> towers;
45  event.getByToken(theCaloTowerCollectionToken,towers);
46 
47  edm::ESHandle<CaloGeometry> caloGeom;
48  eventSetup.get<CaloGeometryRecord>().get(caloGeom);
49 
50  edm::ESHandle<MagneticField> bField;
51  eventSetup.get<IdealMagneticFieldRecord>().get(bField);
52  double bz = bField->inInverseGeV(GlobalPoint(0.,0.,0.)).z();
53 
54  TrackCollection::const_iterator mu;
55  TrackCollection::const_iterator muEnd(muons.end());
56 
57  CaloTowerCollection::const_iterator cal;
58  CaloTowerCollection::const_iterator calEnd(towers->end());
59 
60  for ( mu = muons.begin(); mu != muEnd; ++mu ) {
61  for ( cal = towers->begin(); cal != calEnd; ++cal ) {
63  double dEta = fabs(mu->eta()-cal->eta());
64  if (fabs(dEta) > theDR_Max) continue;
65 
66  double deltar0 = reco::deltaR(*mu,*cal);
67  if (deltar0>theDR_Max) continue;
68 
69  double etecal = cal->emEt();
70  double eecal = cal->emEnergy();
71  bool doEcal = theWeight_E>0 && etecal>theThreshold_E && eecal>3*noiseEcal(*cal);
72  double ethcal = cal->hadEt();
73  double ehcal = cal->hadEnergy();
74  bool doHcal = theWeight_H>0 && ethcal>theThreshold_H && ehcal>3*noiseHcal(*cal);
75  if ((!doEcal) && (!doHcal)) continue;
76 
77  DetId calId = cal->id();
78  GlobalPoint endpos = caloGeom->getPosition(calId);
79  GlobalPoint muatcal = MuonAtCaloPosition(*mu,bz,endpos, vertexConstraintFlag_XY, vertexConstraintFlag_Z);
80  double deltar = reco::deltaR(muatcal,endpos);
81 
82  if (doEcal) {
83  if (deltar<theDR_Veto_E) theVetoCollection.push_back(calId);
84  } else {
85  if (deltar<theDR_Veto_H) theVetoCollection.push_back(calId);
86  }
87  }
88  }
89 
90 }
muonisolation::CaloExtractor::theCaloTowerCollectionToken
edm::EDGetTokenT< CaloTowerCollection > theCaloTowerCollectionToken
Definition: CaloExtractor.h:37
muonisolation::CaloExtractor::noiseEcal
double noiseEcal(const CaloTower &tower) const
Definition: CaloExtractor.cc:293
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:82
reco::deltaR
double deltaR(const T1 &t1, const T2 &t2)
Definition: deltaR.h:48
ecalTB2006H4_GenSimDigiReco_cfg.bField
tuple bField
Definition: ecalTB2006H4_GenSimDigiReco_cfg.py:235
muonisolation::CaloExtractor::noiseHcal
double noiseHcal(const CaloTower &tower) const
Definition: CaloExtractor.cc:300
RPCpg::mu
const int mu
Definition: Constants.h:22
DetId
Definition: DetId.h:18
muonisolation::CaloExtractor::MuonAtCaloPosition
static GlobalPoint MuonAtCaloPosition(const reco::Track &muon, const double bz, const GlobalPoint &endpos, bool fixVxy=false, bool fixVz=false)
Extrapolate muons to calorimeter-object positions.
Definition: CaloExtractor.cc:200
patZpeak.muons
tuple muons
Definition: patZpeak.py:38
muonisolation::CaloExtractor::theVetoCollection
std::vector< DetId > theVetoCollection
Definition: CaloExtractor.h:54
GlobalPoint CaloExtractor::MuonAtCaloPosition ( const reco::Track &  muon,
const double  bz,
const GlobalPoint &  endpos,
bool  fixVxy = false,
bool  fixVz = false 
)
static

Extrapolate muons to calorimeter-object positions.

Definition at line 200 of file CaloExtractor.cc.

References reco::TrackBase::charge(), funct::cos(), reco::TrackBase::covariance(), reco::TrackBase::dsz(), reco::TrackBase::dxy(), reco::TrackBase::dz(), reco::TrackBase::i_dsz, reco::TrackBase::i_dxy, reco::TrackBase::i_lambda, reco::TrackBase::i_phi, reco::TrackBase::i_qoverp, normalizedPhi(), reco::TrackBase::p(), reco::TrackBase::phi(), funct::pow(), reco::TrackBase::pt(), reco::TrackBase::qoverp(), funct::sin(), funct::tan(), theta(), reco::TrackBase::theta(), reco::TrackBase::vx(), reco::TrackBase::vy(), x, PV3DBase< T, PVType, FrameType >::x(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by deposit(), and fillVetos().

200  {
201  double qoverp= muon.qoverp();
202  double cur = bz*muon.charge()/muon.pt();
203  double phi0 = muon.phi();
204  double dca = muon.dxy();
205  double theta = muon.theta();
206  double dz = muon.dz();
207 
208  //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
209  //<< " Pt(GeV): " << muon.pt()
210  //<< ", phi0 " << muon.phi0()
211  //<< ", eta " << muon.eta();
212  //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
213  //<< " d0 " << muon.d0()
214  //<< ", dz " << muon.dz();
215  //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
216  //<< " rhocal " << endpos.perp()
217  //<< ", zcal " << endpos.z();
218 
219  if (fixVxy && fixVz) {
220  // Note that here we assume no correlation between XY and Z projections
221  // This should be a reasonable approximation for our purposes
222  double errd02 = muon.covariance(muon.i_dxy,muon.i_dxy);
223  if (pow(muon.dxy(),2)<4*errd02) {
224  phi0 -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_phi)
225  /errd02;
226  cur -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_qoverp)
227  /errd02 * (cur/qoverp);
228  dca = 0;
229  }
230  double errdsz2 = muon.covariance(muon.i_dsz,muon.i_dsz);
231  if (pow(muon.dsz(),2)<4*errdsz2) {
232  theta += muon.dsz()*muon.covariance(muon.i_dsz,muon.i_lambda)
233  /errdsz2;
234  dz = 0;
235  }
236  } else if (fixVxy) {
237  double errd02 = muon.covariance(muon.i_dxy,muon.i_dxy);
238  if (pow(muon.dxy(),2)<4*errd02) {
239  phi0 -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_phi)
240  /errd02;
241  cur -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_qoverp)
242  /errd02 * (cur/qoverp);
243  theta += muon.dxy()*muon.covariance(muon.i_dxy,muon.i_lambda)
244  /errd02;
245  dz -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_dsz)
246  /errd02 * muon.p()/muon.pt();
247  dca = 0;
248  }
249  } else if (fixVz) {
250  double errdsz2 = muon.covariance(muon.i_dsz,muon.i_dsz);
251  if (pow(muon.dsz(),2)<4*errdsz2) {
252  theta += muon.dsz()*muon.covariance(muon.i_dsz,muon.i_lambda)
253  /errdsz2;
254  phi0 -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_phi)
255  /errdsz2;
256  cur -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_qoverp)
257  /errdsz2 * (cur/qoverp);
258  dca -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_dxy)
259  /errdsz2;
260  dz = 0;
261  }
262  }
263 
264  double sphi0 = sin(phi0);
265  double cphi0 = cos(phi0);
266 
267  double xsin = endpos.x()*sphi0 - endpos.y()*cphi0;
268  double xcos = endpos.x()*cphi0 + endpos.y()*sphi0;
269  double fcdca = fabs(1-cur*dca);
270  double phif = atan2( fcdca*sphi0-cur*endpos.x()
271  , fcdca*cphi0+cur*endpos.y());
272  double tphif2 = tan(0.5*(phif-phi0));
273  double dcaf = dca + xsin + xcos*tphif2;
274 
275  double x = endpos.x() - dcaf*sin(phif);
276  double y = endpos.y() + dcaf*cos(phif);
277 
278  double deltas = (x-muon.vx())*cphi0 + (y-muon.vy())*sphi0;
279  double deltaphi = normalizedPhi(phif-phi0);
280  if (deltaphi!=0) deltas = deltas*deltaphi/sin(deltaphi);
281 
282  double z =dz;
283  double tantheta = tan(theta);
284  if (tantheta!=0) {
285  z += deltas/tan(theta);
286  } else {
287  z = endpos.z();
288  }
289 
290  return GlobalPoint(x,y,z);
291 }
reco::TrackBase::qoverp
double qoverp() const
q / p
Definition: TrackBase.h:621
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:663
reco::TrackBase::theta
double theta() const
polar angle
Definition: TrackBase.h:627
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:170
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:693
reco::TrackBase::dsz
double dsz() const
dsz parameter (THIS IS NOT the SZ impact parameter to (0,0,0) if refPoint is far from (0...
Definition: TrackBase.h:651
detailsBasic3DVector::z
float float float z
Definition: extBasic3DVector.h:15
reco::TrackBase::covariance
CovarianceMatrix covariance() const
return track covariance matrix
Definition: TrackBase.h:779
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:669
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
funct::tan
Tan< T >::type tan(const T &t)
Definition: Tan.h:22
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:657
reco::TrackBase::vy
double vy() const
y coordinate of the reference point on track
Definition: TrackBase.h:711
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:615
x
Definition: DDAxes.h:10
normalizedPhi
double normalizedPhi(double phi)
Definition: normalizedPhi.cc:5
reco::TrackBase::dxy
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
Definition: TrackBase.h:639
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
reco::TrackBase::vx
double vx() const
x coordinate of the reference point on track
Definition: TrackBase.h:705
double CaloExtractor::noiseEcal ( const CaloTower &  tower) const
private

Definition at line 293 of file CaloExtractor.cc.

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

Referenced by deposit(), and fillVetos().

293  {
294  double noise = 0.04;
295  double eta = tower.eta();
296  if (fabs(eta)>1.479) noise = 0.15;
297  return noise;
298 }
eta
T eta() const
Definition: Basic3DVectorLD.h:177
reco::LeafCandidate::eta
virtual float eta() const
momentum pseudorapidity
Definition: LeafCandidate.h:184
double CaloExtractor::noiseHcal ( const CaloTower &  tower) const
private

Definition at line 300 of file CaloExtractor.cc.

Referenced by deposit(), and fillVetos().

300  {
301  double noise = 0.2;
302  return noise;
303 }

Member Data Documentation

edm::EDGetTokenT<CaloTowerCollection> muonisolation::CaloExtractor::theCaloTowerCollectionToken
private

Definition at line 37 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

std::string muonisolation::CaloExtractor::theDepositLabel
private

Definition at line 40 of file CaloExtractor.h.

double muonisolation::CaloExtractor::theDR_Max
private

Definition at line 49 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theDR_Veto_E
private

Definition at line 47 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theDR_Veto_H
private

Definition at line 48 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theThreshold_E
private

Definition at line 45 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theThreshold_H
private

Definition at line 46 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

std::vector<DetId> muonisolation::CaloExtractor::theVetoCollection
private

Definition at line 54 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theWeight_E
private

Definition at line 43 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

double muonisolation::CaloExtractor::theWeight_H
private

Definition at line 44 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

bool muonisolation::CaloExtractor::vertexConstraintFlag_XY
private

Definition at line 50 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

bool muonisolation::CaloExtractor::vertexConstraintFlag_Z
private

Definition at line 51 of file CaloExtractor.h.

Referenced by deposit(), and fillVetos().

Generated for CMSSW Reference Manual by doxygen 1.8.5