CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Private Member Functions | Private Attributes
HcalDetDiagNoiseMonitor Class Reference

#include <HcalDetDiagNoiseMonitor.h>

Inheritance diagram for HcalDetDiagNoiseMonitor:
HcalBaseDQMonitor edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

void analyze (edm::Event const &e, edm::EventSetup const &s)
 
void beginRun (const edm::Run &run, const edm::EventSetup &c)
 
void cleanup ()
 
void done ()
 
void endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
 
int GetStatistics ()
 
 HcalDetDiagNoiseMonitor (const edm::ParameterSet &ps)
 
void reset ()
 
void setup ()
 
void UpdateHistos ()
 
 ~HcalDetDiagNoiseMonitor ()
 
- Public Member Functions inherited from HcalBaseDQMonitor
 HcalBaseDQMonitor (const edm::ParameterSet &ps)
 
 HcalBaseDQMonitor ()
 
virtual ~HcalBaseDQMonitor ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
 
std::vector
< ProductHolderIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
 
bool registeredToConsume (ProductHolderIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Private Member Functions

void SaveRates ()
 

Private Attributes

int dataset_seq_number
 
edm::InputTag digiLabel_
 
int FirstOrbit
 
int FirstOrbitLS
 
MonitorElement * HB_RBXmapSpikeCnt
 
MonitorElement * HBM_Rate300
 
MonitorElement * HBM_Rate50
 
MonitorElement * HBP_Rate300
 
MonitorElement * HBP_Rate50
 
MonitorElement * HE_RBXmapSpikeCnt
 
MonitorElement * HEM_Rate300
 
MonitorElement * HEM_Rate50
 
MonitorElement * HEP_Rate300
 
MonitorElement * HEP_Rate50
 
MonitorElement * HO0_Rate300
 
MonitorElement * HO0_Rate50
 
MonitorElement * HO1M_Rate300
 
MonitorElement * HO1M_Rate50
 
MonitorElement * HO1P_Rate300
 
MonitorElement * HO1P_Rate50
 
MonitorElement * HO_RBXmapSpikeCnt
 
MonitorElement * HPDEnergy
 
double HPDthresholdHi
 
double HPDthresholdLo
 
int ievt_
 
int LastOrbit
 
int LastOrbitLS
 
bool LocalRun
 
MonitorElement * meEVT_
 
bool newLS
 
int NoisyEvents
 
MonitorElement * NZeroes
 
std::string OutputFilePath
 
bool Overwrite
 
MonitorElement * PixelMult
 
MonitorElement * RBXEnergy
 
HcalDetDiagNoiseRMSummary * RMSummary
 
int run_number
 
bool setupDone_
 
double SpikeThreshold
 
edm::EDGetTokenT
< HBHEDigiCollection > 		tok_hbhe_
 
edm::EDGetTokenT
< HODigiCollection > 		tok_ho_
 
edm::EDGetTokenT
< L1GlobalTriggerReadoutRecord > 		tok_l1_
 
edm::EDGetTokenT
< FEDRawDataCollection > 		tok_raw_
 
edm::EDGetTokenT
< HcalTBTriggerData > 		tok_tb_
 
MonitorElement * TriggerBx11
 
MonitorElement * TriggerBx12
 
bool UseDB
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- Protected Member Functions inherited from HcalBaseDQMonitor
virtual void beginJob ()
 
virtual void beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
 
virtual void endJob (void)
 
virtual void endRun (const edm::Run &run, const edm::EventSetup &c)
 
void getLogicalMap (const edm::EventSetup &c)
 
bool IsAllowedCalibType ()
 
bool LumiInOrder (int lumisec)
 
void SetupEtaPhiHists (EtaPhiHists &hh, std::string Name, std::string Units)
 
- 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)
 
- Protected Attributes inherited from HcalBaseDQMonitor
std::vector< int > AllowedCalibTypes_
 
int badChannelStatusMask_
 
int currentLS
 
int currenttype_
 
DQMStore * dbe_
 
int debug_
 
bool enableCleanup_
 
bool eventAllowed_
 
bool HBpresent_
 
bool HEpresent_
 
bool HFpresent_
 
bool HOpresent_
 
int ievt_
 
std::map< unsigned int, int > KnownBadCells_
 
int levt_
 
HcalLogicalMap * logicalMap_
 
bool makeDiagnostics_
 
MonitorElement * meIevt_
 
MonitorElement * meLevt_
 
bool mergeRuns_
 
MonitorElement * meTevt_
 
MonitorElement * meTevtHist_
 
bool needLogicalMap_
 
int NLumiBlocks_
 
bool Online_
 
std::string prefixME_
 
MonitorElement * ProblemsCurrentLB
 
MonitorElement * ProblemsVsLB
 
MonitorElement * ProblemsVsLB_HB
 
MonitorElement * ProblemsVsLB_HBHEHF
 
MonitorElement * ProblemsVsLB_HE
 
MonitorElement * ProblemsVsLB_HF
 
MonitorElement * ProblemsVsLB_HO
 
bool skipOutOfOrderLS_
 
std::string subdir_
 
int tevt_
 

Detailed Description

Author
D. Vishnevskiy

Definition at line 36 of file HcalDetDiagNoiseMonitor.h.

Constructor & Destructor Documentation

HcalDetDiagNoiseMonitor::HcalDetDiagNoiseMonitor ( const edm::ParameterSet &  ps)

Definition at line 155 of file HcalDetDiagNoiseMonitor.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), reco_calib_source_client_cfg::OutputFilePath, cmsHarvester::run_number, and AlCaHLTBitMon_QueryRunRegistry::string.

