Multi5x5ClusterProducer Class Reference

Inheritance diagram for Multi5x5ClusterProducer:

Public Member Functions
	Multi5x5ClusterProducer (const edm::ParameterSet &ps)
void	produce (edm::Event &, const edm::EventSetup &) override
	~Multi5x5ClusterProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	clusterizeECALPart (edm::Event &evt, const edm::EventSetup &es, const edm::EDGetTokenT< EcalRecHitCollection > &token, const std::string &clusterCollection, const reco::CaloID::Detectors detector)
bool	counterExceeded () const
const EcalRecHitCollection *	getCollection (edm::Event &evt, const edm::EDGetTokenT< EcalRecHitCollection > &token)
void	outputValidationInfo (reco::CaloClusterPtrVector &clusterPtrVector)

Private Attributes
std::string	barrelClusterCollection_
edm::EDGetTokenT< EcalRecHitCollection >	barrelHitToken_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
bool	doBarrel_
bool	doEndcap_
std::string	endcapClusterCollection_
edm::EDGetTokenT< EcalRecHitCollection >	endcapHitToken_
Multi5x5ClusterAlgo *	island_p
int	nEvt_
int	nMaxPrintout_
PositionCalc	posCalculator_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 33 of file Multi5x5ClusterProducer.cc.

Constructor & Destructor Documentation

Multi5x5ClusterProducer()

Multi5x5ClusterProducer::Multi5x5ClusterProducer

(

const edm::ParameterSet &

ps

)

Definition at line 75 of file Multi5x5ClusterProducer.cc.

References barrelClusterCollection_, barrelHitToken_, caloGeometryToken_, doBarrel_, doEndcap_, endcapClusterCollection_, endcapHitToken_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), island_p, nEvt_, ecalPiZeroTask_cfi::posCalcParameters, posCalculator_, multi5x5BasicClusters_cfi::reassignSeedCrysToClusterItSeeds, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                          {
  // Parameters to identify the hit collections
  barrelHitToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelHitTag"));

  endcapHitToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapHitTag"));

  //EventSetup Token for CaloGeometry
  caloGeometryToken_ = esConsumes<CaloGeometry, CaloGeometryRecord>();

  // should cluster algo be run in barrel and endcap?
  doEndcap_ = ps.getParameter<bool>("doEndcap");
  doBarrel_ = ps.getParameter<bool>("doBarrel");

  // The names of the produced cluster collections
  barrelClusterCollection_ = ps.getParameter<std::string>("barrelClusterCollection");
  endcapClusterCollection_ = ps.getParameter<std::string>("endcapClusterCollection");

  // Island algorithm parameters
  double barrelSeedThreshold = ps.getParameter<double>("IslandBarrelSeedThr");
  double endcapSeedThreshold = ps.getParameter<double>("IslandEndcapSeedThr");

  const std::vector<std::string> flagnames = ps.getParameter<std::vector<std::string> >("RecHitFlagToBeExcluded");

  const std::vector<int> v_chstatus = StringToEnumValue<EcalRecHit::Flags>(flagnames);

  // Parameters for the position calculation:
  edm::ParameterSet posCalcParameters = ps.getParameter<edm::ParameterSet>("posCalcParameters");
  posCalculator_ = PositionCalc(posCalcParameters);

  // Produces a collection of barrel and a collection of endcap clusters
  produces<reco::BasicClusterCollection>(endcapClusterCollection_);
  produces<reco::BasicClusterCollection>(barrelClusterCollection_);

  bool reassignSeedCrysToClusterItSeeds = false;
  if (ps.exists("reassignSeedCrysToClusterItSeeds"))
    reassignSeedCrysToClusterItSeeds = ps.getParameter<bool>("reassignSeedCrysToClusterItSeeds");

  island_p = new Multi5x5ClusterAlgo(
      barrelSeedThreshold, endcapSeedThreshold, v_chstatus, posCalculator_, reassignSeedCrysToClusterItSeeds);

  nEvt_ = 0;
}

~Multi5x5ClusterProducer()

Multi5x5ClusterProducer::~Multi5x5ClusterProducer ( )

override

Definition at line 118 of file Multi5x5ClusterProducer.cc.

References island_p.

{ delete island_p; }

Member Function Documentation

clusterizeECALPart()

void Multi5x5ClusterProducer::clusterizeECALPart ( edm::Event &  evt, const edm::EventSetup &  es, const edm::EDGetTokenT< EcalRecHitCollection > &  token, const std::string &  clusterCollection, const reco::CaloID::Detectors  detector  )

private

Definition at line 138 of file Multi5x5ClusterProducer.cc.

References barrelClusterCollection_, caloGeometryToken_, bsc_activity_cfg::clusters, reco::CaloID::DET_ECAL_BARREL, hgcalTestNeighbor_cfi::detector, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, endcapClusterCollection_, getCollection(), edm::EventSetup::getHandle(), CaloGeometry::getSubdetectorGeometry(), island_p, Multi5x5ClusterAlgo::makeClusters(), eostools::move(), edm::Event::put(), and unpackBuffers-CaloStage2::token.

