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

22  :
23  theCaloTowerCollectionToken(iC.consumes<CaloTowerCollection>(par.getParameter<edm::InputTag>("CaloTowerCollectionLabel"))),
24  theDepositLabel(par.getUntrackedParameter<string>("DepositLabel")),
25  theWeight_E(par.getParameter<double>("Weight_E")),
26  theWeight_H(par.getParameter<double>("Weight_H")),
27  theThreshold_E(par.getParameter<double>("Threshold_E")),
28  theThreshold_H(par.getParameter<double>("Threshold_H")),
29  theDR_Veto_E(par.getParameter<double>("DR_Veto_E")),
30  theDR_Veto_H(par.getParameter<double>("DR_Veto_H")),
31  theDR_Max(par.getParameter<double>("DR_Max")),
32  vertexConstraintFlag_XY(par.getParameter<bool>("Vertex_Constraint_XY")),
33  vertexConstraintFlag_Z(par.getParameter<bool>("Vertex_Constraint_Z"))
34 {
35 }
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:52
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 89 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, HLT_FULL_cff::towers, vertexConstraintFlag_XY, and vertexConstraintFlag_Z.

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

38 {
39  theVetoCollection.clear();
40 
41  Handle<CaloTowerCollection> towers;
42  event.getByToken(theCaloTowerCollectionToken,towers);
43 
44  edm::ESHandle<CaloGeometry> caloGeom;
45  eventSetup.get<CaloGeometryRecord>().get(caloGeom);
46 
47  edm::ESHandle<MagneticField> bField;
48  eventSetup.get<IdealMagneticFieldRecord>().get(bField);
49  double bz = bField->inInverseGeV(GlobalPoint(0.,0.,0.)).z();
50 
51  TrackCollection::const_iterator mu;
52  TrackCollection::const_iterator muEnd(muons.end());
53 
54  CaloTowerCollection::const_iterator cal;
55  CaloTowerCollection::const_iterator calEnd(towers->end());
56 
57  for ( mu = muons.begin(); mu != muEnd; ++mu ) {
58  for ( cal = towers->begin(); cal != calEnd; ++cal ) {
60  double dEta = fabs(mu->eta()-cal->eta());
61  if (fabs(dEta) > theDR_Max) continue;
62 
63  double deltar0 = reco::deltaR(*mu,*cal);
64  if (deltar0>theDR_Max) continue;
65 
66  double etecal = cal->emEt();
67  double eecal = cal->emEnergy();
68  bool doEcal = theWeight_E>0 && etecal>theThreshold_E && eecal>3*noiseEcal(*cal);
69  double ethcal = cal->hadEt();
70  double ehcal = cal->hadEnergy();
71  bool doHcal = theWeight_H>0 && ethcal>theThreshold_H && ehcal>3*noiseHcal(*cal);
72  if ((!doEcal) && (!doHcal)) continue;
73 
74  DetId calId = cal->id();
75  GlobalPoint endpos = caloGeom->getPosition(calId);
76  GlobalPoint muatcal = MuonAtCaloPosition(*mu,bz,endpos, vertexConstraintFlag_XY, vertexConstraintFlag_Z);
77  double deltar = reco::deltaR(muatcal,endpos);
78 
79  if (doEcal) {
80  if (deltar<theDR_Veto_E) theVetoCollection.push_back(calId);
81  } else {
82  if (deltar<theDR_Veto_H) theVetoCollection.push_back(calId);
83  }
84  }
85  }
86 
87 }
muonisolation::CaloExtractor::theCaloTowerCollectionToken
edm::EDGetTokenT< CaloTowerCollection > theCaloTowerCollectionToken
Definition: CaloExtractor.h:37
muonisolation::CaloExtractor::noiseEcal
double noiseEcal(const CaloTower &tower) const
Definition: CaloExtractor.cc:290
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
edm::SortedCollection< CaloTower >::const_iterator
std::vector< CaloTower >::const_iterator const_iterator
Definition: SortedCollection.h:82
HLT_FULL_cff.towers
tuple towers
Definition: HLT_FULL_cff.py:12739
reco::deltaR
auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
Definition: deltaR.h:28
ecalTB2006H4_GenSimDigiReco_cfg.bField
tuple bField
Definition: ecalTB2006H4_GenSimDigiReco_cfg.py:235
muonisolation::CaloExtractor::noiseHcal
double noiseHcal(const CaloTower &tower) const
Definition: CaloExtractor.cc:297
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:197
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 197 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(), PV3DBase< T, PVType, FrameType >::x(), x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by deposit(), and fillVetos().

