CosmicClusterProducer Class Reference

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
	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 > &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<>
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 36 of file CosmicClusterProducer.cc.

Constructor & Destructor Documentation

CosmicClusterProducer::CosmicClusterProducer

(

const edm::ParameterSet &

ps

)

Definition at line 87 of file CosmicClusterProducer.cc.

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

                                                                      {
  // 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;
}

CosmicClusterProducer::~CosmicClusterProducer ( )

override

Definition at line 158 of file CosmicClusterProducer.cc.

References island_p.

{ delete island_p; }

Member Function Documentation

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 178 of file CosmicClusterProducer.cc.

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

Referenced by produce().

                                                                                          {
  // 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);
}

bool CosmicClusterProducer::counterExceeded ( ) const

inlineprivate

Definition at line 71 of file CosmicClusterProducer.cc.

References nEvt_, and nMaxPrintout_.

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

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

private

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

override

Definition at line 160 of file CosmicClusterProducer.cc.

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

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

Member Data Documentation

std::string CosmicClusterProducer::barrelClusterCollection_

private

Definition at line 57 of file CosmicClusterProducer.cc.

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

std::string CosmicClusterProducer::barrelClusterShapeAssociation_

private

Definition at line 64 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer(), and produce().

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

private

Definition at line 55 of file CosmicClusterProducer.cc.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

std::string CosmicClusterProducer::clustershapecollectionEB_

private

Definition at line 60 of file CosmicClusterProducer.cc.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

std::string CosmicClusterProducer::clustershapecollectionEE_

private

Definition at line 61 of file CosmicClusterProducer.cc.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

edm::EDGetTokenT<EcalRecHitCollection> CosmicClusterProducer::ebHitsToken_

private

Definition at line 50 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer().

edm::EDGetTokenT<EcalUncalibratedRecHitCollection> CosmicClusterProducer::ebUHitsToken_

private

Definition at line 53 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer().

edm::EDGetTokenT<EcalRecHitCollection> CosmicClusterProducer::eeHitsToken_

private

Definition at line 51 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer(), and produce().

edm::EDGetTokenT<EcalUncalibratedRecHitCollection> CosmicClusterProducer::eeUHitsToken_

private

Definition at line 54 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer(), and produce().

std::string CosmicClusterProducer::endcapClusterCollection_

private

Definition at line 58 of file CosmicClusterProducer.cc.

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

std::string CosmicClusterProducer::endcapClusterShapeAssociation_

private

Definition at line 65 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer(), and produce().

CosmicClusterAlgo* CosmicClusterProducer::island_p

private

Definition at line 69 of file CosmicClusterProducer.cc.

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

int CosmicClusterProducer::nEvt_

private

Definition at line 46 of file CosmicClusterProducer.cc.

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

int CosmicClusterProducer::nMaxPrintout_

private

Definition at line 45 of file CosmicClusterProducer.cc.

Referenced by counterExceeded().

PositionCalc CosmicClusterProducer::posCalculator_

private

Definition at line 67 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer().

ClusterShapeAlgo CosmicClusterProducer::shapeAlgo_

private

Definition at line 68 of file CosmicClusterProducer.cc.

Referenced by clusterizeECALPart(), and CosmicClusterProducer().

CosmicClusterAlgo::VerbosityLevel CosmicClusterProducer::verbosity

private

Definition at line 48 of file CosmicClusterProducer.cc.

Referenced by CosmicClusterProducer().