CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
reco::CentralityProducer Class Reference
Inheritance diagram for reco::CentralityProducer:
edm::EDProducer edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 CentralityProducer (const edm::ParameterSet &)
 
 ~CentralityProducer ()
 
- Public Member Functions inherited from edm::EDProducer
 EDProducer ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducer ()
 
- Public Member Functions inherited from edm::ProducerBase
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription
const &)> 
registrationCallback () const
 used by the fwk to register list of products More...
 
virtual ~ProducerBase ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 
itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

virtual void beginJob () override
 
virtual void endJob () override
 
virtual void produce (edm::Event &, const edm::EventSetup &) override
 

Private Attributes

edm::ESHandle< CaloGeometrycGeo
 
bool doPixelCut_
 
double hfEtaCut_
 
bool lowGainZDC_
 
double midRapidityRange_
 
bool produceEcalhits_
 
bool produceETmidRap_
 
bool produceHFhits_
 
bool produceHFtowers_
 
bool producePixelhits_
 
bool producePixelTracks_
 
bool produceTracks_
 
bool produceZDChits_
 
bool recoLevel_
 
bool reuseAny_
 
edm::EDGetTokenT< CentralityreuseTag_
 
edm::EDGetTokenT
< EcalRecHitCollection
srcEBhits_
 
edm::EDGetTokenT
< EcalRecHitCollection
srcEEhits_
 
edm::EDGetTokenT
< HFRecHitCollection
srcHFhits_
 
edm::EDGetTokenT
< SiPixelRecHitCollection
srcPixelhits_
 
edm::EDGetTokenT< TrackCollectionsrcPixelTracks_
 
edm::EDGetTokenT
< CaloTowerCollection
srcTowers_
 
edm::EDGetTokenT< TrackCollectionsrcTracks_
 
edm::EDGetTokenT
< VertexCollection
srcVertex_
 
edm::EDGetTokenT
< ZDCRecHitCollection
srcZDChits_
 
edm::ESHandle< TrackerGeometrytGeo
 
double trackEtaCut_
 
double trackPtCut_
 
reco::TrackBase::TrackQuality trackQuality_
 
bool useQuality_
 

Additional Inherited Members

- Public Types inherited from edm::EDProducer
typedef EDProducer ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef
ProductRegistryHelper::TypeLabelList 
TypeLabelList
 
- Static Public Member Functions inherited from edm::EDProducer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 54 of file CentralityProducer.cc.

Constructor & Destructor Documentation

reco::CentralityProducer::CentralityProducer ( const edm::ParameterSet iConfig)
explicit

Definition at line 118 of file CentralityProducer.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

118  {
119  //register your products
120  produceHFhits_ = iConfig.getParameter<bool>("produceHFhits");
121  produceHFtowers_ = iConfig.getParameter<bool>("produceHFtowers");
122  produceEcalhits_ = iConfig.getParameter<bool>("produceEcalhits");
123  produceZDChits_ = iConfig.getParameter<bool>("produceZDChits");
124  produceETmidRap_ = iConfig.getParameter<bool>("produceETmidRapidity");
125  producePixelhits_ = iConfig.getParameter<bool>("producePixelhits");
126  produceTracks_ = iConfig.getParameter<bool>("produceTracks");
127  producePixelTracks_ = iConfig.getParameter<bool>("producePixelTracks");
128 
129  midRapidityRange_ = iConfig.getParameter<double>("midRapidityRange");
130  trackPtCut_ = iConfig.getParameter<double>("trackPtCut");
131  trackEtaCut_ = iConfig.getParameter<double>("trackEtaCut");
132 
133  hfEtaCut_ = iConfig.getParameter<double>("hfEtaCut");
134 
135  if(produceHFhits_) srcHFhits_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcHFhits"));
136  if(produceHFtowers_ || produceETmidRap_) srcTowers_ = consumes<CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("srcTowers"));
137 
138  if(produceEcalhits_){
139  srcEBhits_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcEBhits"));
140  srcEEhits_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcEEhits"));
141  }
142  if(produceZDChits_){
143  srcZDChits_ = consumes<ZDCRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcZDChits"));
144  lowGainZDC_ = iConfig.getParameter<bool>("lowGainZDC");
145  }
146  if(producePixelhits_){
147  srcPixelhits_ = consumes<SiPixelRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcPixelhits"));
148  doPixelCut_ = iConfig.getParameter<bool>("doPixelCut");
149  srcVertex_ = consumes<VertexCollection>(iConfig.getParameter<edm::InputTag>("srcVertex"));
150  }
151  if(produceTracks_) {
152  srcTracks_ = consumes<TrackCollection>(iConfig.getParameter<edm::InputTag>("srcTracks"));
153  srcVertex_ = consumes<VertexCollection>(iConfig.getParameter<edm::InputTag>("srcVertex"));
154  }
155  if(producePixelTracks_) srcPixelTracks_ = consumes<TrackCollection>(iConfig.getParameter<edm::InputTag>("srcPixelTracks"));
156 
157  reuseAny_ = iConfig.getParameter<bool>("reUseCentrality");
158  if(reuseAny_) reuseTag_ = consumes<Centrality>(iConfig.getParameter<edm::InputTag>("srcReUse"));
159 
160  useQuality_ = iConfig.getParameter<bool>("UseQuality");
162 
163  produces<Centrality>();
164 
165 }
reco::TrackBase::TrackQuality trackQuality_
T getParameter(std::string const &) const
edm::EDGetTokenT< EcalRecHitCollection > srcEBhits_
edm::EDGetTokenT< Centrality > reuseTag_
edm::EDGetTokenT< VertexCollection > srcVertex_
edm::EDGetTokenT< CaloTowerCollection > srcTowers_
edm::EDGetTokenT< HFRecHitCollection > srcHFhits_
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:114
edm::EDGetTokenT< ZDCRecHitCollection > srcZDChits_
edm::EDGetTokenT< SiPixelRecHitCollection > srcPixelhits_
edm::EDGetTokenT< TrackCollection > srcPixelTracks_
edm::EDGetTokenT< TrackCollection > srcTracks_
edm::EDGetTokenT< EcalRecHitCollection > srcEEhits_
reco::CentralityProducer::~CentralityProducer ( )