Referenced by produce().

                                                                                     {
  // get the hit collection from the event:
  const EcalRecHitCollection* hitCollection_p = getCollection(evt, token);

  // get the geometry and topology from the event setup:
  edm::ESHandle<CaloGeometry> geoHandle = es.getHandle(caloGeometryToken_);

  const CaloSubdetectorGeometry* geometry_p;
  std::unique_ptr<CaloSubdetectorTopology> topology_p;

  if (detector == reco::CaloID::DET_ECAL_BARREL) {
    geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
    topology_p = std::make_unique<EcalBarrelTopology>(*geoHandle);
  } else {
    geometry_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalEndcap);
    topology_p = std::make_unique<EcalEndcapTopology>(*geoHandle);
  }

  const CaloSubdetectorGeometry* geometryES_p;
  geometryES_p = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);

  // Run the clusterization algorithm:
  reco::BasicClusterCollection clusters;
  clusters = island_p->makeClusters(hitCollection_p, geometry_p, topology_p.get(), geometryES_p, detector);

  // create a unique_ptr to a BasicClusterCollection, copy the barrel clusters into it and put in the Event:
  auto clusters_p = std::make_unique<reco::BasicClusterCollection>();
  clusters_p->assign(clusters.begin(), clusters.end());
  edm::OrphanHandle<reco::BasicClusterCollection> bccHandle;
  if (detector == reco::CaloID::DET_ECAL_BARREL)
    bccHandle = evt.put(std::move(clusters_p), barrelClusterCollection_);
  else
    bccHandle = evt.put(std::move(clusters_p), endcapClusterCollection_);
}

counterExceeded()

bool Multi5x5ClusterProducer::counterExceeded ( ) const

inlineprivate

Definition at line 59 of file Multi5x5ClusterProducer.cc.

References nEvt_, and nMaxPrintout_.

{ return ((nEvt_ > nMaxPrintout_) || (nMaxPrintout_ < 0)); }

getCollection()

const EcalRecHitCollection * Multi5x5ClusterProducer::getCollection ( edm::Event &  evt, const edm::EDGetTokenT< EcalRecHitCollection > &  token  )

private

Definition at line 131 of file Multi5x5ClusterProducer.cc.

References edm::Event::getByToken(), edm::Handle< T >::product(), and unpackBuffers-CaloStage2::token.

Referenced by clusterizeECALPart().

                                                                      {
  edm::Handle<EcalRecHitCollection> rhcHandle;
  evt.getByToken(token, rhcHandle);
  return rhcHandle.product();
}

outputValidationInfo()

void Multi5x5ClusterProducer::outputValidationInfo ( reco::CaloClusterPtrVector &  clusterPtrVector )

private

produce()

void Multi5x5ClusterProducer::produce ( edm::Event &  evt, const edm::EventSetup &  es  )

override

Definition at line 120 of file Multi5x5ClusterProducer.cc.

References barrelClusterCollection_, barrelHitToken_, clusterizeECALPart(), reco::CaloID::DET_ECAL_BARREL, reco::CaloID::DET_ECAL_ENDCAP, doBarrel_, doEndcap_, endcapClusterCollection_, endcapHitToken_, and nEvt_.

                                                                            {
  if (doEndcap_) {
    clusterizeECALPart(evt, es, endcapHitToken_, endcapClusterCollection_, reco::CaloID::DET_ECAL_ENDCAP);
  }
  if (doBarrel_) {
    clusterizeECALPart(evt, es, barrelHitToken_, barrelClusterCollection_, reco::CaloID::DET_ECAL_BARREL);
  }

  nEvt_++;
}

Member Data Documentation

barrelClusterCollection_

std::string Multi5x5ClusterProducer::barrelClusterCollection_

private

Definition at line 53 of file Multi5x5ClusterProducer.cc.

Referenced by clusterizeECALPart(), Multi5x5ClusterProducer(), and produce().

barrelHitToken_

edm::EDGetTokenT<EcalRecHitCollection> Multi5x5ClusterProducer::barrelHitToken_

private

Definition at line 49 of file Multi5x5ClusterProducer.cc.

Referenced by Multi5x5ClusterProducer(), and produce().

caloGeometryToken_

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> Multi5x5ClusterProducer::caloGeometryToken_

private

Definition at line 51 of file Multi5x5ClusterProducer.cc.

Referenced by clusterizeECALPart(), and Multi5x5ClusterProducer().

doBarrel_

bool Multi5x5ClusterProducer::doBarrel_

private

Definition at line 46 of file Multi5x5ClusterProducer.cc.

Referenced by Multi5x5ClusterProducer(), and produce().

doEndcap_

bool Multi5x5ClusterProducer::doEndcap_

private

Definition at line 47 of file Multi5x5ClusterProducer.cc.

Referenced by Multi5x5ClusterProducer(), and produce().

endcapClusterCollection_

std::string Multi5x5ClusterProducer::endcapClusterCollection_

private

Definition at line 54 of file Multi5x5ClusterProducer.cc.

Referenced by clusterizeECALPart(), Multi5x5ClusterProducer(), and produce().

endcapHitToken_

edm::EDGetTokenT<EcalRecHitCollection> Multi5x5ClusterProducer::endcapHitToken_

private

Definition at line 50 of file Multi5x5ClusterProducer.cc.

Referenced by Multi5x5ClusterProducer(), and produce().

island_p

Multi5x5ClusterAlgo* Multi5x5ClusterProducer::island_p

private

Definition at line 57 of file Multi5x5ClusterProducer.cc.

Referenced by clusterizeECALPart(), Multi5x5ClusterProducer(), and ~Multi5x5ClusterProducer().

nEvt_

int Multi5x5ClusterProducer::nEvt_

private

Definition at line 43 of file Multi5x5ClusterProducer.cc.

Referenced by counterExceeded(), Multi5x5ClusterProducer(), and produce().

nMaxPrintout_

int Multi5x5ClusterProducer::nMaxPrintout_

private

Definition at line 42 of file Multi5x5ClusterProducer.cc.

Referenced by counterExceeded().

posCalculator_

PositionCalc Multi5x5ClusterProducer::posCalculator_

private

Definition at line 56 of file Multi5x5ClusterProducer.cc.

Referenced by Multi5x5ClusterProducer().