CMS 3D CMS Logo

List of all members | Public Member Functions | Protected Member Functions | Protected Attributes
PFEcalBarrelRecHitCreator Class Reference

#include <PFEcalBarrelRecHitCreator.h>

Inheritance diagram for PFEcalBarrelRecHitCreator:
PFRecHitCreatorBase

Public Member Functions

void importRecHits (std::unique_ptr< reco::PFRecHitCollection > &out, std::unique_ptr< reco::PFRecHitCollection > &cleaned, const edm::Event &iEvent, const edm::EventSetup &iSetup) override
 
void init (const edm::EventSetup &es) override
 
 PFEcalBarrelRecHitCreator (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
 
- Public Member Functions inherited from PFRecHitCreatorBase
 PFRecHitCreatorBase ()
 
 PFRecHitCreatorBase (const edm::ParameterSet &iConfig, edm::ConsumesCollector &iC)
 
virtual ~PFRecHitCreatorBase ()=default
 

Protected Member Functions

bool isHighInterest (const EBDetId &detid)
 
EcalTrigTowerDetId readOutUnitOf (const EBDetId &detid) const
 
- Protected Member Functions inherited from PFRecHitCreatorBase
void beginEvent (const edm::Event &event, const edm::EventSetup &setup)
 

Protected Attributes

edm::EDGetTokenT< EcalRecHitCollectionrecHitToken_
 
edm::Handle< EBSrFlagCollectionsrFlagHandle_
 
edm::EDGetTokenT< EBSrFlagCollectionsrFlagToken_
 
const EcalTrigTowerConstituentsMaptriggerTowerMap_
 
- Protected Attributes inherited from PFRecHitCreatorBase
std::vector< std::unique_ptr< PFRecHitQTestBase > > qualityTests_
 

Detailed Description

Definition at line 28 of file PFEcalBarrelRecHitCreator.h.

Constructor & Destructor Documentation

PFEcalBarrelRecHitCreator::PFEcalBarrelRecHitCreator ( const edm::ParameterSet iConfig,
edm::ConsumesCollector iC 
)
inline

Definition at line 30 of file PFEcalBarrelRecHitCreator.h.

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), recHitToken_, srFlagToken_, and triggerTowerMap_.

31  : PFRecHitCreatorBase(iConfig, iC) {
33  auto srF = iConfig.getParameter<edm::InputTag>("srFlags");
34  if (not srF.label().empty())
36  triggerTowerMap_ = nullptr;
37  }
edm::EDGetTokenT< EBSrFlagCollection > srFlagToken_
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
T getParameter(std::string const &) const
edm::EDGetTokenT< EcalRecHitCollection > recHitToken_
const EcalTrigTowerConstituentsMap * triggerTowerMap_

Member Function Documentation

void PFEcalBarrelRecHitCreator::importRecHits ( std::unique_ptr< reco::PFRecHitCollection > &  out,
std::unique_ptr< reco::PFRecHitCollection > &  cleaned,
const edm::Event iEvent,
const edm::EventSetup iSetup 
)
inlineoverridevirtual

Implements PFRecHitCreatorBase.

Definition at line 39 of file PFEcalBarrelRecHitCreator.h.

References PFRecHitCreatorBase::beginEvent(), DetId::Ecal, PFLayer::ECAL_BARREL, EcalBarrel, HCALHighEnergyHPDFilter_cfi::energy, Exception, edm::EventSetup::get(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), CaloGeometry::getSubdetectorGeometry(), isHighInterest(), keep, eostools::move(), PFRecHitCreatorBase::qualityTests_, DetId::rawId(), recHitToken_, srFlagHandle_, srFlagToken_, and ntuplemaker::time.