Definition at line 168 of file CentralityProducer.cc.

169 {
170 
171  // do anything here that needs to be done at desctruction time
172  // (e.g. close files, deallocate resources etc.)
173 
174 }

Member Function Documentation

void reco::CentralityProducer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 467 of file CentralityProducer.cc.

468 {
469 }
void reco::CentralityProducer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 473 of file CentralityProducer.cc.

474 {
475 }
void reco::CentralityProducer::produce ( edm::Event iEvent,
const edm::EventSetup iSetup 
)
overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 183 of file CentralityProducer.cc.

References funct::abs(), edmNew::DetSet< T >::begin(), edmNew::DetSetVector< T >::begin(), SiPixelRecHit::cluster(), edmNew::DetSet< T >::detId(), reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), reco::TrackBase::dzError(), edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), CaloRecHit::energy(), EcalRecHit::energy(), reco::LeafCandidate::eta(), eta(), reco::TrackBase::eta(), edmNew::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByToken(), i, ZDCRecHit::id(), HFRecHit::id(), EcalRecHit::id(), HcalDetId::ieta(), CaloTower::ietaAbs(), isHF(), BaseTrackerRecHit::localPosition(), ZDCRecHit::lowGainEnergy(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), reco::TrackBase::pt(), reco::TrackBase::ptError(), edm::Event::put(), reco::TrackBase::quality(), HI_PhotonSkim_cff::rechits, GeomDet::toGlobal(), HLT_25ns14e33_v1_cff::towers, testEve_cfg::tracks, and HcalZDCDetId::zside().

