CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes | Static Private Attributes
SuperClusterImporter Class Reference
Inheritance diagram for SuperClusterImporter:
BlockElementImporterBase

Public Member Functions

void importToBlock (const edm::Event &, ElementList &) const override
 
 SuperClusterImporter (const edm::ParameterSet &, edm::ConsumesCollector &)
 
void updateEventSetup (const edm::EventSetup &es) override
 
- Public Member Functions inherited from BlockElementImporterBase
 BlockElementImporterBase (const BlockElementImporterBase &)=delete
 
 BlockElementImporterBase (const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)
 
const std::string & name () const
 
BlockElementImporterBaseoperator= (const BlockElementImporterBase &)=delete
 
virtual ~BlockElementImporterBase ()=default
 

Private Attributes

std::unique_ptr< EgammaHadTower_hadTower
 
const double _maxHoverE
 
const double _minSCPt
 
const double _pTbyPass
 
edm::EDGetTokenT< reco::SuperClusterCollection_srcEB
 
edm::EDGetTokenT< reco::SuperClusterCollection_srcEE
 
edm::EDGetTokenT< CaloTowerCollection_srcTowers
 
bool _superClustersArePF
 

Static Private Attributes

static const math::XYZPoint _zero = math::XYZPoint(0, 0, 0)
 

Additional Inherited Members

- Public Types inherited from BlockElementImporterBase
typedef std::vector< std::unique_ptr< reco::PFBlockElement > > ElementList
 

Detailed Description

Definition at line 19 of file SuperClusterImporter.cc.

Constructor & Destructor Documentation

◆ SuperClusterImporter()

SuperClusterImporter::SuperClusterImporter ( const edm::ParameterSet conf,
edm::ConsumesCollector sumes 
)

Definition at line 40 of file SuperClusterImporter.cc.

41  : BlockElementImporterBase(conf, sumes),
44  _srcTowers(sumes.consumes<CaloTowerCollection>(conf.getParameter<edm::InputTag>("source_towers"))),
45  _maxHoverE(conf.getParameter<double>("maximumHoverE")),
46  _pTbyPass(conf.getParameter<double>("minPTforBypass")),
47  _minSCPt(conf.getParameter<double>("minSuperClusterPt")),
48  _hadTower(nullptr),
49  _superClustersArePF(conf.getParameter<bool>("superClustersArePF")) {}

Member Function Documentation

◆ importToBlock()

void SuperClusterImporter::importToBlock ( const edm::Event e,
BlockElementImporterBase::ElementList elems 
) const
overridevirtual

Implements BlockElementImporterBase.

Definition at line 55 of file SuperClusterImporter.cc.

55  {
56  auto eb_scs = e.getHandle(_srcEB);
57  auto ee_scs = e.getHandle(_srcEE);
58  auto const& towers = e.get(_srcTowers);
59  elems.reserve(elems.size() + eb_scs->size() + ee_scs->size());
60  // setup our elements so that all the SCs are grouped together
61  auto SCs_end =
62  std::partition(elems.begin(), elems.end(), [](auto const& a) { return a->type() == reco::PFBlockElement::SC; });
63  // add eb superclusters
64  auto bsc = eb_scs->cbegin();
65  auto esc = eb_scs->cend();
66  reco::PFBlockElementSuperCluster* scbe = nullptr;
68  for (auto sc = bsc; sc != esc; ++sc) {
69  scref = reco::SuperClusterRef(eb_scs, std::distance(bsc, sc));
70  PFBlockElementSCEqual myEqual(scref);
71  auto sc_elem = std::find_if(elems.begin(), SCs_end, myEqual);
72  const double scpT = ptFast(sc->energy(), sc->position(), _zero);
73  const auto towersBehindCluster = _hadTower->towersOf(*sc);
74  const double H_tower = (_hadTower->getDepth1HcalESum(towersBehindCluster, towers) +
75  _hadTower->getDepth2HcalESum(towersBehindCluster, towers));
76  const double HoverE = H_tower / sc->energy();
77  if (sc_elem == SCs_end && scpT > _minSCPt && (scpT > _pTbyPass || HoverE < _maxHoverE)) {
78  scbe = new reco::PFBlockElementSuperCluster(scref);
80  SCs_end = elems.emplace(SCs_end, scbe);
81  ++SCs_end; // point to element *after* the new one
82  }
83  } // loop on eb superclusters
84  // add ee superclusters
85  bsc = ee_scs->cbegin();
86  esc = ee_scs->cend();
87  for (auto sc = bsc; sc != esc; ++sc) {
88  scref = reco::SuperClusterRef(ee_scs, std::distance(bsc, sc));
89  PFBlockElementSCEqual myEqual(scref);
90  auto sc_elem = std::find_if(elems.begin(), SCs_end, myEqual);
91  const double scpT = ptFast(sc->energy(), sc->position(), _zero);
92  const auto towersBehindCluster = _hadTower->towersOf(*sc);
93  const double H_tower = (_hadTower->getDepth1HcalESum(towersBehindCluster, towers) +
94  _hadTower->getDepth2HcalESum(towersBehindCluster, towers));
95  const double HoverE = H_tower / sc->energy();
96  if (sc_elem == SCs_end && scpT > _minSCPt && (scpT > _pTbyPass || HoverE < _maxHoverE)) {
97  scbe = new reco::PFBlockElementSuperCluster(scref);
99  SCs_end = elems.emplace(SCs_end, scbe);
100  ++SCs_end; // point to element *after* the new one
101  }
102  } // loop on ee superclusters
103  elems.shrink_to_fit();
104 }

