CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
MaterialBudgetHcalAnalysis Class Reference
Inheritance diagram for MaterialBudgetHcalAnalysis:
edm::one::EDAnalyzer< edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 MaterialBudgetHcalAnalysis (const edm::ParameterSet &p)
 
 ~MaterialBudgetHcalAnalysis () override=default
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueueglobalLuminosityBlocksQueue () final
 
SerialTaskQueueglobalRunsQueue () final
 
const EDAnalyzeroperator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () 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
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void analyze (edm::Event const &, edm::EventSetup const &) override
 
void beginJob () override
 

Private Attributes

const int binEta_
 
const int binPhi_
 
const double etaHigh_
 
const double etaHighMax_
 
const double etaHighMin_
 
const double etaLow_
 
const double etaLowMax_
 
const double etaLowMin_
 
const double etaMaxP_
 
const double etaMidMax_
 
const double etaMidMin_
 
const double etaMinP_
 
const edm::InputTag labelMBCalo_
 
const double maxEta_
 
TProfile2D * me1000_ [maxSet_]
 
TProfile * me100_ [maxSet_]
 
TProfile2D * me1100_ [maxSet_]
 
TH2F * me1200_ [maxSet_]
 
TH1F * me1300_ [maxSet2_]
 
TH2F * me1400_ [maxSet2_]
 
TProfile * me1500_ [maxSet2_]
 
TProfile * me1600_ [maxSet_]
 
TProfile * me1700_ [maxSet_]
 
TProfile * me1800_ [maxSet_]
 
TProfile * me1900_ [maxSet_]
 
TProfile * me2000_ [maxSet_]
 
TProfile * me200_ [maxSet_]
 
TProfile * me2100_ [maxSet_]
 
TProfile * me2200_ [maxSet_]
 
TProfile * me2300_ [maxSet_]
 
TProfile * me2400_ [maxSet_]
 
TProfile * me300_ [maxSet_]
 
TH1F * me400_ [maxSet_]
 
TProfile * me500_ [maxSet_]
 
TProfile * me600_ [maxSet_]
 
TProfile * me700_ [maxSet_]
 
TH1F * me800_ [maxSet_]
 
TProfile2D * me900_ [maxSet_]
 
const edm::EDGetTokenT< MaterialAccountingCaloCollectiontokMBCalo_
 

Static Private Attributes

static const uint32_t maxSet2_ = 9
 
static const uint32_t maxSet_ = 25
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- 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)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< Bconsumes (edm::InputTag tag) noexcept
 
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
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)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Definition at line 24 of file MaterialBudgetHcalAnalysis.cc.

Constructor & Destructor Documentation

◆ MaterialBudgetHcalAnalysis()

MaterialBudgetHcalAnalysis::MaterialBudgetHcalAnalysis ( const edm::ParameterSet p)

Definition at line 52 of file MaterialBudgetHcalAnalysis.cc.

References binEta_, binPhi_, etaHigh_, etaHighMax_, etaHighMin_, etaLow_, etaLowMax_, etaLowMin_, etaMaxP_, etaMidMax_, etaMidMin_, etaMinP_, TFileService::kSharedResource, and maxEta_.

