CMS 3D CMS Logo
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
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

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

Private Attributes

edm::ESHandle< CaloGeometry > cGeo
 
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< Centrality > reuseTag_
 
edm::EDGetTokenT< EcalRecHitCollection > srcEBhits_
 
edm::EDGetTokenT< EcalRecHitCollection > srcEEhits_
 
edm::EDGetTokenT< HFRecHitCollection > srcHFhits_
 
edm::EDGetTokenT< SiPixelRecHitCollection > srcPixelhits_
 
edm::EDGetTokenT< TrackCollection > srcPixelTracks_
 
edm::EDGetTokenT< CaloTowerCollection > srcTowers_
 
edm::EDGetTokenT< TrackCollection > srcTracks_
 
edm::EDGetTokenT< VertexCollection > srcVertex_
 
edm::EDGetTokenT< ZDCRecHitCollection > srcZDChits_
 
edm::ESHandle< TrackerGeometry > tGeo
 
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
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- 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::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");
161  trackQuality_ = TrackBase::qualityByName(iConfig.getParameter<std::string>("TrackQuality"));
162 
163  produces<Centrality>();
164 
165 }
reco::CentralityProducer::trackQuality_
reco::TrackBase::TrackQuality trackQuality_
Definition: CentralityProducer.cc:99
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
reco::CentralityProducer::srcEBhits_
edm::EDGetTokenT< EcalRecHitCollection > srcEBhits_
Definition: CentralityProducer.cc:90
reco::CentralityProducer::reuseTag_
edm::EDGetTokenT< Centrality > reuseTag_
Definition: CentralityProducer.cc:96
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
reco::CentralityProducer::srcVertex_
edm::EDGetTokenT< VertexCollection > srcVertex_
Definition: CentralityProducer.cc:95
reco::CentralityProducer::srcTowers_
edm::EDGetTokenT< CaloTowerCollection > srcTowers_
Definition: CentralityProducer.cc:88
reco::CentralityProducer::srcHFhits_
edm::EDGetTokenT< HFRecHitCollection > srcHFhits_
Definition: CentralityProducer.cc:87
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:125
reco::CentralityProducer::srcZDChits_
edm::EDGetTokenT< ZDCRecHitCollection > srcZDChits_
Definition: CentralityProducer.cc:91
reco::CentralityProducer::srcPixelhits_
edm::EDGetTokenT< SiPixelRecHitCollection > srcPixelhits_
Definition: CentralityProducer.cc:92
reco::CentralityProducer::srcPixelTracks_
edm::EDGetTokenT< TrackCollection > srcPixelTracks_
Definition: CentralityProducer.cc:94
reco::CentralityProducer::srcTracks_
edm::EDGetTokenT< TrackCollection > srcTracks_
Definition: CentralityProducer.cc:93
reco::CentralityProducer::srcEEhits_
edm::EDGetTokenT< EcalRecHitCollection > srcEEhits_
Definition: CentralityProducer.cc:89
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.

References DEFINE_FWK_MODULE.

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

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(), particleFlowDisplacedVertexCandidate_cfi::dxy, reco::TrackBase::dxy(), reco::TrackBase::dxyError(), allConversions_cfi::dz, reco::TrackBase::dz(), reco::TrackBase::dzError(), edmNew::DetSet< T >::end(), edmNew::DetSetVector< T >::end(), CaloRecHit::energy(), EcalRecHit::energy(), stringResolutionProvider_cfi::et, stringResolutionProvider_cfi::eta, reco::LeafCandidate::eta(), reco::TrackBase::eta(), edmNew::DetSetVector< T >::find(), TrackingRecHit::geographicalId(), edm::EventSetup::get(), edm::Event::getByToken(), hfClusterShapes_cfi::hits, mps_fire::i, ZDCRecHit::id(), HFRecHit::id(), EcalRecHit::id(), HcalDetId::ieta(), CaloTower::ietaAbs(), HcalObjRepresent::isHF(), BaseTrackerRecHit::localPosition(), ZDCRecHit::lowGainEnergy(), eostools::move(), nTracks(), edm::Handle< T >::product(), reco::LeafCandidate::pt(), reco::TrackBase::pt(), reco::TrackBase::ptError(), edm::Event::put(), reco::TrackBase::quality(), rpcPointValidation_cfi::recHit, TrackInfoProducer_cfi::rechits, edm::SortedCollection< T, SORT >::size(), GeomDet::toGlobal(), HiIsolationCommonParameters_cff::track, l1t::tracks, and HcalZDCDetId::zside().

