CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Public Attributes | Static Private Member Functions | Private Attributes
CaloTowersCreator Class Reference

#include <CaloTowersCreator.h>

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

Public Member Functions

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

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Public Attributes

double EBEScale
 
double EEEScale
 
double HBEScale
 
double HEDEScale
 
double HESEScale
 
double HF1EScale
 
double HF2EScale
 
double HOEScale
 

Static Private Member Functions

static const std::vector< double > & getGridValues ()
 

Private Attributes

CaloTowersCreationAlgo algo_
 
bool allowMissingInputs_
 
edm::ESWatcher< IdealGeometryRecord > caloTowerConstituentsWatcher_
 
std::vector< edm::InputTag > ecalLabels_
 
edm::ESWatcher< EcalSeverityLevelAlgoRcd > ecalSevLevelWatcher_
 
EScales eScales_
 
edm::ESWatcher< HcalChannelQualityRcd > hcalChStatusWatcher_
 
edm::ESWatcher< HcalSeverityLevelComputerRcd > hcalSevLevelWatcher_
 
std::vector< int > theEcalSeveritiesToBeExcluded_
 
std::vector< int > theEcalSeveritiesToBeUsedInBadTowers_
 
unsigned int theHcalAcceptSeverityLevel_
 
unsigned int theHcalAcceptSeverityLevelForRejectedHit_
 
bool theRecoveredEcalHitsAreUsed_
 
bool theRecoveredHcalHitsAreUsed_
 
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
 
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
 
edm::EDGetTokenT< HORecHitCollection > tok_ho_
 
std::vector< edm::EDGetTokenT< EcalRecHitCollection > > toks_ecal_
 
bool useRejectedHitsOnly_
 
bool useRejectedRecoveredEcalHits_
 
bool useRejectedRecoveredHcalHits_
 

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

Original author: J. Mans - Minnesota

Definition at line 25 of file CaloTowersCreator.h.

Constructor & Destructor Documentation

◆ CaloTowersCreator()

CaloTowersCreator::CaloTowersCreator ( const edm::ParameterSet &  ps)
explicit

Definition at line 13 of file CaloTowersCreator.cc.