53  : binEta_(p.getParameter<int>("nBinEta")),
54  binPhi_(p.getParameter<int>("nBinPhi")),
55  maxEta_(p.getParameter<double>("maxEta")),
56  etaLow_(p.getParameter<double>("etaLow")),
57  etaHigh_(p.getParameter<double>("etaHigh")),
58  etaLowMin_(p.getParameter<double>("etaLowMin")),
59  etaLowMax_(p.getParameter<double>("etaLowMax")),
60  etaMidMin_(p.getParameter<double>("etaMidMin")),
61  etaMidMax_(p.getParameter<double>("etaMidMax")),
62  etaHighMin_(p.getParameter<double>("etaHighMin")),
63  etaHighMax_(p.getParameter<double>("etaHighMax")),
64  etaMinP_(p.getParameter<double>("etaMinP")),
65  etaMaxP_(p.getParameter<double>("etaMaxP")),
66  labelMBCalo_(p.getParameter<edm::InputTag>("labelMBCaloLabel")),
67  tokMBCalo_(consumes<MaterialAccountingCaloCollection>(labelMBCalo_)) {
68  usesResource(TFileService::kSharedResource);
69  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalAnalysis: == Eta plot: NX " << binEta_ << " Range "
70  << -maxEta_ << ":" << maxEta_ << " Phi plot: NX " << binPhi_ << " Range " << -1._pi
71  << ":" << 1._pi << " (Eta limit " << etaLow_ << ":" << etaHigh_ << ")"
72  << " Eta range (" << etaLowMin_ << ":" << etaLowMax_ << "), (" << etaMidMin_ << ":"
73  << etaMidMax_ << "), (" << etaHighMin_ << ":" << etaHighMax_
74  << ") Debug for eta range " << etaMinP_ << ":" << etaMaxP_;
75 }
static const std::string kSharedResource
Definition: TFileService.h:76
Log< level::Info, true > LogVerbatim
const edm::EDGetTokenT< MaterialAccountingCaloCollection > tokMBCalo_

◆ ~MaterialBudgetHcalAnalysis()

MaterialBudgetHcalAnalysis::~MaterialBudgetHcalAnalysis ( )
overridedefault

Member Function Documentation

◆ analyze()

void MaterialBudgetHcalAnalysis::analyze ( edm::Event const &  iEvent,
edm::EventSetup const &  iSetup 
)
overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 235 of file MaterialBudgetHcalAnalysis.cc.

References funct::abs(), etaHigh_, etaHighMax_, etaHighMin_, etaLow_, etaLowMax_, etaLowMin_, etaMaxP_, etaMidMax_, etaMidMin_, etaMinP_, iEvent, cuy::ii, maxSet_, me1000_, me100_, me1100_, me1200_, me1300_, me1400_, me1500_, me1600_, me1700_, me1800_, me1900_, me2000_, me200_, me2100_, me2200_, me2300_, me2400_, me300_, me400_, me500_, me600_, me700_, me800_, me900_, and tokMBCalo_.

