Inheritance diagram for reco::CentralityProducer:

Public Member Functions
	CentralityProducer (const edm::ParameterSet &)

	~CentralityProducer () override

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

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 ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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 &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override

void	endJob () override

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

edm::ESHandle< TrackerTopology >	topo_

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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

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)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 51 of file CentralityProducer.cc.

Constructor & Destructor Documentation

◆ CentralityProducer()

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

explicit

Definition at line 114 of file CentralityProducer.cc.

                                                                        {
     //register your products
     produceHFhits_ = iConfig.getParameter<bool>("produceHFhits");
     produceHFtowers_ = iConfig.getParameter<bool>("produceHFtowers");
     produceEcalhits_ = iConfig.getParameter<bool>("produceEcalhits");
     produceZDChits_ = iConfig.getParameter<bool>("produceZDChits");
     produceETmidRap_ = iConfig.getParameter<bool>("produceETmidRapidity");
     producePixelhits_ = iConfig.getParameter<bool>("producePixelhits");
     produceTracks_ = iConfig.getParameter<bool>("produceTracks");
     producePixelTracks_ = iConfig.getParameter<bool>("producePixelTracks");
  
     midRapidityRange_ = iConfig.getParameter<double>("midRapidityRange");
     trackPtCut_ = iConfig.getParameter<double>("trackPtCut");
     trackEtaCut_ = iConfig.getParameter<double>("trackEtaCut");
  
     hfEtaCut_ = iConfig.getParameter<double>("hfEtaCut");
  
     if (produceHFhits_)
       srcHFhits_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcHFhits"));
     if (produceHFtowers_ || produceETmidRap_)
       srcTowers_ = consumes<CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("srcTowers"));
  
     if (produceEcalhits_) {
       srcEBhits_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcEBhits"));
       srcEEhits_ = consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcEEhits"));
     }
     if (produceZDChits_) {
       srcZDChits_ = consumes<ZDCRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcZDChits"));
       lowGainZDC_ = iConfig.getParameter<bool>("lowGainZDC");
     }
     if (producePixelhits_) {
       srcPixelhits_ = consumes<SiPixelRecHitCollection>(iConfig.getParameter<edm::InputTag>("srcPixelhits"));
       doPixelCut_ = iConfig.getParameter<bool>("doPixelCut");
       srcVertex_ = consumes<VertexCollection>(iConfig.getParameter<edm::InputTag>("srcVertex"));
     }
     if (produceTracks_) {
       srcTracks_ = consumes<TrackCollection>(iConfig.getParameter<edm::InputTag>("srcTracks"));
       srcVertex_ = consumes<VertexCollection>(iConfig.getParameter<edm::InputTag>("srcVertex"));
     }
     if (producePixelTracks_)
       srcPixelTracks_ = consumes<TrackCollection>(iConfig.getParameter<edm::InputTag>("srcPixelTracks"));
  
     reuseAny_ = iConfig.getParameter<bool>("reUseCentrality");
     if (reuseAny_)
       reuseTag_ = consumes<Centrality>(iConfig.getParameter<edm::InputTag>("srcReUse"));
  
     useQuality_ = iConfig.getParameter<bool>("UseQuality");
     trackQuality_ = TrackBase::qualityByName(iConfig.getParameter<std::string>("TrackQuality"));
  
     produces<Centrality>();
   }

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

◆ ~CentralityProducer()

reco::CentralityProducer::~CentralityProducer ( )

override

Definition at line 166 of file CentralityProducer.cc.

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

Member Function Documentation

◆ beginJob()

void reco::CentralityProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 503 of file CentralityProducer.cc.

503 {}

◆ endJob()

void reco::CentralityProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 506 of file CentralityProducer.cc.

506 {}

◆ produce()

void reco::CentralityProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDProducer.

