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<> edm::stream::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 CaloTowersCreator (const edm::ParameterSet &ps)
 
virtual void produce (edm::Event &e, const edm::EventSetup &c)
 
virtual ~CaloTowersCreator ()
 
- Public Member Functions inherited from edm::stream::EDProducer<>
 EDProducer ()=default
 
- Public Member Functions inherited from edm::stream::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducerBase ()
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (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
 
- Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Original author: J. Mans - Minnesota

Definition at line 25 of file CaloTowersCreator.h.

Constructor & Destructor Documentation

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

Definition at line 12 of file CaloTowersCreator.cc.

References 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, theEcalSeveritiesToBeExcluded_, theEcalSeveritiesToBeUsedInBadTowers_, tok_hbhe_, tok_hf_, tok_ho_, and toks_ecal_.

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

Definition at line 28 of file CaloTowersCreator.h.

References EnergyCorrector::c, MillePedeFileConverter_cfg::e, fillDescriptions(), and produce().

28 { }

Member Function Documentation

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

Definition at line 296 of file CaloTowersCreator.cc.

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

Referenced by ~CaloTowersCreator().

296  {
297  edm::ParameterSetDescription desc;
298  desc.add<double>("EBSumThreshold", 0.2);
299  desc.add<double>("HF2Weight", 1.0);
300  desc.add<double>("EBWeight", 1.0);
301  desc.add<double>("EESumThreshold", 0.45);
302  desc.add<double>("HOThreshold0", 1.1);
303  desc.add<double>("HOThresholdPlus1", 3.5);
304  desc.add<double>("HOThresholdMinus1", 3.5);
305  desc.add<double>("HOThresholdPlus2", 3.5);
306  desc.add<double>("HOThresholdMinus2", 3.5);
307  desc.add<double>("HBThreshold", 0.7);
308  desc.add<double>("HF1Threshold", 0.5);
309  desc.add<double>("HEDWeight", 1.0);
310  desc.add<double>("EEWeight", 1.0);
311  desc.add<double>("HESWeight", 1.0);
312  desc.add<double>("HF1Weight", 1.0);
313  desc.add<double>("HOWeight", 1.0);
314  desc.add<double>("EBThreshold", 0.07);
315  desc.add<double>("EEThreshold", 0.3);
316  desc.add<double>("HcalThreshold", -1000.0);
317  desc.add<double>("HF2Threshold", 0.85);
318  desc.add<double>("HESThreshold", 0.8);
319  desc.add<double>("HEDThreshold", 0.8);
320  desc.add<double>("EcutTower", -1000.0);
321  desc.add<double>("HBWeight", 1.0);
322  desc.add<double>("MomHBDepth", 0.2);
323  desc.add<double>("MomHEDepth", 0.4);
324  desc.add<double>("MomEBDepth", 0.3);
325  desc.add<double>("MomEEDepth", 0.0);
326  desc.add<bool>("UseHO", true);
327  desc.add<bool>("UseEtEBTreshold", false);
328  desc.add<bool>("UseSymEBTreshold", true);
329  desc.add<bool>("UseEtEETreshold", false);
330  desc.add<bool>("UseSymEETreshold", true);
331  desc.add<bool>("UseHcalRecoveredHits", true);
332  desc.add<bool>("UseEcalRecoveredHits", false);
333  desc.add<bool>("UseRejectedHitsOnly", false);
334  desc.add<bool>("UseRejectedRecoveredHcalHits", true);
335  desc.add<bool>("UseRejectedRecoveredEcalHits", false);
336  desc.add<bool>("AllowMissingInputs", false);
337  desc.add<std::vector<double> >("HBGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
338  desc.add<std::vector<double> >("EEWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
339  desc.add<std::vector<double> >("HF2Weights", {1.0, 1.0, 1.0, 1.0, 1.0});
340  desc.add<std::vector<double> >("HOWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
341  desc.add<std::vector<double> >("EEGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
342  desc.add<std::vector<double> >("HBWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
343  desc.add<std::vector<double> >("HF2Grid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
344  desc.add<std::vector<double> >("HEDWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
345  desc.add<std::vector<double> >("HF1Grid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
346  desc.add<std::vector<double> >("EBWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
347  desc.add<std::vector<double> >("HF1Weights", {1.0, 1.0, 1.0, 1.0, 1.0});
348  desc.add<std::vector<double> >("HESGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
349  desc.add<std::vector<double> >("HESWeights", {1.0, 1.0, 1.0, 1.0, 1.0});
350  desc.add<std::vector<double> >("HEDGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
351  desc.add<std::vector<double> >("HOGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
352  desc.add<std::vector<double> >("EBGrid", {-1.0, 1.0, 10.0, 100.0, 1000.0});
353  desc.add<edm::InputTag>("hfInput", edm::InputTag("hfreco"));
354  desc.add<edm::InputTag>("hbheInput", edm::InputTag("hbhereco"));
355  desc.add<edm::InputTag>("hoInput", edm::InputTag("horeco"));
356  desc.add<std::vector<edm::InputTag> >("ecalInputs", {edm::InputTag("ecalRecHit","EcalRecHitsEB"), edm::InputTag("ecalRecHit","EcalRecHitsEE")});
357  desc.add<int>("MomConstrMethod", 1);
358  desc.add<unsigned int>("HcalAcceptSeverityLevel", 9);
359  desc.add<std::vector<std::string> >("EcalRecHitSeveritiesToBeExcluded", {"kTime","kWeird","kBad"});
360  desc.add<unsigned int>("HcalAcceptSeverityLevelForRejectedHit", 9999);
361  desc.add<std::vector<std::string> >("EcalSeveritiesToBeUsedInBadTowers", {});
362  desc.add<int>("HcalPhase", 0);
363 
364  descriptions.addDefault(desc);
365 }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89
static const std::vector<double>& CaloTowersCreator::getGridValues ( )
staticprivate
void CaloTowersCreator::produce ( edm::Event &  e,
const edm::EventSetup &  c 
)
virtual

Definition at line 133 of file CaloTowersCreator.cc.

References algo_, allowMissingInputs_, CaloTowersCreationAlgo::begin(), edm::SortedCollection< T, SORT >::begin(), caloTowerConstituentsWatcher_, edm::ESWatcher< T >::check(), EBEScale, RecoEcal_EventContent_cff::ec, EcalBarrel, EcalEndcap, ecalSevLevelWatcher_, EEEScale, eScales_, CaloTowersCreationAlgo::finish(), edm::EventSetup::get(), edm::Event::getByToken(), HBEScale, photonIsolationHIProducer_cfi::hbhe, hcalChStatusWatcher_, hcalSevLevelWatcher_, HEDEScale, HESEScale, photonIsolationHIProducer_cfi::hf, HF1EScale, HF2EScale, photonIsolationHIProducer_cfi::ho, HOEScale, mps_fire::i, EScales::instanceLabel, CaloTowersCreationAlgo::makeEcalBadChs(), CaloTowersCreationAlgo::makeHcalDropChMap(), eostools::move(), CaloTowersCreationAlgo::process(), parseEventContent::prod, edm::ESHandle< T >::product(), edm::Event::put(), 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(), edm::SortedCollection< T, SORT >::size(), theEcalSeveritiesToBeExcluded_, theEcalSeveritiesToBeUsedInBadTowers_, theHcalAcceptSeverityLevel_, theHcalAcceptSeverityLevelForRejectedHit_, theRecoveredEcalHitsAreUsed_, theRecoveredHcalHitsAreUsed_, tok_hbhe_, tok_hf_, tok_ho_, toks_ecal_, useRejectedHitsOnly_, useRejectedRecoveredEcalHits_, and useRejectedRecoveredHcalHits_.

Referenced by ~CaloTowersCreator().

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 
165  algo_.setEBEScale(EBEScale);
166  algo_.setEEEScale(EEEScale);
167  algo_.setHBEScale(HBEScale);
168  algo_.setHESEScale(HESEScale);
169  algo_.setHEDEScale(HEDEScale);
170  algo_.setHOEScale(HOEScale);
171  algo_.setHF1EScale(HF1EScale);
172  algo_.setHF2EScale(HF2EScale);
173  algo_.setGeometry(cttopo.product(),ctmap.product(),htopo.product(),pG.product());
174 
175  // for treatment of problematic and anomalous cells
176 
177  algo_.setHcalChStatusFromDB(dbHcalChStatus);
178  algo_.setEcalChStatusFromDB(dbEcalChStatus);
179 
180  algo_.setHcalAcceptSeverityLevel(theHcalAcceptSeverityLevel_);
181  algo_.setEcalSeveritiesToBeExcluded(theEcalSeveritiesToBeExcluded_);
182 
183  algo_.setRecoveredHcalHitsAreUsed(theRecoveredHcalHitsAreUsed_);
184  algo_.setRecoveredEcalHitsAreUsed(theRecoveredEcalHitsAreUsed_);
185 
186  algo_.setHcalSevLvlComputer(hcalSevLvlComputer);
187  algo_.setEcalSevLvlAlgo(ecalSevLvlAlgo);
188 
189 
190  algo_.setUseRejectedHitsOnly(useRejectedHitsOnly_);
191 
192  algo_.setHcalAcceptSeverityLevelForRejectedHit(theHcalAcceptSeverityLevelForRejectedHit_);
193  algo_.SetEcalSeveritiesToBeUsedInBadTowers (theEcalSeveritiesToBeUsedInBadTowers_);
194 
195  algo_.setUseRejectedRecoveredHcalHits(useRejectedRecoveredHcalHits_);
196  algo_.setUseRejectedRecoveredEcalHits(useRejectedRecoveredEcalHits_);
197 
198  /*
199  std::cout << "VI Produce: "
200  << (useRejectedHitsOnly_ ? "use rejectOnly " : " ")
201  << (allowMissingInputs_ ? "allowMissing " : " " )
202  << (theRecoveredEcalHitsAreUsed_ ? "use RecoveredEcal ": " " )
203  << toks_ecal_.size()
204  << ' ' << theEcalSeveritiesToBeExcluded_.size()
205  << ' ' << theEcalSeveritiesToBeUsedInBadTowers_.size()
206  << std::endl;
207  */
208 
209  algo_.begin(); // clear the internal buffer
210 
211  // can't chain these in a big OR statement, or else it'll
212  // get triggered for each of the first three events
213  bool check1 = hcalSevLevelWatcher_.check(c);
214  bool check2 = hcalChStatusWatcher_.check(c);
215  bool check3 = caloTowerConstituentsWatcher_.check(c);
216  if(check1 || check2 || check3)
217  {
218  algo_.makeHcalDropChMap();
219  }
220 
221  // check ecal SevLev
222  if (ecalSevLevelWatcher_.check(c)) algo_.makeEcalBadChs();
223 
224  // ----------------------------------------------------------
225  // For ecal error handling need to
226  // have access to the EB and EE collections at the end of
227  // tower reconstruction.
228 
229  edm::Handle<EcalRecHitCollection> ebHandle;
230  edm::Handle<EcalRecHitCollection> eeHandle;
231 
232  for (std::vector<edm::EDGetTokenT<EcalRecHitCollection> >::const_iterator i=toks_ecal_.begin();
233  i!=toks_ecal_.end(); i++) {
234 
235  edm::Handle<EcalRecHitCollection> ec_tmp;
236 
237  if (! e.getByToken(*i,ec_tmp) ) continue;
238  if (ec_tmp->size()==0) continue;
239 
240  // check if this is EB or EE
241  if ( (ec_tmp->begin()->detid()).subdetId() == EcalBarrel ) {
242  ebHandle = ec_tmp;
243  }
244  else if ((ec_tmp->begin()->detid()).subdetId() == EcalEndcap) {
245  eeHandle = ec_tmp;
246  }
247  }
248 
249  algo_.setEbHandle(ebHandle);
250  algo_.setEeHandle(eeHandle);
251 
252  //-----------------------------------------------------------
253 
254  bool present;
255 
256  // Step A/C: Get Inputs and process (repeatedly)
257  edm::Handle<HBHERecHitCollection> hbhe;
258  present=e.getByToken(tok_hbhe_,hbhe);
259  if (present || !allowMissingInputs_) algo_.process(*hbhe);
260 
261  edm::Handle<HORecHitCollection> ho;
262  present=e.getByToken(tok_ho_,ho);
263  if (present || !allowMissingInputs_) algo_.process(*ho);
264 
265  edm::Handle<HFRecHitCollection> hf;
266  present=e.getByToken(tok_hf_,hf);
267  if (present || !allowMissingInputs_) algo_.process(*hf);
268 
269  std::vector<edm::EDGetTokenT<EcalRecHitCollection> >::const_iterator i;
270  for (i=toks_ecal_.begin(); i!=toks_ecal_.end(); i++) {
271  edm::Handle<EcalRecHitCollection> ec;
272  present=e.getByToken(*i,ec);
273  if (present || !allowMissingInputs_) 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  /*
283  int totc=0; float totE=0;
284  reco::LeafCandidate::LorentzVector totP4;
285  for (auto const & tw : (*prod) ) { totc += tw.constituents().size(); totE+=tw.energy(); totP4+=tw.p4();}
286  std::cout << "VI " << (*prod).size() << " " << totc << " " << totE << " " << totP4 << std::endl;
287  */
288 
289  // Step D: Put into the event
290  if (eScales_.instanceLabel=="") e.put(std::move(prod));
291  else e.put(std::move(prod),eScales_.instanceLabel);
292 
293 
294 }
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
CaloTowersCreationAlgo::setEcalChStatusFromDB
void setEcalChStatusFromDB(const EcalChannelStatus *s)
Definition: CaloTowersCreationAlgo.h:115
mps_fire.i
i
Definition: mps_fire.py:156
CaloTowersCreator::algo_
CaloTowersCreationAlgo algo_
Definition: CaloTowersCreator.h:38
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:460
CaloTowersCreationAlgo::setHcalChStatusFromDB
void setHcalChStatusFromDB(const HcalChannelQuality *s)
Definition: CaloTowersCreationAlgo.h:114
CaloTowersCreationAlgo::finish
void finish(CaloTowerCollection &destCollection)
Definition: CaloTowersCreationAlgo.cc:417
CaloTowersCreationAlgo::setEeHandle
void setEeHandle(const edm::Handle< EcalRecHitCollection > ee)
Definition: CaloTowersCreationAlgo.h:193
photonIsolationHIProducer_cfi.ho
ho
Definition: photonIsolationHIProducer_cfi.py:10
CaloTowersCreator::theEcalSeveritiesToBeExcluded_
std::vector< int > theEcalSeveritiesToBeExcluded_
Definition: CaloTowersCreator.h:50
CaloTowersCreator::tok_ho_
edm::EDGetTokenT< HORecHitCollection > tok_ho_
Definition: CaloTowersCreator.h:40
CaloTowersCreationAlgo::setGeometry
void setGeometry(const CaloTowerTopology *cttopo, const CaloTowerConstituentsMap *ctmap, const HcalTopology *htopo, const CaloGeometry *geo)
Definition: CaloTowersCreationAlgo.cc:305
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
CaloTowersCreationAlgo::setEbHandle
void setEbHandle(const edm::Handle< EcalRecHitCollection > eb)
Definition: CaloTowersCreationAlgo.h:192
CaloTowersCreationAlgo::setUseRejectedRecoveredEcalHits
void setUseRejectedRecoveredEcalHits(bool flag)
Definition: CaloTowersCreationAlgo.h:184
RecoEcal_EventContent_cff.ec
ec
Definition: RecoEcal_EventContent_cff.py:106
CaloTowersCreationAlgo::setEcalSeveritiesToBeExcluded
void setEcalSeveritiesToBeExcluded(const std::vector< int > &ecalSev)
Definition: CaloTowersCreationAlgo.h:159
CaloTowersCreationAlgo::setHcalAcceptSeverityLevelForRejectedHit
void setHcalAcceptSeverityLevelForRejectedHit(unsigned int level)
Definition: CaloTowersCreationAlgo.h:178
CaloTowersCreationAlgo::setHcalSevLvlComputer
void setHcalSevLvlComputer(const HcalSeverityLevelComputer *c)
Definition: CaloTowersCreationAlgo.h:166
CaloTowersCreator::toks_ecal_
std::vector< edm::EDGetTokenT< EcalRecHitCollection > > toks_ecal_
Definition: CaloTowersCreator.h:43
CaloTowersCreator::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: CaloTowersCreator.h:39
CaloTowersCreationAlgo::SetEcalSeveritiesToBeUsedInBadTowers
void SetEcalSeveritiesToBeUsedInBadTowers(const std::vector< int > &ecalSev)
Definition: CaloTowersCreationAlgo.h:180
CaloTowersCreationAlgo::setRecoveredEcalHitsAreUsed
void setRecoveredEcalHitsAreUsed(bool flag)
Definition: CaloTowersCreationAlgo.h:163
CaloTowersCreator::hcalChStatusWatcher_
edm::ESWatcher< HcalChannelQualityRcd > hcalChStatusWatcher_
Definition: CaloTowersCreator.h:75
CaloTowersCreationAlgo::process
void process(const HBHERecHitCollection &hbhe)
Definition: CaloTowersCreationAlgo.cc:380
photonIsolationHIProducer_cfi.hf
hf
Definition: photonIsolationHIProducer_cfi.py:9
EScales::instanceLabel
std::string instanceLabel
Definition: EScales.h:21
CaloTowersCreator::theHcalAcceptSeverityLevelForRejectedHit_
unsigned int theHcalAcceptSeverityLevelForRejectedHit_
Definition: CaloTowersCreator.h:60
CaloTowersCreationAlgo::setRecoveredHcalHitsAreUsed
void setRecoveredHcalHitsAreUsed(bool flag)
Definition: CaloTowersCreationAlgo.h:162
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
CaloTowersCreationAlgo::setUseRejectedHitsOnly
void setUseRejectedHitsOnly(bool flag)
Definition: CaloTowersCreationAlgo.h:176
edm::ESWatcher::check
bool check(const edm::EventSetup &iSetup)
Definition: ESWatcher.h:57
CaloTowersCreator::theHcalAcceptSeverityLevel_
unsigned int theHcalAcceptSeverityLevel_
Definition: CaloTowersCreator.h:49
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:236
CaloTowersCreationAlgo::setHcalAcceptSeverityLevel
void setHcalAcceptSeverityLevel(unsigned int level)
Definition: CaloTowersCreationAlgo.h:158
CaloTowersCreator::hcalSevLevelWatcher_
edm::ESWatcher< HcalSeverityLevelComputerRcd > hcalSevLevelWatcher_
Definition: CaloTowersCreator.h:74
CaloTowersCreationAlgo::setEcalSevLvlAlgo
void setEcalSevLvlAlgo(const EcalSeverityLevelAlgo *a)
Definition: CaloTowersCreationAlgo.h:169
CaloTowersCreationAlgo::setUseRejectedRecoveredHcalHits
void setUseRejectedRecoveredHcalHits(bool flag)
Definition: CaloTowersCreationAlgo.h:183
CaloTowersCreator::caloTowerConstituentsWatcher_
edm::ESWatcher< IdealGeometryRecord > caloTowerConstituentsWatcher_
Definition: CaloTowersCreator.h:76
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
CaloTowersCreator::theEcalSeveritiesToBeUsedInBadTowers_
std::vector< int > theEcalSeveritiesToBeUsedInBadTowers_
Definition: CaloTowersCreator.h:62
eostools.move
def move(src, dest)
Definition: eostools.py:510
CaloTowersCreator::tok_hf_
edm::EDGetTokenT< HFRecHitCollection > tok_hf_
Definition: CaloTowersCreator.h:41
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
CaloTowersCreator::ecalSevLevelWatcher_
edm::ESWatcher< EcalSeverityLevelAlgoRcd > ecalSevLevelWatcher_
Definition: CaloTowersCreator.h:77

Member Data Documentation

CaloTowersCreationAlgo CaloTowersCreator::algo_
private

Definition at line 38 of file CaloTowersCreator.h.

Referenced by produce().

bool CaloTowersCreator::allowMissingInputs_
private

Definition at line 44 of file CaloTowersCreator.h.

Referenced by produce().

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

Definition at line 76 of file CaloTowersCreator.h.

Referenced by produce().

double CaloTowersCreator::EBEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 42 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator().

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

Definition at line 77 of file CaloTowersCreator.h.

Referenced by produce().

double CaloTowersCreator::EEEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

EScales CaloTowersCreator::eScales_
private

Definition at line 78 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

double CaloTowersCreator::HBEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 75 of file CaloTowersCreator.h.

Referenced by produce().

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

Definition at line 74 of file CaloTowersCreator.h.

Referenced by produce().

double CaloTowersCreator::HEDEScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

double CaloTowersCreator::HESEScale

Definition at line 31 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

double CaloTowersCreator::HF1EScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

double CaloTowersCreator::HF2EScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

double CaloTowersCreator::HOEScale

Definition at line 32 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 50 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 62 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

unsigned int CaloTowersCreator::theHcalAcceptSeverityLevel_
private

Definition at line 49 of file CaloTowersCreator.h.

Referenced by produce().

unsigned int CaloTowersCreator::theHcalAcceptSeverityLevelForRejectedHit_
private

Definition at line 60 of file CaloTowersCreator.h.

Referenced by produce().

bool CaloTowersCreator::theRecoveredEcalHitsAreUsed_
private

Definition at line 54 of file CaloTowersCreator.h.

Referenced by produce().

bool CaloTowersCreator::theRecoveredHcalHitsAreUsed_
private

Definition at line 53 of file CaloTowersCreator.h.

Referenced by produce().

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

Definition at line 39 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 41 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 40 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

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

Definition at line 43 of file CaloTowersCreator.h.

Referenced by CaloTowersCreator(), and produce().

bool CaloTowersCreator::useRejectedHitsOnly_
private

Definition at line 59 of file CaloTowersCreator.h.

Referenced by produce().

bool CaloTowersCreator::useRejectedRecoveredEcalHits_
private

Definition at line 72 of file CaloTowersCreator.h.

Referenced by produce().

bool CaloTowersCreator::useRejectedRecoveredHcalHits_
private

Definition at line 71 of file CaloTowersCreator.h.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11