155  {
156 
157  tok_tb_ = consumes<HcalTBTriggerData>(ps.getParameter<edm::InputTag>("hcalTBTriggerDataTag"));
158 
159  ievt_=0;
160  run_number=-1;
161  NoisyEvents=0;
162  LocalRun=false;
163  dataset_seq_number=1;
164  FirstOrbit=FirstOrbitLS=0xFFFFFFFF;
165  LastOrbit=LastOrbitLS=0;
166 
167  Online_ = ps.getUntrackedParameter<bool>("online",false);
168  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns",false);
169  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
170  debug_ = ps.getUntrackedParameter<int>("debug",0);
171  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
172  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
173  prefixME_.append("/");
174  subdir_ = ps.getUntrackedParameter<std::string>("TaskFolder","DetDiagNoiseMonitor_Hcal");
175  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
176  subdir_.append("/");
177  subdir_=prefixME_+subdir_;
178  AllowedCalibTypes_ = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
179  skipOutOfOrderLS_ = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",false);
180  NLumiBlocks_ = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
181  makeDiagnostics_ = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
182 
183  UseDB = ps.getUntrackedParameter<bool> ("UseDB" , false);
184  OutputFilePath = ps.getUntrackedParameter<std::string>("OutputFilePath", "");
185  HPDthresholdHi = ps.getUntrackedParameter<double>("NoiseThresholdHPDhi",49.0);
186  HPDthresholdLo = ps.getUntrackedParameter<double>("NoiseThresholdHPDlo",10.0);
187  SpikeThreshold = ps.getUntrackedParameter<double>("NoiseSpikeThreshold",0.5);
188  Overwrite = ps.getUntrackedParameter<bool> ("Overwrite",true);
189 
190  tok_raw_ = consumes<FEDRawDataCollection>(ps.getUntrackedParameter<edm::InputTag>("RawDataLabel",edm::InputTag("source","")));
191  digiLabel_ = ps.getUntrackedParameter<edm::InputTag>("digiLabel",edm::InputTag("hcalDigis"));
192  tok_l1_ = consumes<L1GlobalTriggerReadoutRecord>(ps.getUntrackedParameter<edm::InputTag>("gtLabel"));
193 
194 
195  tok_hbhe_ = consumes<HBHEDigiCollection>(digiLabel_);
196  tok_ho_ = consumes<HODigiCollection>(digiLabel_);
197 
198  RMSummary = 0;
199  needLogicalMap_=true;
200  setupDone_ = false;
201 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalBaseDQMonitor::AllowedCalibTypes_
std::vector< int > AllowedCalibTypes_
Definition: HcalBaseDQMonitor.h:86
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:114
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HcalDetDiagNoiseMonitor::tok_l1_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_l1_
Definition: HcalDetDiagNoiseMonitor.h:56
HcalDetDiagNoiseMonitor::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalDetDiagNoiseMonitor.h:57
HcalDetDiagNoiseMonitor::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalDetDiagNoiseMonitor.h:58
HcalDetDiagNoiseMonitor::tok_raw_
edm::EDGetTokenT< FEDRawDataCollection > tok_raw_
Definition: HcalDetDiagNoiseMonitor.h:54
HcalDetDiagNoiseMonitor::tok_tb_
edm::EDGetTokenT< HcalTBTriggerData > tok_tb_
Definition: HcalDetDiagNoiseMonitor.h:55
HcalDetDiagNoiseMonitor::~HcalDetDiagNoiseMonitor ( )

Definition at line 636 of file HcalDetDiagNoiseMonitor.cc.

636 {if(LocalRun) UpdateHistos(); SaveRates(); }

Member Function Documentation

void HcalDetDiagNoiseMonitor::analyze ( edm::Event const &  e,
edm::EventSetup const &  s 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 289 of file HcalDetDiagNoiseMonitor.cc.

References adc2fC, HcalBaseDQMonitor::analyze(), gather_cfg::cout, FEDRawData::data(), dbe_, relval_parameters_module::energy, HcalDetDiagNoiseRMEvent::energy, reco::tau::disc::Eta(), edm::Event::getByToken(), hc_Null, i, edm::EventBase::id(), cmsHarvester::index, edm::HandleBase::isValid(), edm::EventBase::luminosityBlock(), bookConverter::max, FEDNumbering::MAXHCALFEDID, HcalFrontEndId::maxRmIndex, FEDNumbering::MINHCALFEDID, HcalDetDiagNoiseRMEvent::n_pix_hi, HcalDetDiagNoiseRMEvent::n_pix_lo, HcalDetDiagNoiseRMEvent::n_zero, colinearityKinematic::Phi, lumiPlot::rawdata, submit::rm, HcalFrontEndId::rmIndex(), edm::EventID::run(), cmsHarvester::run_number, sd, FEDRawData::size(), AlCaHLTBitMon_QueryRunRegistry::string, and subdets.

289  {
290  getLogicalMap(iSetup);
291  if (!IsAllowedCalibType()) return;
292  if (LumiInOrder(iEvent.luminosityBlock())==false) return;
293  HcalBaseDQMonitor::analyze(iEvent, iSetup);
294  bool isNoiseEvent=false;
295  if(!dbe_) return;
296  int orbit=-1111;
297  int bx=-1111;
298 
299  // for local runs
300  edm::Handle<HcalTBTriggerData> trigger_data;
301  iEvent.getByToken(tok_tb_, trigger_data);
302  if(trigger_data.isValid()){
303  if(trigger_data->triggerWord()>1000) isNoiseEvent=true;
304  LocalRun=true;
305  }
306 
307  // We do not want to look at Abort Gap events
308  edm::Handle<FEDRawDataCollection> rawdata;
309  iEvent.getByToken(tok_raw_,rawdata);
310  //checking FEDs for calibration information
311  for(int i=FEDNumbering::MINHCALFEDID;i<=FEDNumbering::MAXHCALFEDID; i++) {
312  const FEDRawData& fedData = rawdata->FEDData(i) ;
313  if ( fedData.size() < 24 ) continue ;
314  orbit= ((const HcalDCCHeader*)(fedData.data()))->getOrbitNumber();
315  bx=((const HcalDCCHeader*)(fedData.data()))->getBunchId();
316  if(((const HcalDCCHeader*)(fedData.data()))->getCalibType()!=hc_Null) return;
317  }
318 
319  // Check GCT trigger bits
320  edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
321  iEvent.getByToken(tok_l1_, gtRecord);
322  if(gtRecord.isValid()){
323  const TechnicalTriggerWord tWord = gtRecord->technicalTriggerWord();
324  if(tWord.at(11) || tWord.at(12)) isNoiseEvent=true;
325  if(tWord.at(11)){ TriggerBx11->Fill(bx);}
326  if(tWord.at(12)){ TriggerBx12->Fill(bx);}
327  }
328 
329  if(!isNoiseEvent) return;
330  if(ievt_==0){ FirstOrbit=orbit; FirstOrbitLS=orbit; newLS=true;}
331  if(LastOrbit <orbit) LastOrbit=orbit;
332  if(FirstOrbit>orbit) FirstOrbit=orbit;
333  if(LastOrbitLS <orbit) LastOrbitLS=orbit;
334  if(FirstOrbitLS>orbit) FirstOrbitLS=orbit;
335  if(newLS){
336  FirstOrbitLS=orbit;
337  newLS=false;
338  }
339 
340  if(!LocalRun){
341  double TIME=(double)(LastOrbit-FirstOrbit)/11223.0;
342  if(TIME>1800.0){
343  UpdateHistos();
344  SaveRates();
345  RMSummary->reset();
346  FirstOrbit=orbit;
347  }
348  }
349 
350  meEVT_->Fill(++ievt_);
351 
352  run_number=iEvent.id().run();
353 
354  HcalDetDiagNoiseRMEvent RMs[HcalFrontEndId::maxRmIndex];
355 
356  edm::Handle<HBHEDigiCollection> hbhe;
357  iEvent.getByToken(tok_hbhe_,hbhe);
358  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
359  double max=-100/*,sum*/,energy=0; int n_zero=0;
360  for(int i=0;i<digi->size();i++){
361  // sum=adc2fC[digi->sample(i).adc()&0xff];
362  if(max<adc2fC[digi->sample(i).adc()&0xff]) max=adc2fC[digi->sample(i).adc()&0xff];
363  if(adc2fC[digi->sample(i).adc()&0xff]==0) n_zero++;
364  }
365  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(digi->id());
366  int index=lmap_entry.rmIndex(); if(index>=HcalFrontEndId::maxRmIndex) continue;
367  RMs[index].n_zero++;
368  if(max>HPDthresholdLo){
369  for(int i=0;i<digi->size();i++) energy+=adc2fC[digi->sample(i).adc()&0xff]-2.5;
370  RMs[index].n_pix_lo++;
371  if(max>HPDthresholdHi){ RMs[index].n_pix_hi++; isNoiseEvent=true;}
372  RMs[index].energy+=energy;
373  }
374  }
375 
376  edm::Handle<HODigiCollection> ho;
377  iEvent.getByToken(tok_ho_,ho);
378  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
379  double max=-100,energy=0; int Eta=digi->id().ieta(); int Phi=digi->id().iphi(); int n_zero=0;
380  for(int i=0;i<digi->size()-1;i++){
381  if(max<adc2fC[digi->sample(i).adc()&0xff]) max=adc2fC[digi->sample(i).adc()&0xff];
382  if(adc2fC[digi->sample(i).adc()&0xff]==0) n_zero++;
383  }
384  if((Eta>=11 && Eta<=15 && Phi>=59 && Phi<=70) || (Eta>=5 && Eta<=10 && Phi>=47 && Phi<=58)){
385  continue; // ignory SiPMs
386  }else{
387  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(digi->id());
388  int index=lmap_entry.rmIndex(); if(index>=HcalFrontEndId::maxRmIndex) continue;
389  RMs[index].n_zero++;
390  if(max>HPDthresholdLo){
391  for(int i=0;i<digi->size();i++) energy+=adc2fC[digi->sample(i).adc()&0xff]-2.5;
392  RMs[index].n_pix_lo++;
393  if(max>HPDthresholdHi){ RMs[index].n_pix_hi++; isNoiseEvent=true;}
394  RMs[index].energy+=energy;
395  }
396  }
397  }
398 
399  NoisyEvents++;
400 
401  // RMs loop
402  for(int i=0;i<HcalFrontEndId::maxRmIndex;i++){
403  if(RMs[i].n_pix_hi>0){
404  HPDEnergy->Fill(RMs[i].energy);
405  RMSummary->rm[i].n_th_hi++;
406  RMSummary->rm[i].n_th_hi_LS++;
407  if(RMs[i].energy>300) RMSummary->rm[i].n_th_300++;
408  if(RMs[i].energy>300) RMSummary->rm[i].n_th_300_LS++;
409  if(RMs[i].n_pix_lo>1) RMSummary->rm[i].n_pix_1++;
410  if(RMs[i].n_pix_lo>8) RMSummary->rm[i].n_pix_8++;
411  RMSummary->rm[i].pix+=RMs[i].n_pix_lo;
412  RMSummary->rm[i].n_pix++;
413  }
414  }
415 
416  // RBX loop
417  for(int sd=0;sd<7;sd++) for(int sect=1;sect<=18;sect++){
418  std::stringstream tempss;
419  tempss << std::setw(2) << std::setfill('0') << sect;
420  std::string rbx= subdets[sd]+tempss.str();
421 
422  double rbx_energy=0;int pix_mult=0; int n_zero=0; bool isValidRBX=false;
423  for(int rm=1;rm<=4;rm++){
424  int index=RMSummary->GetRMindex(rbx,rm);
425  if(index>0 && index<HcalFrontEndId::maxRmIndex){
426  rbx_energy+=RMs[index].energy;
427  pix_mult+=RMs[index].n_pix_lo;
428  n_zero+=RMs[index].n_zero;
429  isValidRBX=true;
430  }
431  }
432  if(isValidRBX){
433  PixelMult->Fill(pix_mult);
434  RBXEnergy->Fill(rbx_energy);
435  NZeroes->Fill(n_zero);
436  }
437  }
438 
439  if((ievt_%100)==0 && debug_>0)
440  std::cout <<ievt_<<"\t"<<NoisyEvents<<std::endl;
441  return;
442 }
i
int i
Definition: DBlmapReader.cc:9
HcalBaseDQMonitor::LumiInOrder
bool LumiInOrder(int lumisec)
Definition: HcalBaseDQMonitor.cc:183
HcalBaseDQMonitor::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:229
HcalDetDiagNoiseRMSummary::rm
HcalDetDiagNoiseRMData rm[HcalFrontEndId::maxRmIndex]
Definition: HcalDetDiagNoiseMonitor.cc:152
edm::SortedCollection< HBHEDataFrame >::const_iterator
std::vector< HBHEDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:82
HcalBaseDQMonitor::getLogicalMap
void getLogicalMap(const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:220
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:114
FEDRawData::size
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
iEvent
int iEvent
Definition: GenABIO.cc:230
HcalDetDiagNoiseRMSummary::GetRMindex
int GetRMindex(const std::string &rbx, int rm)
Definition: HcalDetDiagNoiseMonitor.cc:136
HcalDetDiagNoiseMonitor::tok_l1_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_l1_
Definition: HcalDetDiagNoiseMonitor.h:56
HcalFrontEndId::rmIndex
int rmIndex() const
Definition: HcalFrontEndId.h:20
HcalBaseDQMonitor::logicalMap_
HcalLogicalMap * logicalMap_
Definition: HcalBaseDQMonitor.h:123
HcalDetDiagNoiseMonitor::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalDetDiagNoiseMonitor.h:57
submit.rm
string rm
Definition: submit.py:76
HcalFrontEndId::maxRmIndex
static const int maxRmIndex
Definition: HcalFrontEndId.h:21
TechnicalTriggerWord
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
Definition: L1GlobalTriggerReadoutSetupFwd.h:41
HcalDetDiagNoiseMonitor::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalDetDiagNoiseMonitor.h:58
subdets
char const * subdets[11]
Definition: HcalDetDiagNoiseMonitor.cc:76
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
HcalLogicalMap::getHcalFrontEndId
const HcalFrontEndId getHcalFrontEndId(const DetId &)
Definition: HcalLogicalMap.cc:263
sd
double sd
Definition: CascadeWrapper.h:137
FEDRawData::data
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
HcalDetDiagNoiseMonitor::tok_raw_
edm::EDGetTokenT< FEDRawDataCollection > tok_raw_
Definition: HcalDetDiagNoiseMonitor.h:54
lumiPlot.rawdata
dictionary rawdata
Definition: lumiPlot.py:393
adc2fC
float adc2fC[128]
Definition: HcalDetDiagNoiseMonitor.cc:65
HcalDetDiagNoiseMonitor::tok_tb_
edm::EDGetTokenT< HcalTBTriggerData > tok_tb_
Definition: HcalDetDiagNoiseMonitor.h:55
void HcalDetDiagNoiseMonitor::beginRun ( const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 213 of file HcalDetDiagNoiseMonitor.cc.

References HcalBaseDQMonitor::beginRun(), gather_cfg::cout, reset(), and HcalObjRepresent::setup().

214 {
215  if (debug_>1) std::cout <<"HcalDetDiagNoiseMonitor::beginRun"<<std::endl;
216  HcalBaseDQMonitor::beginRun(run,c);
217 
218  if (tevt_==0) this->setup(); // set up histograms if they have not been created before
219  if (mergeRuns_==false)
220  this->reset();
221 
222  return;
223 
224 }
HcalBaseDQMonitor::beginRun
virtual void beginRun(const edm::Run &run, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:87
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
void HcalDetDiagNoiseMonitor::cleanup ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 203 of file HcalDetDiagNoiseMonitor.cc.

References dbe_, DQMStore::removeContents(), and DQMStore::setCurrentFolder().

203  {
204  if(dbe_){
205  dbe_->setCurrentFolder(subdir_);
206  dbe_->removeContents();
207  dbe_ = 0;
208  }
209 }
DQMStore::removeContents
void removeContents(void)
erase all monitoring elements in current directory (not including subfolders);
Definition: DQMStore.cc:3135
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:667
void HcalDetDiagNoiseMonitor::done ( )

Definition at line 634 of file HcalDetDiagNoiseMonitor.cc.

634 {}
void HcalDetDiagNoiseMonitor::endLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 444 of file HcalDetDiagNoiseMonitor.cc.

References HB_RBX, HE_RBX, HO_RBX, cmsHarvester::index, HcalFrontEndId::maxRmIndex, submit::rm, and sd.

444  {
445 int first_rbx=0,last_rbx=0;
446  //double TIME=(double)(LastOrbitLS-FirstOrbitLS)/11223.0;
447  double TIME=23.0;
448  newLS=true;
449  if(TIME==0) return;
450 
451  for(int sd=0;sd<9;sd++){
452  if(sd==0){ first_rbx=0; last_rbx=18;} //HBM
453  if(sd==1){ first_rbx=18; last_rbx=36;} //HBP
454  if(sd==0 || sd==1){ // update HB plots
455  for(int rbx=first_rbx;rbx<last_rbx;rbx++)for(int rm=1;rm<=4;rm++){
456  int index=RMSummary->GetRMindex(HB_RBX[rbx],rm);
457  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
458  double val=RMSummary->rm[index].n_th_hi_LS/TIME;
459  if(val>SpikeThreshold){
460  HB_RBXmapSpikeCnt->Fill(rm,rbx+1,1);
461  //printf("%s %i %f (%f)\n",HO_RBX[rbx].c_str(),rm,RMSummary->rm[index].n_th_hi_LS/TIME,TIME);
462  }
463  RMSummary->rm[index].reset_LS();
464  }
465  }
466  if(sd==2){ first_rbx=0; last_rbx=18;} //HEM
467  if(sd==3){ first_rbx=18; last_rbx=36;} //HEP
468  if(sd==2 || sd==3){ // update HB plots
469  for(int rbx=first_rbx;rbx<last_rbx;rbx++)for(int rm=1;rm<=4;rm++){
470  int index=RMSummary->GetRMindex(HE_RBX[rbx],rm);
471  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
472  double val=RMSummary->rm[index].n_th_hi_LS/TIME;
473  if(val>SpikeThreshold){
474  HE_RBXmapSpikeCnt->Fill(rm,rbx+1,1);
475  //printf("%s %i %f (%f)\n",HO_RBX[rbx].c_str(),rm,RMSummary->rm[index].n_th_hi_LS/TIME,TIME);
476  }
477  RMSummary->rm[index].reset_LS();
478  }
479  }
480  if(sd==4){ first_rbx=6; last_rbx=12;} //HO1M
481  if(sd==5){ first_rbx=12; last_rbx=24;} //HO0
482  if(sd==6){ first_rbx=24; last_rbx=30;} //HO1P
483  if(sd>3){
484  for(int rbx=first_rbx;rbx<last_rbx;rbx++)for(int rm=1;rm<=4;rm++){
485  int index=RMSummary->GetRMindex(HO_RBX[rbx],rm);
486  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
487  double val=RMSummary->rm[index].n_th_hi_LS/TIME;
488  if(val>SpikeThreshold){
489  HO_RBXmapSpikeCnt->Fill(rm,rbx+1,1);
490  //printf("%s %i %f (%f)\n",HO_RBX[rbx].c_str(),rm,RMSummary->rm[index].n_th_hi_LS/TIME,TIME);
491  }
492  RMSummary->rm[index].reset_LS();
493  }
494  }
495  } //sd=0;sd<9
496 }
HcalDetDiagNoiseRMSummary::rm
HcalDetDiagNoiseRMData rm[HcalFrontEndId::maxRmIndex]
Definition: HcalDetDiagNoiseMonitor.cc:152
HO_RBX
char const * HO_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:83
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:114
HB_RBX
char const * HB_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:77
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
HcalDetDiagNoiseRMSummary::GetRMindex
int GetRMindex(const std::string &rbx, int rm)
Definition: HcalDetDiagNoiseMonitor.cc:136
submit.rm
string rm
Definition: submit.py:76
HcalFrontEndId::maxRmIndex
static const int maxRmIndex
Definition: HcalFrontEndId.h:21
HE_RBX
char const * HE_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:80
sd
double sd
Definition: CascadeWrapper.h:137
int HcalDetDiagNoiseMonitor::GetStatistics ( )
inline

Definition at line 50 of file HcalDetDiagNoiseMonitor.h.

References ievt_.

50 { return ievt_; }
void HcalDetDiagNoiseMonitor::reset ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 210 of file HcalDetDiagNoiseMonitor.cc.

210 {}
void HcalDetDiagNoiseMonitor::SaveRates ( )
private

Definition at line 560 of file HcalDetDiagNoiseMonitor.cc.

References patZpeak::events, HB_RBX, HE_RBX, HO_RBX, cmsHarvester::index, HcalFrontEndId::maxRmIndex, reco_calib_source_client_cfg::OutputFilePath, submit::rm, DTTTrigCorrFirst::run, cmsHarvester::run_number, edmStreamStallGrapher::t, interactiveExample::theFile, cond::rpcobgas::time, and MainPageGenerator::tree.

560  {
561 char RBX[20];
562 int RM;
563 float VAL1,VAL2,VAL3,VAL4,VAL5;
564 char str[500];
565  double TIME=(double)(LastOrbit-FirstOrbit)/11223.0;
566  if(TIME==0) return;
567  if(OutputFilePath.size()>0){
568  if(!Overwrite){
569  sprintf(str,"%sHcalDetDiagNoiseData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
570  }else{
571  sprintf(str,"%sHcalDetDiagNoiseData.root",OutputFilePath.c_str());
572  }
573  TFile *theFile = new TFile(str, "RECREATE");
574  if(!theFile->IsOpen()) return;
575  theFile->cd();
576  sprintf(str,"%d",run_number); TObjString run(str); run.Write("run number");
577  sprintf(str,"%d",ievt_); TObjString events(str); events.Write("Total events processed");
578  sprintf(str,"%d",dataset_seq_number); TObjString dsnum(str); dsnum.Write("Dataset number");
579  Long_t t; t=time(0); strftime(str,30,"%F %T",localtime(&t)); TObjString tm(str); tm.Write("Dataset creation time");
580 
581  TTree *tree =new TTree("HCAL Noise data","HCAL Noise data");
582  if(tree==0) return;
583  tree->Branch("RBX", &RBX, "RBX/C");
584  tree->Branch("rm", &RM, "rm/I");
585  tree->Branch("RATE_50", &VAL1, "RATE_50");
586  tree->Branch("RATE_300", &VAL2, "RATE_300");
587  tree->Branch("RATE_PIX1", &VAL3, "RATE_PIX1");
588  tree->Branch("RATE_PIX8", &VAL4, "RATE_PIX8");
589  tree->Branch("RATE_PIXMEAN", &VAL5, "RATE_PIXMEAN");
590  for(int rbx=0;rbx<36;rbx++) for(int rm=1;rm<=4;rm++){
591  int index=RMSummary->GetRMindex(HB_RBX[rbx],rm);
592  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
593  sprintf(RBX,"%s",HB_RBX[rbx]);
594  RM=rm;
595  VAL1=RMSummary->rm[index].n_th_hi/TIME;
596  VAL2=RMSummary->rm[index].n_th_300/TIME;
597  VAL3=RMSummary->rm[index].n_pix_1/TIME;
598  VAL4=RMSummary->rm[index].n_pix_8/TIME;
599  if(RMSummary->rm[index].n_pix>0)VAL5=RMSummary->rm[index].pix/RMSummary->rm[index].n_pix; else VAL5=0;
600  tree->Fill();
601  }
602  for(int rbx=0;rbx<36;rbx++) for(int rm=1;rm<=4;rm++){
603  int index=RMSummary->GetRMindex(HE_RBX[rbx],rm);
604  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
605  sprintf(RBX,"%s",HE_RBX[rbx]);
606  RM=rm;
607  VAL1=RMSummary->rm[index].n_th_hi/TIME;
608  VAL2=RMSummary->rm[index].n_th_300/TIME;
609  VAL3=RMSummary->rm[index].n_pix_1/TIME;
610  VAL4=RMSummary->rm[index].n_pix_8/TIME;
611  if(RMSummary->rm[index].n_pix>0)VAL5=RMSummary->rm[index].pix/RMSummary->rm[index].n_pix; else VAL5=0;
612  tree->Fill();
613  }
614  for(int rbx=0;rbx<36;rbx++) for(int rm=1;rm<=4;rm++){
615  int index=RMSummary->GetRMindex(HO_RBX[rbx],rm);
616  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
617  sprintf(RBX,"%s",HO_RBX[rbx]);
618  RM=rm;
619  VAL1=RMSummary->rm[index].n_th_hi/TIME;
620  VAL2=RMSummary->rm[index].n_th_300/TIME;
621  VAL3=RMSummary->rm[index].n_pix_1/TIME;
622  VAL4=RMSummary->rm[index].n_pix_8/TIME;
623  if(RMSummary->rm[index].n_pix>0)VAL5=RMSummary->rm[index].pix/RMSummary->rm[index].n_pix; else VAL5=0;
624  tree->Fill();
625  }
626  theFile->Write();
627  theFile->Close();
628  dataset_seq_number++;
629 
630  }
631 }
HcalDetDiagNoiseRMSummary::rm
HcalDetDiagNoiseRMData rm[HcalFrontEndId::maxRmIndex]
Definition: HcalDetDiagNoiseMonitor.cc:152
HO_RBX
char const * HO_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:83
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:114
HB_RBX
char const * HB_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:77
HcalDetDiagNoiseRMSummary::GetRMindex
int GetRMindex(const std::string &rbx, int rm)
Definition: HcalDetDiagNoiseMonitor.cc:136
RBX_struct
Definition: myFilter.cc:51
submit.rm
string rm
Definition: submit.py:76
HcalFrontEndId::maxRmIndex
static const int maxRmIndex
Definition: HcalFrontEndId.h:21
HE_RBX
char const * HE_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:80
patZpeak.events
tuple events
Definition: patZpeak.py:19
void HcalDetDiagNoiseMonitor::setup ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 226 of file HcalDetDiagNoiseMonitor.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookInt(), dbe_, HB_RBX, HE_RBX, HO_RBX, i, mergeVDriftHistosByStation::name, NULL, DQMStore::setCurrentFolder(), HcalBaseDQMonitor::setup(), AlCaHLTBitMon_QueryRunRegistry::string, and indexGen::title.

226  {
227  if (setupDone_)
228  return;
229  setupDone_ = true;
230  // Call base class setup
231  HcalBaseDQMonitor::setup();
232  if (!dbe_) return;
233  RMSummary = new HcalDetDiagNoiseRMSummary();
234 
235  std::string name;
236  if(dbe_!=NULL){
237  dbe_->setCurrentFolder(subdir_);
238  meEVT_ = dbe_->bookInt("HcalNoiseMonitor Event Number");
239  dbe_->setCurrentFolder(subdir_+"Common Plots");
240 
241  name="RBX Pixel multiplicity"; PixelMult = dbe_->book1D(name,name,73,0,73);
242  name="HPD energy"; HPDEnergy = dbe_->book1D(name,name,200,0,2500);
243  name="RBX energy"; RBXEnergy = dbe_->book1D(name,name,200,0,3500);
244  name="Number of zero TS per RBX"; NZeroes = dbe_->book1D(name,name,100,0,100);
245  name="Trigger BX Tbit11"; TriggerBx11 = dbe_->book1D(name,name,4000,0,4000);
246  name="Trigger BX Tbit12"; TriggerBx12 = dbe_->book1D(name,name,4000,0,4000);
247 
248  dbe_->setCurrentFolder(subdir_+"HBHE Plots");
249  name="HBP HPD Noise Rate Pixel above 50fC"; HBP_Rate50 = dbe_->book1D(name,name,73,0,73);
250  name="HBM HPD Noise Rate Pixel above 50fC"; HBM_Rate50 = dbe_->book1D(name,name,73,0,73);
251  name="HEP HPD Noise Rate Pixel above 50fC"; HEP_Rate50 = dbe_->book1D(name,name,73,0,73);
252  name="HEM HPD Noise Rate Pixel above 50fC"; HEM_Rate50 = dbe_->book1D(name,name,73,0,73);
253  name="HBP HPD Noise Rate HPD above 300fC"; HBP_Rate300 = dbe_->book1D(name,name,73,0,73);
254  name="HBM HPD Noise Rate HPD above 300fC"; HBM_Rate300 = dbe_->book1D(name,name,73,0,73);
255  name="HEP HPD Noise Rate HPD above 300fC"; HEP_Rate300 = dbe_->book1D(name,name,73,0,73);
256  name="HEM HPD Noise Rate HPD above 300fC"; HEM_Rate300 = dbe_->book1D(name,name,73,0,73);
257 
258  dbe_->setCurrentFolder(subdir_+"HO Plots");
259  name="HO0 HPD Noise Rate Pixel above 50fC"; HO0_Rate50 = dbe_->book1D(name,name,49,0,49);
260  name="HO1P HPD Noise Rate Pixel above 50fC"; HO1P_Rate50 = dbe_->book1D(name,name,48,0,48);
261  name="HO1M HPD Noise Rate Pixel above 50fC"; HO1M_Rate50 = dbe_->book1D(name,name,48,0,48);
262  name="HO0 HPD Noise Rate HPD above 300fC"; HO0_Rate300 = dbe_->book1D(name,name,48,0,48);
263  name="HO1P HPD Noise Rate HPD abGetRMindexove 300fC"; HO1P_Rate300 = dbe_->book1D(name,name,48,0,48);
264  name="HO1M HPD Noise Rate HPD above 300fC"; HO1M_Rate300 = dbe_->book1D(name,name,48,0,48);
265 
266 
267  dbe_->setCurrentFolder(subdir_+"Noise Spike Plots");
268 
269  name="HB RM Spike Map"; HB_RBXmapSpikeCnt= dbe_->book2D(name,name,4,0.5,4.5,36,0.5,36.5);
270  name="HE RM Spike Map"; HE_RBXmapSpikeCnt= dbe_->book2D(name,name,4,0.5,4.5,36,0.5,36.5);
271  name="HO RM Spike Map"; HO_RBXmapSpikeCnt= dbe_->book2D(name,name,4,0.5,4.5,36,0.5,36.5);
272 
273  std::string title="RM";
274  HB_RBXmapSpikeCnt->setAxisTitle(title);
275  HE_RBXmapSpikeCnt->setAxisTitle(title);
276  HO_RBXmapSpikeCnt->setAxisTitle(title);
277 
278  for(int i=0;i<36;i++){
279  HB_RBXmapSpikeCnt->setBinLabel(i+1,HB_RBX[i],2);
280  HE_RBXmapSpikeCnt->setBinLabel(i+1,HE_RBX[i],2);
281  HO_RBXmapSpikeCnt->setBinLabel(i+1,HO_RBX[i],2);
282  }
283  }
284 
285 
286  return;
287 }
i
int i
Definition: DBlmapReader.cc:9
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:954
HO_RBX
char const * HO_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:83
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:942
NULL
#define NULL
Definition: scimark2.h:8
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:114
HB_RBX
char const * HB_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:77
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HE_RBX
char const * HE_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:80
DQMStore::bookInt
MonitorElement * bookInt(const char *name)
Book int.
Definition: DQMStore.cc:861
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1082
HcalBaseDQMonitor::setup
virtual void setup(void)
Definition: HcalBaseDQMonitor.cc:131
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:969
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:667
void HcalDetDiagNoiseMonitor::UpdateHistos ( )

Definition at line 498 of file HcalDetDiagNoiseMonitor.cc.

References HB_RBX, HE_RBX, HO_RBX, cmsHarvester::index, HcalFrontEndId::maxRmIndex, n, submit::rm, and sd.

498  {
499  int first_rbx=0;
500  double TIME=(double)(LastOrbit-FirstOrbit)/11223.0;
501  if(TIME==0) return;
502  for(int sd=0;sd<9;sd++){
503  if(sd==0){ first_rbx=0; } //HBM
504  if(sd==1){ first_rbx=18;} //HBP
505  if(sd==0 || sd==1){ // update HB plots
506  for(int rbx=0;rbx<18;rbx++)for(int rm=1;rm<=4;rm++){
507  int index=RMSummary->GetRMindex(HB_RBX[rbx+first_rbx],rm);
508  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
509  if(sd==0){
510  HBM_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
511  HBM_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
512  }
513  if(sd==1){
514  HBP_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
515  HBP_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
516  }
517  }
518  }
519  if(sd==2){ first_rbx=0;} //HEM
520  if(sd==3){ first_rbx=18;} //HEP
521  if(sd==2 || sd==3){ // update HB plots
522  for(int rbx=0;rbx<18;rbx++)for(int rm=1;rm<=4;rm++){
523  int index=RMSummary->GetRMindex(HE_RBX[rbx+first_rbx],rm);
524  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
525  if(sd==2){
526  HEM_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
527  HEM_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
528  }
529  if(sd==3){
530  HEP_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
531  HEP_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
532  }
533  }
534  }
535  int n=0;
536  if(sd==4){ first_rbx=6; n=6;} //HO1M
537  if(sd==5){ first_rbx=12;n=12;} //HO0
538  if(sd==6){ first_rbx=24;n=6;} //HO1P
539  if(sd>3){ // update HO plots
540  for(int rbx=0;rbx<n;rbx++)for(int rm=1;rm<=4;rm++){
541  int index=RMSummary->GetRMindex(HO_RBX[rbx+first_rbx],rm);
542  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
543  if(sd==4){
544  HO1M_Rate50->setBinContent(rbx*4*2+rm,RMSummary->rm[index].n_th_hi/TIME);
545  HO1M_Rate300->setBinContent(rbx*4*2+rm,RMSummary->rm[index].n_th_300/TIME);
546  }
547  if(sd==5){
548  HO0_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
549  HO0_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
550  }
551  if(sd==5){
552  HO1P_Rate50->setBinContent(rbx*4*2+rm,RMSummary->rm[index].n_th_hi/TIME);
553  HO1P_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
554  }
555  }
556  }
557  } //sd=0;sd<9
558 }
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:866
HcalDetDiagNoiseRMSummary::rm
HcalDetDiagNoiseRMData rm[HcalFrontEndId::maxRmIndex]
Definition: HcalDetDiagNoiseMonitor.cc:152
HO_RBX
char const * HO_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:83
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:114
HB_RBX
char const * HB_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:77
HcalDetDiagNoiseRMSummary::GetRMindex
int GetRMindex(const std::string &rbx, int rm)
Definition: HcalDetDiagNoiseMonitor.cc:136
submit.rm
string rm
Definition: submit.py:76
HcalFrontEndId::maxRmIndex
static const int maxRmIndex
Definition: HcalFrontEndId.h:21
HE_RBX
char const * HE_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:80
sd
double sd
Definition: CascadeWrapper.h:137

Member Data Documentation

int HcalDetDiagNoiseMonitor::dataset_seq_number
private

Definition at line 63 of file HcalDetDiagNoiseMonitor.h.

edm::InputTag HcalDetDiagNoiseMonitor::digiLabel_
private

Definition at line 52 of file HcalDetDiagNoiseMonitor.h.

int HcalDetDiagNoiseMonitor::FirstOrbit
private

Definition at line 76 of file HcalDetDiagNoiseMonitor.h.

int HcalDetDiagNoiseMonitor::FirstOrbitLS
private

Definition at line 78 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HB_RBXmapSpikeCnt
private

Definition at line 102 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HBM_Rate300
private

Definition at line 91 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HBM_Rate50
private

Definition at line 87 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HBP_Rate300
private

Definition at line 90 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HBP_Rate50
private

Definition at line 86 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HE_RBXmapSpikeCnt
private

Definition at line 103 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HEM_Rate300
private

Definition at line 93 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HEM_Rate50
private

Definition at line 89 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HEP_Rate300
private

Definition at line 92 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HEP_Rate50
private

Definition at line 88 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO0_Rate300
private

Definition at line 98 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO0_Rate50
private

Definition at line 95 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO1M_Rate300
private

Definition at line 100 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO1M_Rate50
private

Definition at line 97 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO1P_Rate300
private

Definition at line 99 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO1P_Rate50
private

Definition at line 96 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HO_RBXmapSpikeCnt
private

Definition at line 104 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::HPDEnergy
private

Definition at line 107 of file HcalDetDiagNoiseMonitor.h.

double HcalDetDiagNoiseMonitor::HPDthresholdHi
private

Definition at line 68 of file HcalDetDiagNoiseMonitor.h.

double HcalDetDiagNoiseMonitor::HPDthresholdLo
private

Definition at line 69 of file HcalDetDiagNoiseMonitor.h.

int HcalDetDiagNoiseMonitor::ievt_
private

Definition at line 74 of file HcalDetDiagNoiseMonitor.h.

Referenced by GetStatistics().

int HcalDetDiagNoiseMonitor::LastOrbit
private

Definition at line 77 of file HcalDetDiagNoiseMonitor.h.

int HcalDetDiagNoiseMonitor::LastOrbitLS
private

Definition at line 79 of file HcalDetDiagNoiseMonitor.h.

bool HcalDetDiagNoiseMonitor::LocalRun
private

Definition at line 82 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::meEVT_
private

Definition at line 84 of file HcalDetDiagNoiseMonitor.h.

bool HcalDetDiagNoiseMonitor::newLS
private

Definition at line 75 of file HcalDetDiagNoiseMonitor.h.

int HcalDetDiagNoiseMonitor::NoisyEvents
private

Definition at line 81 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::NZeroes
private

Definition at line 109 of file HcalDetDiagNoiseMonitor.h.

std::string HcalDetDiagNoiseMonitor::OutputFilePath
private

Definition at line 61 of file HcalDetDiagNoiseMonitor.h.

bool HcalDetDiagNoiseMonitor::Overwrite
private

Definition at line 62 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::PixelMult
private

Definition at line 106 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::RBXEnergy
private

Definition at line 108 of file HcalDetDiagNoiseMonitor.h.

HcalDetDiagNoiseRMSummary* HcalDetDiagNoiseMonitor::RMSummary
private

Definition at line 114 of file HcalDetDiagNoiseMonitor.h.

int HcalDetDiagNoiseMonitor::run_number
private

Definition at line 80 of file HcalDetDiagNoiseMonitor.h.

bool HcalDetDiagNoiseMonitor::setupDone_
private

Definition at line 66 of file HcalDetDiagNoiseMonitor.h.

double HcalDetDiagNoiseMonitor::SpikeThreshold
private

Definition at line 70 of file HcalDetDiagNoiseMonitor.h.

edm::EDGetTokenT<HBHEDigiCollection> HcalDetDiagNoiseMonitor::tok_hbhe_
private

Definition at line 57 of file HcalDetDiagNoiseMonitor.h.

edm::EDGetTokenT<HODigiCollection> HcalDetDiagNoiseMonitor::tok_ho_
private

Definition at line 58 of file HcalDetDiagNoiseMonitor.h.

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> HcalDetDiagNoiseMonitor::tok_l1_
private

Definition at line 56 of file HcalDetDiagNoiseMonitor.h.

edm::EDGetTokenT<FEDRawDataCollection> HcalDetDiagNoiseMonitor::tok_raw_
private

Definition at line 54 of file HcalDetDiagNoiseMonitor.h.

edm::EDGetTokenT<HcalTBTriggerData> HcalDetDiagNoiseMonitor::tok_tb_
private

Definition at line 55 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::TriggerBx11
private

Definition at line 110 of file HcalDetDiagNoiseMonitor.h.

MonitorElement* HcalDetDiagNoiseMonitor::TriggerBx12
private

Definition at line 111 of file HcalDetDiagNoiseMonitor.h.

bool HcalDetDiagNoiseMonitor::UseDB
private

Definition at line 65 of file HcalDetDiagNoiseMonitor.h.

Generated for CMSSW Reference Manual by doxygen 1.8.5