Definition at line 176 of file CentralityProducer.cc.

                                                                                 {
     if (producePixelhits_)
       iSetup.get<TrackerDigiGeometryRecord>().get(tGeo);
     if (produceEcalhits_)
       iSetup.get<CaloGeometryRecord>().get(cGeo);
     if (producePixelhits_)
       iSetup.get<TrackerTopologyRcd>().get(topo_);
  
     auto creco = std::make_unique<Centrality>();
     Handle<Centrality> inputCentrality;
  
     if (reuseAny_)
       iEvent.getByToken(reuseTag_, inputCentrality);
  
     if (produceHFhits_) {
       creco->etHFhitSumPlus_ = 0;
       creco->etHFhitSumMinus_ = 0;
  
       Handle<HFRecHitCollection> hits;
       iEvent.getByToken(srcHFhits_, hits);
       for (size_t ihit = 0; ihit < hits->size(); ++ihit) {
         const HFRecHit& rechit = (*hits)[ihit];
         if (rechit.id().ieta() > 0)
           creco->etHFhitSumPlus_ += rechit.energy();
         if (rechit.id().ieta() < 0)
           creco->etHFhitSumMinus_ += rechit.energy();
       }
     } else {
       if (reuseAny_) {
         creco->etHFhitSumMinus_ = inputCentrality->EtHFhitSumMinus();
         creco->etHFhitSumPlus_ = inputCentrality->EtHFhitSumPlus();
       }
     }
  
     if (produceHFtowers_ || produceETmidRap_) {
       creco->etHFtowerSumPlus_ = 0;
       creco->etHFtowerSumMinus_ = 0;
       creco->etHFtowerSumECutPlus_ = 0;
       creco->etHFtowerSumECutMinus_ = 0;
       creco->etMidRapiditySum_ = 0;
  
       Handle<CaloTowerCollection> towers;
  
       iEvent.getByToken(srcTowers_, towers);
  
       for (size_t i = 0; i < towers->size(); ++i) {
         const CaloTower& tower = (*towers)[i];
         double eta = tower.eta();
         bool isHF = tower.ietaAbs() > 29;
         if (produceHFtowers_) {
           if (isHF && eta > 0) {
             creco->etHFtowerSumPlus_ += tower.pt();
             if (tower.energy() > 1.5)
               creco->etHFtowerSumECutPlus_ += tower.pt();
             if (eta > hfEtaCut_)
               creco->etHFtruncatedPlus_ += tower.pt();
           }
           if (isHF && eta < 0) {
             creco->etHFtowerSumMinus_ += tower.pt();
             if (tower.energy() > 1.5)
               creco->etHFtowerSumECutMinus_ += tower.pt();
             if (eta < -hfEtaCut_)
               creco->etHFtruncatedMinus_ += tower.pt();
           }
         } else {
           if (reuseAny_) {
             creco->etHFtowerSumMinus_ = inputCentrality->EtHFtowerSumMinus();
             creco->etHFtowerSumPlus_ = inputCentrality->EtHFtowerSumPlus();
             creco->etHFtowerSumECutMinus_ = inputCentrality->EtHFtowerSumECutMinus();
             creco->etHFtowerSumECutPlus_ = inputCentrality->EtHFtowerSumECutPlus();
             creco->etHFtruncatedMinus_ = inputCentrality->EtHFtruncatedMinus();
             creco->etHFtruncatedPlus_ = inputCentrality->EtHFtruncatedPlus();
           }
         }
         if (produceETmidRap_) {
           if (std::abs(eta) < midRapidityRange_)
             creco->etMidRapiditySum_ += tower.pt() / (midRapidityRange_ * 2.);
         } else if (reuseAny_)
           creco->etMidRapiditySum_ = inputCentrality->EtMidRapiditySum();
       }
     } else {
       if (reuseAny_) {
         creco->etHFtowerSumMinus_ = inputCentrality->EtHFtowerSumMinus();
         creco->etHFtowerSumPlus_ = inputCentrality->EtHFtowerSumPlus();
         creco->etHFtowerSumECutMinus_ = inputCentrality->EtHFtowerSumECutMinus();
         creco->etHFtowerSumECutPlus_ = inputCentrality->EtHFtowerSumECutPlus();
         creco->etMidRapiditySum_ = inputCentrality->EtMidRapiditySum();
       }
     }
  
     if (produceEcalhits_) {
       creco->etEESumPlus_ = 0;
       creco->etEESumMinus_ = 0;
       creco->etEBSum_ = 0;
  
       Handle<EcalRecHitCollection> ebHits;
       Handle<EcalRecHitCollection> eeHits;
  
       iEvent.getByToken(srcEBhits_, ebHits);
       iEvent.getByToken(srcEEhits_, eeHits);
  
       for (unsigned int i = 0; i < ebHits->size(); ++i) {
         const EcalRecHit& hit = (*ebHits)[i];
         const GlobalPoint& pos = cGeo->getPosition(hit.id());
         double et = hit.energy() * (pos.perp() / pos.mag());
         creco->etEBSum_ += et;
       }
  
       for (unsigned int i = 0; i < eeHits->size(); ++i) {
         const EcalRecHit& hit = (*eeHits)[i];
         const GlobalPoint& pos = cGeo->getPosition(hit.id());
         double et = hit.energy() * (pos.perp() / pos.mag());
         if (pos.z() > 0) {
           creco->etEESumPlus_ += et;
         } else {
           creco->etEESumMinus_ += et;
         }
       }
     } else {
       if (reuseAny_) {
         creco->etEESumMinus_ = inputCentrality->EtEESumMinus();
         creco->etEESumPlus_ = inputCentrality->EtEESumPlus();
         creco->etEBSum_ = inputCentrality->EtEBSum();
       }
     }
  
     if (producePixelhits_) {
       creco->pixelMultiplicity_ = 0;
       const SiPixelRecHitCollection* rechits;
       Handle<SiPixelRecHitCollection> rchts;
       iEvent.getByToken(srcPixelhits_, rchts);
       rechits = rchts.product();
       int nPixel = 0;
       int nPixel_plus = 0;
       int nPixel_minus = 0;
       for (SiPixelRecHitCollection::const_iterator it = rechits->begin(); it != rechits->end(); it++) {
         SiPixelRecHitCollection::DetSet hits = *it;
         DetId detId = DetId(hits.detId());
         SiPixelRecHitCollection::const_iterator recHitMatch = rechits->find(detId);
         const SiPixelRecHitCollection::DetSet recHitRange = *recHitMatch;
         for (SiPixelRecHitCollection::DetSet::const_iterator recHitIterator = recHitRange.begin();
              recHitIterator != recHitRange.end();
              ++recHitIterator) {
           // add selection if needed, now all hits.
           const SiPixelRecHit* recHit = &(*recHitIterator);
           const PixelGeomDetUnit* pixelLayer =
               dynamic_cast<const PixelGeomDetUnit*>(tGeo->idToDet(recHit->geographicalId()));
           GlobalPoint gpos = pixelLayer->toGlobal(recHit->localPosition());
           math::XYZVector rechitPos(gpos.x(), gpos.y(), gpos.z());
           double eta = rechitPos.eta();
           double abeta = std::abs(rechitPos.eta());
           int clusterSize = recHit->cluster()->size();
           unsigned layer = topo_->layer(detId);
           if (doPixelCut_) {
             if (detId.det() == DetId::Tracker && detId.subdetId() == PixelSubdetector::PixelBarrel) {
               if (layer == 1 && 18 * abeta - 40 > clusterSize)
                 continue;
               if (layer == 2 && 6 * abeta - 7.2 > clusterSize)
                 continue;
               if ((layer == 3 || layer == 4) && 4 * abeta - 2.4 > clusterSize)
                 continue;
             }
           }
           nPixel++;
           if (eta >= 0)
             nPixel_plus++;
           if (eta < 0)
             nPixel_minus++;
         }
       }
       creco->pixelMultiplicity_ = nPixel;
       creco->pixelMultiplicityPlus_ = nPixel_plus;
       creco->pixelMultiplicityMinus_ = nPixel_minus;
     } else {
       if (reuseAny_) {
         creco->pixelMultiplicity_ = inputCentrality->multiplicityPixel();
         creco->pixelMultiplicityPlus_ = inputCentrality->multiplicityPixelPlus();
         creco->pixelMultiplicityMinus_ = inputCentrality->multiplicityPixelMinus();
       }
     }
  
     if (produceTracks_) {
       double vx = -999.;
       double vy = -999.;
       double vz = -999.;
       double vxError = -999.;
       double vyError = -999.;
       double vzError = -999.;
       math::XYZVector vtxPos(0, 0, 0);
  
       Handle<VertexCollection> recoVertices;
       iEvent.getByToken(srcVertex_, recoVertices);
       unsigned int daughter = 0;
       int greatestvtx = 0;
  
       for (unsigned int i = 0; i < recoVertices->size(); ++i) {
         daughter = (*recoVertices)[i].tracksSize();
         if (daughter > (*recoVertices)[greatestvtx].tracksSize())
           greatestvtx = i;
       }
  
       if (!recoVertices->empty()) {
         vx = (*recoVertices)[greatestvtx].position().x();
         vy = (*recoVertices)[greatestvtx].position().y();
         vz = (*recoVertices)[greatestvtx].position().z();
         vxError = (*recoVertices)[greatestvtx].xError();
         vyError = (*recoVertices)[greatestvtx].yError();
         vzError = (*recoVertices)[greatestvtx].zError();
       }
  
       vtxPos = math::XYZVector(vx, vy, vz);
  
       Handle<TrackCollection> tracks;
       iEvent.getByToken(srcTracks_, tracks);
       int nTracks = 0;
  
       double trackCounter = 0;
       double trackCounterEta = 0;
       double trackCounterEtaPt = 0;
  
       for (unsigned int i = 0; i < tracks->size(); ++i) {
         const Track& track = (*tracks)[i];
         if (useQuality_ && !track.quality(trackQuality_))
           continue;
  
         if (track.pt() > trackPtCut_)
           trackCounter++;
         if (std::abs(track.eta()) < trackEtaCut_) {
           trackCounterEta++;
           if (track.pt() > trackPtCut_)
             trackCounterEtaPt++;
         }
  
         math::XYZPoint v1(vx, vy, vz);
         double dz = track.dz(v1);
         double dzsigma2 = track.dzError() * track.dzError() + vzError * vzError;
         double dxy = track.dxy(v1);
         double dxysigma2 = track.dxyError() * track.dxyError() + vxError * vyError;
  
         const double pterrcut = 0.1;
         const double dzrelcut = 3.0;
         const double dxyrelcut = 3.0;
  
         if (track.quality(trackQuality_) && track.pt() > 0.4 && std::abs(track.eta()) < 2.4 &&
             track.ptError() / track.pt() < pterrcut && dz * dz < dzrelcut * dzrelcut * dzsigma2 &&
             dxy * dxy < dxyrelcut * dxyrelcut * dxysigma2) {
           nTracks++;
         }
       }
  
       creco->trackMultiplicity_ = nTracks;
       creco->ntracksPtCut_ = trackCounter;
       creco->ntracksEtaCut_ = trackCounterEta;
       creco->ntracksEtaPtCut_ = trackCounterEtaPt;
  
     } else {
       if (reuseAny_) {
         creco->trackMultiplicity_ = inputCentrality->Ntracks();
         creco->ntracksPtCut_ = inputCentrality->NtracksPtCut();
         creco->ntracksEtaCut_ = inputCentrality->NtracksEtaCut();
         creco->ntracksEtaPtCut_ = inputCentrality->NtracksEtaPtCut();
       }
     }
  
     if (producePixelTracks_) {
       Handle<TrackCollection> pixeltracks;
       iEvent.getByToken(srcPixelTracks_, pixeltracks);
       int nPixelTracks = pixeltracks->size();
       int nPixelTracksPlus = 0;
       int nPixelTracksMinus = 0;
  
       for (auto const& track : *pixeltracks) {
         if (track.eta() < 0)
           nPixelTracksMinus++;
         else
           nPixelTracksPlus++;
       }
       creco->nPixelTracks_ = nPixelTracks;
       creco->nPixelTracksPlus_ = nPixelTracksPlus;
       creco->nPixelTracksMinus_ = nPixelTracksMinus;
     } else {
       if (reuseAny_) {
         creco->nPixelTracks_ = inputCentrality->NpixelTracks();
         creco->nPixelTracksPlus_ = inputCentrality->NpixelTracksPlus();
         creco->nPixelTracksMinus_ = inputCentrality->NpixelTracksMinus();
       }
     }
  
     if (produceZDChits_) {
       creco->zdcSumPlus_ = 0;
       creco->zdcSumMinus_ = 0;
  
       Handle<ZDCRecHitCollection> hits;
       bool zdcAvailable = iEvent.getByToken(srcZDChits_, hits);
       if (zdcAvailable) {
         for (size_t ihit = 0; ihit < hits->size(); ++ihit) {
           const ZDCRecHit& rechit = (*hits)[ihit];
           if (rechit.id().zside() > 0) {
             if (lowGainZDC_) {
               creco->zdcSumPlus_ += rechit.lowGainEnergy();
             } else {
               creco->zdcSumPlus_ += rechit.energy();
             }
           }
           if (rechit.id().zside() < 0) {
             if (lowGainZDC_) {
               creco->zdcSumMinus_ += rechit.lowGainEnergy();
             } else {
               creco->zdcSumMinus_ += rechit.energy();
             }
           }
         }
       } else {
         creco->zdcSumPlus_ = -9;
         creco->zdcSumMinus_ = -9;
       }
     } else {
       if (reuseAny_) {
         creco->zdcSumMinus_ = inputCentrality->zdcSumMinus();
         creco->zdcSumPlus_ = inputCentrality->zdcSumPlus();
       }
     }
  
     iEvent.put(std::move(creco));
   }