197  {
198  double qoverp= muon.qoverp();
199  double cur = bz*muon.charge()/muon.pt();
200  double phi0 = muon.phi();
201  double dca = muon.dxy();
202  double theta = muon.theta();
203  double dz = muon.dz();
204 
205  //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
206  //<< " Pt(GeV): " << muon.pt()
207  //<< ", phi0 " << muon.phi0()
208  //<< ", eta " << muon.eta();
209  //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
210  //<< " d0 " << muon.d0()
211  //<< ", dz " << muon.dz();
212  //LogDebug("Muon|RecoMuon|L2MuonIsolationProducer")
213  //<< " rhocal " << endpos.perp()
214  //<< ", zcal " << endpos.z();
215 
216  if (fixVxy && fixVz) {
217  // Note that here we assume no correlation between XY and Z projections
218  // This should be a reasonable approximation for our purposes
219  double errd02 = muon.covariance(muon.i_dxy,muon.i_dxy);
220  if (pow(muon.dxy(),2)<4*errd02) {
221  phi0 -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_phi)
222  /errd02;
223  cur -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_qoverp)
224  /errd02 * (cur/qoverp);
225  dca = 0;
226  }
227  double errdsz2 = muon.covariance(muon.i_dsz,muon.i_dsz);
228  if (pow(muon.dsz(),2)<4*errdsz2) {
229  theta += muon.dsz()*muon.covariance(muon.i_dsz,muon.i_lambda)
230  /errdsz2;
231  dz = 0;
232  }
233  } else if (fixVxy) {
234  double errd02 = muon.covariance(muon.i_dxy,muon.i_dxy);
235  if (pow(muon.dxy(),2)<4*errd02) {
236  phi0 -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_phi)
237  /errd02;
238  cur -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_qoverp)
239  /errd02 * (cur/qoverp);
240  theta += muon.dxy()*muon.covariance(muon.i_dxy,muon.i_lambda)
241  /errd02;
242  dz -= muon.dxy()*muon.covariance(muon.i_dxy,muon.i_dsz)
243  /errd02 * muon.p()/muon.pt();
244  dca = 0;
245  }
246  } else if (fixVz) {
247  double errdsz2 = muon.covariance(muon.i_dsz,muon.i_dsz);
248  if (pow(muon.dsz(),2)<4*errdsz2) {
249  theta += muon.dsz()*muon.covariance(muon.i_dsz,muon.i_lambda)
250  /errdsz2;
251  phi0 -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_phi)
252  /errdsz2;
253  cur -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_qoverp)
254  /errdsz2 * (cur/qoverp);
255  dca -= muon.dsz()*muon.covariance(muon.i_dsz,muon.i_dxy)
256  /errdsz2;
257  dz = 0;
258  }
259  }
260 
261  double sphi0 = sin(phi0);
262  double cphi0 = cos(phi0);
263 
264  double xsin = endpos.x()*sphi0 - endpos.y()*cphi0;
265  double xcos = endpos.x()*cphi0 + endpos.y()*sphi0;
266  double fcdca = fabs(1-cur*dca);
267  double phif = atan2( fcdca*sphi0-cur*endpos.x()
268  , fcdca*cphi0+cur*endpos.y());
269  double tphif2 = tan(0.5*(phif-phi0));
270  double dcaf = dca + xsin + xcos*tphif2;
271 
272  double x = endpos.x() - dcaf*sin(phif);
273  double y = endpos.y() + dcaf*cos(phif);
274 
275  double deltas = (x-muon.vx())*cphi0 + (y-muon.vy())*sphi0;
276  double deltaphi = normalizedPhi(phif-phi0);
277  if (deltaphi!=0) deltas = deltas*deltaphi/sin(deltaphi);
278 
279  double z =dz;
280  double tantheta = tan(theta);
281  if (tantheta!=0) {
282  z += deltas/tan(theta);
283  } else {
284  z = endpos.z();
285  }
286 
287  return GlobalPoint(x,y,z);
288 }
reco::TrackBase::qoverp
double qoverp() const
q / p
Definition: TrackBase.h:568
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:610
reco::TrackBase::theta
double theta() const
polar angle
Definition: TrackBase.h:574
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:640
deltas
Definition: deltas.h:18
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:598
reco::TrackBase::covariance
CovarianceMatrix covariance() const
return track covariance matrix
Definition: TrackBase.h:726
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:616
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:604
reco::TrackBase::vy
double vy() const
y coordinate of the reference point on track
Definition: TrackBase.h:658
reco::TrackBase::charge
int charge() const
track electric charge
Definition: TrackBase.h:562
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:586
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
normalizedPhi
T normalizedPhi(T phi)
Definition: normalizedPhi.h:8
reco::TrackBase::vx
double vx() const
x coordinate of the reference point on track
Definition: TrackBase.h:652
double CaloExtractor::noiseEcal ( const CaloTower &  tower) const
private

Definition at line 290 of file CaloExtractor.cc.

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

Referenced by deposit(), and fillVetos().

290  {
291  double noise = 0.04;
292  double eta = tower.eta();
293  if (fabs(eta)>1.479) noise = 0.15;
294  return noise;
295 }
HLT_FULL_cff.noise
tuple noise
Definition: HLT_FULL_cff.py:4549
reco::LeafCandidate::eta
virtual double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:137
double CaloExtractor::noiseHcal ( const CaloTower &  tower) const
private

Definition at line 297 of file CaloExtractor.cc.

References HLT_FULL_cff::noise.

Referenced by deposit(), and fillVetos().

297  {
298  double noise = 0.2;
299  return noise;
300 }
HLT_FULL_cff.noise
tuple noise
Definition: HLT_FULL_cff.py:4549

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