184 {
185 
187  if(produceEcalhits_) iSetup.get<CaloGeometryRecord>().get(cGeo);
188 
189  std::auto_ptr<Centrality> creco(new Centrality());
190  Handle<Centrality> inputCentrality;
191 
192  if(reuseAny_) iEvent.getByToken(reuseTag_,inputCentrality);
193 
194  if(produceHFhits_){
195  creco->etHFhitSumPlus_ = 0;
196  creco->etHFhitSumMinus_ = 0;
197 
199  iEvent.getByToken(srcHFhits_,hits);
200  for( size_t ihit = 0; ihit<hits->size(); ++ ihit){
201  const HFRecHit & rechit = (*hits)[ ihit ];
202  if(rechit.id().ieta() > 0 )
203  creco->etHFhitSumPlus_ += rechit.energy();
204  if(rechit.id().ieta() < 0)
205  creco->etHFhitSumMinus_ += rechit.energy();
206  }
207  }else{
208  if(reuseAny_){
209  creco->etHFhitSumMinus_ = inputCentrality->EtHFhitSumMinus();
210  creco->etHFhitSumPlus_ = inputCentrality->EtHFhitSumPlus();
211  }
212  }
213 
215  creco->etHFtowerSumPlus_ = 0;
216  creco->etHFtowerSumMinus_ = 0;
217  creco->etMidRapiditySum_ = 0;
218 
220 
221  iEvent.getByToken(srcTowers_,towers);
222 
223  for( size_t i = 0; i<towers->size(); ++ i){
224  const CaloTower & tower = (*towers)[ i ];
225  double eta = tower.eta();
226  bool isHF = tower.ietaAbs() > 29;
227  if(produceHFtowers_){
228  if(isHF && eta > 0){
229  creco->etHFtowerSumPlus_ += tower.pt();
230  if(eta > hfEtaCut_) creco->etHFtruncatedPlus_ += tower.pt();
231  }
232  if(isHF && eta < 0){
233  creco->etHFtowerSumMinus_ += tower.pt();
234  if(eta < -hfEtaCut_) creco->etHFtruncatedMinus_ += tower.pt();
235  }
236  }else{
237  if(reuseAny_){
238  creco->etHFtowerSumMinus_ = inputCentrality->EtHFtowerSumMinus();
239  creco->etHFtowerSumPlus_ = inputCentrality->EtHFtowerSumPlus();
240  creco->etHFtruncatedMinus_ = inputCentrality->EtHFtruncatedMinus();
241  creco->etHFtruncatedPlus_ = inputCentrality->EtHFtruncatedPlus();
242  }
243  }
244  if(produceETmidRap_){
245  if(std::abs(eta) < midRapidityRange_) creco->etMidRapiditySum_ += tower.pt()/(midRapidityRange_*2.);
246  }else if(reuseAny_) creco->etMidRapiditySum_ = inputCentrality->EtMidRapiditySum();
247  }
248  }else{
249  if(reuseAny_){
250  creco->etHFtowerSumMinus_ = inputCentrality->EtHFtowerSumMinus();
251  creco->etHFtowerSumPlus_ = inputCentrality->EtHFtowerSumPlus();
252  creco->etMidRapiditySum_ = inputCentrality->EtMidRapiditySum();
253  }
254  }
255 
256  if(produceEcalhits_){
257  creco->etEESumPlus_ = 0;
258  creco->etEESumMinus_ = 0;
259  creco->etEBSum_ = 0;
260 
263 
264  iEvent.getByToken(srcEBhits_,ebHits);
265  iEvent.getByToken(srcEEhits_,eeHits);
266 
267  for(unsigned int i = 0; i < ebHits->size(); ++i){
268  const EcalRecHit & hit= (*ebHits)[i];
269  const GlobalPoint& pos=cGeo->getPosition(hit.id());
270  double et = hit.energy()*(pos.perp()/pos.mag());
271  creco->etEBSum_ += et;
272  }
273 
274  for(unsigned int i = 0; i < eeHits->size(); ++i){
275  const EcalRecHit & hit= (*eeHits)[i];
276  const GlobalPoint& pos=cGeo->getPosition(hit.id());
277  double et = hit.energy()*(pos.perp()/pos.mag());
278  if(pos.z() > 0){
279  creco->etEESumPlus_ += et;
280  }else{
281  creco->etEESumMinus_ += et;
282  }
283  }
284  }else{
285  if(reuseAny_){
286  creco->etEESumMinus_ = inputCentrality->EtEESumMinus();
287  creco->etEESumPlus_ = inputCentrality->EtEESumPlus();
288  creco->etEBSum_ = inputCentrality->EtEBSum();
289  }
290  }
291 
292  if(producePixelhits_){
293  creco->pixelMultiplicity_ = 0;
296  iEvent.getByToken(srcPixelhits_,rchts);
297  rechits = rchts.product();
298  int nPixel =0 ;
299  for (SiPixelRecHitCollection::const_iterator it = rechits->begin(); it!=rechits->end();it++)
300  {
302  DetId detId = DetId(hits.detId());
303  SiPixelRecHitCollection::const_iterator recHitMatch = rechits->find(detId);
304  const SiPixelRecHitCollection::DetSet recHitRange = *recHitMatch;
305  for ( SiPixelRecHitCollection::DetSet::const_iterator recHitIterator = recHitRange.begin();
306  recHitIterator != recHitRange.end(); ++recHitIterator) {
307  // add selection if needed, now all hits.
308  if(doPixelCut_){
309  const SiPixelRecHit * recHit = &(*recHitIterator);
310  const PixelGeomDetUnit* pixelLayer = dynamic_cast<const PixelGeomDetUnit*> (tGeo->idToDet(recHit->geographicalId()));
311  GlobalPoint gpos = pixelLayer->toGlobal(recHit->localPosition());
312  math::XYZVector rechitPos(gpos.x(),gpos.y(),gpos.z());
313  double abeta = std::abs(rechitPos.eta());
314  int clusterSize = recHit->cluster()->size();
315  if ( abeta < 0.5 && clusterSize < 1) continue;
316  if ( abeta > 0.5 && abeta < 1 && clusterSize < 2) continue;
317  if ( abeta > 1. && abeta < 1.5 && clusterSize < 3) continue;
318  if ( abeta > 1.5 && abeta < 2. && clusterSize < 4) continue;
319  if ( abeta > 2. && abeta < 2.5 && clusterSize < 6) continue;
320  if ( abeta > 2.5 && abeta < 5 && clusterSize < 9) continue;
321  }
322  nPixel++;
323  }
324  }
325  creco->pixelMultiplicity_ = nPixel;
326  }else{
327  if(reuseAny_){
328  creco->pixelMultiplicity_ = inputCentrality->multiplicityPixel();
329  }
330  }
331 
332  if(produceTracks_){
333 
334  double vx=-999.;
335  double vy=-999.;
336  double vz=-999.;
337  double vxError=-999.;
338  double vyError=-999.;
339  double vzError=-999.;
340  math::XYZVector vtxPos(0,0,0);
341 
342  Handle<VertexCollection> recoVertices;
343  iEvent.getByToken(srcVertex_,recoVertices);
344  unsigned int daughter = 0;
345  int greatestvtx = 0;
346 
347  for (unsigned int i = 0 ; i< recoVertices->size(); ++i){
348  daughter = (*recoVertices)[i].tracksSize();
349  if( daughter > (*recoVertices)[greatestvtx].tracksSize()) greatestvtx = i;
350  }
351 
352  if(recoVertices->size()>0){
353  vx = (*recoVertices)[greatestvtx].position().x();
354  vy = (*recoVertices)[greatestvtx].position().y();
355  vz = (*recoVertices)[greatestvtx].position().z();
356  vxError = (*recoVertices)[greatestvtx].xError();
357  vyError = (*recoVertices)[greatestvtx].yError();
358  vzError = (*recoVertices)[greatestvtx].zError();
359  }
360 
361  vtxPos = math::XYZVector(vx,vy,vz);
362 
364  iEvent.getByToken(srcTracks_,tracks);
365  int nTracks = 0;
366 
367  double trackCounter = 0;
368  double trackCounterEta = 0;
369  double trackCounterEtaPt = 0;
370 
371  for(unsigned int i = 0 ; i < tracks->size(); ++i){
372  const Track& track = (*tracks)[i];
373  if(useQuality_ && !track.quality(trackQuality_)) continue;
374 
375  if( track.pt() > trackPtCut_) trackCounter++;
376  if(std::abs(track.eta())<trackEtaCut_) {
377  trackCounterEta++;
378  if (track.pt() > trackPtCut_) trackCounterEtaPt++;
379  }
380 
381  math::XYZPoint v1(vx,vy, vz);
382  double dz= track.dz(v1);
383  double dzsigma2 = track.dzError()*track.dzError()+vzError*vzError;
384  double dxy= track.dxy(v1);
385  double dxysigma2 = track.dxyError()*track.dxyError()+vxError*vyError;
386 
387  const double pterrcut = 0.1;
388  const double dzrelcut = 3.0;
389  const double dxyrelcut = 3.0;
390 
391  if( track.quality(trackQuality_) &&
392  track.pt()>0.4 && std::abs(track.eta())<2.4 &&
393  track.ptError()/track.pt() < pterrcut &&
394  dz*dz < dzrelcut*dzrelcut * dzsigma2 &&
395  dxy*dxy < dxyrelcut*dxyrelcut * dxysigma2 ){
396  nTracks++;
397  }
398  }
399 
400  creco->trackMultiplicity_ = nTracks;
401  creco->ntracksPtCut_ = trackCounter;
402  creco->ntracksEtaCut_ = trackCounterEta;
403  creco->ntracksEtaPtCut_ = trackCounterEtaPt;
404 
405  }else{
406  if(reuseAny_){
407  creco->trackMultiplicity_ = inputCentrality->Ntracks();
408  creco->ntracksPtCut_= inputCentrality->NtracksPtCut();
409  creco->ntracksEtaCut_= inputCentrality->NtracksEtaCut();
410  creco->ntracksEtaPtCut_= inputCentrality->NtracksEtaPtCut();
411  }
412  }
413 
415  Handle<TrackCollection> pixeltracks;
416  iEvent.getByToken(srcPixelTracks_,pixeltracks);
417  int nPixelTracks = pixeltracks->size();
418  creco->nPixelTracks_ = nPixelTracks;
419  }
420  else{
421  if(reuseAny_){
422  creco->nPixelTracks_ = inputCentrality->NpixelTracks();
423  }
424  }
425 
426  if(produceZDChits_){
427  creco->zdcSumPlus_ = 0;
428  creco->zdcSumMinus_ = 0;
429 
431  bool zdcAvailable = iEvent.getByToken(srcZDChits_,hits);
432  if(zdcAvailable){
433  for( size_t ihit = 0; ihit<hits->size(); ++ ihit){
434  const ZDCRecHit & rechit = (*hits)[ ihit ];
435  if(rechit.id().zside() > 0 ) {
436  if(lowGainZDC_){
437  creco->zdcSumPlus_ += rechit.lowGainEnergy();
438  }else{
439  creco->zdcSumPlus_ += rechit.energy();
440  }
441  }
442  if(rechit.id().zside() < 0) {
443  if(lowGainZDC_){
444  creco->zdcSumMinus_ += rechit.lowGainEnergy();
445  }else{
446  creco->zdcSumMinus_ += rechit.energy();
447  }
448  }
449  }
450  }else{
451  creco->zdcSumPlus_ = -9;
452  creco->zdcSumMinus_ = -9;
453  }
454  }else{
455  if(reuseAny_){
456  creco->zdcSumMinus_ = inputCentrality->zdcSumMinus();
457  creco->zdcSumPlus_ = inputCentrality->zdcSumPlus();
458  }
459  }
460 
461  iEvent.put(creco);
462 
463 }
reco::TrackBase::TrackQuality trackQuality_
int i
Definition: DBlmapReader.cc:9
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
const_iterator end(bool update=false) const
edm::EDGetTokenT< EcalRecHitCollection > srcEBhits_
T perp() const
Definition: PV3DBase.h:72
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:449
edm::EDGetTokenT< Centrality > reuseTag_
double dxyError() const
error on dxy
Definition: TrackBase.h:749
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:52
int zside() const
get the z-side of the cell (1/-1)
Definition: HcalZDCDetId.h:34
T y() const
Definition: PV3DBase.h:63
edm::EDGetTokenT< VertexCollection > srcVertex_
T eta() const
virtual double eta() const
momentum pseudorapidity
virtual double pt() const
transverse momentum
edm::EDGetTokenT< CaloTowerCollection > srcTowers_
T mag() const
Definition: PV3DBase.h:67
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:604
OrphanHandle< PROD > put(std::auto_ptr< PROD > product)
Put a new product.
Definition: Event.h:113
float energy() const
Definition: CaloRecHit.h:17
double pt() const
track transverse momentum
Definition: TrackBase.h:574
T z() const
Definition: PV3DBase.h:64
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:716
int ieta() const
get the cell ieta
Definition: HcalDetId.h:36
edm::EDGetTokenT< HFRecHitCollection > srcHFhits_
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
float energy() const
Definition: EcalRecHit.h:68
float lowGainEnergy() const
Definition: ZDCRecHit.h:21
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:562
double dzError() const
error on dz
Definition: TrackBase.h:767
Definition: DetId.h:18
bool isHF(int etabin, int depth)
DetId id() const
get the id
Definition: EcalRecHit.h:76
T const * product() const
Definition: Handle.h:81
tuple tracks
Definition: testEve_cfg.py:39
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
edm::EDGetTokenT< ZDCRecHitCollection > srcZDChits_
ClusterRef cluster() const
Definition: SiPixelRecHit.h:49
const T & get() const
Definition: EventSetup.h:55
id_type detId() const
Definition: DetSetNew.h:83
edm::ESHandle< CaloGeometry > cGeo
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:463
const_iterator find(id_type i, bool update=false) const
edm::EDGetTokenT< SiPixelRecHitCollection > srcPixelhits_
int ietaAbs() const
Definition: CaloTower.h:170
iterator end()
Definition: DetSetNew.h:70
edm::EDGetTokenT< TrackCollection > srcPixelTracks_
edm::EDGetTokenT< TrackCollection > srcTracks_
HcalZDCDetId id() const
get the id
Definition: ZDCRecHit.h:19
edm::ESHandle< TrackerGeometry > tGeo
HcalDetId id() const
Definition: HFRecHit.h:23
DetId geographicalId() const
edm::EDGetTokenT< EcalRecHitCollection > srcEEhits_
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:544
T x() const
Definition: PV3DBase.h:62
virtual LocalPoint localPosition() const
const_iterator begin(bool update=false) const
Our base class.
Definition: SiPixelRecHit.h:23
iterator begin()
Definition: DetSetNew.h:67