14  : algo_(conf.getParameter<double>("EBThreshold"),
15  conf.getParameter<double>("EEThreshold"),
16 
17  conf.getParameter<bool>("UseEtEBTreshold"),
18  conf.getParameter<bool>("UseEtEETreshold"),
19  conf.getParameter<bool>("UseSymEBTreshold"),
20  conf.getParameter<bool>("UseSymEETreshold"),
21 
22  conf.getParameter<double>("HcalThreshold"),
23  conf.getParameter<double>("HBThreshold"),
24  conf.getParameter<double>("HBThreshold1"),
25  conf.getParameter<double>("HBThreshold2"),
26  conf.getParameter<double>("HESThreshold"),
27  conf.getParameter<double>("HESThreshold1"),
28  conf.getParameter<double>("HEDThreshold"),
29  conf.getParameter<double>("HEDThreshold1"),
30  conf.getParameter<double>("HOThreshold0"),
31  conf.getParameter<double>("HOThresholdPlus1"),
32  conf.getParameter<double>("HOThresholdMinus1"),
33  conf.getParameter<double>("HOThresholdPlus2"),
34  conf.getParameter<double>("HOThresholdMinus2"),
35  conf.getParameter<double>("HF1Threshold"),
36  conf.getParameter<double>("HF2Threshold"),
37  conf.getParameter<std::vector<double> >("EBGrid"),
38  conf.getParameter<std::vector<double> >("EBWeights"),
39  conf.getParameter<std::vector<double> >("EEGrid"),
40  conf.getParameter<std::vector<double> >("EEWeights"),
41  conf.getParameter<std::vector<double> >("HBGrid"),
42  conf.getParameter<std::vector<double> >("HBWeights"),
43  conf.getParameter<std::vector<double> >("HESGrid"),
44  conf.getParameter<std::vector<double> >("HESWeights"),
45  conf.getParameter<std::vector<double> >("HEDGrid"),
46  conf.getParameter<std::vector<double> >("HEDWeights"),
47  conf.getParameter<std::vector<double> >("HOGrid"),
48  conf.getParameter<std::vector<double> >("HOWeights"),
49  conf.getParameter<std::vector<double> >("HF1Grid"),
50  conf.getParameter<std::vector<double> >("HF1Weights"),
51  conf.getParameter<std::vector<double> >("HF2Grid"),
52  conf.getParameter<std::vector<double> >("HF2Weights"),
53  conf.getParameter<double>("EBWeight"),
54  conf.getParameter<double>("EEWeight"),
55  conf.getParameter<double>("HBWeight"),
56  conf.getParameter<double>("HESWeight"),
57  conf.getParameter<double>("HEDWeight"),
58  conf.getParameter<double>("HOWeight"),
59  conf.getParameter<double>("HF1Weight"),
60  conf.getParameter<double>("HF2Weight"),
61  conf.getParameter<double>("EcutTower"),
62  conf.getParameter<double>("EBSumThreshold"),
63  conf.getParameter<double>("EESumThreshold"),
64  conf.getParameter<bool>("UseHO"),
65  // (for momentum reconstruction algorithm)
66  conf.getParameter<int>("MomConstrMethod"),
67  conf.getParameter<double>("MomHBDepth"),
68  conf.getParameter<double>("MomHEDepth"),
69  conf.getParameter<double>("MomEBDepth"),
70  conf.getParameter<double>("MomEEDepth"),
71  conf.getParameter<int>("HcalPhase")),
72 
73  ecalLabels_(conf.getParameter<std::vector<edm::InputTag> >("ecalInputs")),
74  allowMissingInputs_(conf.getParameter<bool>("AllowMissingInputs")),
75 
76  theHcalAcceptSeverityLevel_(conf.getParameter<unsigned int>("HcalAcceptSeverityLevel")),
77 
78  theRecoveredHcalHitsAreUsed_(conf.getParameter<bool>("UseHcalRecoveredHits")),
79  theRecoveredEcalHitsAreUsed_(conf.getParameter<bool>("UseEcalRecoveredHits")),
80 
81  // paramaters controlling the use of rejected hits
82 
83  useRejectedHitsOnly_(conf.getParameter<bool>("UseRejectedHitsOnly")),
84 
85  theHcalAcceptSeverityLevelForRejectedHit_(
86  conf.getParameter<unsigned int>("HcalAcceptSeverityLevelForRejectedHit")),
87 
88  useRejectedRecoveredHcalHits_(conf.getParameter<bool>("UseRejectedRecoveredHcalHits")),
89  useRejectedRecoveredEcalHits_(conf.getParameter<bool>("UseRejectedRecoveredEcalHits"))
90 
91 {
92  algo_.setMissingHcalRescaleFactorForEcal(conf.getParameter<double>("missingHcalRescaleFactorForEcal"));
93 
94  // register for data access
95  tok_hbhe_ = consumes<HBHERecHitCollection>(conf.getParameter<edm::InputTag>("hbheInput"));
96  tok_ho_ = consumes<HORecHitCollection>(conf.getParameter<edm::InputTag>("hoInput"));
97  tok_hf_ = consumes<HFRecHitCollection>(conf.getParameter<edm::InputTag>("hfInput"));
98 
99  const unsigned nLabels = ecalLabels_.size();
100  for (unsigned i = 0; i != nLabels; i++)
101  toks_ecal_.push_back(consumes<EcalRecHitCollection>(ecalLabels_[i]));
102 
103  EBEScale = eScales_.EBScale;
104  EEEScale = eScales_.EEScale;
105  HBEScale = eScales_.HBScale;
106  HESEScale = eScales_.HESScale;
107  HEDEScale = eScales_.HEDScale;
108  HOEScale = eScales_.HOScale;
109  HF1EScale = eScales_.HF1Scale;
110  HF2EScale = eScales_.HF2Scale;
111 
112  // get the Ecal severities to be excluded
113  const std::vector<std::string> severitynames =
114  conf.getParameter<std::vector<std::string> >("EcalRecHitSeveritiesToBeExcluded");
115 
116  theEcalSeveritiesToBeExcluded_ = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynames);
117 
118  // get the Ecal severities to be used for bad towers
119  theEcalSeveritiesToBeUsedInBadTowers_ = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(
120  conf.getParameter<std::vector<std::string> >("EcalSeveritiesToBeUsedInBadTowers"));
121 
122  if (eScales_.instanceLabel.empty())
123  produces<CaloTowerCollection>();
124  else
125  produces<CaloTowerCollection>(eScales_.instanceLabel);
126 
127 #ifdef EDM_ML_DEBUG
128  std::cout << "VI Producer " << (useRejectedHitsOnly_ ? "use rejectOnly " : " ")
129  << (allowMissingInputs_ ? "allowMissing " : " ") << nLabels << ' ' << severitynames.size() << std::endl;
130 #endif
131 }

References algo_, allowMissingInputs_, gather_cfg::cout, EBEScale, EScales::EBScale, ecalLabels_, EEEScale, EScales::EEScale, eScales_, edm::ParameterSet::getParameter(), HBEScale, EScales::HBScale, HEDEScale, EScales::HEDScale, HESEScale, EScales::HESScale, HF1EScale, EScales::HF1Scale, HF2EScale, EScales::HF2Scale, HOEScale, EScales::HOScale, mps_fire::i, EScales::instanceLabel, CaloTowersCreationAlgo::setMissingHcalRescaleFactorForEcal(), theEcalSeveritiesToBeExcluded_, theEcalSeveritiesToBeUsedInBadTowers_, tok_hbhe_, tok_hf_, tok_ho_, toks_ecal_, and useRejectedHitsOnly_.

◆ ~CaloTowersCreator()

CaloTowersCreator::~CaloTowersCreator ( )
inlineoverride

Definition at line 28 of file CaloTowersCreator.h.

28 {}

Member Function Documentation

◆ fillDescriptions()

void CaloTowersCreator::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 303 of file CaloTowersCreator.cc.