Member Data Documentation

◆ cGeo

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

private

Definition at line 99 of file CentralityProducer.cc.

◆ doPixelCut_

bool reco::CentralityProducer::doPixelCut_

private

Definition at line 75 of file CentralityProducer.cc.

◆ hfEtaCut_

double reco::CentralityProducer::hfEtaCut_

private

Definition at line 80 of file CentralityProducer.cc.

◆ lowGainZDC_

bool reco::CentralityProducer::lowGainZDC_

private

Definition at line 82 of file CentralityProducer.cc.

◆ midRapidityRange_

double reco::CentralityProducer::midRapidityRange_

private

Definition at line 77 of file CentralityProducer.cc.

◆ produceEcalhits_

bool reco::CentralityProducer::produceEcalhits_

private

Definition at line 67 of file CentralityProducer.cc.

◆ produceETmidRap_

bool reco::CentralityProducer::produceETmidRap_

private

Definition at line 69 of file CentralityProducer.cc.

◆ produceHFhits_

bool reco::CentralityProducer::produceHFhits_

private

Definition at line 65 of file CentralityProducer.cc.

◆ produceHFtowers_

bool reco::CentralityProducer::produceHFtowers_

private

Definition at line 66 of file CentralityProducer.cc.

◆ producePixelhits_

