CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
HybridClusterProducer Class Reference

#include <HybridClusterProducer.h>

Inheritance diagram for HybridClusterProducer:
edm::stream::EDProducer<>

Public Member Functions

 HybridClusterProducer (const edm::ParameterSet &ps)
 
void produce (edm::Event &, const edm::EventSetup &) override
 
 ~HybridClusterProducer () 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 Attributes

std::string basicclusterCollection_
 
edm::EDGetTokenT< EcalRecHitCollectionhitsToken_
 
HybridClusterAlgohybrid_p
 
int nEvt_
 
PositionCalc posCalculator_
 
std::string superclusterCollection_
 

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, GlobalCacheLuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCacheRunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 

Detailed Description

Definition at line 18 of file HybridClusterProducer.h.

Constructor & Destructor Documentation

HybridClusterProducer::HybridClusterProducer ( const edm::ParameterSet ps)

Definition at line 38 of file HybridClusterProducer.cc.

References basicclusterCollection_, cosmicSuperClusters_cfi::bremRecoveryPset, cosmicSuperClusters_cfi::dynamicPhiRoad, edm::ParameterSet::getParameter(), hitsToken_, hybrid_p, nEvt_, MonitorAlCaEcalPi0_cfi::posCalcParameters, posCalculator_, AlCaHLTBitMon_QueryRunRegistry::string, and superclusterCollection_.

38  {
39  basicclusterCollection_ = ps.getParameter<std::string>("basicclusterCollection");
40  superclusterCollection_ = ps.getParameter<std::string>("superclusterCollection");
41  hitsToken_ = consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("recHitsCollection"));
42 
43  //Setup for core tools objects.
45 
46  posCalculator_ = PositionCalc(posCalcParameters);
47 
48  const std::vector<std::string> flagnames = ps.getParameter<std::vector<std::string> >("RecHitFlagToBeExcluded");
49 
50  const std::vector<int> flagsexcl = StringToEnumValue<EcalRecHit::Flags>(flagnames);
51 
52  const std::vector<std::string> severitynames =
53  ps.getParameter<std::vector<std::string> >("RecHitSeverityToBeExcluded");
54 
55  const std::vector<int> severitiesexcl = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynames);
56 
57  hybrid_p = new HybridClusterAlgo(ps.getParameter<double>("HybridBarrelSeedThr"),
58  ps.getParameter<int>("step"),
59  ps.getParameter<double>("ethresh"),
60  ps.getParameter<double>("eseed"),
61  ps.getParameter<double>("xi"),
62  ps.getParameter<bool>("useEtForXi"),
63  ps.getParameter<double>("ewing"),
64  flagsexcl,
65  posCalculator_,
66  ps.getParameter<bool>("dynamicEThresh"),
67  ps.getParameter<double>("eThreshA"),
68  ps.getParameter<double>("eThreshB"),
69  severitiesexcl,
70  ps.getParameter<bool>("excludeFlagged"));
71  //bremRecoveryPset,
72 
73  // get brem recovery parameters
74  bool dynamicPhiRoad = ps.getParameter<bool>("dynamicPhiRoad");
75  if (dynamicPhiRoad) {
77  hybrid_p->setDynamicPhiRoad(bremRecoveryPset);
78  }
79 
80  produces<reco::BasicClusterCollection>(basicclusterCollection_);
81  produces<reco::SuperClusterCollection>(superclusterCollection_);
82  nEvt_ = 0;
83 }
T getParameter(std::string const &) const
edm::EDGetTokenT< EcalRecHitCollection > hitsToken_
void setDynamicPhiRoad(const edm::ParameterSet &bremRecoveryPset)
HybridClusterAlgo * hybrid_p
HybridClusterProducer::~HybridClusterProducer ( )
override

Definition at line 85 of file HybridClusterProducer.cc.

References hybrid_p.

85 { delete hybrid_p; }
HybridClusterAlgo * hybrid_p

Member Function Documentation

void HybridClusterProducer::produce ( edm::Event evt,
const edm::EventSetup es 
)
override

Definition at line 87 of file HybridClusterProducer.cc.

References basicclusterCollection_, DetId::Ecal, EcalBarrel, edm::EventSetup::get(), edm::Event::getByToken(), CaloGeometry::getSubdetectorGeometry(), hitsToken_, hybrid_p, mps_fire::i, edm::OrphanHandleBase::isValid(), edm::HandleBase::isValid(), LogTrace, HybridClusterAlgo::makeClusters(), HybridClusterAlgo::makeSuperClusters(), eostools::move(), nEvt_, edm::Handle< T >::product(), edm::ESHandle< T >::product(), edm::PtrVector< T >::push_back(), edm::Event::put(), superclusterCollection_, HLT_2018_cff::superClusters, and ecaldqm::topology().

