CosmicClusterProducer Class Reference

#include <CosmicClusterProducer.h>

Inheritance diagram for CosmicClusterProducer:

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

Private Member Functions
void	clusterizeECALPart (edm::Event &evt, const edm::EventSetup &es, const edm::EDGetTokenT< EcalRecHitCollection > &hitsToken, const edm::EDGetTokenT< EcalUncalibratedRecHitCollection > &uhitsToken, const std::string &clusterCollection, const std::string &clusterShapeAssociation, const CosmicClusterAlgo::EcalPart &ecalPart)
bool	counterExceeded () const
void	outputValidationInfo (reco::CaloClusterPtrVector &clusterPtrVector)

Private Attributes
std::string	barrelClusterCollection_
std::string	barrelClusterShapeAssociation_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
std::string	clustershapecollectionEB_
std::string	clustershapecollectionEE_
edm::EDGetTokenT< EcalRecHitCollection >	ebHitsToken_
edm::EDGetTokenT< EcalUncalibratedRecHitCollection >	ebUHitsToken_
edm::EDGetTokenT< EcalRecHitCollection >	eeHitsToken_
edm::EDGetTokenT< EcalUncalibratedRecHitCollection >	eeUHitsToken_
std::string	endcapClusterCollection_
std::string	endcapClusterShapeAssociation_
CosmicClusterAlgo *	island_p
int	nEvt_
int	nMaxPrintout_
PositionCalc	posCalculator_
ClusterShapeAlgo	shapeAlgo_
CosmicClusterAlgo::VerbosityLevel	verbosity

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

Detailed Description

Definition at line 29 of file CosmicClusterProducer.h.

Constructor & Destructor Documentation

CosmicClusterProducer()

CosmicClusterProducer::CosmicClusterProducer

(

const edm::ParameterSet &

ps

)

Definition at line 37 of file CosmicClusterProducer.cc.

                                                                      {
  // The verbosity level
  std::string verbosityString = ps.getParameter<std::string>("VerbosityLevel");
  if (verbosityString == "DEBUG")
    verbosity = CosmicClusterAlgo::pDEBUG;
  else if (verbosityString == "WARNING")
    verbosity = CosmicClusterAlgo::pWARNING;
  else if (verbosityString == "INFO")
    verbosity = CosmicClusterAlgo::pINFO;
  else
    verbosity = CosmicClusterAlgo::pERROR;
 
  // Parameters to identify the hit collections
  ebHitsToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("barrelHits"));
  eeHitsToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("endcapHits"));
 
  ebUHitsToken_ = consumes<EcalUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("barrelUncalibHits"));
  eeUHitsToken_ = consumes<EcalUncalibratedRecHitCollection>(ps.getParameter<edm::InputTag>("endcapUncalibHits"));
 
  caloGeometryToken_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
 
  // 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>("BarrelSeedThr");
  double barrelSingleThreshold = ps.getParameter<double>("BarrelSingleThr");
  double barrelSecondThreshold = ps.getParameter<double>("BarrelSecondThr");
  double barrelSupThreshold = ps.getParameter<double>("BarrelSupThr");
  double endcapSeedThreshold = ps.getParameter<double>("EndcapSeedThr");
  double endcapSingleThreshold = ps.getParameter<double>("EndcapSingleThr");
  double endcapSecondThreshold = ps.getParameter<double>("EndcapSecondThr");
  double endcapSupThreshold = ps.getParameter<double>("EndcapSupThr");
 
  // Parameters for the position calculation:
  edm::ParameterSet posCalcParameters = ps.getParameter<edm::ParameterSet>("posCalcParameters");
 
  posCalculator_ = PositionCalc(posCalcParameters);
  shapeAlgo_ = ClusterShapeAlgo(posCalcParameters);
 
  clustershapecollectionEB_ = ps.getParameter<std::string>("clustershapecollectionEB");
  clustershapecollectionEE_ = ps.getParameter<std::string>("clustershapecollectionEE");
 
  //AssociationMap
  barrelClusterShapeAssociation_ = ps.getParameter<std::string>("barrelShapeAssociation");
  endcapClusterShapeAssociation_ = ps.getParameter<std::string>("endcapShapeAssociation");
 
  // Produces a collection of barrel and a collection of endcap clusters
 
  produces<reco::ClusterShapeCollection>(clustershapecollectionEE_);
  produces<reco::BasicClusterCollection>(endcapClusterCollection_);
  produces<reco::ClusterShapeCollection>(clustershapecollectionEB_);
  produces<reco::BasicClusterCollection>(barrelClusterCollection_);
  produces<reco::BasicClusterShapeAssociationCollection>(barrelClusterShapeAssociation_);
  produces<reco::BasicClusterShapeAssociationCollection>(endcapClusterShapeAssociation_);
 
  island_p = new CosmicClusterAlgo(barrelSeedThreshold,
                                   barrelSingleThreshold,
                                   barrelSecondThreshold,
                                   barrelSupThreshold,
                                   endcapSeedThreshold,
                                   endcapSingleThreshold,
                                   endcapSecondThreshold,
                                   endcapSupThreshold,
                                   posCalculator_,
                                   verbosity);
 
  nEvt_ = 0;
}