bool reco::CentralityProducer::producePixelhits_

private

Definition at line 70 of file CentralityProducer.cc.

◆ producePixelTracks_

bool reco::CentralityProducer::producePixelTracks_

private

Definition at line 73 of file CentralityProducer.cc.

◆ produceTracks_

bool reco::CentralityProducer::produceTracks_

private

Definition at line 71 of file CentralityProducer.cc.

◆ produceZDChits_

bool reco::CentralityProducer::produceZDChits_

private

Definition at line 68 of file CentralityProducer.cc.

◆ recoLevel_

bool reco::CentralityProducer::recoLevel_

private

Definition at line 63 of file CentralityProducer.cc.

◆ reuseAny_

bool reco::CentralityProducer::reuseAny_

private

Definition at line 72 of file CentralityProducer.cc.

◆ reuseTag_

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

private

Definition at line 93 of file CentralityProducer.cc.

◆ srcEBhits_

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

private

Definition at line 87 of file CentralityProducer.cc.

◆ srcEEhits_

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

private

Definition at line 86 of file CentralityProducer.cc.

◆ srcHFhits_

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

private

Definition at line 84 of file CentralityProducer.cc.

◆ srcPixelhits_

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

private

Definition at line 89 of file CentralityProducer.cc.

◆ srcPixelTracks_

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