Member Data Documentation

edm::ESHandle<CaloGeometry> reco::CentralityProducer::cGeo
private

Definition at line 102 of file CentralityProducer.cc.

bool reco::CentralityProducer::doPixelCut_
private

Definition at line 78 of file CentralityProducer.cc.

double reco::CentralityProducer::hfEtaCut_
private

Definition at line 83 of file CentralityProducer.cc.

bool reco::CentralityProducer::lowGainZDC_
private

Definition at line 85 of file CentralityProducer.cc.

double reco::CentralityProducer::midRapidityRange_
private

Definition at line 80 of file CentralityProducer.cc.

bool reco::CentralityProducer::produceEcalhits_
private

Definition at line 70 of file CentralityProducer.cc.

bool reco::CentralityProducer::produceETmidRap_
private

Definition at line 72 of file CentralityProducer.cc.

bool reco::CentralityProducer::produceHFhits_
private

Definition at line 68 of file CentralityProducer.cc.

bool reco::CentralityProducer::produceHFtowers_
private

Definition at line 69 of file CentralityProducer.cc.

bool reco::CentralityProducer::producePixelhits_
private

Definition at line 73 of file CentralityProducer.cc.

bool reco::CentralityProducer::producePixelTracks_
private

Definition at line 76 of file CentralityProducer.cc.