235  {
236 #ifdef EDM_ML_DEBUG
237  edm::LogVerbatim("MaterialBudgetFull") << "Run " << iEvent.id().run() << " Event " << iEvent.id().event()
238  << " Luminosity " << iEvent.luminosityBlock() << " Bunch "
239  << iEvent.bunchCrossing();
240 #endif
241 
242  // Fill from the MB collection
243  auto const &hcalMBColl = iEvent.getHandle(tokMBCalo_);
244  if (hcalMBColl.isValid()) {
245  auto hcalMB = hcalMBColl.product();
246 #ifdef EDM_ML_DEBUG
247  edm::LogVerbatim("MaterialBudgetFull")
248  << "Finds HcalMaterialBudgetCollection with " << hcalMB->size() << " entries";
249 #endif
250 
251  for (auto itr = hcalMB->begin(); itr != hcalMB->end(); ++itr) {
252  for (uint32_t ii = 0; ii < itr->m_stepLen.size(); ++ii) {
253 #ifdef EDM_ML_DEBUG
254  if ((std::abs(itr->m_eta) >= etaMinP_) && (std::abs(itr->m_eta) <= etaMaxP_))
255  edm::LogVerbatim("MaterialBudget")
256  << "MaterialBudgetHcalAnalysis:FillHisto called with index " << ii << " integrated step "
257  << itr->m_stepLen[ii] << " X0 " << itr->m_radLen[ii] << " Lamda " << itr->m_intLen[ii];
258 #endif
259  if (ii < maxSet_) {
260  me100_[ii]->Fill(itr->m_eta, itr->m_radLen[ii]);
261  me200_[ii]->Fill(itr->m_eta, itr->m_intLen[ii]);
262  me300_[ii]->Fill(itr->m_eta, itr->m_stepLen[ii]);
263  me400_[ii]->Fill(itr->m_eta);
264 
265  if (itr->m_eta >= etaLow_ && itr->m_eta <= etaHigh_) {
266  me500_[ii]->Fill(itr->m_phi, itr->m_radLen[ii]);
267  me600_[ii]->Fill(itr->m_phi, itr->m_intLen[ii]);
268  me700_[ii]->Fill(itr->m_phi, itr->m_stepLen[ii]);
269  me800_[ii]->Fill(itr->m_phi);
270  }
271 
272  me900_[ii]->Fill(itr->m_eta, itr->m_phi, itr->m_radLen[ii]);
273  me1000_[ii]->Fill(itr->m_eta, itr->m_phi, itr->m_intLen[ii]);
274  me1100_[ii]->Fill(itr->m_eta, itr->m_phi, itr->m_stepLen[ii]);
275  me1200_[ii]->Fill(itr->m_eta, itr->m_phi);
276 
277  if ((std::abs(itr->m_eta) >= etaMidMin_) && (std::abs(itr->m_eta) <= etaMidMax_)) {
278  me1600_[ii]->Fill(itr->m_phi, itr->m_radLen[ii]);
279  me1700_[ii]->Fill(itr->m_phi, itr->m_intLen[ii]);
280  me1800_[ii]->Fill(itr->m_phi, itr->m_stepLen[ii]);
281  }
282 
283  if ((std::abs(itr->m_eta) >= etaHighMin_) && (std::abs(itr->m_eta) <= etaHighMax_)) {
284  me1900_[ii]->Fill(itr->m_phi, itr->m_radLen[ii]);
285  me2000_[ii]->Fill(itr->m_phi, itr->m_intLen[ii]);
286  me2100_[ii]->Fill(itr->m_phi, itr->m_stepLen[ii]);
287  }
288 
289  if ((std::abs(itr->m_eta) >= etaLowMin_) && (std::abs(itr->m_eta) <= etaLowMax_)) {
290  me2200_[ii]->Fill(itr->m_phi, itr->m_radLen[ii]);
291  me2300_[ii]->Fill(itr->m_phi, itr->m_intLen[ii]);
292  me2400_[ii]->Fill(itr->m_phi, itr->m_stepLen[ii]);
293  }
294  }
295  }
296 
297  me1300_[0]->Fill(itr->m_eta);
298  me1400_[0]->Fill(itr->m_eta, itr->m_phi);
299  if (itr->m_layers[0] > 0) {
300  me1300_[1]->Fill(itr->m_eta);
301  me1400_[1]->Fill(itr->m_eta, itr->m_phi);
302  }
303  if (itr->m_layers[0] >= 16) {
304  me1300_[2]->Fill(itr->m_eta);
305  me1400_[2]->Fill(itr->m_eta, itr->m_phi);
306  }
307  if (itr->m_layers[1] > 0) {
308  me1300_[3]->Fill(itr->m_eta);
309  me1400_[3]->Fill(itr->m_eta, itr->m_phi);
310  }
311  if (itr->m_layers[1] >= 16) {
312  me1300_[4]->Fill(itr->m_eta);
313  me1400_[4]->Fill(itr->m_eta, itr->m_phi);
314  }
315  if (itr->m_layers[2] > 0) {
316  me1300_[5]->Fill(itr->m_eta);
317  me1400_[5]->Fill(itr->m_eta, itr->m_phi);
318  }
319  if (itr->m_layers[2] >= 2) {
320  me1300_[6]->Fill(itr->m_eta);
321  me1400_[6]->Fill(itr->m_eta, itr->m_phi);
322  }
323  if (itr->m_layers[3] > 0) {
324  me1300_[7]->Fill(itr->m_eta);
325  me1400_[7]->Fill(itr->m_eta, itr->m_phi);
326  }
327  if (itr->m_layers[0] > 0 || itr->m_layers[1] > 0 || (itr->m_layers[3] > 0 && std::abs(itr->m_eta) > 3.0)) {
328  me1300_[8]->Fill(itr->m_eta);
329  me1400_[8]->Fill(itr->m_eta, itr->m_phi);
330  }
331  me1500_[0]->Fill(itr->m_eta, (double)(itr->m_layers[0] + itr->m_layers[1] + itr->m_layers[2] + itr->m_layers[3]));
332  me1500_[1]->Fill(itr->m_eta, (double)(itr->m_layers[0]));
333  me1500_[2]->Fill(itr->m_eta, (double)(itr->m_layers[1]));
334  me1500_[4]->Fill(itr->m_eta, (double)(itr->m_layers[3]));
335  }
336  }
337 }
Log< level::Info, true > LogVerbatim
const edm::EDGetTokenT< MaterialAccountingCaloCollection > tokMBCalo_
int iEvent
Definition: GenABIO.cc:224
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
ii
Definition: cuy.py:589