References barrelClusterCollection_, barrelClusterShapeAssociation_, caloGeometryToken_, clustershapecollectionEB_, clustershapecollectionEE_, ebHitsToken_, ebUHitsToken_, eeHitsToken_, eeUHitsToken_, endcapClusterCollection_, endcapClusterShapeAssociation_, edm::ParameterSet::getParameter(), island_p, nEvt_, CosmicClusterAlgo::pDEBUG, CosmicClusterAlgo::pERROR, CosmicClusterAlgo::pINFO, MonitorAlCaEcalPi0_cfi::posCalcParameters, posCalculator_, CosmicClusterAlgo::pWARNING, shapeAlgo_, AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

~CosmicClusterProducer()

CosmicClusterProducer::~CosmicClusterProducer ( )

override

Definition at line 108 of file CosmicClusterProducer.cc.

{ delete island_p; }

References island_p.

Member Function Documentation

clusterizeECALPart()

void CosmicClusterProducer::clusterizeECALPart ( edm::Event &  evt, const edm::EventSetup &  es, const edm::EDGetTokenT< EcalRecHitCollection > &  hitsToken, const edm::EDGetTokenT< EcalUncalibratedRecHitCollection > &  uhitsToken, const std::string &  clusterCollection, const std::string &  clusterShapeAssociation, const CosmicClusterAlgo::EcalPart &  ecalPart  )

private

Definition at line 128 of file CosmicClusterProducer.cc.

                                                                                          {
  // get the hit collection from the event:
 
  edm::Handle<EcalRecHitCollection> hits_h;
  edm::Handle<EcalUncalibratedRecHitCollection> uhits_h;
 
  evt.getByToken(hitsToken, hits_h);
  evt.getByToken(uhitsToken, uhits_h);
 
  const EcalRecHitCollection* hitCollection_p = hits_h.product();
  const EcalUncalibratedRecHitCollection* uhitCollection_p = uhits_h.product();
 
  // 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;
 
  std::string clustershapetag;
  if (ecalPart == CosmicClusterAlgo::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, uhitCollection_p, geometry_p, topology_p.get(), geometryES_p, ecalPart);
 
  //Create associated ClusterShape objects.
  std::vector<reco::ClusterShape> ClusVec;
 
  for (int erg = 0; erg < int(clusters.size()); ++erg) {
    reco::ClusterShape TestShape = shapeAlgo_.Calculate(clusters[erg], hitCollection_p, geometry_p, topology_p.get());
    ClusVec.push_back(TestShape);
  }
 
  //Put clustershapes in event, but retain a Handle on them.
  auto clustersshapes_p = std::make_unique<reco::ClusterShapeCollection>();
  clustersshapes_p->assign(ClusVec.begin(), ClusVec.end());
  edm::OrphanHandle<reco::ClusterShapeCollection> clusHandle;
  if (ecalPart == CosmicClusterAlgo::barrel)
    clusHandle = evt.put(std::move(clustersshapes_p), clustershapecollectionEB_);
  else
    clusHandle = evt.put(std::move(clustersshapes_p), clustershapecollectionEE_);
 
  // 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 (ecalPart == CosmicClusterAlgo::barrel)
    bccHandle = evt.put(std::move(clusters_p), barrelClusterCollection_);
  else
    bccHandle = evt.put(std::move(clusters_p), endcapClusterCollection_);
 
  // BasicClusterShapeAssociationMap
  auto shapeAssocs_p = std::make_unique<reco::BasicClusterShapeAssociationCollection>(bccHandle, clusHandle);
 
  for (unsigned int i = 0; i < clusHandle->size(); i++) {
    shapeAssocs_p->insert(edm::Ref<reco::BasicClusterCollection>(bccHandle, i),
                          edm::Ref<reco::ClusterShapeCollection>(clusHandle, i));
  }
  evt.put(std::move(shapeAssocs_p), clusterShapeAssociation);
}