bool reco::CentralityProducer::produceTracks_
private

Definition at line 74 of file CentralityProducer.cc.

bool reco::CentralityProducer::produceZDChits_
private

Definition at line 71 of file CentralityProducer.cc.

bool reco::CentralityProducer::recoLevel_
private

Definition at line 66 of file CentralityProducer.cc.

bool reco::CentralityProducer::reuseAny_
private

Definition at line 75 of file CentralityProducer.cc.

edm::EDGetTokenT<Centrality> reco::CentralityProducer::reuseTag_
private

Definition at line 96 of file CentralityProducer.cc.

edm::EDGetTokenT<EcalRecHitCollection> reco::CentralityProducer::srcEBhits_
private

Definition at line 90 of file CentralityProducer.cc.

edm::EDGetTokenT<EcalRecHitCollection> reco::CentralityProducer::srcEEhits_
private

Definition at line 89 of file CentralityProducer.cc.

edm::EDGetTokenT<HFRecHitCollection> reco::CentralityProducer::srcHFhits_
private

Definition at line 87 of file CentralityProducer.cc.

edm::EDGetTokenT<SiPixelRecHitCollection> reco::CentralityProducer::srcPixelhits_
private

Definition at line 92 of file CentralityProducer.cc.

edm::EDGetTokenT<TrackCollection> reco::CentralityProducer::srcPixelTracks_
private