87  {
88  // get the hit collection from the event:
90 
91  evt.getByToken(hitsToken_, rhcHandle);
92  if (!(rhcHandle.isValid())) {
93  edm::LogError("MissingProduct") << "could not get a handle on the EcalRecHitCollection!";
94  return;
95  }
96  const EcalRecHitCollection* hit_collection = rhcHandle.product();
97 
98  // get the collection geometry:
100  es.get<CaloGeometryRecord>().get(geoHandle);
101  const CaloGeometry& geometry = *geoHandle;
102  const CaloSubdetectorGeometry* geometry_p;
103  std::unique_ptr<const CaloSubdetectorTopology> topology;
104 
106  es.get<EcalSeverityLevelAlgoRcd>().get(sevLv);
107 
108  geometry_p = geometry.getSubdetectorGeometry(DetId::Ecal, EcalBarrel);
109  topology = std::make_unique<EcalBarrelTopology>(*geoHandle);
110 
111  // make the Basic clusters!
112  reco::BasicClusterCollection basicClusters;
114  hit_collection, geometry_p, basicClusters, sevLv.product(), false, std::vector<RectangularEtaPhiRegion>());
115 
116  LogTrace("EcalClusters") << "Finished clustering - BasicClusterCollection returned to producer...";
117 
118  // create a unique_ptr to a BasicClusterCollection, copy the clusters into it and put in the Event:
119  auto basicclusters_p = std::make_unique<reco::BasicClusterCollection>();
120  basicclusters_p->assign(basicClusters.begin(), basicClusters.end());
122  evt.put(std::move(basicclusters_p), basicclusterCollection_);
123 
124  //Basic clusters now in the event.
125  LogTrace("EcalClusters") << "Basic Clusters now put into event.";
126 
127  //Weird though it is, get the BasicClusters back out of the event. We need the
128  //edm::Ref to these guys to make our superclusters for Hybrid.
129 
130  if (!(bccHandle.isValid())) {
131  edm::LogError("Missing Product") << "could not get a handle on the BasicClusterCollection!";
132  return;
133  }
134 
135  reco::BasicClusterCollection clusterCollection = *bccHandle;
136 
137  LogTrace("EcalClusters") << "Got the BasicClusterCollection" << std::endl;
138 
139  reco::CaloClusterPtrVector clusterPtrVector;
140  for (unsigned int i = 0; i < clusterCollection.size(); i++) {
141  clusterPtrVector.push_back(reco::CaloClusterPtr(bccHandle, i));
142  }
143 
145  LogTrace("EcalClusters") << "Found: " << superClusters.size() << " superclusters.";
146 
147  auto superclusters_p = std::make_unique<reco::SuperClusterCollection>();
148  superclusters_p->assign(superClusters.begin(), superClusters.end());
149 
150  evt.put(std::move(superclusters_p), superclusterCollection_);
151  LogTrace("EcalClusters") << "Hybrid Clusters (Basic/Super) added to the Event! :-)";
152 
153  nEvt_++;
154 }
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:34
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:131
CaloTopology const * topology(0)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
void push_back(Ptr< T > const &iPtr)
Definition: PtrVector.h:149
edm::EDGetTokenT< EcalRecHitCollection > hitsToken_
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
reco::SuperClusterCollection makeSuperClusters(const reco::CaloClusterPtrVector &)
bool isValid() const
Definition: HandleBase.h:70
#define LogTrace(id)
T const * product() const
Definition: Handle.h:69
std::vector< BasicCluster > BasicClusterCollection
collection of BasicCluster objects
T get() const
Definition: EventSetup.h:73
HybridClusterAlgo * hybrid_p
void makeClusters(const EcalRecHitCollection *, const CaloSubdetectorGeometry *geometry, reco::BasicClusterCollection &basicClusters, const EcalSeverityLevelAlgo *sevLv, bool regional=false, const std::vector< RectangularEtaPhiRegion > &regions=std::vector< RectangularEtaPhiRegion >())
T const * product() const
Definition: ESHandle.h:86
def move(src, dest)
Definition: eostools.py:511

Member Data Documentation

std::string HybridClusterProducer::basicclusterCollection_
private

Definition at line 29 of file HybridClusterProducer.h.

Referenced by HybridClusterProducer(), and produce().

edm::EDGetTokenT<EcalRecHitCollection> HybridClusterProducer::hitsToken_
private

Definition at line 32 of file HybridClusterProducer.h.

Referenced by HybridClusterProducer(), and produce().

HybridClusterAlgo* HybridClusterProducer::hybrid_p
private

Definition at line 34 of file HybridClusterProducer.h.

Referenced by HybridClusterProducer(), produce(), and ~HybridClusterProducer().

int HybridClusterProducer::nEvt_
private

Definition at line 27 of file HybridClusterProducer.h.

Referenced by HybridClusterProducer(), and produce().

PositionCalc HybridClusterProducer::posCalculator_
private

Definition at line 35 of file HybridClusterProducer.h.

Referenced by HybridClusterProducer().

std::string HybridClusterProducer::superclusterCollection_
private

Definition at line 30 of file HybridClusterProducer.h.

Referenced by HybridClusterProducer(), and produce().