CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
CaloTowersDQMClient Class Reference

#include <CaloTowersDQMClient.h>

Inheritance diagram for CaloTowersDQMClient:
DQMEDHarvester edm::one::EDProducer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources, edm::EndLuminosityBlockProducer > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void beginJob (void) override
 
void beginRun (const edm::Run &run, const edm::EventSetup &c) override
 
 CaloTowersDQMClient (const edm::ParameterSet &)
 
int CaloTowersEndjob (const std::vector< MonitorElement * > &hcalMEs)
 
void dqmEndJob (DQMStore::IBooker &, DQMStore::IGetter &) override
 
 ~CaloTowersDQMClient () override
 
- Public Member Functions inherited from DQMEDHarvester
void beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void beginRun (edm::Run const &, edm::EventSetup const &) override
 
 DQMEDHarvester ()
 
virtual void dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
 
void endJob () final
 
void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
void produce (edm::Event &, edm::EventSetup const &) final
 
 ~DQMEDHarvester () override=default
 
- Public Member Functions inherited from edm::one::EDProducer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources, edm::EndLuminosityBlockProducer >
 EDProducer ()=default
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDProducerBase () override
 
- 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 ()
 
std::vector< edm::ProductResolverIndex > const & putTokenIndexToProductResolverIndex () const
 
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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- 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
 
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
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

edm::ParameterSet conf_
 
bool debug_
 
std::string dirName_
 
std::string dirNameJet_
 
std::string dirNameMET_
 
std::string outputFile_
 
bool verbose_
 

Additional Inherited Members

- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
using ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
 
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- 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

Definition at line 40 of file CaloTowersDQMClient.h.

Constructor & Destructor Documentation

CaloTowersDQMClient::CaloTowersDQMClient ( const edm::ParameterSet iConfig)
explicit

Definition at line 12 of file CaloTowersDQMClient.cc.

References debug_, dirName_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, and verbose_.

12  :conf_(iConfig)
13 {
14 
15  outputFile_ = iConfig.getUntrackedParameter<std::string>("outputFile", "myfile.root");
16  debug_ = false;
17  verbose_ = false;
18  dirName_=iConfig.getParameter<std::string>("DQMDirName");
19 
20  }
T getParameter(std::string const &) const
T getUntrackedParameter(std::string const &, T const &) const
edm::ParameterSet conf_
CaloTowersDQMClient::~CaloTowersDQMClient ( )
override

Definition at line 23 of file CaloTowersDQMClient.cc.

24 {
25 
26 }

Member Function Documentation

void CaloTowersDQMClient::beginJob ( void  )
overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 28 of file CaloTowersDQMClient.cc.

29 {
30 
31 
32 }
void CaloTowersDQMClient::beginRun ( const edm::Run run,
const edm::EventSetup c 
)
override

Definition at line 35 of file CaloTowersDQMClient.cc.

36 {
37 
38 }
int CaloTowersDQMClient::CaloTowersEndjob ( const std::vector< MonitorElement * > &  hcalMEs)

Definition at line 43 of file CaloTowersDQMClient.cc.

References gather_cfg::cout, MillePedeFileConverter_cfg::e, objects.autophobj::float, MonitorElement::getEntries(), MonitorElement::getNbinsX(), MonitorElement::getNbinsY(), mps_fire::i, nevent, funct::pow(), and verbose_.

Referenced by dqmEndJob().