◆ beginJob()

void MaterialBudgetHcalAnalysis::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 96 of file MaterialBudgetHcalAnalysis.cc.

References binEta_, binPhi_, etaHighMax_, etaHighMin_, etaLowMax_, etaLowMin_, etaMidMax_, etaMidMin_, mps_fire::i, maxEta_, HLT_2023v12_cff::maxPhi, maxSet2_, maxSet_, me1000_, me100_, me1100_, me1200_, me1300_, me1400_, me1500_, me1600_, me1700_, me1800_, me1900_, me2000_, me200_, me2100_, me2200_, me2300_, me2400_, me300_, me400_, me500_, me600_, me700_, me800_, me900_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, and to_string().

96  {
97  // Book histograms
99 
100  if (!tfile.isAvailable())
101  throw cms::Exception("BadConfig") << "TFileService unavailable: "
102  << "please add it to config file";
103 
104  double maxPhi = 1._pi;
105  edm::LogVerbatim("MaterialBudgetFull") << "MaterialBudgetHcalAnalysis: Booking user histos === with " << binEta_
106  << " bins in eta from " << -maxEta_ << " to " << maxEta_ << " and " << binPhi_
107  << " bins in phi from " << -maxPhi << " to " << maxPhi;
108 
109  std::string iter;
110  std::string range0 = "(" + std::to_string(etaMidMin_) + ":" + std::to_string(etaMidMax_) + ") ";
111  std::string range1 = "(" + std::to_string(etaHighMin_) + ":" + std::to_string(etaHighMax_) + ") ";
112  std::string range2 = "(" + std::to_string(etaLowMin_) + ":" + std::to_string(etaLowMax_) + ") ";
113  // total X0
114  for (uint32_t i = 0; i < maxSet_; i++) {
115  iter = std::to_string(i);
116  me100_[i] = tfile->make<TProfile>(
117  std::to_string(i + 100).c_str(), ("MB(X0) prof Eta in region " + iter).c_str(), binEta_, -maxEta_, maxEta_);
118  me200_[i] = tfile->make<TProfile>(
119  std::to_string(i + 200).c_str(), ("MB(L0) prof Eta in region " + iter).c_str(), binEta_, -maxEta_, maxEta_);
120  me300_[i] = tfile->make<TProfile>(
121  std::to_string(i + 300).c_str(), ("MB(Step) prof Eta in region " + iter).c_str(), binEta_, -maxEta_, maxEta_);
122  me400_[i] = tfile->make<TH1F>(
123  std::to_string(i + 400).c_str(), ("Eta in region " + iter).c_str(), binEta_, -maxEta_, maxEta_);
124  me500_[i] = tfile->make<TProfile>(
125  std::to_string(i + 500).c_str(), ("MB(X0) prof Ph in region " + iter).c_str(), binPhi_, -maxPhi, maxPhi);
126  me600_[i] = tfile->make<TProfile>(
127  std::to_string(i + 600).c_str(), ("MB(L0) prof Ph in region " + iter).c_str(), binPhi_, -maxPhi, maxPhi);
128  me700_[i] = tfile->make<TProfile>(
129  std::to_string(i + 700).c_str(), ("MB(Step) prof Ph in region " + iter).c_str(), binPhi_, -maxPhi, maxPhi);
130  me800_[i] =
131  tfile->make<TH1F>(std::to_string(i + 800).c_str(), ("Phi in region " + iter).c_str(), binPhi_, -maxPhi, maxPhi);
132  me900_[i] = tfile->make<TProfile2D>(std::to_string(i + 900).c_str(),
133  ("MB(X0) prof Eta Phi in region " + iter).c_str(),
134  binEta_ / 2,
135  -maxEta_,
136  maxEta_,
137  binPhi_ / 2,
138  -maxPhi,
139  maxPhi);
140  me1000_[i] = tfile->make<TProfile2D>(std::to_string(i + 1000).c_str(),
141  ("MB(L0) prof Eta Phi in region " + iter).c_str(),
142  binEta_ / 2,
143  -maxEta_,
144  maxEta_,
145  binPhi_ / 2,
146  -maxPhi,
147  maxPhi);
148  me1100_[i] = tfile->make<TProfile2D>(std::to_string(i + 1100).c_str(),
149  ("MB(Step) prof Eta Phi in region " + iter).c_str(),
150  binEta_ / 2,
151  -maxEta_,
152  maxEta_,
153  binPhi_ / 2,
154  -maxPhi,
155  maxPhi);
156  me1200_[i] = tfile->make<TH2F>(std::to_string(i + 1200).c_str(),
157  ("Eta vs Phi in region " + iter).c_str(),
158  binEta_ / 2,
159  -maxEta_,
160  maxEta_,
161  binPhi_ / 2,
162  -maxPhi,
163  maxPhi);
164  me1600_[i] = tfile->make<TProfile>(std::to_string(i + 1600).c_str(),
165  ("MB(X0) prof Ph in region " + range0 + iter).c_str(),
166  binPhi_,
167  -maxPhi,
168  maxPhi);
169  me1700_[i] = tfile->make<TProfile>(std::to_string(i + 1700).c_str(),
170  ("MB(L0) prof Ph in region " + range0 + iter).c_str(),
171  binPhi_,
172  -maxPhi,
173  maxPhi);
174  me1800_[i] = tfile->make<TProfile>(std::to_string(i + 1800).c_str(),
175  ("MB(Step) prof Ph in region " + range0 + iter).c_str(),
176  binPhi_,
177  -maxPhi,
178  maxPhi);
179  me1900_[i] = tfile->make<TProfile>(std::to_string(i + 1900).c_str(),
180  ("MB(X0) prof Ph in region " + range1 + iter).c_str(),
181  binPhi_,
182  -maxPhi,
183  maxPhi);
184  me2000_[i] = tfile->make<TProfile>(std::to_string(i + 2000).c_str(),
185  ("MB(L0) prof Ph in region " + range1 + iter).c_str(),
186  binPhi_,
187  -maxPhi,
188  maxPhi);
189  me2100_[i] = tfile->make<TProfile>(std::to_string(i + 2100).c_str(),
190  ("MB(Step) prof Ph in region " + range1 + iter).c_str(),
191  binPhi_,
192  -maxPhi,
193  maxPhi);
194  me2200_[i] = tfile->make<TProfile>(std::to_string(i + 2200).c_str(),
195  ("MB(X0) prof Ph in region " + range2 + iter).c_str(),
196  binPhi_,
197  -maxPhi,
198  maxPhi);
199  me2300_[i] = tfile->make<TProfile>(std::to_string(i + 2300).c_str(),
200  ("MB(L0) prof Ph in region " + range2 + iter).c_str(),
201  binPhi_,
202  -maxPhi,
203  maxPhi);
204  me2400_[i] = tfile->make<TProfile>(std::to_string(i + 2400).c_str(),
205  ("MB(Step) prof Ph in region " + range2 + iter).c_str(),
206  binPhi_,
207  -maxPhi,
208  maxPhi);
209  }
210  for (uint32_t i = 0; i < maxSet2_; i++) {
211  iter = std::to_string(i);
212  me1300_[i] = tfile->make<TH1F>(std::to_string(i + 1300).c_str(),
213  ("Events with layers Hit (0 all, 1 HB, ..) for " + iter).c_str(),
214  binEta_,
215  -maxEta_,
216  maxEta_);
217  me1400_[i] = tfile->make<TH2F>(std::to_string(i + 1400).c_str(),
218  ("Eta vs Phi for layers hit in " + iter).c_str(),
219  binEta_ / 2,
220  -maxEta_,
221  maxEta_,
222  binPhi_ / 2,
223  -maxPhi,
224  maxPhi);
225  me1500_[i] = tfile->make<TProfile>(std::to_string(i + 1500).c_str(),
226  ("Number of layers crossed (0 all, 1 HB, ..) for " + iter).c_str(),
227  binEta_,
228  -maxEta_,
229  maxEta_);
230  }
231 
232  edm::LogVerbatim("MaterialBudget") << "MaterialBudgetHcalAnalysis: Booking user histos done ===";
233 }
Log< level::Info, true > LogVerbatim
static std::string to_string(const XMLCh *ch)
Definition: tfile.py:1