References _hadTower, _maxHoverE, _minSCPt, _pTbyPass, _srcEB, _srcEE, _srcTowers, _superClustersArePF, _zero, a, bsc_activity_cfg::bsc, HLT_2018_cff::distance, MillePedeFileConverter_cfg::e, ptFast(), reco::PFBlockElement::SC, reco::PFBlockElementSuperCluster::setFromPFSuperCluster(), and HLT_2018_cff::towers.

◆ updateEventSetup()

void SuperClusterImporter::updateEventSetup ( const edm::EventSetup es)
overridevirtual

Reimplemented from BlockElementImporterBase.

Definition at line 51 of file SuperClusterImporter.cc.

51  {
53 }

References _hadTower, and EgammaHadTower::SingleTower.

Member Data Documentation

◆ _hadTower

std::unique_ptr<EgammaHadTower> SuperClusterImporter::_hadTower
private

Definition at line 31 of file SuperClusterImporter.cc.

Referenced by importToBlock(), and updateEventSetup().

◆ _maxHoverE

const double SuperClusterImporter::_maxHoverE
private

Definition at line 30 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _minSCPt

const double SuperClusterImporter::_minSCPt
private

Definition at line 30 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _pTbyPass

const double SuperClusterImporter::_pTbyPass
private

Definition at line 30 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _srcEB

edm::EDGetTokenT<reco::SuperClusterCollection> SuperClusterImporter::_srcEB
private

Definition at line 28 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _srcEE

edm::EDGetTokenT<reco::SuperClusterCollection> SuperClusterImporter::_srcEE
private

Definition at line 28 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _srcTowers

edm::EDGetTokenT<CaloTowerCollection> SuperClusterImporter::_srcTowers
private

Definition at line 29 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _superClustersArePF

bool SuperClusterImporter::_superClustersArePF
private

Definition at line 32 of file SuperClusterImporter.cc.

Referenced by importToBlock().

◆ _zero

const math::XYZPoint SuperClusterImporter::_zero = math::XYZPoint(0, 0, 0)
staticprivate

Definition at line 33 of file SuperClusterImporter.cc.

Referenced by importToBlock().

EgammaHadTower::SingleTower
Definition: EgammaHadTower.h:19
reco::PFBlockElementSuperCluster
Cluster Element.
Definition: PFBlockElementSuperCluster.h:15
reco::PFBlockElement::SC
Definition: PFBlockElement.h:41
HLT_2018_cff.distance
distance
Definition: HLT_2018_cff.py:6417
edm::SortedCollection< CaloTower >
BlockElementImporterBase::BlockElementImporterBase
BlockElementImporterBase(const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)
Definition: BlockElementImporterBase.h:17
edm::Ref< SuperClusterCollection >
reco::SuperClusterCollection
std::vector< SuperCluster > SuperClusterCollection
collection of SuperCluser objectr
Definition: SuperClusterFwd.h:9
bsc_activity_cfg.bsc
bsc
Definition: bsc_activity_cfg.py:31
PFBlockElementSCEqual
Definition: PFBlockElementSCEqual.h:7
reco::SuperClusterRef
edm::Ref< SuperClusterCollection > SuperClusterRef
reference to an object in a collection of SuperCluster objects
Definition: SuperClusterFwd.h:15
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:49
SuperClusterImporter::_hadTower
std::unique_ptr< EgammaHadTower > _hadTower
Definition: SuperClusterImporter.cc:31
SuperClusterImporter::_srcEB
edm::EDGetTokenT< reco::SuperClusterCollection > _srcEB
Definition: SuperClusterImporter.cc:28
reco::PFBlockElementSuperCluster::setFromPFSuperCluster
void setFromPFSuperCluster(bool val)
Definition: PFBlockElementSuperCluster.h:61
a
double a
Definition: hdecay.h:119
SuperClusterImporter::_srcTowers
edm::EDGetTokenT< CaloTowerCollection > _srcTowers
Definition: SuperClusterImporter.cc:29
SuperClusterImporter::_minSCPt
const double _minSCPt
Definition: SuperClusterImporter.cc:30
SuperClusterImporter::_maxHoverE
const double _maxHoverE
Definition: SuperClusterImporter.cc:30
SuperClusterImporter::_srcEE
edm::EDGetTokenT< reco::SuperClusterCollection > _srcEE
Definition: SuperClusterImporter.cc:28
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
HLT_2018_cff.towers
towers
Definition: HLT_2018_cff.py:35030
SuperClusterImporter::_zero
static const math::XYZPoint _zero
Definition: SuperClusterImporter.cc:33
ptFast
double ptFast(const double energy, const math::XYZPoint &position, const math::XYZPoint &origin)
Definition: SuperClusterImporter.cc:12
SuperClusterImporter::_superClustersArePF
bool _superClustersArePF
Definition: SuperClusterImporter.cc:32
EgammaHadTower
Definition: EgammaHadTower.h:17
SuperClusterImporter::_pTbyPass
const double _pTbyPass
Definition: SuperClusterImporter.cc:30
edm::InputTag
Definition: InputTag.h:15
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37