43  {
44 
45  int useAllHistos = 0;
46  MonitorElement* Ntowers_vs_ieta =nullptr;
47  MonitorElement* mapEnergy_N =nullptr, *mapEnergy_E =nullptr, *mapEnergy_H =nullptr, *mapEnergy_EH =nullptr;
48  MonitorElement* occupancy_map =nullptr, *occupancy_vs_ieta =nullptr;
49  for(unsigned int ih=0; ih<hcalMEs.size(); ih++){
50  if( strcmp(hcalMEs[ih]->getName().c_str(), "Ntowers_per_event_vs_ieta") ==0 ){
51  Ntowers_vs_ieta = hcalMEs[ih];
52  }
53  if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_Nentries") ==0 ){
54  mapEnergy_N = hcalMEs[ih];
55  }
56  if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_H") ==0 ){
57  useAllHistos++; mapEnergy_H = hcalMEs[ih];
58  }
59  if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_E") ==0 ){
60  useAllHistos++; mapEnergy_E = hcalMEs[ih];
61  }
62  if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_energy_EH") ==0 ){
63  useAllHistos++; mapEnergy_EH = hcalMEs[ih];
64  }
65  if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_map_occupancy") ==0 ){
66  occupancy_map = hcalMEs[ih];
67  }
68  if( strcmp(hcalMEs[ih]->getName().c_str(), "CaloTowersTask_occupancy_vs_ieta") ==0 ){
69  occupancy_vs_ieta = hcalMEs[ih];
70  }
71 
72  }
73  if( useAllHistos !=0 && useAllHistos !=3 ) return 0;
74 
75  double nevent = mapEnergy_N->getEntries();
76  if(verbose_) std::cout<<"nevent : "<<nevent<<std::endl;
77 
78  // mean number of towers per ieta
79  int nx = Ntowers_vs_ieta->getNbinsX();
80  float cont;
81  float conte;
82  float fev = float(nevent);
83 
84  for (int i = 1; i <= nx; i++) {
85  cont = Ntowers_vs_ieta -> getBinContent(i) / fev ;
86  conte = pow(Ntowers_vs_ieta -> getBinContent(i),0.5) / fev ;
87  Ntowers_vs_ieta -> setBinContent(i,cont);
88  Ntowers_vs_ieta -> setBinError(i,conte);
89  }
90 
91  // mean energies & occupancies evaluation
92 
93  nx = mapEnergy_N->getNbinsX();
94  int ny = mapEnergy_N->getNbinsY();
95  float cnorm;
96  float cnorme;
97  float phi_factor;
98 
99  for (int i = 1; i <= nx; i++) {
100  float sumphi = 0.;
101 
102  for (int j = 1; j <= ny; j++) {
103 
104  // Emean
105  cnorm = mapEnergy_N -> getBinContent(i,j);
106  //Phi histos are not used in the macros
107  if(cnorm > 0.000001 && useAllHistos) {
108 
109  cont = mapEnergy_E -> getBinContent(i,j) / cnorm ;
110  conte = mapEnergy_E -> getBinError(i,j) / cnorm ;
111  mapEnergy_E -> setBinContent(i,j,cont);
112  mapEnergy_E -> setBinError(i,j,conte);
113 
114  cont = mapEnergy_H -> getBinContent(i,j) / cnorm ;
115  conte = mapEnergy_H -> getBinError(i,j) / cnorm ;
116  mapEnergy_H -> setBinContent(i,j,cont);
117  mapEnergy_H -> setBinError(i,j,conte);
118 
119  cont = mapEnergy_EH -> getBinContent(i,j) / cnorm ;
120  conte = mapEnergy_EH -> getBinError(i,j) / cnorm ;
121  mapEnergy_EH -> setBinContent(i,j,cont);
122  mapEnergy_EH -> setBinError(i,j,conte);
123  }
124 
125  // Occupancy (needed for occupancy vs ieta)
126  cont = occupancy_map -> getBinContent(i,j);
127  conte = occupancy_map -> getBinError(i,j);
128  if(fev>0. && cnorm>1.e-30){
129  occupancy_map -> setBinContent(i,j,cont/fev);
130  occupancy_map -> setBinError(i,j,conte/fev);
131  }
132 
133  sumphi += cont;
134 
135  } // end of iphy cycle (j)
136 
137  //Occupancy vs ieta histo is drawn
138  // phi-factor evaluation for occupancy_vs_ieta calculation
139  int ieta = i - 43; //should be the same as int ieta = int(occupancy_vs_ieta->getBinCenter(i));
140 
141  if(ieta >= -20 && ieta <= 20 )
142  {phi_factor = 72.;}
143  else {
144  if(ieta >= 40 || ieta <= -40 ) {phi_factor = 18.;}
145  else
146  phi_factor = 36.;
147  }
148 
149  cnorm = sumphi / phi_factor;
150  cnorme = pow(sumphi,0.5) / phi_factor;
151  if(fev>0. && cnorm>1.e-30){
152  occupancy_vs_ieta->setBinContent(i, cnorm/fev);
153  occupancy_vs_ieta->setBinError(i, cnorme/fev);
154  }
155 
156  } // end of ieta cycle (i)
157 
158  return 1;
159 
160 }
int nevent
Definition: AMPTWrapper.h:74
double getEntries() const
get # of entries
int getNbinsY() const
get # of bins in Y-axis
int getNbinsX() const
get # of bins in X-axis
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
void CaloTowersDQMClient::dqmEndJob ( DQMStore::IBooker ibooker,
DQMStore::IGetter igetter 
)
overridevirtual