◆ fillDescriptions()

void MaterialBudgetHcalAnalysis::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 77 of file MaterialBudgetHcalAnalysis.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

77  {
79  desc.add<int>("nBinEta", 260);
80  desc.add<int>("nBinPhi", 180);
81  desc.add<double>("maxEta", 5.2);
82  desc.add<double>("etaLow", -5.2);
83  desc.add<double>("etaHigh", 5.2);
84  desc.add<double>("etaMinP", 5.2);
85  desc.add<double>("etaMaxP", 0.0);
86  desc.add<double>("etaLowMin", 0.783);
87  desc.add<double>("etaLowMax", 0.870);
88  desc.add<double>("etaMidMin", 2.650);
89  desc.add<double>("etaMidMax", 2.868);
90  desc.add<double>("etaHighMin", 2.868);
91  desc.add<double>("etaHighMax", 3.000);
92  desc.add<edm::InputTag>("labelMBCaloLabel", edm::InputTag("g4SimHits", "HcalMatBCalo"));
93  descriptions.add("materialBudgetHcalAnalysis", desc);
94 }
void add(std::string const &label, ParameterSetDescription const &psetDescription)

Member Data Documentation

◆ binEta_

const int MaterialBudgetHcalAnalysis::binEta_
private

Definition at line 36 of file MaterialBudgetHcalAnalysis.cc.