184 {
185 
186  if(producePixelhits_) iSetup.get<TrackerDigiGeometryRecord>().get(tGeo);
187  if(produceEcalhits_) iSetup.get<CaloGeometryRecord>().get(cGeo);
188 
189  auto creco = std::make_unique<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 
198  Handle<HFRecHitCollection> hits;
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 
214  if(produceHFtowers_ || produceETmidRap_){
215  creco->etHFtowerSumPlus_ = 0;
216  creco->etHFtowerSumMinus_ = 0;
217  creco->etMidRapiditySum_ = 0;
218 
219  Handle<CaloTowerCollection> towers;
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 
261  Handle<EcalRecHitCollection> ebHits;
262  Handle<EcalRecHitCollection> eeHits;
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;
294  const SiPixelRecHitCollection* rechits;
295  Handle<SiPixelRecHitCollection> rchts;
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  {
301  SiPixelRecHitCollection::DetSet hits = *it;
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 
363  Handle<TrackCollection> tracks;
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 
414  if(producePixelTracks_){
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 
430  Handle<ZDCRecHitCollection> hits;
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(std::move(creco));
462 
463 }
rpcPointValidation_cfi.recHit
recHit
Definition: rpcPointValidation_cfi.py:6
reco::CentralityProducer::trackQuality_
reco::TrackBase::TrackQuality trackQuality_
Definition: CentralityProducer.cc:99
reco::LeafCandidate::pt
virtual double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:188
nTracks
const unsigned int nTracks(const reco::Vertex &sv)
Definition: TemplatedVertexArbitrator.h:48
edmNew::DetSetVector::end
const_iterator end(bool update=false) const
Definition: DetSetVectorNew.h:553
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
mps_fire.i
i
Definition: mps_fire.py:156
reco::CentralityProducer::srcEBhits_
edm::EDGetTokenT< EcalRecHitCollection > srcEBhits_
Definition: CentralityProducer.cc:90
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
stringResolutionProvider_cfi.eta
eta
Definition: stringResolutionProvider_cfi.py:14
reco::LeafCandidate::eta
virtual double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:137
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
reco::CentralityProducer::reuseTag_
edm::EDGetTokenT< Centrality > reuseTag_
Definition: CentralityProducer.cc:96
reco::TrackBase::dxyError
double dxyError() const
error on dxy
Definition: TrackBase.h:791
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:54
HcalZDCDetId::zside
int zside() const
get the z-side of the cell (1/-1)
Definition: HcalZDCDetId.h:39
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
reco::CentralityProducer::srcVertex_
edm::EDGetTokenT< VertexCollection > srcVertex_
Definition: CentralityProducer.cc:95
allConversions_cfi.dz
dz
Definition: allConversions_cfi.py:41
reco::CentralityProducer::srcTowers_
edm::EDGetTokenT< CaloTowerCollection > srcTowers_
Definition: CentralityProducer.cc:88
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:67
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:646
CaloRecHit::energy
float energy() const
Definition: CaloRecHit.h:17
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:616
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
particleFlowDisplacedVertexCandidate_cfi.dxy
dxy
Definition: particleFlowDisplacedVertexCandidate_cfi.py:41
reco::TrackBase::ptError
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:758
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
reco::CentralityProducer::srcHFhits_
edm::EDGetTokenT< HFRecHitCollection > srcHFhits_
Definition: CentralityProducer.cc:87
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
TrackInfoProducer_cfi.rechits
rechits
Definition: TrackInfoProducer_cfi.py:9
EcalRecHit::energy
float energy() const
Definition: EcalRecHit.h:68
CaloGeometry::getPosition
const GlobalPoint & getPosition(const DetId &id) const
Get the position of a given detector id.
Definition: CaloGeometry.cc:70
BaseTrackerRecHit::localPosition
virtual LocalPoint localPosition() const final
Definition: BaseTrackerRecHit.h:76
ZDCRecHit::lowGainEnergy
float lowGainEnergy() const
Definition: ZDCRecHit.h:21
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::dzError
double dzError() const
error on dz
Definition: TrackBase.h:809
DetId
Definition: DetId.h:18
EcalRecHit::id
DetId id() const
get the id
Definition: EcalRecHit.h:77
edm::Handle::product
T const * product() const
Definition: Handle.h:81
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
math::XYZVector
XYZVectorD XYZVector
spatial vector with cartesian internal representation
Definition: Vector3D.h:30
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
reco::CentralityProducer::srcZDChits_
edm::EDGetTokenT< ZDCRecHitCollection > srcZDChits_
Definition: CentralityProducer.cc:91
SiPixelRecHit::cluster
ClusterRef cluster() const
Definition: SiPixelRecHit.h:49
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
edmNew::DetSet::detId
id_type detId() const
Definition: DetSetNew.h:84
HcalObjRepresent::isHF
bool isHF(int etabin, int depth)
Definition: HcalObjRepresent.h:397
reco::CentralityProducer::cGeo
edm::ESHandle< CaloGeometry > cGeo
Definition: CentralityProducer.cc:102
stringResolutionProvider_cfi.et
et
define resolution functions of each parameter
Definition: stringResolutionProvider_cfi.py:13
reco::TrackBase::quality
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:505
edmNew::DetSetVector::find
const_iterator find(id_type i, bool update=false) const
Definition: DetSetVectorNew.h:532
reco::CentralityProducer::srcPixelhits_
edm::EDGetTokenT< SiPixelRecHitCollection > srcPixelhits_
Definition: CentralityProducer.cc:92
CaloTower::ietaAbs
int ietaAbs() const
Definition: CaloTower.h:186
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:70
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:236
reco::CentralityProducer::srcPixelTracks_
edm::EDGetTokenT< TrackCollection > srcPixelTracks_
Definition: CentralityProducer.cc:94
reco::CentralityProducer::srcTracks_
edm::EDGetTokenT< TrackCollection > srcTracks_
Definition: CentralityProducer.cc:93
ZDCRecHit::id
HcalZDCDetId id() const
get the id
Definition: ZDCRecHit.h:19
reco::CentralityProducer::tGeo
edm::ESHandle< TrackerGeometry > tGeo
Definition: CentralityProducer.cc:101
HFRecHit::id
HcalDetId id() const
Definition: HFRecHit.h:23
TrackingRecHit::geographicalId
DetId geographicalId() const
Definition: TrackingRecHit.h:116
reco::CentralityProducer::srcEEhits_
edm::EDGetTokenT< EcalRecHitCollection > srcEEhits_
Definition: CentralityProducer.cc:89
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
eostools.move
def move(src, dest)
Definition: eostools.py:510
edmNew::DetSetVector::begin
const_iterator begin(bool update=false) const
Definition: DetSetVectorNew.h:547
TrackerGeometry::idToDet
const TrackerGeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:214
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
edmNew::DetSet::begin
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.

Generated for CMSSW Reference Manual by doxygen 1.8.11