private

Definition at line 91 of file CentralityProducer.cc.

◆ srcTowers_

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

private

Definition at line 85 of file CentralityProducer.cc.

◆ srcTracks_

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

private

Definition at line 90 of file CentralityProducer.cc.

◆ srcVertex_

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

private

Definition at line 92 of file CentralityProducer.cc.

◆ srcZDChits_

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

private

Definition at line 88 of file CentralityProducer.cc.

◆ tGeo

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

private

Definition at line 98 of file CentralityProducer.cc.

◆ topo_

edm::ESHandle<TrackerTopology> reco::CentralityProducer::topo_

private

Definition at line 100 of file CentralityProducer.cc.

◆ trackEtaCut_

double reco::CentralityProducer::trackEtaCut_

private

Definition at line 79 of file CentralityProducer.cc.

◆ trackPtCut_

double reco::CentralityProducer::trackPtCut_

private

Definition at line 78 of file CentralityProducer.cc.

◆ trackQuality_

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

private

Definition at line 96 of file CentralityProducer.cc.

◆ useQuality_

bool reco::CentralityProducer::useQuality_

private

Definition at line 95 of file CentralityProducer.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CentralityProducer()

◆ ~CentralityProducer()

Member Function Documentation

◆ beginJob()

◆ endJob()

◆ produce()

Member Data Documentation

◆ cGeo

◆ doPixelCut_

◆ hfEtaCut_

◆ lowGainZDC_

◆ midRapidityRange_

◆ produceEcalhits_

◆ produceETmidRap_

◆ produceHFhits_

◆ produceHFtowers_

◆ producePixelhits_

◆ producePixelTracks_

◆ produceTracks_

◆ produceZDChits_

◆ recoLevel_

◆ reuseAny_

◆ reuseTag_

◆ srcEBhits_

◆ srcEEhits_

◆ srcHFhits_

◆ srcPixelhits_

◆ srcPixelTracks_

◆ srcTowers_

◆ srcTracks_

◆ srcVertex_

◆ srcZDChits_

◆ tGeo

◆ topo_

◆ trackEtaCut_

◆ trackPtCut_

◆ trackQuality_

◆ useQuality_