42  {
43  beginEvent(iEvent, iSetup);
44 
46 
48  iSetup.get<CaloGeometryRecord>().get(geoHandle);
49 
50  bool useSrF = false;
51  if (not srFlagToken_.isUninitialized()) {
53  useSrF = true;
54  }
55 
56  // get the ecal geometry
58 
59  const EcalBarrelGeometry* ecalGeo = dynamic_cast<const EcalBarrelGeometry*>(gTmp);
60 
61  iEvent.getByToken(recHitToken_, recHitHandle);
62  for (const auto& erh : *recHitHandle) {
63  const DetId& detid = erh.detid();
64  auto energy = erh.energy();
65  auto time = erh.time();
66  bool hi = (useSrF ? isHighInterest(detid) : true);
67 
68  const auto thisCell = ecalGeo->getGeometry(detid);
69 
70  // find rechit geometry
71  if (!thisCell) {
72  throw cms::Exception("PFEcalBarrelRecHitCreator") << "detid " << detid.rawId() << "not found in geometry";
73  }
74 
75  out->emplace_back(thisCell, detid.rawId(), PFLayer::ECAL_BARREL, energy);
76 
77  auto& rh = out->back();
78 
79  bool rcleaned = false;
80  bool keep = true;
81 
82  //Apply Q tests
83  for (const auto& qtest : qualityTests_) {
84  if (!qtest->test(rh, erh, rcleaned, hi)) {
85  keep = false;
86  }
87  }
88 
89  if (keep) {
90  rh.setTime(time);
91  rh.setDepth(1);
92  } else {
93  if (rcleaned)
94  cleaned->push_back(std::move(out->back()));
95  out->pop_back();
96  }
97  }
98  }
edm::EDGetTokenT< EBSrFlagCollection > srFlagToken_
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:34
std::vector< std::unique_ptr< PFRecHitQTestBase > > qualityTests_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:525
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
edm::EDGetTokenT< EcalRecHitCollection > recHitToken_
const int keep
void beginEvent(const edm::Event &event, const edm::EventSetup &setup)
Definition: DetId.h:17
edm::Handle< EBSrFlagCollection > srFlagHandle_
virtual std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
T get() const
Definition: EventSetup.h:73
bool isHighInterest(const EBDetId &detid)
def move(src, dest)
Definition: eostools.py:511
void PFEcalBarrelRecHitCreator::init ( const edm::EventSetup es)
inlineoverridevirtual

Reimplemented from PFRecHitCreatorBase.

Definition at line 100 of file PFEcalBarrelRecHitCreator.h.

References edm::EventSetup::get(), edm::ESHandle< T >::product(), and triggerTowerMap_.

100  {
102  es.get<IdealGeometryRecord>().get(hTriggerTowerMap);
103  triggerTowerMap_ = hTriggerTowerMap.product();
104  }
const EcalTrigTowerConstituentsMap * triggerTowerMap_
T get() const
Definition: EventSetup.h:73
T const * product() const
Definition: ESHandle.h:86
bool PFEcalBarrelRecHitCreator::isHighInterest ( const EBDetId detid)
inlineprotected

Definition at line 107 of file PFEcalBarrelRecHitCreator.h.

References readOutUnitOf(), mps_fire::result, EcalSrFlag::SRF_FULL, and srFlagHandle_.

Referenced by importRecHits().

107  {
108  bool result = false;
109  auto srf = srFlagHandle_->find(readOutUnitOf(detid));
110  if (srf == srFlagHandle_->end())
111  return false;
112  else
113  result = ((srf->value() & ~EcalSrFlag::SRF_FORCED_MASK) == EcalSrFlag::SRF_FULL);
114  return result;
115  }
static const int SRF_FULL
Definition: EcalSrFlag.h:24
edm::Handle< EBSrFlagCollection > srFlagHandle_
EcalTrigTowerDetId readOutUnitOf(const EBDetId &detid) const
EcalTrigTowerDetId PFEcalBarrelRecHitCreator::readOutUnitOf ( const EBDetId detid) const
inlineprotected

Definition at line 117 of file PFEcalBarrelRecHitCreator.h.

References EcalTrigTowerConstituentsMap::towerOf(), and triggerTowerMap_.

Referenced by isHighInterest().

117 { return triggerTowerMap_->towerOf(detid); }
EcalTrigTowerDetId towerOf(const DetId &id) const
Get the tower id for this det id (or null if not known)
const EcalTrigTowerConstituentsMap * triggerTowerMap_

Member Data Documentation

edm::EDGetTokenT<EcalRecHitCollection> PFEcalBarrelRecHitCreator::recHitToken_
protected

Definition at line 119 of file PFEcalBarrelRecHitCreator.h.

Referenced by importRecHits(), and PFEcalBarrelRecHitCreator().

edm::Handle<EBSrFlagCollection> PFEcalBarrelRecHitCreator::srFlagHandle_
protected

Definition at line 125 of file PFEcalBarrelRecHitCreator.h.

Referenced by importRecHits(), and isHighInterest().

edm::EDGetTokenT<EBSrFlagCollection> PFEcalBarrelRecHitCreator::srFlagToken_
protected

Definition at line 120 of file PFEcalBarrelRecHitCreator.h.

Referenced by importRecHits(), and PFEcalBarrelRecHitCreator().

const EcalTrigTowerConstituentsMap* PFEcalBarrelRecHitCreator::triggerTowerMap_
protected

Definition at line 123 of file PFEcalBarrelRecHitCreator.h.

Referenced by init(), PFEcalBarrelRecHitCreator(), and readOutUnitOf().