303  {
304  edm::ParameterSetDescription desc;
305  desc.add<double>("EBSumThreshold", 0.2);
306  desc.add<double>("HF2Weight", 1.0);
307  desc.add<double>("EBWeight", 1.0);
308  desc.add<double>("EESumThreshold", 0.45);
309  desc.add<double>("HOThreshold0", 1.1);
310  desc.add<double>("HOThresholdPlus1", 3.5);
311  desc.add<double>("HOThresholdMinus1", 3.5);
312  desc.add<double>("HOThresholdPlus2", 3.5);
313  desc.add<double>("HOThresholdMinus2", 3.5);
314  desc.add<double>("HBThreshold", 0.7);
315  desc.add<double>("HBThreshold1", 0.7);
316  desc.add<double>("HBThreshold2", 0.7);
317  desc.add<double>("HF1Threshold", 0.5);
318  desc.add<double>("HEDWeight", 1.0);
319  desc.add<double>("EEWeight", 1.0);
320  desc.add<double>("HESWeight", 1.0);
321  desc.add<double>("HF1Weight", 1.0);
322  desc.add<double>("HOWeight", 1.0);
323  desc.add<double>("EBThreshold", 0.07);
324  desc.add<double>("EEThreshold", 0.3);
325  desc.add<double>("HcalThreshold", -1000.0);
326  desc.add<double>("HF2Threshold", 0.85);
327  desc.add<double>("HESThreshold", 0.8);
328  desc.add<double>("HESThreshold1", 0.8);
329  desc.add<double>("HEDThreshold", 0.8);
330  desc.add<double>("HEDThreshold1", 0.8);
331  desc.add<double>("EcutTower", -1000.0);
332  desc.add<double>("HBWeight", 1.0);
333  desc.add<double>("MomHBDepth", 0.2);
334  desc.add<double>("MomHEDepth", 0.4);
335  desc.add<double>("MomEBDepth", 0.3);
336  desc.add<double>("MomEEDepth", 0.0);
337  desc.add<bool>("UseHO", true);
338  desc.add<bool>("UseEtEBTreshold", false);
339  desc.add<bool>("UseSymEBTreshold", true);
340  desc.add<bool>("UseEtEETreshold", false);
341  desc.add<bool>("UseSymEETreshold", true);
342  desc.add<bool>("UseHcalRecoveredHits", true);
343  desc.add<bool>("UseEcalRecoveredHits", false);
344  desc.add<bool>("UseRejectedHitsOnly", false);
345  desc.add<bool>("UseRejectedRecoveredHcalHits", true);
346  desc.add<bool>("UseRejectedRecoveredEcalHits", false);
347  desc.add<double>("missingHcalRescaleFactorForEcal", 0.0);
348  desc.add<bool>("AllowMissingInputs", false);
349  desc.add<std::vector<double> >("HBGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
350  desc.add<std::vector<double> >("EEWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
351  desc.add<std::vector<double> >("HF2Weights", {1.0, 1.0, 1.0, 1.0, 1.0});
352  desc.add<std::vector<double> >("HOWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
353  desc.add<std::vector<double> >("EEGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
354  desc.add<std::vector<double> >("HBWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
355  desc.add<std::vector<double> >("HF2Grid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
356  desc.add<std::vector<double> >("HEDWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
357  desc.add<std::vector<double> >("HF1Grid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
358  desc.add<std::vector<double> >("EBWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
359  desc.add<std::vector<double> >("HF1Weights", {1.0, 1.0, 1.0, 1.0, 1.0});
360  desc.add<std::vector<double> >("HESGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
361  desc.add<std::vector<double> >("HESWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
362  desc.add<std::vector<double> >("HEDGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
363  desc.add<std::vector<double> >("HOGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
364  desc.add<std::vector<double> >("EBGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
365  desc.add<edm::InputTag>("hfInput", edm::InputTag("hfreco"));
366  desc.add<edm::InputTag>("hbheInput", edm::InputTag("hbhereco"));
367  desc.add<edm::InputTag>("hoInput", edm::InputTag("horeco"));
368  desc.add<std::vector<edm::InputTag> >(
369  "ecalInputs", {edm::InputTag("ecalRecHit", "EcalRecHitsEB"), edm::InputTag("ecalRecHit", "EcalRecHitsEE")});
370  desc.add<int>("MomConstrMethod", 1);
371  desc.add<unsigned int>("HcalAcceptSeverityLevel", 9);
372  desc.add<std::vector<std::string> >("EcalRecHitSeveritiesToBeExcluded", {"kTime", "kWeird", "kBad"});
373  desc.add<unsigned int>("HcalAcceptSeverityLevelForRejectedHit", 9999);
374  desc.add<std::vector<std::string> >("EcalSeveritiesToBeUsedInBadTowers", {});
375  desc.add<int>("HcalPhase", 0);
376 
377  descriptions.addDefault(desc);
378 }

References edm::ParameterSetDescription::add(), edm::ConfigurationDescriptions::addDefault(), and HLT_2018_cff::InputTag.

◆ getGridValues()

static const std::vector<double>& CaloTowersCreator::getGridValues ( )
staticprivate

◆ produce()

void CaloTowersCreator::produce ( edm::Event &  e,
const edm::EventSetup &  c 
)
override

Definition at line 133 of file CaloTowersCreator.cc.

133  {
134  // get the necessary event setup objects...
135  edm::ESHandle<CaloGeometry> pG;
136  edm::ESHandle<HcalTopology> htopo;
137  edm::ESHandle<CaloTowerTopology> cttopo;
138  edm::ESHandle<CaloTowerConstituentsMap> ctmap;
139  c.get<CaloGeometryRecord>().get(pG);
140  c.get<HcalRecNumberingRecord>().get(htopo);
141  c.get<HcalRecNumberingRecord>().get(cttopo);
142  c.get<CaloGeometryRecord>().get(ctmap);
143 
144  // ECAL channel status map ****************************************
145  edm::ESHandle<EcalChannelStatus> ecalChStatus;
146  c.get<EcalChannelStatusRcd>().get(ecalChStatus);
147  const EcalChannelStatus* dbEcalChStatus = ecalChStatus.product();
148 
149  // HCAL channel status map ****************************************
150  edm::ESHandle<HcalChannelQuality> hcalChStatus;
151  c.get<HcalChannelQualityRcd>().get("withTopo", hcalChStatus);
152 
153  const HcalChannelQuality* dbHcalChStatus = hcalChStatus.product();
154 
155  // Assignment of severity levels **********************************
156  edm::ESHandle<HcalSeverityLevelComputer> hcalSevLvlComputerHndl;
157  c.get<HcalSeverityLevelComputerRcd>().get(hcalSevLvlComputerHndl);
158  const HcalSeverityLevelComputer* hcalSevLvlComputer = hcalSevLvlComputerHndl.product();
159 
160  edm::ESHandle<EcalSeverityLevelAlgo> ecalSevLvlAlgoHndl;
161  c.get<EcalSeverityLevelAlgoRcd>().get(ecalSevLvlAlgoHndl);
162  const EcalSeverityLevelAlgo* ecalSevLvlAlgo = ecalSevLvlAlgoHndl.product();
163 
164  algo_.setEBEScale(EBEScale);
165  algo_.setEEEScale(EEEScale);
166  algo_.setHBEScale(HBEScale);
167  algo_.setHESEScale(HESEScale);
168  algo_.setHEDEScale(HEDEScale);
169  algo_.setHOEScale(HOEScale);
170  algo_.setHF1EScale(HF1EScale);
171  algo_.setHF2EScale(HF2EScale);
172  algo_.setGeometry(cttopo.product(), ctmap.product(), htopo.product(), pG.product());
173 
174  // for treatment of problematic and anomalous cells
175 
176  algo_.setHcalChStatusFromDB(dbHcalChStatus);
177  algo_.setEcalChStatusFromDB(dbEcalChStatus);
178 
179  algo_.setHcalAcceptSeverityLevel(theHcalAcceptSeverityLevel_);
180  algo_.setEcalSeveritiesToBeExcluded(theEcalSeveritiesToBeExcluded_);
181 
182  algo_.setRecoveredHcalHitsAreUsed(theRecoveredHcalHitsAreUsed_);
183  algo_.setRecoveredEcalHitsAreUsed(theRecoveredEcalHitsAreUsed_);
184 
185  algo_.setHcalSevLvlComputer(hcalSevLvlComputer);
186  algo_.setEcalSevLvlAlgo(ecalSevLvlAlgo);
187 
188  algo_.setUseRejectedHitsOnly(useRejectedHitsOnly_);
189 
190  algo_.setHcalAcceptSeverityLevelForRejectedHit(theHcalAcceptSeverityLevelForRejectedHit_);
191  algo_.SetEcalSeveritiesToBeUsedInBadTowers(theEcalSeveritiesToBeUsedInBadTowers_);
192 
193  algo_.setUseRejectedRecoveredHcalHits(useRejectedRecoveredHcalHits_);
194  algo_.setUseRejectedRecoveredEcalHits(useRejectedRecoveredEcalHits_);
195 
196 #ifdef EDM_ML_DEBUG
197  std::cout << "VI Produce: " << (useRejectedHitsOnly_ ? "use rejectOnly " : " ")
198  << (allowMissingInputs_ ? "allowMissing " : " ")
199  << (theRecoveredEcalHitsAreUsed_ ? "use RecoveredEcal " : " ") << toks_ecal_.size() << ' '
200  << theEcalSeveritiesToBeExcluded_.size() << ' ' << theEcalSeveritiesToBeUsedInBadTowers_.size()
201  << std::endl;
202 #endif
203 
204  algo_.begin(); // clear the internal buffer
205 
206  // can't chain these in a big OR statement, or else it'll
207  // get triggered for each of the first three events
208  bool check1 = hcalSevLevelWatcher_.check(c);
209  bool check2 = hcalChStatusWatcher_.check(c);
210  bool check3 = caloTowerConstituentsWatcher_.check(c);
211  if (check1 || check2 || check3) {
212  algo_.makeHcalDropChMap();
213  }
214 
215  // check ecal SevLev
216  if (ecalSevLevelWatcher_.check(c))
217  algo_.makeEcalBadChs();
218 
219  // ----------------------------------------------------------
220  // For ecal error handling need to
221  // have access to the EB and EE collections at the end of
222  // tower reconstruction.
223 
224  edm::Handle<EcalRecHitCollection> ebHandle;
225  edm::Handle<EcalRecHitCollection> eeHandle;
226 
227  for (std::vector<edm::EDGetTokenT<EcalRecHitCollection> >::const_iterator i = toks_ecal_.begin();
228  i != toks_ecal_.end();
229  i++) {
230  edm::Handle<EcalRecHitCollection> ec_tmp;
231 
232  if (!e.getByToken(*i, ec_tmp))
233  continue;
234  if (ec_tmp->empty())
235  continue;
236 
237  // check if this is EB or EE
238  if ((ec_tmp->begin()->detid()).subdetId() == EcalBarrel) {
239  ebHandle = ec_tmp;
240  } else if ((ec_tmp->begin()->detid()).subdetId() == EcalEndcap) {
241  eeHandle = ec_tmp;
242  }
243  }
244 
245  algo_.setEbHandle(ebHandle);
246  algo_.setEeHandle(eeHandle);
247 
248  //-----------------------------------------------------------
249 
250  bool present;
251 
252  // Step A/C: Get Inputs and process (repeatedly)
253  edm::Handle<HBHERecHitCollection> hbhe;
254  present = e.getByToken(tok_hbhe_, hbhe);
255  if (present || !allowMissingInputs_)
256  algo_.process(*hbhe);
257 
258  edm::Handle<HORecHitCollection> ho;
259  present = e.getByToken(tok_ho_, ho);
260  if (present || !allowMissingInputs_)
261  algo_.process(*ho);
262 
263  edm::Handle<HFRecHitCollection> hf;
264  present = e.getByToken(tok_hf_, hf);
265  if (present || !allowMissingInputs_)
266  algo_.process(*hf);
267 
268  std::vector<edm::EDGetTokenT<EcalRecHitCollection> >::const_iterator i;
269  for (i = toks_ecal_.begin(); i != toks_ecal_.end(); i++) {
270  edm::Handle<EcalRecHitCollection> ec;
271  present = e.getByToken(*i, ec);
272  if (present || !allowMissingInputs_)
273  algo_.process(*ec);
274  }
275 
276  // Step B: Create empty output
277  auto prod = std::make_unique<CaloTowerCollection>();
278 
279  // Step C: Process
280  algo_.finish(*prod);
281 
282 #ifdef EDM_ML_DEBUG
283  int totc = 0;
284  float totE = 0;
285  reco::LeafCandidate::LorentzVector totP4;
286  for (auto const& tw : (*prod)) {
287  totc += tw.constituents().size();
288  totE += tw.energy();
289  totP4 += tw.p4();
290  std::cout << "CaloTowerCreator: " << tw.id() << " with E " << tw.energy() << " and " << tw.constituents().size()
291  << " constituents\n";
292  }
293  std::cout << "VI " << (*prod).size() << " " << totc << " " << totE << " " << totP4 << std::endl;
294 #endif
295 
296  // Step D: Put into the event
297  if (eScales_.instanceLabel.empty())
298  e.put(std::move(prod));
299  else
300  e.put(std::move(prod), eScales_.instanceLabel);
301 }

References algo_, allowMissingInputs_, edm::SortedCollection< T, SORT >::begin(), CaloTowersCreationAlgo::begin(), HltBtagPostValidation_cff::c, caloTowerConstituentsWatcher_, edm::ESWatcher< T >::check(), gather_cfg::cout, MillePedeFileConverter_cfg::e, EBEScale, EcalBarrel, EcalEndcap, ecalSevLevelWatcher_, EEEScale, edm::SortedCollection< T, SORT >::empty(), eScales_, CaloTowersCreationAlgo::finish(), get, HBEScale, photonIsolationHIProducer_cfi::hbhe, hcalChStatusWatcher_, hcalSevLevelWatcher_, HEDEScale, HESEScale, es_hardcode_cfi::hf, HF1EScale, HF2EScale, es_hardcode_cfi::ho, HOEScale, mps_fire::i, EScales::instanceLabel, CaloTowersCreationAlgo::makeEcalBadChs(), CaloTowersCreationAlgo::makeHcalDropChMap(), eostools::move(), CaloTowersCreationAlgo::process(), dumpMFGeometry_cfg::prod, edm::ESHandle< T >::product(), CaloTowersCreationAlgo::setEBEScale(), CaloTowersCreationAlgo::setEbHandle(), CaloTowersCreationAlgo::setEcalChStatusFromDB(), CaloTowersCreationAlgo::setEcalSeveritiesToBeExcluded(), CaloTowersCreationAlgo::SetEcalSeveritiesToBeUsedInBadTowers(), CaloTowersCreationAlgo::setEcalSevLvlAlgo(), CaloTowersCreationAlgo::setEEEScale(), CaloTowersCreationAlgo::setEeHandle(), CaloTowersCreationAlgo::setGeometry(), CaloTowersCreationAlgo::setHBEScale(), CaloTowersCreationAlgo::setHcalAcceptSeverityLevel(), CaloTowersCreationAlgo::setHcalAcceptSeverityLevelForRejectedHit(), CaloTowersCreationAlgo::setHcalChStatusFromDB(), CaloTowersCreationAlgo::setHcalSevLvlComputer(), CaloTowersCreationAlgo::setHEDEScale(), CaloTowersCreationAlgo::setHESEScale(), CaloTowersCreationAlgo::setHF1EScale(), CaloTowersCreationAlgo::setHF2EScale(), CaloTowersCreationAlgo::setHOEScale(), CaloTowersCreationAlgo::setRecoveredEcalHitsAreUsed(), CaloTowersCreationAlgo::setRecoveredHcalHitsAreUsed(), CaloTowersCreationAlgo::setUseRejectedHitsOnly(), CaloTowersCreationAlgo::setUseRejectedRecoveredEcalHits(), CaloTowersCreationAlgo::setUseRejectedRecoveredHcalHits(), theEcalSeveritiesToBeExcluded_, theEcalSeveritiesToBeUsedInBadTowers_, theHcalAcceptSeverityLevel_, theHcalAcceptSeverityLevelForRejectedHit_, theRecoveredEcalHitsAreUsed_, theRecoveredHcalHitsAreUsed_, tok_hbhe_, tok_hf_, tok_ho_, toks_ecal_, useRejectedHitsOnly_, useRejectedRecoveredEcalHits_, and useRejectedRecoveredHcalHits_.

Member Data Documentation

◆ algo_

CaloTowersCreationAlgo CaloTowersCreator::algo_
private

Definition at line 37 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ allowMissingInputs_

bool CaloTowersCreator::allowMissingInputs_
private

Definition at line 43 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ caloTowerConstituentsWatcher_

edm::ESWatcher<IdealGeometryRecord> CaloTowersCreator::caloTowerConstituentsWatcher_
private

Definition at line 71 of file CaloTowersCreator.h.

Referenced by produce().

◆ EBEScale

double CaloTowersCreator::EBEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ ecalLabels_

std::vector<edm::InputTag> CaloTowersCreator::ecalLabels_
private

Definition at line 41 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator().

◆ ecalSevLevelWatcher_

edm::ESWatcher<EcalSeverityLevelAlgoRcd> CaloTowersCreator::ecalSevLevelWatcher_
private

Definition at line 72 of file CaloTowersCreator.h.

Referenced by produce().

◆ EEEScale

double CaloTowersCreator::EEEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ eScales_

EScales CaloTowersCreator::eScales_
private

Definition at line 73 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ HBEScale

double CaloTowersCreator::HBEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ hcalChStatusWatcher_

edm::ESWatcher<HcalChannelQualityRcd> CaloTowersCreator::hcalChStatusWatcher_
private

Definition at line 70 of file CaloTowersCreator.h.

Referenced by produce().

◆ hcalSevLevelWatcher_

edm::ESWatcher<HcalSeverityLevelComputerRcd> CaloTowersCreator::hcalSevLevelWatcher_
private

Definition at line 69 of file CaloTowersCreator.h.

Referenced by produce().

◆ HEDEScale

double CaloTowersCreator::HEDEScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ HESEScale

double CaloTowersCreator::HESEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ HF1EScale

double CaloTowersCreator::HF1EScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ HF2EScale

double CaloTowersCreator::HF2EScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ HOEScale

double CaloTowersCreator::HOEScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ theEcalSeveritiesToBeExcluded_

std::vector<int> CaloTowersCreator::theEcalSeveritiesToBeExcluded_
private

Definition at line 48 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ theEcalSeveritiesToBeUsedInBadTowers_

std::vector<int> CaloTowersCreator::theEcalSeveritiesToBeUsedInBadTowers_
private

Definition at line 59 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ theHcalAcceptSeverityLevel_

unsigned int CaloTowersCreator::theHcalAcceptSeverityLevel_
private

Definition at line 47 of file CaloTowersCreator.h.

Referenced by produce().

◆ theHcalAcceptSeverityLevelForRejectedHit_

unsigned int CaloTowersCreator::theHcalAcceptSeverityLevelForRejectedHit_
private

Definition at line 57 of file CaloTowersCreator.h.

Referenced by produce().

◆ theRecoveredEcalHitsAreUsed_

bool CaloTowersCreator::theRecoveredEcalHitsAreUsed_
private

Definition at line 52 of file CaloTowersCreator.h.

Referenced by produce().

◆ theRecoveredHcalHitsAreUsed_

bool CaloTowersCreator::theRecoveredHcalHitsAreUsed_
private

Definition at line 51 of file CaloTowersCreator.h.

Referenced by produce().

◆ tok_hbhe_

edm::EDGetTokenT<HBHERecHitCollection> CaloTowersCreator::tok_hbhe_
private

Definition at line 38 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ tok_hf_

edm::EDGetTokenT<HFRecHitCollection> CaloTowersCreator::tok_hf_
private

Definition at line 40 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ tok_ho_

edm::EDGetTokenT<HORecHitCollection> CaloTowersCreator::tok_ho_
private

Definition at line 39 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ toks_ecal_

std::vector<edm::EDGetTokenT<EcalRecHitCollection> > CaloTowersCreator::toks_ecal_
private

Definition at line 42 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ useRejectedHitsOnly_

bool CaloTowersCreator::useRejectedHitsOnly_
private

Definition at line 56 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

◆ useRejectedRecoveredEcalHits_

bool CaloTowersCreator::useRejectedRecoveredEcalHits_
private

Definition at line 67 of file CaloTowersCreator.h.

Referenced by produce().

◆ useRejectedRecoveredHcalHits_

bool CaloTowersCreator::useRejectedRecoveredHcalHits_
private

Definition at line 66 of file CaloTowersCreator.h.

Referenced by produce().

edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:52
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
EcalSeverityLevelAlgo
Definition: EcalSeverityLevelAlgo.h:33
CaloTowersCreator::hcalSevLevelWatcher_
edm::ESWatcher< HcalSeverityLevelComputerRcd > hcalSevLevelWatcher_
Definition: CaloTowersCreator.h:69
CaloTowersCreationAlgo::setEcalChStatusFromDB
void setEcalChStatusFromDB(const EcalChannelStatus *s)
Definition: CaloTowersCreationAlgo.h:164
mps_fire.i
i
Definition: mps_fire.py:355
CaloTowersCreationAlgo::setHF1EScale
void setHF1EScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1383
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:95
CaloTowersCreator::theRecoveredHcalHitsAreUsed_
bool theRecoveredHcalHitsAreUsed_
Definition: CaloTowersCreator.h:51
CaloTowersCreationAlgo::setRecoveredHcalHitsAreUsed
void setRecoveredHcalHitsAreUsed(bool flag)
Definition: CaloTowersCreationAlgo.h:209
CaloTowersCreator::EEEScale
double EEEScale
Definition: CaloTowersCreator.h:31
CaloTowersCreationAlgo::setGeometry
void setGeometry(const CaloTowerTopology *cttopo, const CaloTowerConstituentsMap *ctmap, const HcalTopology *htopo, const CaloGeometry *geo)
Definition: CaloTowersCreationAlgo.cc:354
CaloTowersCreationAlgo::finish
void finish(CaloTowerCollection &destCollection)
Definition: CaloTowersCreationAlgo.cc:404
CaloTowersCreationAlgo::setHOEScale
void setHOEScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1376
HcalChannelQualityRcd
Definition: HcalChannelQualityRcd.h:8
CaloTowersCreator::ecalSevLevelWatcher_
edm::ESWatcher< EcalSeverityLevelAlgoRcd > ecalSevLevelWatcher_
Definition: CaloTowersCreator.h:72
edm::EDGetTokenT
Definition: EDGetToken.h:33
CaloGeometryRecord
Definition: CaloGeometryRecord.h:30
CaloTowersCreator::eScales_
EScales eScales_
Definition: CaloTowersCreator.h:73
CaloTowersCreator::useRejectedRecoveredEcalHits_
bool useRejectedRecoveredEcalHits_
Definition: CaloTowersCreator.h:67
CaloTowersCreationAlgo::setEcalSevLvlAlgo
void setEcalSevLvlAlgo(const EcalSeverityLevelAlgo *a)
Definition: CaloTowersCreationAlgo.h:216
EScales::instanceLabel
std::string instanceLabel
Definition: EScales.h:21
CaloTowersCreator::EBEScale
double EBEScale
Definition: CaloTowersCreator.h:31
gather_cfg.cout
cout
Definition: gather_cfg.py:144
es_hardcode_cfi.hf
hf
Definition: es_hardcode_cfi.py:161
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
CaloTowersCreationAlgo::setHF2EScale
void setHF2EScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1390
HcalChannelQuality
Definition: HcalChannelQuality.h:17
CaloTowersCreator::algo_
CaloTowersCreationAlgo algo_
Definition: CaloTowersCreator.h:37
EcalCondObjectContainer< EcalChannelStatusCode >
edm::Handle
Definition: AssociativeIterator.h:50
CaloTowersCreationAlgo::setUseRejectedRecoveredEcalHits
void setUseRejectedRecoveredEcalHits(bool flag)
Definition: CaloTowersCreationAlgo.h:233
EcalBarrel
Definition: EcalSubdetector.h:10
EScales::HBScale
double HBScale
Definition: EScales.h:15
CaloTowersCreationAlgo::setHBEScale
void setHBEScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1355
CaloTowersCreator::HEDEScale
double HEDEScale
Definition: CaloTowersCreator.h:32
CaloTowersCreator::HESEScale
double HESEScale
Definition: CaloTowersCreator.h:31
CaloTowersCreationAlgo::begin
void begin()
Definition: CaloTowersCreationAlgo.cc:368
CaloTowersCreator::theEcalSeveritiesToBeExcluded_
std::vector< int > theEcalSeveritiesToBeExcluded_
Definition: CaloTowersCreator.h:48
EScales::HESScale
double HESScale
Definition: EScales.h:16
EScales::EBScale
double EBScale
Definition: EScales.h:13
CaloTowersCreationAlgo::makeEcalBadChs
void makeEcalBadChs()
Definition: CaloTowersCreationAlgo.cc:1710
CaloTowersCreator::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: CaloTowersCreator.h:38
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:262
CaloTowersCreationAlgo::setEeHandle
void setEeHandle(const edm::Handle< EcalRecHitCollection > ee)
Definition: CaloTowersCreationAlgo.h:241
CaloTowersCreator::tok_ho_
edm::EDGetTokenT< HORecHitCollection > tok_ho_
Definition: CaloTowersCreator.h:39
HcalSeverityLevelComputer
Definition: HcalSeverityLevelComputer.h:24
dumpMFGeometry_cfg.prod
prod
Definition: dumpMFGeometry_cfg.py:24
edm::ESHandle< CaloGeometry >
EcalSeverityLevelAlgoRcd
Definition: EcalSeverityLevelAlgoRcd.h:12
CaloTowersCreationAlgo::SetEcalSeveritiesToBeUsedInBadTowers
void SetEcalSeveritiesToBeUsedInBadTowers(const std::vector< int > &ecalSev)
Definition: CaloTowersCreationAlgo.h:228
CaloTowersCreator::HF1EScale
double HF1EScale
Definition: CaloTowersCreator.h:32
CaloTowersCreationAlgo::setEcalSeveritiesToBeExcluded
void setEcalSeveritiesToBeExcluded(const std::vector< int > &ecalSev)
Definition: CaloTowersCreationAlgo.h:206
CaloTowersCreationAlgo::setHcalChStatusFromDB
void setHcalChStatusFromDB(const HcalChannelQuality *s)
Definition: CaloTowersCreationAlgo.h:163
HcalRecNumberingRecord
Definition: HcalRecNumberingRecord.h:23
CaloTowersCreator::toks_ecal_
std::vector< edm::EDGetTokenT< EcalRecHitCollection > > toks_ecal_
Definition: CaloTowersCreator.h:42
CaloTowersCreator::useRejectedHitsOnly_
bool useRejectedHitsOnly_
Definition: CaloTowersCreator.h:56
EcalEndcap
Definition: EcalSubdetector.h:10
CaloTowersCreationAlgo::setHEDEScale
void setHEDEScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1369
CaloTowersCreationAlgo::setHcalSevLvlComputer
void setHcalSevLvlComputer(const HcalSeverityLevelComputer *c)
Definition: CaloTowersCreationAlgo.h:213
EScales::HF2Scale
double HF2Scale
Definition: EScales.h:20
CaloTowersCreationAlgo::setRecoveredEcalHitsAreUsed
void setRecoveredEcalHitsAreUsed(bool flag)
Definition: CaloTowersCreationAlgo.h:210
CaloTowersCreationAlgo::setEEEScale
void setEEEScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1348
HLT_2018_cff.InputTag
InputTag
Definition: HLT_2018_cff.py:79016
EScales::EEScale
double EEScale
Definition: EScales.h:14
CaloTowersCreationAlgo::setHcalAcceptSeverityLevel
void setHcalAcceptSeverityLevel(unsigned int level)
Definition: CaloTowersCreationAlgo.h:205
EcalChannelStatusRcd
Definition: EcalChannelStatusRcd.h:5
CaloTowersCreationAlgo::process
void process(const HBHERecHitCollection &hbhe)
Definition: CaloTowersCreationAlgo.cc:374
CaloTowersCreationAlgo::makeHcalDropChMap
void makeHcalDropChMap()
Definition: CaloTowersCreationAlgo.cc:1649
CaloTowersCreationAlgo::setHcalAcceptSeverityLevelForRejectedHit
void setHcalAcceptSeverityLevelForRejectedHit(unsigned int level)
Definition: CaloTowersCreationAlgo.h:224
CaloTowersCreator::HBEScale
double HBEScale
Definition: CaloTowersCreator.h:31
CaloTowersCreator::hcalChStatusWatcher_
edm::ESWatcher< HcalChannelQualityRcd > hcalChStatusWatcher_
Definition: CaloTowersCreator.h:70
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
get
#define get
CaloTowersCreator::useRejectedRecoveredHcalHits_
bool useRejectedRecoveredHcalHits_
Definition: CaloTowersCreator.h:66
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
CaloTowersCreationAlgo::setMissingHcalRescaleFactorForEcal
void setMissingHcalRescaleFactorForEcal(float factor)
Definition: CaloTowersCreationAlgo.h:234
CaloTowersCreationAlgo::setEbHandle
void setEbHandle(const edm::Handle< EcalRecHitCollection > eb)
Definition: CaloTowersCreationAlgo.h:240
CaloTowersCreator::theHcalAcceptSeverityLevel_
unsigned int theHcalAcceptSeverityLevel_
Definition: CaloTowersCreator.h:47
CaloTowersCreator::caloTowerConstituentsWatcher_
edm::ESWatcher< IdealGeometryRecord > caloTowerConstituentsWatcher_
Definition: CaloTowersCreator.h:71
eostools.move
def move(src, dest)
Definition: eostools.py:511
CaloTowersCreator::HF2EScale
double HF2EScale
Definition: CaloTowersCreator.h:32
CaloTowersCreator::theHcalAcceptSeverityLevelForRejectedHit_
unsigned int theHcalAcceptSeverityLevelForRejectedHit_
Definition: CaloTowersCreator.h:57
EScales::HF1Scale
double HF1Scale
Definition: EScales.h:19
CaloTowersCreator::allowMissingInputs_
bool allowMissingInputs_
Definition: CaloTowersCreator.h:43
CaloTowersCreationAlgo::setEBEScale
void setEBEScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1341
EScales::HEDScale
double HEDScale
Definition: EScales.h:17
CaloTowersCreator::tok_hf_
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
Definition: CaloTowersCreator.h:40
CaloTowersCreationAlgo::setUseRejectedHitsOnly
void setUseRejectedHitsOnly(bool flag)
Definition: CaloTowersCreationAlgo.h:222
CaloTowersCreator::theEcalSeveritiesToBeUsedInBadTowers_
std::vector< int > theEcalSeveritiesToBeUsedInBadTowers_
Definition: CaloTowersCreator.h:59
es_hardcode_cfi.ho
ho
Definition: es_hardcode_cfi.py:191
CaloTowersCreationAlgo::setHESEScale
void setHESEScale(double scale)
Definition: CaloTowersCreationAlgo.cc:1362
EScales::HOScale
double HOScale
Definition: EScales.h:18
CaloTowersCreationAlgo::setUseRejectedRecoveredHcalHits
void setUseRejectedRecoveredHcalHits(bool flag)
Definition: CaloTowersCreationAlgo.h:232
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99
CaloTowersCreator::ecalLabels_
std::vector< edm::InputTag > ecalLabels_
Definition: CaloTowersCreator.h:41
edm::InputTag
Definition: InputTag.h:15
HcalSeverityLevelComputerRcd
Definition: HcalSeverityLevelComputerRcd.h:23
CaloTowersCreator::HOEScale
double HOEScale
Definition: CaloTowersCreator.h:32
CaloTowersCreator::theRecoveredEcalHitsAreUsed_
bool theRecoveredEcalHitsAreUsed_
Definition: CaloTowersCreator.h:52
reco::LeafCandidate::LorentzVector
math::XYZTLorentzVector LorentzVector
Lorentz vector.
Definition: LeafCandidate.h:23
edm::SortedCollection::empty
bool empty() const
Definition: SortedCollection.h:210
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
Generated for CMSSW Reference Manual by doxygen 1.8.16