Referenced by beginJob(), and MaterialBudgetHcalAnalysis().

◆ binPhi_

const int MaterialBudgetHcalAnalysis::binPhi_
private

Definition at line 36 of file MaterialBudgetHcalAnalysis.cc.

Referenced by beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaHigh_

const double MaterialBudgetHcalAnalysis::etaHigh_
private

Definition at line 37 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and MaterialBudgetHcalAnalysis().

◆ etaHighMax_

const double MaterialBudgetHcalAnalysis::etaHighMax_
private

Definition at line 38 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaHighMin_

const double MaterialBudgetHcalAnalysis::etaHighMin_
private

Definition at line 38 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaLow_

const double MaterialBudgetHcalAnalysis::etaLow_
private

Definition at line 37 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and MaterialBudgetHcalAnalysis().

◆ etaLowMax_

const double MaterialBudgetHcalAnalysis::etaLowMax_
private

Definition at line 37 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaLowMin_

const double MaterialBudgetHcalAnalysis::etaLowMin_
private

Definition at line 37 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaMaxP_

const double MaterialBudgetHcalAnalysis::etaMaxP_
private

Definition at line 38 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and MaterialBudgetHcalAnalysis().

◆ etaMidMax_

const double MaterialBudgetHcalAnalysis::etaMidMax_
private

Definition at line 38 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaMidMin_

const double MaterialBudgetHcalAnalysis::etaMidMin_
private

Definition at line 37 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), beginJob(), and MaterialBudgetHcalAnalysis().

◆ etaMinP_

const double MaterialBudgetHcalAnalysis::etaMinP_
private

Definition at line 38 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and MaterialBudgetHcalAnalysis().

◆ labelMBCalo_