Implements DQMEDHarvester.

Definition at line 162 of file CaloTowersDQMClient.cc.

References CaloTowersEndjob(), gather_cfg::cout, DEFINE_FWK_MODULE, dirName_, DQMStore::IGetter::getContents(), DQMStore::IGetter::getSubdirs(), mps_fire::i, DQMStore::IGetter::setCurrentFolder(), and verbose_.

163 {
164 
165  igetter.setCurrentFolder(dirName_);
166  if (verbose_) std::cout << "\nrunClient" << std::endl;
167 
168  std::vector<MonitorElement*> hcalMEs;
169 
170  // Since out folders are fixed to three, we can just go over these three folders
171  // i.e., CaloTowersD/CaloTowersTask, HcalRecHitsD/HcalRecHitTask, NoiseRatesV/NoiseRatesTask.
172  std::vector<std::string> fullPathHLTFolders = igetter.getSubdirs();
173  for(unsigned int i=0;i<fullPathHLTFolders.size();i++) {
174 
175  if (verbose_) std::cout <<"\nfullPath: "<< fullPathHLTFolders[i] << std::endl;
176  igetter.setCurrentFolder(fullPathHLTFolders[i]);
177 
178  std::vector<std::string> fullSubPathHLTFolders = igetter.getSubdirs();
179  for(unsigned int j=0;j<fullSubPathHLTFolders.size();j++) {
180 
181  if (verbose_) std::cout <<"fullSub: "<<fullSubPathHLTFolders[j] << std::endl;
182 
183  if( strcmp(fullSubPathHLTFolders[j].c_str(), "CaloTowersD/CaloTowersTask") ==0 ){
184  hcalMEs = igetter.getContents(fullSubPathHLTFolders[j]);
185  if (verbose_) std::cout <<"hltMES size : "<<hcalMEs.size()<<std::endl;
186  if( !CaloTowersEndjob(hcalMEs) ) std::cout<<"\nError in CaloTowersEndjob!"<<std::endl<<std::endl;
187  }
188 
189  }
190 
191  }
192 
193 }
std::vector< MonitorElement * > getContents(Args &&...args)
Definition: DQMStore.h:305
int CaloTowersEndjob(const std::vector< MonitorElement * > &hcalMEs)
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:349
std::vector< std::string > getSubdirs()
Definition: DQMStore.cc:325

Member Data Documentation

edm::ParameterSet CaloTowersDQMClient::conf_
private

Definition at line 44 of file CaloTowersDQMClient.h.

bool CaloTowersDQMClient::debug_
private

Definition at line 47 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient().

std::string CaloTowersDQMClient::dirName_
private

Definition at line 49 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient(), and dqmEndJob().

std::string CaloTowersDQMClient::dirNameJet_
private

Definition at line 50 of file CaloTowersDQMClient.h.

std::string CaloTowersDQMClient::dirNameMET_
private

Definition at line 51 of file CaloTowersDQMClient.h.

std::string CaloTowersDQMClient::outputFile_
private

Definition at line 43 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient().

bool CaloTowersDQMClient::verbose_
private

Definition at line 46 of file CaloTowersDQMClient.h.

Referenced by CaloTowersDQMClient(), CaloTowersEndjob(), and dqmEndJob().