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	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () 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_
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 26 of file CosmicClusterProducer.h.

Constructor & Destructor Documentation

CosmicClusterProducer()

CosmicClusterProducer::CosmicClusterProducer

(

const edm::ParameterSet &

ps

)

Definition at line 39 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"));
 
  // 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_, 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.get<CaloGeometryRecord>().get(geoHandle);
 
  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(), bsc_activity_cfg::clusters, clustershapecollectionEB_, clustershapecollectionEE_, DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, endcapClusterCollection_, edm::EventSetup::get(), get, edm::Event::getByToken(), 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 60 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 46 of file CosmicClusterProducer.h.

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

barrelClusterShapeAssociation_

std::string CosmicClusterProducer::barrelClusterShapeAssociation_

private

Definition at line 53 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

clustershapecollectionEB_

std::string CosmicClusterProducer::clustershapecollectionEB_

private

Definition at line 49 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

clustershapecollectionEE_

std::string CosmicClusterProducer::clustershapecollectionEE_

private

Definition at line 50 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

ebHitsToken_

edm::EDGetTokenT<EcalRecHitCollection> CosmicClusterProducer::ebHitsToken_

private

Definition at line 40 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().

ebUHitsToken_

edm::EDGetTokenT<EcalUncalibratedRecHitCollection> CosmicClusterProducer::ebUHitsToken_

private

Definition at line 43 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().

eeHitsToken_

edm::EDGetTokenT<EcalRecHitCollection> CosmicClusterProducer::eeHitsToken_

private

Definition at line 41 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

eeUHitsToken_

edm::EDGetTokenT<EcalUncalibratedRecHitCollection> CosmicClusterProducer::eeUHitsToken_

private

Definition at line 44 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

endcapClusterCollection_

std::string CosmicClusterProducer::endcapClusterCollection_

private

Definition at line 47 of file CosmicClusterProducer.h.

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

endcapClusterShapeAssociation_

std::string CosmicClusterProducer::endcapClusterShapeAssociation_

private

Definition at line 54 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer(), and produce().

island_p

CosmicClusterAlgo* CosmicClusterProducer::island_p

private

Definition at line 58 of file CosmicClusterProducer.h.

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

nEvt_

int CosmicClusterProducer::nEvt_

private

Definition at line 36 of file CosmicClusterProducer.h.

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

nMaxPrintout_

int CosmicClusterProducer::nMaxPrintout_

private

Definition at line 35 of file CosmicClusterProducer.h.

Referenced by counterExceeded().

posCalculator_

PositionCalc CosmicClusterProducer::posCalculator_

private

Definition at line 56 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().

shapeAlgo_

ClusterShapeAlgo CosmicClusterProducer::shapeAlgo_

private

Definition at line 57 of file CosmicClusterProducer.h.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

verbosity

CosmicClusterAlgo::VerbosityLevel CosmicClusterProducer::verbosity

private

Definition at line 38 of file CosmicClusterProducer.h.

Referenced by CosmicClusterProducer().