const edm::InputTag MaterialBudgetHcalAnalysis::labelMBCalo_
private

Definition at line 39 of file MaterialBudgetHcalAnalysis.cc.

◆ maxEta_

const double MaterialBudgetHcalAnalysis::maxEta_
private

Definition at line 37 of file MaterialBudgetHcalAnalysis.cc.

Referenced by beginJob(), and MaterialBudgetHcalAnalysis().

◆ maxSet2_

const uint32_t MaterialBudgetHcalAnalysis::maxSet2_ = 9
staticprivate

Definition at line 35 of file MaterialBudgetHcalAnalysis.cc.

Referenced by beginJob().

◆ maxSet_

const uint32_t MaterialBudgetHcalAnalysis::maxSet_ = 25
staticprivate

Definition at line 35 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1000_

TProfile2D * MaterialBudgetHcalAnalysis::me1000_[maxSet_]
private

Definition at line 49 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me100_

TProfile* MaterialBudgetHcalAnalysis::me100_[maxSet_]
private

Definition at line 43 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1100_

TProfile2D * MaterialBudgetHcalAnalysis::me1100_[maxSet_]
private

Definition at line 49 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1200_

TH2F* MaterialBudgetHcalAnalysis::me1200_[maxSet_]
private

Definition at line 42 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1300_

TH1F * MaterialBudgetHcalAnalysis::me1300_[maxSet2_]
private

Definition at line 41 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1400_

TH2F * MaterialBudgetHcalAnalysis::me1400_[maxSet2_]
private

Definition at line 42 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1500_

TProfile* MaterialBudgetHcalAnalysis::me1500_[maxSet2_]
private

Definition at line 45 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1600_

TProfile* MaterialBudgetHcalAnalysis::me1600_[maxSet_]
private

Definition at line 46 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1700_

TProfile * MaterialBudgetHcalAnalysis::me1700_[maxSet_]
private

Definition at line 46 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1800_

TProfile * MaterialBudgetHcalAnalysis::me1800_[maxSet_]
private

Definition at line 46 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me1900_

TProfile* MaterialBudgetHcalAnalysis::me1900_[maxSet_]
private

Definition at line 47 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me2000_

TProfile * MaterialBudgetHcalAnalysis::me2000_[maxSet_]
private

Definition at line 47 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me200_

TProfile * MaterialBudgetHcalAnalysis::me200_[maxSet_]
private

Definition at line 43 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me2100_

TProfile * MaterialBudgetHcalAnalysis::me2100_[maxSet_]
private

Definition at line 47 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me2200_

TProfile* MaterialBudgetHcalAnalysis::me2200_[maxSet_]
private

Definition at line 48 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me2300_

TProfile * MaterialBudgetHcalAnalysis::me2300_[maxSet_]
private

Definition at line 48 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me2400_

TProfile * MaterialBudgetHcalAnalysis::me2400_[maxSet_]
private

Definition at line 48 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me300_

TProfile * MaterialBudgetHcalAnalysis::me300_[maxSet_]
private

Definition at line 43 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me400_

TH1F* MaterialBudgetHcalAnalysis::me400_[maxSet_]
private

Definition at line 41 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me500_

TProfile* MaterialBudgetHcalAnalysis::me500_[maxSet_]
private

Definition at line 44 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me600_

TProfile * MaterialBudgetHcalAnalysis::me600_[maxSet_]
private

Definition at line 44 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me700_

TProfile * MaterialBudgetHcalAnalysis::me700_[maxSet_]
private

Definition at line 44 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me800_

TH1F * MaterialBudgetHcalAnalysis::me800_[maxSet_]
private

Definition at line 41 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ me900_

TProfile2D* MaterialBudgetHcalAnalysis::me900_[maxSet_]
private

Definition at line 49 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze(), and beginJob().

◆ tokMBCalo_

const edm::EDGetTokenT<MaterialAccountingCaloCollection> MaterialBudgetHcalAnalysis::tokMBCalo_
private

Definition at line 40 of file MaterialBudgetHcalAnalysis.cc.

Referenced by analyze().