Definition at line 94 of file CentralityProducer.cc.

edm::EDGetTokenT<CaloTowerCollection> reco::CentralityProducer::srcTowers_
private

Definition at line 88 of file CentralityProducer.cc.

edm::EDGetTokenT<TrackCollection> reco::CentralityProducer::srcTracks_
private

Definition at line 93 of file CentralityProducer.cc.

edm::EDGetTokenT<VertexCollection> reco::CentralityProducer::srcVertex_
private

Definition at line 95 of file CentralityProducer.cc.

edm::EDGetTokenT<ZDCRecHitCollection> reco::CentralityProducer::srcZDChits_
private

Definition at line 91 of file CentralityProducer.cc.

edm::ESHandle<TrackerGeometry> reco::CentralityProducer::tGeo
private

Definition at line 101 of file CentralityProducer.cc.

double reco::CentralityProducer::trackEtaCut_
private

Definition at line 82 of file CentralityProducer.cc.

double reco::CentralityProducer::trackPtCut_
private

Definition at line 81 of file CentralityProducer.cc.

reco::TrackBase::TrackQuality reco::CentralityProducer::trackQuality_
private

Definition at line 99 of file CentralityProducer.cc.

bool reco::CentralityProducer::useQuality_
private

Definition at line 98 of file CentralityProducer.cc.