References CosmicClusterAlgo::barrel, barrelClusterCollection_, ClusterShapeAlgo::Calculate(), caloGeometryToken_, bsc_activity_cfg::clusters, clustershapecollectionEB_, clustershapecollectionEE_, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, endcapClusterCollection_, edm::Event::getByToken(), edm::EventSetup::getHandle(), CaloGeometry::getSubdetectorGeometry(), mps_fire::i, createfilelist::int, island_p, CosmicClusterAlgo::makeClusters(), eostools::move(), edm::Handle< T >::product(), edm::Event::put(), shapeAlgo_, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

counterExceeded()

bool CosmicClusterProducer::counterExceeded ( ) const

inlineprivate

Definition at line 64 of file CosmicClusterProducer.h.

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

References nEvt_, and nMaxPrintout_.

outputValidationInfo()

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

private

produce()

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

override

Definition at line 110 of file CosmicClusterProducer.cc.

                                                                          {
  clusterizeECALPart(evt,
                     es,
                     eeHitsToken_,
                     eeUHitsToken_,
                     endcapClusterCollection_,
                     endcapClusterShapeAssociation_,
                     CosmicClusterAlgo::endcap);
  clusterizeECALPart(evt,
                     es,
                     eeHitsToken_,
                     eeUHitsToken_,
                     barrelClusterCollection_,
                     barrelClusterShapeAssociation_,
                     CosmicClusterAlgo::barrel);
  nEvt_++;
}

References CosmicClusterAlgo::barrel, barrelClusterCollection_, barrelClusterShapeAssociation_, clusterizeECALPart(), eeHitsToken_, eeUHitsToken_, CosmicClusterAlgo::endcap, endcapClusterCollection_, endcapClusterShapeAssociation_, and nEvt_.

Member Data Documentation

barrelClusterCollection_

std::string CosmicClusterProducer::barrelClusterCollection_

private

Definition at line 50 of file CosmicClusterProducer.h.

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

barrelClusterShapeAssociation_

std::string CosmicClusterProducer::barrelClusterShapeAssociation_

private

Definition at line 57 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

caloGeometryToken_

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

private

Definition at line 48 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

clustershapecollectionEB_

std::string CosmicClusterProducer::clustershapecollectionEB_

private

Definition at line 53 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

clustershapecollectionEE_

std::string CosmicClusterProducer::clustershapecollectionEE_

private

Definition at line 54 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

ebHitsToken_

edm::EDGetTokenT<EcalRecHitCollection> CosmicClusterProducer::ebHitsToken_

private

Definition at line 43 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().

ebUHitsToken_

edm::EDGetTokenT<EcalUncalibratedRecHitCollection> CosmicClusterProducer::ebUHitsToken_

private

Definition at line 46 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().

eeHitsToken_

edm::EDGetTokenT<EcalRecHitCollection> CosmicClusterProducer::eeHitsToken_

private

Definition at line 44 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

eeUHitsToken_

edm::EDGetTokenT<EcalUncalibratedRecHitCollection> CosmicClusterProducer::eeUHitsToken_

private

Definition at line 47 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

endcapClusterCollection_

std::string CosmicClusterProducer::endcapClusterCollection_

private

Definition at line 51 of file CosmicClusterProducer.h.

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

endcapClusterShapeAssociation_

std::string CosmicClusterProducer::endcapClusterShapeAssociation_

private

Definition at line 58 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

island_p

CosmicClusterAlgo* CosmicClusterProducer::island_p

private

Definition at line 62 of file CosmicClusterProducer.h.

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

nEvt_

int CosmicClusterProducer::nEvt_

private

Definition at line 39 of file CosmicClusterProducer.h.

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

nMaxPrintout_

int CosmicClusterProducer::nMaxPrintout_

private

Definition at line 38 of file CosmicClusterProducer.h.

Referenced by counterExceeded().

posCalculator_

PositionCalc CosmicClusterProducer::posCalculator_

private

Definition at line 60 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().

shapeAlgo_

ClusterShapeAlgo CosmicClusterProducer::shapeAlgo_

private

Definition at line 61 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

verbosity

CosmicClusterAlgo::VerbosityLevel CosmicClusterProducer::verbosity

private

Definition at line 41 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().