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
HcalDetDiagLEDMonitor Class Reference

#include <HcalDetDiagLEDMonitor.h>

Inheritance diagram for HcalDetDiagLEDMonitor:
HcalBaseDQMonitor HcalBaseDQMonitor edm::EDAnalyzer edm::EDAnalyzer

Public Member Functions

void analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)
 
void analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)
 
void beginRun (const edm::Run &run, const edm::EventSetup &c)
 
void beginRun (const edm::Run &run, const edm::EventSetup &c)
 
void cleanup ()
 
void cleanup ()
 
void done ()
 
void endRun (const edm::Run &run, const edm::EventSetup &c)
 
void fillHistos ()
 
void fillHistos ()
 
int GetStatistics ()
 
int GetStatistics ()
 
 HcalDetDiagLEDMonitor (const edm::ParameterSet &ps)
 
 HcalDetDiagLEDMonitor (const edm::ParameterSet &ps)
 
void reset ()
 
void reset ()
 
void setup ()
 
void setup ()
 
 ~HcalDetDiagLEDMonitor ()
 
 ~HcalDetDiagLEDMonitor ()
 
- Public Member Functions inherited from HcalBaseDQMonitor
 HcalBaseDQMonitor (const edm::ParameterSet &ps)
 
 HcalBaseDQMonitor ()
 
 ~HcalBaseDQMonitor ()
 
- Public Member Functions inherited from edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 

Private Member Functions

void CheckStatus ()
 
void CheckStatus ()
 
void fill_channel_status (std::string subdet, int eta, int phi, int depth, int type, double status)
 
void fill_channel_status (std::string subdet, int eta, int phi, int depth, int type, double status)
 
void fill_energy (std::string subdet, int eta, int phi, int depth, double e, int type)
 
void fill_energy (std::string subdet, int eta, int phi, int depth, double e, int type)
 
double get_energy (std::string subdet, int eta, int phi, int depth, int type)
 
double get_energy (std::string subdet, int eta, int phi, int depth, int type)
 
HcalDetDiagLEDData * GetCalib (std::string sd, int eta, int phi)
 
HcalDetDiagLEDData * GetCalib (std::string sd, int eta, int phi)
 
void LoadReference ()
 
void LoadReference ()
 
void SaveReference ()
 
void SaveReference ()
 

Private Attributes

HcalDetDiagLEDData * calib_data [5][5][72]
 
HcalDetDiagLEDData calib_data [5][5][72]
 
edm::InputTag calibDigiLabel_
 
EtaPhiHists ChannelsLEDEnergy
 
EtaPhiHists * ChannelsLEDEnergy
 
EtaPhiHists ChannelsLEDEnergyRef
 
EtaPhiHists * ChannelsLEDEnergyRef
 
EtaPhiHists ChannelStatusLEDMean
 
EtaPhiHists * ChannelStatusLEDMean
 
EtaPhiHists ChannelStatusLEDRMS
 
EtaPhiHists * ChannelStatusLEDRMS
 
EtaPhiHists ChannelStatusMissingChannels
 
EtaPhiHists * ChannelStatusMissingChannels
 
EtaPhiHists ChannelStatusTimeMean
 
EtaPhiHists * ChannelStatusTimeMean
 
EtaPhiHists ChannelStatusTimeRMS
 
EtaPhiHists * ChannelStatusTimeRMS
 
EtaPhiHists ChannelStatusUnstableChannels
 
EtaPhiHists * ChannelStatusUnstableChannels
 
EtaPhiHists ChannelStatusUnstableLEDsignal
 
EtaPhiHists * ChannelStatusUnstableLEDsignal
 
int dataset_seq_number
 
edm::InputTag digiLabel_
 
HcalElectronicsMap emap
 
MonitorElement * Energy
 
MonitorElement * Energy2Dhbhehf
 
MonitorElement * Energy2Dho
 
MonitorElement * EnergyCorr
 
MonitorElement * EnergyHF
 
MonitorElement * EnergyRMS
 
MonitorElement * EnergyRMSHF
 
HcalLogicalMapGenerator * gen
 
HcalDetDiagLEDData * hb_data [85][72][4]
 
HcalDetDiagLEDData hb_data [85][72][4]
 
MonitorElement * HBMphi
 
MonitorElement * HBPphi
 
HcalDetDiagLEDData * he_data [85][72][4]
 
HcalDetDiagLEDData he_data [85][72][4]
 
MonitorElement * HEMphi
 
MonitorElement * HEPphi
 
HcalDetDiagLEDData * hf_data [85][72][4]
 
HcalDetDiagLEDData hf_data [85][72][4]
 
MonitorElement * HFMphi
 
MonitorElement * HFPphi
 
MonitorElement * HO0phi
 
MonitorElement * HO1Mphi
 
MonitorElement * HO1Pphi
 
MonitorElement * HO2Mphi
 
MonitorElement * HO2Pphi
 
HcalDetDiagLEDData * ho_data [85][72][4]
 
HcalDetDiagLEDData ho_data [85][72][4]
 
int ievt_
 
bool IsReference
 
std::map< unsigned int, int > KnownBadCells_
 
double LEDMeanTreshold
 
double LEDRmsTreshold
 
HcalLogicalMap * lmap
 
MonitorElement * meEVT_
 
MonitorElement * meRUN_
 
std::string OutputFilePath
 
std::string ReferenceData
 
std::string ReferenceRun
 
MonitorElement * RefRun_
 
int run_number
 
MonitorElement * Time
 
MonitorElement * Time2Dhbhehf
 
MonitorElement * Time2Dho
 
MonitorElement * TimeHF
 
MonitorElement * TimeRMS
 
MonitorElement * TimeRMSHF
 
edm::InputTag triggerLabel_
 
bool UseDB
 
std::string XmlFilePath
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzer > WorkerType
 
- 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 endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, 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::EDAnalyzer
CurrentProcessingContext const * currentContext () const
 
- 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_
 
bool makeDiagnostics_
 
MonitorElement * meIevt_
 
MonitorElement * meLevt_
 
bool mergeRuns_
 
MonitorElement * meTevt_
 
MonitorElement * meTevtHist_
 
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

Date:
2010/03/25 11:00:57
Revision:
1.6
Author
D. Vishnevskiy

Definition at line 29 of file HcalDetDiagLEDMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 239 of file HcalDetDiagLEDMonitor.cc.

References HcalBaseDQMonitor::AllowedCalibTypes_, calibDigiLabel_, dataset_seq_number, HcalBaseDQMonitor::debug_, digiLabel_, HcalBaseDQMonitor::enableCleanup_, edm::ParameterSet::getUntrackedParameter(), ievt_, IsReference, LEDMeanTreshold, LEDRmsTreshold, HcalBaseDQMonitor::makeDiagnostics_, HcalBaseDQMonitor::mergeRuns_, HcalBaseDQMonitor::NLumiBlocks_, HcalBaseDQMonitor::Online_, OutputFilePath, HcalBaseDQMonitor::prefixME_, ReferenceData, run_number, HcalBaseDQMonitor::skipOutOfOrderLS_, HcalBaseDQMonitor::subdir_, and XmlFilePath.

239  {
240  ievt_=0;
241  dataset_seq_number=1;
242  run_number=-1;
243  IsReference=false;
244 
245  Online_ = ps.getUntrackedParameter<bool>("online",false);
246  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns",false);
247  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
248  debug_ = ps.getUntrackedParameter<int>("debug",0);
249  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
250  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
251  prefixME_.append("/");
252  subdir_ = ps.getUntrackedParameter<std::string>("TaskFolder","DetDiagLEDMonitor_Hcal");
253  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
254  subdir_.append("/");
255  subdir_=prefixME_+subdir_;
256  AllowedCalibTypes_ = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
257  skipOutOfOrderLS_ = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",false);
258  NLumiBlocks_ = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
259  makeDiagnostics_ = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
260 
261  LEDMeanTreshold = ps.getUntrackedParameter<double>("LEDMeanTreshold" , 0.1);
262  LEDRmsTreshold = ps.getUntrackedParameter<double>("LEDRmsTreshold" , 0.1);
263 
264  ReferenceData = ps.getUntrackedParameter<std::string>("LEDReferenceData" ,"");
265  OutputFilePath = ps.getUntrackedParameter<std::string>("OutputFilePath", "");
266  XmlFilePath = ps.getUntrackedParameter<std::string>("XmlFilePath", "");
267 
268  digiLabel_ = ps.getUntrackedParameter<edm::InputTag>("digiLabel", edm::InputTag("hcalDigis"));
269  calibDigiLabel_ = ps.getUntrackedParameter<edm::InputTag>("calibDigiLabel",edm::InputTag("hcalDigis"));
270 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalBaseDQMonitor::AllowedCalibTypes_
std::vector< int > AllowedCalibTypes_
Definition: HcalBaseDQMonitor.h:82
HcalDetDiagLEDMonitor::~HcalDetDiagLEDMonitor ( )

Definition at line 272 of file HcalDetDiagLEDMonitor.cc.

272 {}
HcalDetDiagLEDMonitor::HcalDetDiagLEDMonitor ( const edm::ParameterSet &  ps)
HcalDetDiagLEDMonitor::~HcalDetDiagLEDMonitor ( )

Member Function Documentation

void HcalDetDiagLEDMonitor::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 377 of file HcalDetDiagLEDMonitor.cc.

References adc2fC, HcalDetDiagLEDData::add_statistics(), HcalBaseDQMonitor::analyze(), calib, calib_data, calibDigiLabel_, CheckStatus(), data, HcalBaseDQMonitor::dbe_, digiLabel_, alignCSCRings::e, eta(), MonitorElement::Fill(), fillHistos(), edm::Event::getByLabel(), edm::Event::getByType(), hb_data, HcalBarrel, HcalEndcap, he_data, hf_data, ho_data, i, edm::EventBase::id(), ievt_, edm::HandleBase::isValid(), meEVT_, meRUN_, phi, edm::EventID::run(), and run_number.

377  {
378 int eta,phi,depth,nTS;
379  if(!dbe_) return;
380  bool LEDEvent=false;
381  bool LocalRun=false;
382  // for local runs
383 
384  edm::Handle<HcalTBTriggerData> trigger_data;
385  iEvent.getByType(trigger_data);
386  if(trigger_data.isValid()){
387  if(trigger_data->triggerWord()==6){ LEDEvent=true;LocalRun=true;}
388  }
389  if(!LocalRun) return;
390  if(!LEDEvent) return;
391 
392  HcalBaseDQMonitor::analyze(iEvent, iSetup);
393  meEVT_->Fill(++ievt_);
394  run_number=iEvent.id().run();
395  meRUN_->Fill(iEvent.id().run());
396 
397  double data[20];
398 
399  edm::Handle<HBHEDigiCollection> hbhe;
400  iEvent.getByLabel(digiLabel_, hbhe);
401  if(hbhe.isValid()) for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
402  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
403  if(digi->id().subdet()==HcalBarrel){
404  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-2.5;
405  hb_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
406  }
407  if(digi->id().subdet()==HcalEndcap){
408  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-2.5;
409  he_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
410  }
411  }
412 
413  edm::Handle<HODigiCollection> ho;
414  iEvent.getByLabel(digiLabel_,ho);
415  if(ho.isValid()) for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
416  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
417  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-2.5;
418  ho_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
419  }
420 
421  edm::Handle<HFDigiCollection> hf;
422  iEvent.getByLabel(digiLabel_,hf);
423  if(hf.isValid()) for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
424  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
425  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-2.5;
426  hf_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
427  }
428 
429  edm::Handle<HcalCalibDigiCollection> calib;
430  iEvent.getByLabel(calibDigiLabel_, calib);
431  if(calib.isValid())for(HcalCalibDigiCollection::const_iterator digi=calib->begin();digi!=calib->end();digi++){
432  if(digi->id().cboxChannel()!=0 || digi->id().hcalSubdet()==0) continue;
433  nTS=digi->size();
434  double e=0;
435  for(int i=0;i<nTS;i++){ data[i]=adc2fC[digi->sample(i).adc()&0xff]; e+=data[i];}
436  if(e<15000) calib_data[digi->id().hcalSubdet()][digi->id().ieta()+2][digi->id().iphi()-1].add_statistics(data,nTS);
437  }
438 
439  if(((ievt_)%500)==0){
440  fillHistos();
441  CheckStatus();
442  }
443  return;
444 }
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
i
int i
Definition: DBlmapReader.cc:9
HcalBaseDQMonitor::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:213
edm::SortedCollection::const_iterator
std::vector< T >::const_iterator const_iterator
Definition: SortedCollection.h:80
HcalDetDiagLEDMonitor::calib_data
HcalDetDiagLEDData * calib_data[5][5][72]
Definition: HcalDetDiagLEDMonitor.h:123
edm::Event::getByType
bool getByType(Handle< PROD > &result) const
Definition: Event.h:398
eta
T eta() const
Definition: Basic3DVectorLD.h:172
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
calib
MVATrainerComputer * calib
Definition: MVATrainer.cc:64
HcalDetDiagLEDMonitor::hb_data
HcalDetDiagLEDData * hb_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:119
HcalDetDiagLEDMonitor::hf_data
HcalDetDiagLEDData * hf_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:122
HcalDetDiagLEDData::add_statistics
void add_statistics(double *data, int nTS)
Definition: HcalDetDiagLEDMonitor.cc:40
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
adc2fC
static const float adc2fC[128]
Definition: HcalDetDiagLEDMonitor.cc:224
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
HcalDetDiagLEDMonitor::ho_data
HcalDetDiagLEDData * ho_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:121
HcalDetDiagLEDMonitor::he_data
HcalDetDiagLEDData * he_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:120
phi
Definition: DDAxes.h:10
void HcalDetDiagLEDMonitor::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
virtual

Reimplemented from HcalBaseDQMonitor.

void HcalDetDiagLEDMonitor::beginRun ( const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 283 of file HcalDetDiagLEDMonitor.cc.

References HcalBaseDQMonitor::beginRun(), gather_cfg::cout, HcalBaseDQMonitor::debug_, edm::EventSetup::get(), HcalCondObjectContainer< Item >::getAllChannels(), reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalChannelStatus::HcalCellMask, HcalChannelStatus::HcalCellOff, i, KnownBadCells_, HcalBaseDQMonitor::mergeRuns_, AlCaHLTBitMon_ParallelJobs::p, edm::ESHandle< class >::product(), reset(), setup(), ntuplemaker::status, and HcalBaseDQMonitor::tevt_.

284 {
285  if (debug_>1) std::cout <<"HcalDetDiagLEDMonitor::beginRun"<<std::endl;
286  HcalBaseDQMonitor::beginRun(run,c);
287 
288  if (tevt_==0) this->setup(); // set up histograms if they have not been created before
289  if (mergeRuns_==false) this->reset();
290 
291  edm::ESHandle<HcalChannelQuality> p;
292  c.get<HcalChannelQualityRcd>().get(p);
293  HcalChannelQuality* chanquality= new HcalChannelQuality(*p.product());
294  std::vector<DetId> mydetids = chanquality->getAllChannels();
295  KnownBadCells_.clear();
296 
297  for (std::vector<DetId>::const_iterator i = mydetids.begin();i!=mydetids.end();++i){
298  if (i->det()!=DetId::Hcal) continue; // not an hcal cell
299  HcalDetId id=HcalDetId(*i);
300  int status=(chanquality->getValues(id))->getValue();
301  if((status & HcalChannelStatus::HcalCellOff) || (status & HcalChannelStatus::HcalCellMask)){
302  KnownBadCells_[id.rawId()]=status;
303  }
304  }
305 
306  return;
307 } // void HcalNDetDiagLEDMonitor::beginRun(...)
i
int i
Definition: DBlmapReader.cc:9
HcalBaseDQMonitor::beginRun
virtual void beginRun(const edm::Run &run, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:83
HcalCondObjectContainer::getAllChannels
std::vector< DetId > getAllChannels() const
Definition: HcalCondObjectContainer.h:300
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
HcalDetDiagLEDMonitor::KnownBadCells_
std::map< unsigned int, int > KnownBadCells_
Definition: HcalDetDiagLEDMonitor.cc:216
reco::JetExtendedAssociation::getValue
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
Definition: JetExtendedAssociation.cc:51
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId) const
Definition: HcalCondObjectContainer.h:113
void HcalDetDiagLEDMonitor::beginRun ( const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

void HcalDetDiagLEDMonitor::CheckStatus ( )
private

Definition at line 1046 of file HcalDetDiagLEDMonitor.cc.

References HcalElectronicsMap::allElectronicsIdPrecision(), HcalDetDiagLEDData::change_status(), ChannelStatusLEDMean, ChannelStatusLEDRMS, ChannelStatusMissingChannels, ChannelStatusTimeMean, ChannelStatusTimeRMS, ChannelStatusUnstableChannels, ChannelStatusUnstableLEDsignal, HcalDetId::depth(), EtaPhiHists::depth, DetId::det(), cond::rpcobgas::detid, emap, eta(), fill_channel_status(), HcalDetDiagLEDData::get_average_time(), HcalDetDiagLEDData::get_overflow(), HcalDetDiagLEDData::get_statistics(), HcalDetDiagLEDData::get_undeflow(), MonitorElement::getMean(), hb_data, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalGenericDetId::HcalGenBarrel, HcalGenericDetId::HcalGenEndcap, HcalGenericDetId::HcalGenForward, HcalGenericDetId::HcalGenOuter, HcalOuter, he_data, hf_data, ho_data, i, HcalDetId::ieta(), ievt_, HcalDetId::iphi(), KnownBadCells_, HcalElectronicsMap::lookup(), phi, plotscripts::rms(), DetId::subdetId(), Time, and TimeHF.

Referenced by analyze(), and endRun().

1046  {
1047  for(int i=0;i<4;i++){
1048  ChannelStatusMissingChannels->depth[i]->Reset();
1049  ChannelStatusUnstableChannels->depth[i]->Reset();
1050  ChannelStatusUnstableLEDsignal->depth[i]->Reset();
1051  ChannelStatusLEDMean->depth[i]->Reset();
1052  ChannelStatusLEDRMS->depth[i]->Reset();
1053  ChannelStatusTimeMean->depth[i]->Reset();
1054  ChannelStatusTimeRMS->depth[i]->Reset();
1055  }
1056 
1057  std::vector <HcalElectronicsId> AllElIds = emap.allElectronicsIdPrecision();
1058  for (std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++) {
1059  DetId detid=emap.lookup(*eid);
1060  if (detid.det()!=DetId::Hcal) continue;
1061  HcalGenericDetId gid(emap.lookup(*eid));
1062  if(!(!(gid.null()) &&
1063  (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
1064  gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap ||
1065  gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
1066  gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
1067  int eta=0,phi=0,depth=0;
1068  if(KnownBadCells_.find(gid.rawId())==KnownBadCells_.end()) continue;
1069 
1070  HcalDetId hid(detid);
1071  eta=hid.ieta();
1072  phi=hid.iphi();
1073  depth=hid.depth();
1074 
1075  double AVE_TIME=Time->getMean();
1076  if(detid.subdetId()==HcalBarrel){
1077  int stat=hb_data[eta+42][phi-1][depth-1].get_statistics()+
1078  hb_data[eta+42][phi-1][depth-1].get_overflow()+hb_data[eta+42][phi-1][depth-1].get_undeflow();
1079  if(stat==0){
1080  fill_channel_status("HB",eta,phi,depth,1,1);
1081  hb_data[eta+42][phi-1][depth-1].change_status(1);
1082  }
1083  if(stat>0 && stat!=(ievt_)){
1084  fill_channel_status("HB",eta,phi,depth,2,(double)stat/(double)(ievt_));
1085  hb_data[eta+42][phi-1][depth-1].change_status(2);
1086  }
1087  if(hb_data[eta+42][phi-1][depth-1].get_statistics()>100){
1088  double ave=0;
1089  double rms=0;
1090  hb_data[eta+42][phi-1][depth-1].get_average_time(&ave,&rms);
1091  if((AVE_TIME-ave)>0.75 || (AVE_TIME-ave)<-0.75){
1092  fill_channel_status("HB",eta,phi,depth,6,AVE_TIME-ave);
1093  hb_data[eta+42][phi-1][depth-1].change_status(8);
1094  }
1095  }
1096  stat=hb_data[eta+42][phi-1][depth-1].get_undeflow();
1097  if(stat>0){
1098  fill_channel_status("HB",eta,phi,depth,3,(double)stat/(double)(ievt_));
1099  hb_data[eta+42][phi-1][depth-1].change_status(4);
1100  }
1101  }
1102  if(detid.subdetId()==HcalEndcap){
1103  int stat=he_data[eta+42][phi-1][depth-1].get_statistics()+
1104  he_data[eta+42][phi-1][depth-1].get_overflow()+he_data[eta+42][phi-1][depth-1].get_undeflow();
1105  if(stat==0){
1106  fill_channel_status("HE",eta,phi,depth,1,1);
1107  he_data[eta+42][phi-1][depth-1].change_status(1);
1108  }
1109  if(stat>0 && stat!=(ievt_)){
1110  fill_channel_status("HE",eta,phi,depth,2,(double)stat/(double)(ievt_));
1111  he_data[eta+42][phi-1][depth-1].change_status(2);
1112  }
1113  if(he_data[eta+42][phi-1][depth-1].get_statistics()>100){
1114  double ave=0;
1115  double rms=0;
1116  he_data[eta+42][phi-1][depth-1].get_average_time(&ave,&rms);
1117  if((AVE_TIME-ave)>0.75 || (AVE_TIME-ave)<-0.75){
1118  fill_channel_status("HE",eta,phi,depth,6,AVE_TIME-ave);
1119  he_data[eta+42][phi-1][depth-1].change_status(8);
1120  }
1121  }
1122  stat=he_data[eta+42][phi-1][depth-1].get_undeflow();
1123  if(stat>0){
1124  fill_channel_status("HE",eta,phi,depth,3,(double)stat/(double)(ievt_));
1125  he_data[eta+42][phi-1][depth-1].change_status(4);
1126  }
1127  }
1128  if(detid.subdetId()==HcalOuter){
1129  int stat=ho_data[eta+42][phi-1][depth-1].get_statistics()+
1130  ho_data[eta+42][phi-1][depth-1].get_overflow()+ho_data[eta+42][phi-1][depth-1].get_undeflow();
1131  if(stat==0){
1132  fill_channel_status("HO",eta,phi,depth,1,1);
1133  ho_data[eta+42][phi-1][depth-1].change_status(1);
1134  }
1135  if(stat>0 && stat!=(ievt_)){
1136  fill_channel_status("HO",eta,phi,depth,2,(double)stat/(double)(ievt_));
1137  ho_data[eta+42][phi-1][depth-1].change_status(2);
1138  }
1139  if(ho_data[eta+42][phi-1][depth-1].get_statistics()>100){
1140  double ave=0;
1141  double rms=0;
1142  ho_data[eta+42][phi-1][depth-1].get_average_time(&ave,&rms);
1143  if((AVE_TIME-ave)>0.75 || (AVE_TIME-ave)<-0.75){
1144  fill_channel_status("HO",eta,phi,depth,6,AVE_TIME-ave);
1145  ho_data[eta+42][phi-1][depth-1].change_status(8);
1146  }
1147  }
1148  stat=ho_data[eta+42][phi-1][depth-1].get_undeflow();
1149  if(stat>0){
1150  fill_channel_status("HO",eta,phi,depth,3,(double)stat/(double)(ievt_));
1151  ho_data[eta+42][phi-1][depth-1].change_status(4);
1152  }
1153  }
1154  if(detid.subdetId()==HcalForward){
1155  AVE_TIME=TimeHF->getMean();
1156  int stat=hf_data[eta+42][phi-1][depth-1].get_statistics()+
1157  hf_data[eta+42][phi-1][depth-1].get_overflow()+hf_data[eta+42][phi-1][depth-1].get_undeflow();
1158  if(stat==0){
1159  fill_channel_status("HF",eta,phi,depth,1,1);
1160  hf_data[eta+42][phi-1][depth-1].change_status(1);
1161  }
1162  if(stat>0 && stat!=(ievt_)){
1163  fill_channel_status("HF",eta,phi,depth,2,(double)stat/(double)(ievt_));
1164  hf_data[eta+42][phi-1][depth-1].change_status(2);
1165  }
1166  if(hf_data[eta+42][phi-1][depth-1].get_statistics()>100){
1167  double ave=0;
1168  double rms=0;
1169  hf_data[eta+42][phi-1][depth-1].get_average_time(&ave,&rms);
1170  if((AVE_TIME-ave)>0.75 || (AVE_TIME-ave)<-0.75){
1171  fill_channel_status("HF",eta,phi,depth,6,AVE_TIME-ave);
1172  hf_data[eta+42][phi-1][depth-1].change_status(8);
1173  }
1174  }
1175  stat=hf_data[eta+42][phi-1][depth-1].get_undeflow();
1176  if(stat>0){
1177  fill_channel_status("HF",eta,phi,depth,3,(double)stat/(double)(ievt_));
1178  hf_data[eta+42][phi-1][depth-1].change_status(4);
1179  }
1180  }
1181  }
1182 }
i
int i
Definition: DBlmapReader.cc:9
HcalElectronicsMap::allElectronicsIdPrecision
std::vector< HcalElectronicsId > allElectronicsIdPrecision() const
Definition: HcalElectronicsMap.cc:122
HcalDetDiagLEDMonitor::fill_channel_status
void fill_channel_status(std::string subdet, int eta, int phi, int depth, int type, double status)
Definition: HcalDetDiagLEDMonitor.cc:1210
eta
T eta() const
Definition: Basic3DVectorLD.h:172
MonitorElement::getMean
double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:669
HcalDetDiagLEDMonitor::hb_data
HcalDetDiagLEDData * hb_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:119
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalDetDiagLEDMonitor::ChannelStatusMissingChannels
EtaPhiHists ChannelStatusMissingChannels
Definition: HcalDetDiagLEDMonitor.h:127
HcalDetDiagLEDMonitor::hf_data
HcalDetDiagLEDData * hf_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:122
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
HcalDetDiagLEDData::get_average_time
bool get_average_time(double *ave, double *rms)
Definition: HcalDetDiagLEDMonitor.cc:65
DetId
Definition: DetId.h:20
HcalDetDiagLEDMonitor::ChannelStatusUnstableLEDsignal
EtaPhiHists ChannelStatusUnstableLEDsignal
Definition: HcalDetDiagLEDMonitor.h:129
HcalDetDiagLEDMonitor::ChannelStatusUnstableChannels
EtaPhiHists ChannelStatusUnstableChannels
Definition: HcalDetDiagLEDMonitor.h:128
HcalDetDiagLEDMonitor::ho_data
HcalDetDiagLEDData * ho_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:121
HcalDetDiagLEDMonitor::KnownBadCells_
std::map< unsigned int, int > KnownBadCells_
Definition: HcalDetDiagLEDMonitor.cc:216
DetId::det
Detector det() const
get the detector field from this detid
Definition: DetId.h:37
HcalElectronicsMap::lookup
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
Definition: HcalElectronicsMap.cc:73
HcalDetDiagLEDMonitor::he_data
HcalDetDiagLEDData * he_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:120
phi
Definition: DDAxes.h:10
void HcalDetDiagLEDMonitor::CheckStatus ( )
private
void HcalDetDiagLEDMonitor::cleanup ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 274 of file HcalDetDiagLEDMonitor.cc.

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

274  {
275  if(dbe_){
276  dbe_->setCurrentFolder(subdir_);
277  dbe_->removeContents();
278  dbe_ = 0;
279  }
280 }
DQMStore::removeContents
void removeContents(void)
erase all monitoring elements in current directory (not including subfolders);
Definition: DQMStore.cc:2564
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
void HcalDetDiagLEDMonitor::cleanup ( )
virtual

Reimplemented from HcalBaseDQMonitor.

void HcalDetDiagLEDMonitor::done ( )
void HcalDetDiagLEDMonitor::endRun ( const edm::Run &  run,
const edm::EventSetup &  c 
)
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 1227 of file HcalDetDiagLEDMonitor.cc.

References CheckStatus(), fillHistos(), ievt_, and SaveReference().

1227  {
1228  if(ievt_>=100){
1229  fillHistos();
1230  CheckStatus();
1231  SaveReference();
1232  }
1233 }
void HcalDetDiagLEDMonitor::fill_channel_status ( std::string  subdet,
int  eta,
int  phi,
int  depth,
int  type,
double  status 
)
private

Definition at line 1210 of file HcalDetDiagLEDMonitor.cc.

References CalcEtaBin(), ChannelStatusLEDMean, ChannelStatusLEDRMS, ChannelStatusMissingChannels, ChannelStatusTimeMean, ChannelStatusTimeRMS, ChannelStatusUnstableChannels, ChannelStatusUnstableLEDsignal, EtaPhiHists::depth, HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

Referenced by CheckStatus().

1210  {
1211  int subdetval=-1;
1212  if (subdet.compare("HB")==0) subdetval=(int)HcalBarrel;
1213  else if (subdet.compare("HE")==0) subdetval=(int)HcalEndcap;
1214  else if (subdet.compare("HO")==0) subdetval=(int)HcalOuter;
1215  else if (subdet.compare("HF")==0) subdetval=(int)HcalForward;
1216  else return;
1217  int ietabin=CalcEtaBin(subdetval, eta, depth)+1;
1218 
1219  if(type==1) ChannelStatusMissingChannels->depth[depth-1] ->setBinContent(ietabin,phi,status);
1220  if(type==2) ChannelStatusUnstableChannels->depth[depth-1] ->setBinContent(ietabin,phi,status);
1221  if(type==3) ChannelStatusUnstableLEDsignal->depth[depth-1]->setBinContent(ietabin,phi,status);
1222  if(type==4) ChannelStatusLEDMean->depth[depth-1] ->setBinContent(ietabin,phi,status);
1223  if(type==5) ChannelStatusLEDRMS->depth[depth-1] ->setBinContent(ietabin,phi,status);
1224  if(type==6) ChannelStatusTimeMean->depth[depth-1] ->setBinContent(ietabin,phi,status);
1225  if(type==7) ChannelStatusTimeRMS->depth[depth-1] ->setBinContent(ietabin,phi,status);
1226 }
type
type
Definition: HCALResponse.h:22
eta
T eta() const
Definition: Basic3DVectorLD.h:172
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalDetDiagLEDMonitor::ChannelStatusMissingChannels
EtaPhiHists ChannelStatusMissingChannels
Definition: HcalDetDiagLEDMonitor.h:127
HcalDetDiagLEDMonitor::ChannelStatusUnstableLEDsignal
EtaPhiHists ChannelStatusUnstableLEDsignal
Definition: HcalDetDiagLEDMonitor.h:129
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
HcalDetDiagLEDMonitor::ChannelStatusUnstableChannels
EtaPhiHists ChannelStatusUnstableChannels
Definition: HcalDetDiagLEDMonitor.h:128
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
phi
Definition: DDAxes.h:10
void HcalDetDiagLEDMonitor::fill_channel_status ( std::string  subdet,
int  eta,
int  phi,
int  depth,
int  type,
double  status 
)
private
void HcalDetDiagLEDMonitor::fill_energy ( std::string  subdet,
int  eta,
int  phi,
int  depth,
double  e,
int  type 
)
private

Definition at line 1183 of file HcalDetDiagLEDMonitor.cc.

References CalcEtaBin(), ChannelsLEDEnergy, ChannelsLEDEnergyRef, EtaPhiHists::depth, HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

Referenced by fillHistos().

1183  {
1184  int subdetval=-1;
1185  if (subdet.compare("HB")==0) subdetval=(int)HcalBarrel;
1186  else if (subdet.compare("HE")==0) subdetval=(int)HcalEndcap;
1187  else if (subdet.compare("HO")==0) subdetval=(int)HcalOuter;
1188  else if (subdet.compare("HF")==0) subdetval=(int)HcalForward;
1189  else return;
1190 
1191  int ietabin=CalcEtaBin(subdetval, eta, depth)+1;
1192  if(type==1) ChannelsLEDEnergy->depth[depth-1] ->setBinContent(ietabin,phi,e);
1193  else if(type==2) ChannelsLEDEnergyRef->depth[depth-1]->setBinContent(ietabin,phi,e);
1194 }
type
type
Definition: HCALResponse.h:22
eta
T eta() const
Definition: Basic3DVectorLD.h:172
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
phi
Definition: DDAxes.h:10
void HcalDetDiagLEDMonitor::fill_energy ( std::string  subdet,
int  eta,
int  phi,
int  depth,
double  e,
int  type 
)
private
void HcalDetDiagLEDMonitor::fillHistos ( )

Definition at line 447 of file HcalDetDiagLEDMonitor.cc.

References abs, HcalElectronicsMap::allElectronicsIdPrecision(), HcalDetId::depth(), DetId::det(), cond::rpcobgas::detid, emap, Energy, Energy2Dhbhehf, Energy2Dho, EnergyCorr, EnergyHF, EnergyRMS, EnergyRMSHF, eta(), MonitorElement::Fill(), fill_energy(), HcalDetDiagLEDData::get_average_led(), HcalDetDiagLEDData::get_average_time(), get_energy(), GetCalib(), HcalLogicalMap::getHcalFrontEndId(), hb_data, HBMphi, HBPphi, DetId::Hcal, HcalGenericDetId::HcalGenBarrel, HcalGenericDetId::HcalGenEndcap, HcalGenericDetId::HcalGenForward, HcalGenericDetId::HcalGenOuter, he_data, HEMphi, HEPphi, hf_data, HFMphi, HFPphi, HO0phi, HO1Mphi, HO1Pphi, HO2Mphi, HO2Pphi, ho_data, i, HcalDetId::ieta(), HcalDetId::iphi(), isSiPM(), j, KnownBadCells_, lmap, HcalElectronicsMap::lookup(), phi, DetId::rawId(), HcalFrontEndId::rbx(), MonitorElement::Reset(), HcalFrontEndId::rm(), plotscripts::rms(), sd, cond::rpcobgas::time, Time, Time2Dhbhehf, Time2Dho, TimeHF, TimeRMS, and TimeRMSHF.

Referenced by analyze(), and endRun().

447  {
448  std::string subdet[4]={"HB","HE","HO","HF"};
449  Energy->Reset();
450  Time->Reset();
451  EnergyRMS->Reset();
452  TimeRMS->Reset();
453  EnergyHF->Reset();
454  TimeHF->Reset();
455  EnergyRMSHF->Reset();
456  TimeRMSHF->Reset();
457  EnergyCorr->Reset();
458  Time2Dhbhehf->Reset();
459  Time2Dho->Reset();
460  Energy2Dhbhehf->Reset();
461  Energy2Dho->Reset();
462  HBPphi->Reset();
463  HBMphi->Reset();
464  HEPphi->Reset();
465  HEMphi->Reset();
466  HFPphi->Reset();
467  HFMphi->Reset();
468  HO0phi->Reset();
469  HO1Pphi->Reset();
470  HO2Pphi->Reset();
471  HO1Mphi->Reset();
472  HO2Mphi->Reset();
473 
474  // HB histograms
475  for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){
476  double T=0,nT=0,E=0,nE=0;
477  for(int depth=1;depth<=2;depth++){
478  if(hb_data[eta+42][phi-1][depth-1].get_statistics()>100){
479  double ave=0;
480  double rms=0;
481  double time=0;
482  double time_rms=0;
483  hb_data[eta+42][phi-1][depth-1].get_average_led(&ave,&rms);
484  hb_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
485  Energy->Fill(ave);
486  if(ave>0)EnergyRMS->Fill(rms/ave);
487  Time->Fill(time);
488  TimeRMS->Fill(time_rms);
489  T+=time; nT++; E+=ave; nE++;
490  if(GetCalib("HB",eta,phi)->get_statistics()>100){
491  double ave_calib=0;
492  double rms_calib=0;
493  GetCalib("HB",eta,phi)->get_average_led(&ave_calib,&rms_calib);
494  fill_energy("HB",eta,phi,depth,ave/ave_calib,1);
495  EnergyCorr->Fill(ave_calib/ave);
496  }
497  }
498  }
499  if(nT>0){Time2Dhbhehf->Fill(eta,phi,T/nT);Energy2Dhbhehf->Fill(eta,phi,E/nE); }
500  }
501  // HE histograms
502  for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
503  double T=0,nT=0,E=0,nE=0;
504  for(int depth=1;depth<=3;depth++){
505  if(he_data[eta+42][phi-1][depth-1].get_statistics()>100){
506  double ave=0;
507  double rms=0;
508  double time=0;
509  double time_rms=0;
510  he_data[eta+42][phi-1][depth-1].get_average_led(&ave,&rms);
511  he_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
512  Energy->Fill(ave);
513  if(ave>0)EnergyRMS->Fill(rms/ave);
514  Time->Fill(time);
515  T+=time; nT++; E+=ave; nE++;
516  TimeRMS->Fill(time_rms);
517  if(GetCalib("HE",eta,phi)->get_statistics()>100){
518  double ave_calib=0;
519  double rms_calib=0;
520  GetCalib("HE",eta,phi)->get_average_led(&ave_calib,&rms_calib);
521  fill_energy("HE",eta,phi,depth,ave/ave_calib,1);
522  EnergyCorr->Fill(ave_calib/ave);
523  }
524  }
525  }
526  if(nT>0 && abs(eta)>16 ){Time2Dhbhehf->Fill(eta,phi,T/nT); Energy2Dhbhehf->Fill(eta,phi,E/nE); }
527  if(nT>0 && abs(eta)>20 ){Time2Dhbhehf->Fill(eta,phi+1,T/nT); Energy2Dhbhehf->Fill(eta,phi+1,E/nE);}
528  }
529  // HF histograms
530  for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
531  double T=0,nT=0,E=0,nE=0;
532  for(int depth=1;depth<=2;depth++){
533  if(hf_data[eta+42][phi-1][depth-1].get_statistics()>100){
534  double ave=0;
535  double rms=0;
536  double time=0;
537  double time_rms=0;
538  hf_data[eta+42][phi-1][depth-1].get_average_led(&ave,&rms);
539  hf_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
540  EnergyHF->Fill(ave);
541  if(ave>0)EnergyRMSHF->Fill(rms/ave);
542  TimeHF->Fill(time);
543  T+=time; nT++; E+=ave; nE++;
544  TimeRMSHF->Fill(time_rms);
545  if(GetCalib("HF",eta,phi)->get_statistics()>100){
546  double ave_calib=0;
547  double rms_calib=0;
548  GetCalib("HF",eta,phi)->get_average_led(&ave_calib,&rms_calib);
549  fill_energy("HF",eta,phi,depth,ave/ave_calib,1);
550  EnergyCorr->Fill(ave_calib/ave);
551  }
552  }
553  }
554  if(nT>0 && abs(eta)>29 ){ Time2Dhbhehf->Fill(eta,phi,T/nT); Time2Dhbhehf->Fill(eta,phi+1,T/nT);}
555  if(nT>0 && abs(eta)>29 ){ Energy2Dhbhehf->Fill(eta,phi,E/nE); Energy2Dhbhehf->Fill(eta,phi+1,E/nE);}
556  }
557  // HO histograms
558  for(int eta=-10;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
559  double T=0,nT=0,E=0,nE=0;
560  if(eta>10 && !isSiPM(eta,phi,4)) continue;
561  for(int depth=4;depth<=4;depth++){
562  if(ho_data[eta+42][phi-1][depth-1].get_statistics()>100){
563  double ave=0;
564  double rms=0;
565  double time=0;
566  double time_rms=0;
567  ho_data[eta+42][phi-1][depth-1].get_average_led(&ave,&rms);
568  ho_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
569  Energy->Fill(ave);
570  if(ave>0)EnergyRMS->Fill(rms/ave);
571  Time->Fill(time);
572  T+=time; nT++; E+=ave; nE++;
573  TimeRMS->Fill(time_rms);
574  if(GetCalib("HO",eta,phi)->get_statistics()>100){
575  double ave_calib=0;
576  double rms_calib=0;
577  GetCalib("HO",eta,phi)->get_average_led(&ave_calib,&rms_calib);
578  fill_energy("HO",eta,phi,depth,ave/ave_calib,1);
579  EnergyCorr->Fill(ave_calib/ave);
580  }
581  }
582  }
583  if(nT>0){ Time2Dho->Fill(eta,phi,T/nT); Energy2Dho->Fill(eta,phi+1,E/nE) ;}
584  }
585 
586  double ave=0.,rms=0.,ave_calib=0.,rms_calib=0.;
587  // HB Ref histograms
588  for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
589  if(hb_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms) && GetCalib("HB",eta,phi)->get_reference(&ave_calib,&rms_calib)){
590  fill_energy("HB",eta,phi,depth,ave/ave_calib,2);
591  }
592  }
593  // HE Ref histograms
594  for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=3;depth++){
595  if(he_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms) && GetCalib("HE",eta,phi)->get_reference(&ave_calib,&rms_calib)){
596  fill_energy("HE",eta,phi,depth,ave/ave_calib,2);
597  }
598  }
599  // HO Ref histograms
600  for(int eta=-10;eta<=15;eta++) for(int phi=1;phi<=72;phi++) for(int depth=4;depth<=4;depth++){
601  if(eta>10 && !isSiPM(eta,phi,4)) continue;
602  if(ho_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms) && GetCalib("HO",eta,phi)->get_reference(&ave_calib,&rms_calib)){
603  fill_energy("HO",eta,phi,depth,ave/ave_calib,2);
604  }
605  }
606  // HF Ref histograms
607  for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
608  if(hf_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms) && GetCalib("HF",eta,phi)->get_reference(&ave_calib,&rms_calib)){
609  fill_energy("HF",eta,phi,depth,ave/ave_calib,2);
610  }
611  }
612 
613  //fill RM histograms: this part is incomplete, will be modefied later
614  double hbp[18][4],nhbp[18][4],hbm[18][4],nhbm[18][4];
615  double hep[18][4],nhep[18][4],hem[18][4],nhem[18][4];
616  double hfp[18][4],nhfp[18][4],hfm[18][4],nhfm[18][4];
617  double ho0[18][4],nho0[18][4];
618  double ho1p[18][4],nho1p[18][4];
619  double ho2p[18][4],nho2p[18][4];
620  double ho1m[18][4],nho1m[18][4];
621  double ho2m[18][4],nho2m[18][4];
622  for(int i=0;i<18;i++) for(int j=0;j<4;j++)
623  hbp[i][j]=nhbp[i][j]=hbm[i][j]=nhbm[i][j]=hep[i][j]=nhep[i][j]=hem[i][j]=nhem[i][j]=hfp[i][j]=nhfp[i][j]=hfm[i][j]=nhfm[i][j]=0;
624  for(int i=0;i<18;i++) for(int j=0;j<4;j++)
625  ho0[i][j]=nho0[i][j]=ho1p[i][j]=nho1p[i][j]=ho2p[i][j]=nho2p[i][j]=ho1m[i][j]=nho1m[i][j]=ho2m[i][j]=nho2m[i][j]=0;
626 
627  std::vector <HcalElectronicsId> AllElIds = emap.allElectronicsIdPrecision();
628  for(std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
629  DetId detid=emap.lookup(*eid);
630  if(detid.det()!=DetId::Hcal) continue;
631  HcalGenericDetId gid(emap.lookup(*eid));
632  if(!(!(gid.null()) &&
633  (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
634  gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap ||
635  gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
636  gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
637  int sd=0,eta=0,phi=0,depth=0;
638  if(gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel) sd=0;
639  else if(gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap) sd=1;
640  else if(gid.genericSubdet()==HcalGenericDetId::HcalGenOuter) sd=2;
641  else if(gid.genericSubdet()==HcalGenericDetId::HcalGenForward)sd=3;
642  HcalDetId hid(detid);
643  if(KnownBadCells_.find(hid.rawId())==KnownBadCells_.end()) continue;
644 
645  eta=hid.ieta();
646  phi=hid.iphi();
647  depth=hid.depth();
648 
649  double ave =get_energy(subdet[sd],eta,phi,depth,1);
650  double ref =get_energy(subdet[sd],eta,phi,depth,2);
651 
652  HcalFrontEndId lmap_entry=lmap->getHcalFrontEndId(hid);
653  int rbx;
654  if(sd==0 || sd==1 || sd==3){
655  sscanf(&(lmap_entry.rbx().c_str())[3],"%d",&rbx);
656  }else{
657  if(abs(eta)<5) sscanf(&(lmap_entry.rbx().c_str())[3],"%d",&rbx);
658  if(abs(eta)>=5) sscanf(&(lmap_entry.rbx().c_str())[4],"%d",&rbx);
659  }
660  if(ave>0 && ref>0){
661  if(sd==0 && eta>0){ hbp[rbx-1][lmap_entry.rm()-1]+=ave/ref; nhbp[rbx-1][lmap_entry.rm()-1]++; }
662  if(sd==0 && eta<0){ hbm[rbx-1][lmap_entry.rm()-1]+=ave/ref; nhbm[rbx-1][lmap_entry.rm()-1]++; }
663  if(sd==1 && eta>0){ hep[rbx-1][lmap_entry.rm()-1]+=ave/ref; nhep[rbx-1][lmap_entry.rm()-1]++; }
664  if(sd==1 && eta<0){ hem[rbx-1][lmap_entry.rm()-1]+=ave/ref; nhem[rbx-1][lmap_entry.rm()-1]++; }
665  if(sd==3 && eta>0){ hfp[rbx-1][lmap_entry.rm()-1]+=ave/ref; nhfp[rbx-1][lmap_entry.rm()-1]++; }
666  if(sd==3 && eta<0){ hfm[rbx-1][lmap_entry.rm()-1]+=ave/ref; nhfm[rbx-1][lmap_entry.rm()-1]++; }
667  if(sd==2 && abs(eta)<5){ ho0[rbx-1][lmap_entry.rm()-1]+=ave/ref; nho0[rbx-1][lmap_entry.rm()-1]++; }
668  if(sd==2 && eta>=5 && eta<=10){ ho1p[rbx-1][lmap_entry.rm()-1]+=ave/ref; nho1p[rbx-1][lmap_entry.rm()-1]++; }
669  if(sd==2 && eta>=11 && eta<=15){ ho2p[rbx-1][lmap_entry.rm()-1]+=ave/ref; nho2p[rbx-1][lmap_entry.rm()-1]++; }
670  if(sd==2 && eta>=-10 && eta<=-5){ ho1m[rbx-1][lmap_entry.rm()-1]+=ave/ref; nho1m[rbx-1][lmap_entry.rm()-1]++; }
671  if(sd==2 && eta>=-15 && eta<=-11){ ho2m[rbx-1][lmap_entry.rm()-1]+=ave/ref; nho2m[rbx-1][lmap_entry.rm()-1]++; }
672  }
673  }
674  for(int i=0;i<18;i++)for(int j=0;j<4;j++){
675  int phi=i*4+j+1;
676  if(nhbp[i][j]>1) HBPphi->Fill(phi+0.5,hbp[i][j]/nhbp[i][j]);
677  if(nhbm[i][j]>1) HBMphi->Fill(phi+0.5,hbm[i][j]/nhbm[i][j]);
678  if(nhep[i][j]>1) HEPphi->Fill(phi+0.5,hep[i][j]/nhep[i][j]);
679  if(nhem[i][j]>1) HEMphi->Fill(phi+0.5,hem[i][j]/nhem[i][j]);
680  }
681  for(int i=0;i<12;i++)for(int j=0;j<3;j++){
682  int phi=i*3+j+1;
683  if(nhfp[i][j]>1) HFPphi->Fill(phi+0.5,hfp[i][j]/nhfp[i][j]);
684  if(nhfm[i][j]>1) HFMphi->Fill(phi+0.5,hfm[i][j]/nhfm[i][j]);
685  }
686  for(int i=0;i<12;i++)for(int j=0;j<4;j++){
687  int phi=i*4+j+1;
688  if(nho0[i][j]>1) HO0phi->Fill(phi+0.5,ho0[i][j]/nho0[i][j]);
689  if(nho1p[i][j]>1) HO1Pphi->Fill(phi+0.5,ho1p[i][j]/nho1p[i][j]);
690  if(nho2p[i][j]>1) HO2Pphi->Fill(phi+0.5,ho2p[i][j]/nho2p[i][j]);
691  if(nho1m[i][j]>1) HO1Mphi->Fill(phi+0.5,ho1m[i][j]/nho1m[i][j]);
692  if(nho2m[i][j]>1) HO2Mphi->Fill(phi+0.5,ho2m[i][j]/nho2m[i][j]);
693  }
694 }
HcalDetDiagLEDMonitor::Time2Dhbhehf
MonitorElement * Time2Dhbhehf
Definition: HcalDetDiagLEDMonitor.h:98
i
int i
Definition: DBlmapReader.cc:9
HcalElectronicsMap::allElectronicsIdPrecision
std::vector< HcalElectronicsId > allElectronicsIdPrecision() const
Definition: HcalElectronicsMap.cc:122
HcalDetDiagLEDData::get_average_led
bool get_average_led(double *ave, double *rms)
Definition: HcalDetDiagLEDMonitor.cc:61
abs
#define abs(x)
Definition: mlp_lapack.h:159
eta
T eta() const
Definition: Basic3DVectorLD.h:172
HcalDetDiagLEDMonitor::GetCalib
HcalDetDiagLEDData * GetCalib(std::string sd, int eta, int phi)
Definition: HcalDetDiagLEDMonitor.h:51
HcalFrontEndId::rbx
std::string rbx() const
Definition: HcalFrontEndId.cc:47
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
isSiPM
bool isSiPM(int ieta, int iphi, int depth)
Definition: HcalEtaPhiHists.h:489
HcalDetDiagLEDMonitor::hb_data
HcalDetDiagLEDData * hb_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:119
HcalFrontEndId::rm
int rm() const
Definition: HcalFrontEndId.h:26
j
int j
Definition: DBlmapReader.cc:9
HcalDetDiagLEDMonitor::hf_data
HcalDetDiagLEDData * hf_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:122
HcalDetDiagLEDData::get_average_time
bool get_average_time(double *ave, double *rms)
Definition: HcalDetDiagLEDMonitor.cc:65
HcalLogicalMap::getHcalFrontEndId
const HcalFrontEndId getHcalFrontEndId(const DetId &)
Definition: HcalLogicalMap.cc:263
DetId
Definition: DetId.h:20
sd
double sd
Definition: CascadeWrapper.h:137
HcalDetDiagLEDMonitor::Energy2Dhbhehf
MonitorElement * Energy2Dhbhehf
Definition: HcalDetDiagLEDMonitor.h:100
HcalDetDiagLEDMonitor::fill_energy
void fill_energy(std::string subdet, int eta, int phi, int depth, double e, int type)
Definition: HcalDetDiagLEDMonitor.cc:1183
HcalDetDiagLEDMonitor::ho_data
HcalDetDiagLEDData * ho_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:121
HcalDetDiagLEDMonitor::KnownBadCells_
std::map< unsigned int, int > KnownBadCells_
Definition: HcalDetDiagLEDMonitor.cc:216
HcalDetDiagLEDMonitor::get_energy
double get_energy(std::string subdet, int eta, int phi, int depth, int type)
Definition: HcalDetDiagLEDMonitor.cc:1196
DetId::det
Detector det() const
get the detector field from this detid
Definition: DetId.h:37
T
long double T
Definition: Basic3DVectorLD.h:59
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:422
HcalElectronicsMap::lookup
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
Definition: HcalElectronicsMap.cc:73
HcalDetDiagLEDMonitor::he_data
HcalDetDiagLEDData * he_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:120
phi
Definition: DDAxes.h:10
void HcalDetDiagLEDMonitor::fillHistos ( )
double HcalDetDiagLEDMonitor::get_energy ( std::string  subdet,
int  eta,
int  phi,
int  depth,
int  type 
)
private

Definition at line 1196 of file HcalDetDiagLEDMonitor.cc.

References CalcEtaBin(), ChannelsLEDEnergy, ChannelsLEDEnergyRef, EtaPhiHists::depth, HcalBarrel, HcalEndcap, HcalForward, and HcalOuter.

Referenced by fillHistos().

1196  {
1197  int subdetval=-1;
1198  if (subdet.compare("HB")==0) subdetval=(int)HcalBarrel;
1199  else if (subdet.compare("HE")==0) subdetval=(int)HcalEndcap;
1200  else if (subdet.compare("HO")==0) subdetval=(int)HcalOuter;
1201  else if (subdet.compare("HF")==0) subdetval=(int)HcalForward;
1202  else return -1.0;
1203 
1204  int ietabin=CalcEtaBin(subdetval, eta, depth)+1;
1205  if(type==1) return ChannelsLEDEnergy->depth[depth-1] ->getBinContent(ietabin, phi);
1206  else if(type==2) return ChannelsLEDEnergyRef->depth[depth-1] ->getBinContent(ietabin,phi);
1207  return -1.0;
1208 }
type
type
Definition: HCALResponse.h:22
eta
T eta() const
Definition: Basic3DVectorLD.h:172
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
phi
Definition: DDAxes.h:10
double HcalDetDiagLEDMonitor::get_energy ( std::string  subdet,
int  eta,
int  phi,
int  depth,
int  type 
)
private
HcalDetDiagLEDData* HcalDetDiagLEDMonitor::GetCalib ( std::string  sd,
int  eta,
int  phi 
)
inlineprivate

Definition at line 51 of file HcalDetDiagLEDMonitor.h.

References abs, calib_data, ETA, and PHI.

Referenced by fillHistos().

51  {
52  int SD=0,ETA=0,PHI=0;
53  if(sd.compare("HB")==0) SD=1;
54  if(sd.compare("HE")==0) SD=2;
55  if(sd.compare("HO")==0) SD=3;
56  if(sd.compare("HF")==0) SD=4;
57  if(SD==1 || SD==2){
58  if(eta>0) ETA=1; else ETA=-1;
59  if(phi==71 ||phi==72 || phi==1 || phi==2) PHI=71; else PHI=((phi-3)/4)*4+3;
60  }else if(SD==3){
61  if(abs(eta)<=4){
62  ETA=0;
63  if(phi==71 ||phi==72 || phi==1 || phi==2 || phi==3 || phi==4) PHI=71; else PHI=((phi-5)/6)*6+5;
64  }else{
65  if(abs(eta)>4 && abs(eta)<=10) ETA=1;
66  if(abs(eta)>10 && abs(eta)<=15) ETA=2;
67  if(eta<0) ETA=-ETA;
68  if(phi==71 ||phi==72 || (phi>=1 && phi<=10)) PHI=71; else PHI=((phi-11)/12)*12+11;
69  }
70  }else if(SD==4){
71  if(eta>0) ETA=1; else ETA=-1;
72  if(phi>=1 && phi<=18) PHI=1;
73  if(phi>=19 && phi<=36) PHI=19;
74  if(phi>=37 && phi<=54) PHI=37;
75  if(phi>=55 && phi<=72) PHI=55;
76  }
77  return calib_data[SD][ETA+2][PHI-1];
78  };
abs
#define abs(x)
Definition: mlp_lapack.h:159
HcalDetDiagLEDMonitor::calib_data
HcalDetDiagLEDData * calib_data[5][5][72]
Definition: HcalDetDiagLEDMonitor.h:123
eta
T eta() const
Definition: Basic3DVectorLD.h:172
ETA
#define ETA
Definition: GenericBenchmark.cc:32
sd
double sd
Definition: CascadeWrapper.h:137
PHI
#define PHI
Definition: GenericBenchmark.cc:33
phi
Definition: DDAxes.h:10
HcalDetDiagLEDData* HcalDetDiagLEDMonitor::GetCalib ( std::string  sd,
int  eta,
int  phi 
)
inlineprivate

Definition at line 129 of file HcalDetDiagLEDMonitor.cc.

References abs, calib_data, ETA, and PHI.

129  {
130  int SD=0,ETA=0,PHI=0;
131  if(sd.compare("HB")==0) SD=1;
132  if(sd.compare("HE")==0) SD=2;
133  if(sd.compare("HO")==0) SD=3;
134  if(sd.compare("HF")==0) SD=4;
135  if(SD==1 || SD==2){
136  if(eta>0) ETA=1; else ETA=-1;
137  if(phi==71 ||phi==72 || phi==1 || phi==2) PHI=71; else PHI=((phi-3)/4)*4+3;
138  }else if(SD==3){
139  if(abs(eta)<=4){
140  ETA=0;
141  if(phi==71 ||phi==72 || phi==1 || phi==2 || phi==3 || phi==4) PHI=71; else PHI=((phi-5)/6)*6+5;
142  }else{
143  if(abs(eta)>4 && abs(eta)<=10) ETA=1;
144  if(abs(eta)>10 && abs(eta)<=15) ETA=2;
145  if(eta<0) ETA=-ETA;
146  if(phi==71 ||phi==72 || (phi>=1 && phi<=10)) PHI=71; else PHI=((phi-11)/12)*12+11;
147  }
148  }else if(SD==4){
149  if(eta>0) ETA=1; else ETA=-1;
150  if(phi>=1 && phi<=18) PHI=1;
151  if(phi>=19 && phi<=36) PHI=19;
152  if(phi>=37 && phi<=54) PHI=37;
153  if(phi>=55 && phi<=72) PHI=55;
154  }
155  return &calib_data[SD][ETA+2][PHI-1];
156  };
abs
#define abs(x)
Definition: mlp_lapack.h:159
HcalDetDiagLEDMonitor::calib_data
HcalDetDiagLEDData * calib_data[5][5][72]
Definition: HcalDetDiagLEDMonitor.h:123
eta
T eta() const
Definition: Basic3DVectorLD.h:172
ETA
#define ETA
Definition: GenericBenchmark.cc:32
sd
double sd
Definition: CascadeWrapper.h:137
PHI
#define PHI
Definition: GenericBenchmark.cc:33
phi
Definition: DDAxes.h:10
int HcalDetDiagLEDMonitor::GetStatistics ( )
inline

Definition at line 41 of file HcalDetDiagLEDMonitor.h.

References ievt_.

41 { return ievt_; }
int HcalDetDiagLEDMonitor::GetStatistics ( )
inline

Definition at line 120 of file HcalDetDiagLEDMonitor.cc.

References ievt_.

120 { return ievt_; }
void HcalDetDiagLEDMonitor::LoadReference ( )
private

Definition at line 1013 of file HcalDetDiagLEDMonitor.cc.

References calib_data, reco::tau::disc::Eta(), f, hb_data, he_data, hf_data, ho_data, IsReference, colinearityKinematic::Phi, ReferenceData, ReferenceRun, plotscripts::rms(), HcalDetDiagLEDData::set_reference(), and lumiQTWidget::t.

Referenced by setup().

1013  {
1014 double led,rms;
1015 int Eta,Phi,Depth;
1016 char subdet[10];
1017 TFile *f;
1018  if(gSystem->AccessPathName(ReferenceData.c_str())) return;
1019  f = new TFile(ReferenceData.c_str(),"READ");
1020  if(!f->IsOpen()) return ;
1021  TObjString *STR=(TObjString *)f->Get("run number");
1022  if(STR){ std::string Ref(STR->String()); ReferenceRun=Ref;}
1023  TTree* t=(TTree*)f->Get("HCAL LED data");
1024  if(!t) return;
1025  t->SetBranchAddress("Subdet", subdet);
1026  t->SetBranchAddress("eta", &Eta);
1027  t->SetBranchAddress("phi", &Phi);
1028  t->SetBranchAddress("depth", &Depth);
1029  t->SetBranchAddress("led", &led);
1030  t->SetBranchAddress("rms", &rms);
1031  for(int ievt=0;ievt<t->GetEntries();ievt++){
1032 
1033  t->GetEntry(ievt);
1034  if(strcmp(subdet,"HB")==0) hb_data[Eta+42][Phi-1][Depth-1].set_reference(led,rms);
1035  if(strcmp(subdet,"HE")==0) he_data[Eta+42][Phi-1][Depth-1].set_reference(led,rms);
1036  if(strcmp(subdet,"HO")==0) ho_data[Eta+42][Phi-1][Depth-1].set_reference(led,rms);
1037  if(strcmp(subdet,"HF")==0) hf_data[Eta+42][Phi-1][Depth-1].set_reference(led,rms);
1038  if(strcmp(subdet,"CALIB_HB")==0) calib_data[1][Eta+2][Phi-1].set_reference(led,rms);
1039  if(strcmp(subdet,"CALIB_HE")==0) calib_data[2][Eta+2][Phi-1].set_reference(led,rms);
1040  if(strcmp(subdet,"CALIB_HO")==0) calib_data[3][Eta+2][Phi-1].set_reference(led,rms);
1041  if(strcmp(subdet,"CALIB_HF")==0) calib_data[4][Eta+2][Phi-1].set_reference(led,rms);
1042  }
1043  f->Close();
1044  IsReference=true;
1045 }
HcalDetDiagLEDMonitor::calib_data
HcalDetDiagLEDData * calib_data[5][5][72]
Definition: HcalDetDiagLEDMonitor.h:123
HcalDetDiagLEDMonitor::hb_data
HcalDetDiagLEDData * hb_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:119
HcalDetDiagLEDMonitor::hf_data
HcalDetDiagLEDData * hf_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:122
f
double f[11][100]
Definition: MuScleFitUtils.cc:80
HcalDetDiagLEDMonitor::ho_data
HcalDetDiagLEDData * ho_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:121
HcalDetDiagLEDMonitor::he_data
HcalDetDiagLEDData * he_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:120
HcalDetDiagLEDData::set_reference
void set_reference(float val, float rms)
Definition: HcalDetDiagLEDMonitor.cc:47
void HcalDetDiagLEDMonitor::LoadReference ( )
private
void HcalDetDiagLEDMonitor::reset ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 281 of file HcalDetDiagLEDMonitor.cc.

Referenced by beginRun().

281 {}
void HcalDetDiagLEDMonitor::reset ( )
virtual

Reimplemented from HcalBaseDQMonitor.

void HcalDetDiagLEDMonitor::SaveReference ( )
private

Definition at line 696 of file HcalDetDiagLEDMonitor.cc.

References HcalElectronicsMap::allElectronicsIdPrecision(), calib_data, dataset_seq_number, HcalDetId::depth(), DetId::det(), cond::rpcobgas::detid, alignCSCRings::e, emap, reco::tau::disc::Eta(), eta(), patZpeak::events, HcalDetDiagLEDData::get_average_led(), HcalDetDiagLEDData::get_average_time(), HcalDetDiagLEDData::get_statistics(), HcalDetDiagLEDData::get_status(), hb_data, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalGenericDetId::HcalGenBarrel, HcalGenericDetId::HcalGenEndcap, HcalGenericDetId::HcalGenForward, HcalGenericDetId::HcalGenOuter, HcalOuter, he_data, hf_data, ho_data, HcalDetId::ieta(), ievt_, HcalDetId::iphi(), HcalElectronicsMap::lookup(), OutputFilePath, phi, colinearityKinematic::Phi, plotscripts::rms(), DTTTrigCorrFirst::run, run_number, sd, step1_ZMM_7Tev::Status, DetId::subdetId(), lumiQTWidget::t, interactiveExample::theFile, cond::rpcobgas::time, Time, MainPageGenerator::tree, xmlFile(), and XmlFilePath.

Referenced by endRun().

696  {
697 double led,rms,Time,time_rms;
698 int Eta,Phi,Depth,Statistic,Status=0;
699 char Subdet[10],str[500];
700  sprintf(str,"%sHcalDetDiagLEDData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
701  TFile *theFile = new TFile(str, "RECREATE");
702  if(!theFile->IsOpen()) return;
703  theFile->cd();
704  sprintf(str,"%d",run_number); TObjString run(str); run.Write("run number");
705  sprintf(str,"%d",ievt_); TObjString events(str); events.Write("Total events processed");
706 
707  TTree *tree =new TTree("HCAL LED data","HCAL LED data");
708  if(tree==0) return;
709  tree->Branch("Subdet", &Subdet, "Subdet/C");
710  tree->Branch("eta", &Eta, "Eta/I");
711  tree->Branch("phi", &Phi, "Phi/I");
712  tree->Branch("depth", &Depth, "Depth/I");
713  tree->Branch("statistic",&Statistic, "Statistic/I");
714  tree->Branch("status", &Status, "Status/I");
715  tree->Branch("led", &led, "led/D");
716  tree->Branch("rms", &rms, "rms/D");
717  tree->Branch("time", &Time, "time/D");
718  tree->Branch("time_rms", &time_rms, "time_rms/D");
719  sprintf(Subdet,"HB");
720  for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
721  if((Statistic=hb_data[eta+42][phi-1][depth-1].get_statistics())>100){
722  Eta=eta; Phi=phi; Depth=depth;
723  Status=hb_data[eta+42][phi-1][depth-1].get_status();
724  hb_data[eta+42][phi-1][depth-1].get_average_led(&led,&rms);
725  hb_data[eta+42][phi-1][depth-1].get_average_time(&Time,&time_rms);
726  tree->Fill();
727  }
728  }
729  sprintf(Subdet,"HE");
730  for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=3;depth++){
731  if((Statistic=he_data[eta+42][phi-1][depth-1].get_statistics())>100){
732  Eta=eta; Phi=phi; Depth=depth;
733  Status=he_data[eta+42][phi-1][depth-1].get_status();
734  he_data[eta+42][phi-1][depth-1].get_average_led(&led,&rms);
735  he_data[eta+42][phi-1][depth-1].get_average_time(&Time,&time_rms);
736  tree->Fill();
737  }
738  }
739  sprintf(Subdet,"HO");
740  for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++) for(int depth=4;depth<=4;depth++){
741  if((Statistic=ho_data[eta+42][phi-1][depth-1].get_statistics())>100){
742  Eta=eta; Phi=phi; Depth=depth;
743  Status=ho_data[eta+42][phi-1][depth-1].get_status();
744  ho_data[eta+42][phi-1][depth-1].get_average_led(&led,&rms);
745  ho_data[eta+42][phi-1][depth-1].get_average_time(&Time,&time_rms);
746  tree->Fill();
747  }
748  }
749  sprintf(Subdet,"HF");
750  for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
751  if((Statistic=hf_data[eta+42][phi-1][depth-1].get_statistics())>100){
752  Eta=eta; Phi=phi; Depth=depth;
753  Status=hf_data[eta+42][phi-1][depth-1].get_status();
754  hf_data[eta+42][phi-1][depth-1].get_average_led(&led,&rms);
755  hf_data[eta+42][phi-1][depth-1].get_average_time(&Time,&time_rms);
756  tree->Fill();
757  }
758  }
759  sprintf(Subdet,"CALIB_HB");
760  for(int eta=-1;eta<=1;eta++) for(int phi=1;phi<=72;phi++){
761  if((calib_data[1][eta+2][phi-1].get_statistics())>100){
762  Eta=eta; Phi=phi; Depth=0;
763  Status=calib_data[1][eta+2][phi-1].get_status();
764  calib_data[1][eta+2][phi-1].get_average_led(&led,&rms);
765  calib_data[1][eta+2][phi-1].get_average_time(&Time,&time_rms);
766  tree->Fill();
767  }
768  }
769  sprintf(Subdet,"CALIB_HE");
770  for(int eta=-1;eta<=1;eta++) for(int phi=1;phi<=72;phi++){
771  if((calib_data[2][eta+2][phi-1].get_statistics())>100){
772  Eta=eta; Phi=phi; Depth=0;
773  Status=calib_data[2][eta+2][phi-1].get_status();
774  calib_data[2][eta+2][phi-1].get_average_led(&led,&rms);
775  calib_data[2][eta+2][phi-1].get_average_time(&Time,&time_rms);
776  tree->Fill();
777  }
778  }
779  sprintf(Subdet,"CALIB_HO");
780  for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
781  if((calib_data[3][eta+2][phi-1].get_statistics())>100){
782  Eta=eta; Phi=phi; Depth=0;
783  Status=calib_data[3][eta+2][phi-1].get_status();
784  calib_data[3][eta+2][phi-1].get_average_led(&led,&rms);
785  calib_data[3][eta+2][phi-1].get_average_time(&Time,&time_rms);
786  tree->Fill();
787  }
788  }
789  sprintf(Subdet,"CALIB_HF");
790  for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
791  if((calib_data[4][eta+2][phi-1].get_statistics())>100){
792  Eta=eta; Phi=phi; Depth=0;
793  Status=calib_data[4][eta+2][phi-1].get_status();
794  calib_data[4][eta+2][phi-1].get_average_led(&led,&rms);
795  calib_data[4][eta+2][phi-1].get_average_time(&Time,&time_rms);
796  tree->Fill();
797  }
798  }
799  theFile->Write();
800  theFile->Close();
801 
802 
803  if(XmlFilePath.size()>0){
804  //create XML file
805  char TIME[40];
806  Long_t t; t=time(0); strftime(TIME,30,"%F %T",localtime(&t));
807 
808  sprintf(str,"HcalDetDiagLED_%i_%i.xml",run_number,dataset_seq_number);
809  std::string xmlName=str;
810  ofstream xmlFile;
811  xmlFile.open(xmlName.c_str());
812 
813  xmlFile<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"yes\"?>\n";
814  xmlFile<<"<ROOT>\n";
815  xmlFile<<" <HEADER>\n";
816  xmlFile<<" <HINTS mode='only-det-root'/>\n";
817  xmlFile<<" <TYPE>\n";
818  xmlFile<<" <EXTENSION_TABLE_NAME>HCAL_DETMON_LED_LASER_V1</EXTENSION_TABLE_NAME>\n";
819  xmlFile<<" <NAME>HCAL LED [local]</NAME>\n";
820  xmlFile<<" </TYPE>\n";
821  xmlFile<<" <!-- run details -->\n";
822  xmlFile<<" <RUN>\n";
823  xmlFile<<" <RUN_TYPE>LOCAL-RUN</RUN_TYPE>\n";
824  xmlFile<<" <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
825  xmlFile<<" <RUN_BEGIN_TIMESTAMP>2009-01-01 00:00:00</RUN_BEGIN_TIMESTAMP>\n";
826  xmlFile<<" <COMMENT_DESCRIPTION>hcal LED data</COMMENT_DESCRIPTION>\n";
827  xmlFile<<" <LOCATION>P5</LOCATION>\n";
828  xmlFile<<" <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
829  xmlFile<<" </RUN>\n";
830  xmlFile<<" </HEADER>\n";
831  xmlFile<<" <DATA_SET>\n";
832  xmlFile<<" <!-- optional dataset metadata -->\n\n";
833  xmlFile<<" <SET_NUMBER>"<<dataset_seq_number<<"</SET_NUMBER>\n";
834  xmlFile<<" <SET_BEGIN_TIMESTAMP>2009-01-01 00:00:00</SET_BEGIN_TIMESTAMP>\n";
835  xmlFile<<" <SET_END_TIMESTAMP>2009-01-01 00:00:00</SET_END_TIMESTAMP>\n";
836  xmlFile<<" <NUMBER_OF_EVENTS_IN_SET>"<<ievt_<<"</NUMBER_OF_EVENTS_IN_SET>\n";
837  xmlFile<<" <COMMENT_DESCRIPTION>Automatic DQM output</COMMENT_DESCRIPTION>\n";
838  xmlFile<<" <DATA_FILE_NAME>"<< xmlName <<"</DATA_FILE_NAME>\n";
839  xmlFile<<" <IMAGE_FILE_NAME>data plot url or file path</IMAGE_FILE_NAME>\n";
840  xmlFile<<" <!-- who and when created this dataset-->\n\n";
841  xmlFile<<" <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
842  xmlFile<<" <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
843  xmlFile<<" <!-- version (string) and subversion (number) -->\n";
844  xmlFile<<" <!-- fields are used to read data back from the database -->\n\n";
845  xmlFile<<" <VERSION>"<<run_number<<dataset_seq_number<<"</VERSION>\n";
846  xmlFile<<" <SUBVERSION>1</SUBVERSION>\n";
847  xmlFile<<" <!-- Assign predefined dataset attributes -->\n\n";
848  xmlFile<<" <PREDEFINED_ATTRIBUTES>\n";
849  xmlFile<<" <ATTRIBUTE>\n";
850  xmlFile<<" <NAME>HCAL Dataset Status</NAME>\n";
851  xmlFile<<" <VALUE>VALID</VALUE>\n";
852  xmlFile<<" </ATTRIBUTE>\n";
853  xmlFile<<" </PREDEFINED_ATTRIBUTES>\n";
854  xmlFile<<" <!-- multiple data block records -->\n\n";
855 
856  std::vector <HcalElectronicsId> AllElIds = emap.allElectronicsIdPrecision();
857  for(std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
858  DetId detid=emap.lookup(*eid);
859  if (detid.det()!=DetId::Hcal) continue;
860  HcalGenericDetId gid(emap.lookup(*eid));
861  if(!(!(gid.null()) &&
862  (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
863  gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap ||
864  gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
865  gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
866  int eta,phi,depth;
867  std::string subdet="";
868  HcalDetId hid(detid);
869  eta=hid.ieta();
870  phi=hid.iphi();
871  depth=hid.depth();
872 
873  double e=0,e_rms=0,t=0,t_rms=0;
874  if(detid.subdetId()==HcalBarrel){
875  subdet="HB";
876  Statistic=hb_data[eta+42][phi-1][depth-1].get_statistics();
877  Status =hb_data[eta+42][phi-1][depth-1].get_status();
878  hb_data[eta+42][phi-1][depth-1].get_average_led(&e,&e_rms);
879  hb_data[eta+42][phi-1][depth-1].get_average_time(&t,&t_rms);
880  }else if(detid.subdetId()==HcalEndcap){
881  subdet="HE";
882  Statistic=he_data[eta+42][phi-1][depth-1].get_statistics();
883  Status =he_data[eta+42][phi-1][depth-1].get_status();
884  he_data[eta+42][phi-1][depth-1].get_average_led(&e,&e_rms);
885  he_data[eta+42][phi-1][depth-1].get_average_time(&t,&t_rms);
886  }else if(detid.subdetId()==HcalForward){
887  subdet="HF";
888  Statistic=hf_data[eta+42][phi-1][depth-1].get_statistics();
889  Status =hf_data[eta+42][phi-1][depth-1].get_status();
890  hf_data[eta+42][phi-1][depth-1].get_average_led(&e,&e_rms);
891  hf_data[eta+42][phi-1][depth-1].get_average_time(&t,&t_rms);
892  }else if(detid.subdetId()==HcalOuter){
893  subdet="HO";
894  Statistic=ho_data[eta+42][phi-1][depth-1].get_statistics();
895  Status =ho_data[eta+42][phi-1][depth-1].get_status();
896  ho_data[eta+42][phi-1][depth-1].get_average_led(&e,&e_rms);
897  ho_data[eta+42][phi-1][depth-1].get_average_time(&t,&t_rms);
898  }else continue;
899  xmlFile<<" <DATA>\n";
900  xmlFile<<" <NUMBER_OF_EVENTS_USED>"<<Statistic<<"</NUMBER_OF_EVENTS_USED>\n";
901  xmlFile<<" <SIGNAL_MEAN>"<<e<<"</SIGNAL_MEAN>\n";
902  xmlFile<<" <SIGNAL_RMS>"<<e_rms<<"</SIGNAL_RMS>\n";
903  xmlFile<<" <TIME_MEAN>"<<t<<"</TIME_MEAN>\n";
904  xmlFile<<" <TIME_RMS>"<<t_rms<<"</TIME_RMS>\n";
905  xmlFile<<" <CHANNEL_STATUS_WORD>"<<Status<<"</CHANNEL_STATUS_WORD>\n";
906  xmlFile<<" <CHANNEL_OBJECTNAME>HcalDetId</CHANNEL_OBJECTNAME>\n";
907  xmlFile<<" <SUBDET>"<<subdet<<"</SUBDET>\n";
908  xmlFile<<" <IETA>"<<eta<<"</IETA>\n";
909  xmlFile<<" <IPHI>"<<phi<<"</IPHI>\n";
910  xmlFile<<" <DEPTH>"<<depth<<"</DEPTH>\n";
911  xmlFile<<" <TYPE>0</TYPE>\n";
912  xmlFile<<" </DATA>\n";
913  }
915  xmlFile<<" </DATA_SET>\n";
916  xmlFile<<"</ROOT>\n";
917  xmlFile.close();
918 
919  //create CALIB XML file
920  sprintf(str,"HcalDetDiagLEDCalib_%i_%i.xml",run_number,dataset_seq_number);
921  std::string xmlNameCalib=str;
922  ofstream xmlFileCalib;
923  xmlFileCalib.open(xmlNameCalib.c_str());
924 
925  xmlFileCalib<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"yes\"?>\n";
926  xmlFileCalib<<"<ROOT>\n";
927  xmlFileCalib<<" <HEADER>\n";
928  xmlFileCalib<<" <HINTS mode='only-det-root'/>\n";
929  xmlFileCalib<<" <TYPE>\n";
930  xmlFileCalib<<" <EXTENSION_TABLE_NAME>HCAL_DETMON_LED_LASER_V1</EXTENSION_TABLE_NAME>\n";
931  xmlFileCalib<<" <NAME>HCAL LED CALIB [local]</NAME>\n";
932  xmlFileCalib<<" </TYPE>\n";
933  xmlFileCalib<<" <!-- run details -->\n";
934  xmlFileCalib<<" <RUN>\n";
935  xmlFileCalib<<" <RUN_TYPE>LOCAL-RUN</RUN_TYPE>\n";
936  xmlFileCalib<<" <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
937  xmlFileCalib<<" <RUN_BEGIN_TIMESTAMP>2009-01-01 00:00:00</RUN_BEGIN_TIMESTAMP>\n";
938  xmlFileCalib<<" <COMMENT_DESCRIPTION>hcal LED CALIB data</COMMENT_DESCRIPTION>\n";
939  xmlFileCalib<<" <LOCATION>P5</LOCATION>\n";
940  xmlFileCalib<<" <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
941  xmlFileCalib<<" </RUN>\n";
942  xmlFileCalib<<" </HEADER>\n";
943  xmlFileCalib<<" <DATA_SET>\n";
944  xmlFileCalib<<" <!-- optional dataset metadata -->\n\n";
945  xmlFileCalib<<" <SET_NUMBER>"<<dataset_seq_number<<"</SET_NUMBER>\n";
946  xmlFileCalib<<" <SET_BEGIN_TIMESTAMP>2009-01-01 00:00:00</SET_BEGIN_TIMESTAMP>\n";
947  xmlFileCalib<<" <SET_END_TIMESTAMP>2009-01-01 00:00:00</SET_END_TIMESTAMP>\n";
948  xmlFileCalib<<" <NUMBER_OF_EVENTS_IN_SET>"<<ievt_<<"</NUMBER_OF_EVENTS_IN_SET>\n";
949  xmlFileCalib<<" <COMMENT_DESCRIPTION>Automatic DQM output</COMMENT_DESCRIPTION>\n";
950  xmlFileCalib<<" <DATA_FILE_NAME>"<< xmlNameCalib <<"</DATA_FILE_NAME>\n";
951  xmlFileCalib<<" <IMAGE_FILE_NAME>data plot url or file path</IMAGE_FILE_NAME>\n";
952  xmlFileCalib<<" <!-- who and when created this dataset-->\n\n";
953  xmlFileCalib<<" <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
954  xmlFileCalib<<" <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
955  xmlFileCalib<<" <!-- version (string) and subversion (number) -->\n";
956  xmlFileCalib<<" <!-- fields are used to read data back from the database -->\n\n";
957  xmlFileCalib<<" <VERSION>"<<run_number<<dataset_seq_number<<"</VERSION>\n";
958  xmlFileCalib<<" <SUBVERSION>1</SUBVERSION>\n";
959  xmlFileCalib<<" <!-- Assign predefined dataset attributes -->\n\n";
960  xmlFileCalib<<" <PREDEFINED_ATTRIBUTES>\n";
961  xmlFileCalib<<" <ATTRIBUTE>\n";
962  xmlFileCalib<<" <NAME>HCAL Dataset Status</NAME>\n";
963  xmlFileCalib<<" <VALUE>VALID</VALUE>\n";
964  xmlFileCalib<<" </ATTRIBUTE>\n";
965  xmlFileCalib<<" </PREDEFINED_ATTRIBUTES>\n";
966  xmlFileCalib<<" <!-- multiple data block records -->\n\n";
967 
968  for(int sd=1;sd<=4;sd++) for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
969  std::string subdet="";
970  if(sd==1) subdet="HB";
971  if(sd==2) subdet="HE";
972  if(sd==3) subdet="HO";
973  if(sd==4) subdet="HF";
974  if((calib_data[sd][eta+2][phi-1].get_statistics())>100){
975  double e=0,e_rms=0,t=0,t_rms=0;
976  Status=calib_data[sd][eta+2][phi-1].get_status();
977  Statistic=calib_data[sd][eta+2][phi-1].get_statistics();
978  calib_data[sd][eta+2][phi-1].get_average_led(&e,&e_rms);
979  calib_data[sd][eta+2][phi-1].get_average_time(&t,&t_rms);
980  xmlFileCalib<<" <DATA>\n";
981  xmlFileCalib<<" <NUMBER_OF_EVENTS_USED>"<<Statistic<<"</NUMBER_OF_EVENTS_USED>\n";
982  xmlFileCalib<<" <SIGNAL_MEAN>"<<e<<"</SIGNAL_MEAN>\n";
983  xmlFileCalib<<" <SIGNAL_RMS>"<<e_rms<<"</SIGNAL_RMS>\n";
984  xmlFileCalib<<" <TIME_MEAN>"<<t<<"</TIME_MEAN>\n";
985  xmlFileCalib<<" <TIME_RMS>"<<t_rms<<"</TIME_RMS>\n";
986  xmlFileCalib<<" <CHANNEL_STATUS_WORD>"<<Status<<"</CHANNEL_STATUS_WORD>\n";
987  xmlFileCalib<<" <CHANNEL_OBJECTNAME>HcalDetId</CHANNEL_OBJECTNAME>\n";
988  xmlFileCalib<<" <SUBDET>"<<subdet<<"</SUBDET>\n";
989  xmlFileCalib<<" <IETA>"<<eta<<"</IETA>\n";
990  xmlFileCalib<<" <IPHI>"<<phi<<"</IPHI>\n";
991  xmlFileCalib<<" <DEPTH>"<<0<<"</DEPTH>\n";
992  xmlFileCalib<<" <TYPE>0</TYPE>\n";
993  xmlFileCalib<<" </DATA>\n";
994  }
995  }
997  xmlFileCalib<<" </DATA_SET>\n";
998  xmlFileCalib<<"</ROOT>\n";
999  xmlFileCalib.close();
1000 
1001  sprintf(str,"zip %s.zip %s %s",xmlName.c_str(),xmlName.c_str(),xmlNameCalib.c_str());
1002  system(str);
1003  sprintf(str,"rm -f %s %s",xmlName.c_str(),xmlNameCalib.c_str());
1004  system(str);
1005  sprintf(str,"mv -f %s.zip %s",xmlName.c_str(),XmlFilePath.c_str());
1006  system(str);
1007  }
1008 
1009 
1010  dataset_seq_number++;
1011 }
HcalElectronicsMap::allElectronicsIdPrecision
std::vector< HcalElectronicsId > allElectronicsIdPrecision() const
Definition: HcalElectronicsMap.cc:122
HcalDetDiagLEDData::get_average_led
bool get_average_led(double *ave, double *rms)
Definition: HcalDetDiagLEDMonitor.cc:61
HcalDetDiagLEDMonitor::calib_data
HcalDetDiagLEDData * calib_data[5][5][72]
Definition: HcalDetDiagLEDMonitor.h:123
eta
T eta() const
Definition: Basic3DVectorLD.h:172
HcalDetDiagLEDMonitor::hb_data
HcalDetDiagLEDData * hb_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:119
HcalDetDiagLEDMonitor::hf_data
HcalDetDiagLEDData * hf_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:122
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
HcalDetDiagLEDData::get_average_time
bool get_average_time(double *ave, double *rms)
Definition: HcalDetDiagLEDMonitor.cc:65
step1_ZMM_7Tev.Status
tuple Status
Definition: step1_ZMM_7Tev.py:112
DetId
Definition: DetId.h:20
sd
double sd
Definition: CascadeWrapper.h:137
Statistic
Statistic
Definition: trackSplitPlot.h:33
xmlFile
bool xmlFile(const std::string fParam)
Definition: hcalCalibrationsCopy.cc:214
patZpeak.events
tuple events
Definition: patZpeak.py:19
HcalDetDiagLEDMonitor::ho_data
HcalDetDiagLEDData * ho_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:121
DetId::det
Detector det() const
get the detector field from this detid
Definition: DetId.h:37
HcalElectronicsMap::lookup
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
Definition: HcalElectronicsMap.cc:73
HcalDetDiagLEDMonitor::he_data
HcalDetDiagLEDData * he_data[85][72][4]
Definition: HcalDetDiagLEDMonitor.h:120
phi
Definition: DDAxes.h:10
void HcalDetDiagLEDMonitor::SaveReference ( )
private
void HcalDetDiagLEDMonitor::setup ( void  )
virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 309 of file HcalDetDiagLEDMonitor.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookInt(), DQMStore::bookString(), ChannelsLEDEnergy, ChannelsLEDEnergyRef, ChannelStatusLEDMean, ChannelStatusLEDRMS, ChannelStatusMissingChannels, ChannelStatusTimeMean, ChannelStatusTimeRMS, ChannelStatusUnstableChannels, ChannelStatusUnstableLEDsignal, HcalLogicalMapGenerator::createMap(), HcalBaseDQMonitor::dbe_, emap, Energy, Energy2Dhbhehf, Energy2Dho, EnergyCorr, EnergyHF, EnergyRMS, EnergyRMSHF, gen, HcalLogicalMap::generateHcalElectronicsMap(), HBMphi, HBPphi, HEMphi, HEPphi, HFMphi, HFPphi, HO0phi, HO1Mphi, HO1Pphi, HO2Mphi, HO2Pphi, lmap, LoadReference(), meEVT_, meRUN_, mergeVDriftHistosByStation::name, ReferenceRun, RefRun_, DQMStore::setCurrentFolder(), EtaPhiHists::setup(), HcalBaseDQMonitor::setup(), HcalBaseDQMonitor::subdir_, Time, Time2Dhbhehf, Time2Dho, TimeHF, TimeRMS, and TimeRMSHF.

Referenced by beginRun().

309  {
310  // Call base class setup
311  HcalBaseDQMonitor::setup();
312  if (!dbe_) return;
313 
314  std::string name;
315  dbe_->setCurrentFolder(subdir_);
316  meEVT_ = dbe_->bookInt("HcalDetDiagLEDMonitor Event Number");
317  meRUN_ = dbe_->bookInt("HcalDetDiagLEDMonitor Run Number");
318  ReferenceRun="UNKNOWN";
319  LoadReference();
320  dbe_->setCurrentFolder(subdir_);
321  RefRun_= dbe_->bookString("HcalDetDiagLEDMonitor Reference Run",ReferenceRun);
322  dbe_->setCurrentFolder(subdir_+"Summary Plots");
323 
324  name="HBHEHO LED Energy Distribution"; Energy = dbe_->book1D(name,name,200,0,3000);
325  name="HBHEHO LED Timing Distribution"; Time = dbe_->book1D(name,name,200,0,10);
326  name="HBHEHO LED Energy RMS_div_Energy Distribution";EnergyRMS = dbe_->book1D(name,name,200,0,0.2);
327  name="HBHEHO LED Timing RMS Distribution"; TimeRMS = dbe_->book1D(name,name,200,0,0.4);
328  name="HF LED Energy Distribution"; EnergyHF = dbe_->book1D(name,name,200,0,3000);
329  name="HF LED Timing Distribution"; TimeHF = dbe_->book1D(name,name,200,0,10);
330  name="HF LED Energy RMS_div_Energy Distribution"; EnergyRMSHF = dbe_->book1D(name,name,200,0,0.5);
331  name="HF LED Timing RMS Distribution"; TimeRMSHF = dbe_->book1D(name,name,200,0,0.4);
332  name="LED Energy Corr(PinDiod) Distribution"; EnergyCorr = dbe_->book1D(name,name,200,0,10);
333  name="LED Timing HBHEHF"; Time2Dhbhehf = dbe_->book2D(name,name,87,-43,43,74,0,73);
334  name="LED Timing HO"; Time2Dho = dbe_->book2D(name,name,33,-16,16,74,0,73);
335  name="LED Energy HBHEHF"; Energy2Dhbhehf = dbe_->book2D(name,name,87,-43,43,74,0,73);
336  name="LED Energy HO"; Energy2Dho = dbe_->book2D(name,name,33,-16,16,74,0,73);
337 
338  name="HBP Average over HPD LED Ref"; HBPphi = dbe_->book2D(name,name,180,1,73,400,0,2);
339  name="HBM Average over HPD LED Ref"; HBMphi = dbe_->book2D(name,name,180,1,73,400,0,2);
340  name="HEP Average over HPD LED Ref"; HEPphi = dbe_->book2D(name,name,180,1,73,400,0,2);
341  name="HEM Average over HPD LED Ref"; HEMphi = dbe_->book2D(name,name,180,1,73,400,0,2);
342  name="HFP Average over RM LED Ref"; HFPphi = dbe_->book2D(name,name,180,1,37,400,0,2);
343  name="HFM Average over RM LED Ref"; HFMphi = dbe_->book2D(name,name,180,1,37,400,0,2);
344  name="HO0 Average over HPD LED Ref"; HO0phi = dbe_->book2D(name,name,180,1,49,400,0,2);
345  name="HO1P Average over HPD LED Ref"; HO1Pphi= dbe_->book2D(name,name,180,1,49,400,0,2);
346  name="HO2P Average over HPD LED Ref"; HO2Pphi= dbe_->book2D(name,name,180,1,49,400,0,2);
347  name="HO1M Average over HPD LED Ref"; HO1Mphi= dbe_->book2D(name,name,180,1,49,400,0,2);
348  name="HO2M Average over HPD LED Ref"; HO2Mphi= dbe_->book2D(name,name,180,1,49,400,0,2);
349 
350  ChannelsLEDEnergy = new EtaPhiHists();
351  ChannelsLEDEnergy->setup(dbe_," Channel LED Energy");
352  ChannelsLEDEnergyRef = new EtaPhiHists();
353  ChannelsLEDEnergyRef->setup(dbe_," Channel LED Energy Reference");
354 
355  dbe_->setCurrentFolder(subdir_+"channel status");
356  ChannelStatusMissingChannels = new EtaPhiHists();
357  ChannelStatusMissingChannels->setup(dbe_," Missing Channels");
358  ChannelStatusUnstableChannels = new EtaPhiHists();
359  ChannelStatusUnstableChannels->setup(dbe_," Unstable Channels");
360  ChannelStatusUnstableLEDsignal = new EtaPhiHists();
361  ChannelStatusUnstableLEDsignal->setup(dbe_," Unstable LED");
362  ChannelStatusLEDMean = new EtaPhiHists();
363  ChannelStatusLEDMean->setup(dbe_," LED Mean");
364  ChannelStatusLEDRMS = new EtaPhiHists();
365  ChannelStatusLEDRMS->setup(dbe_," LED RMS");
366  ChannelStatusTimeMean = new EtaPhiHists();
367  ChannelStatusTimeMean->setup(dbe_," Time Mean");
368  ChannelStatusTimeRMS = new EtaPhiHists();
369  ChannelStatusTimeRMS->setup(dbe_," Time RMS");
370 
371  gen=new HcalLogicalMapGenerator();
372  lmap =new HcalLogicalMap(gen->createMap());
373  emap=lmap->generateHcalElectronicsMap();
374  return;
375 }
HcalDetDiagLEDMonitor::Time2Dhbhehf
MonitorElement * Time2Dhbhehf
Definition: HcalDetDiagLEDMonitor.h:98
HcalDetDiagLEDMonitor::gen
HcalLogicalMapGenerator * gen
Definition: HcalDetDiagLEDMonitor.h:43
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:717
HcalLogicalMapGenerator::createMap
HcalLogicalMap createMap(unsigned int mapIOV=5)
Definition: HcalLogicalMapGenerator.cc:38
EtaPhiHists::setup
void setup(DQMStore *&m_dbe, std::string Name, std::string Units="")
Definition: HcalEtaPhiHists.h:24
HcalLogicalMap::generateHcalElectronicsMap
HcalElectronicsMap generateHcalElectronicsMap()
Definition: HcalLogicalMap.cc:156
HcalDetDiagLEDMonitor::ChannelStatusMissingChannels
EtaPhiHists ChannelStatusMissingChannels
Definition: HcalDetDiagLEDMonitor.h:127
DQMStore::bookString
MonitorElement * bookString(const char *name, const char *value)
Book string.
Definition: DQMStore.cc:683
HcalDetDiagLEDMonitor::ChannelStatusUnstableLEDsignal
EtaPhiHists ChannelStatusUnstableLEDsignal
Definition: HcalDetDiagLEDMonitor.h:129
HcalDetDiagLEDMonitor::Energy2Dhbhehf
MonitorElement * Energy2Dhbhehf
Definition: HcalDetDiagLEDMonitor.h:100
HcalDetDiagLEDMonitor::ChannelStatusUnstableChannels
EtaPhiHists ChannelStatusUnstableChannels
Definition: HcalDetDiagLEDMonitor.h:128
DQMStore::bookInt
MonitorElement * bookInt(const char *name)
Book int.
Definition: DQMStore.cc:624
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:845
HcalBaseDQMonitor::setup
virtual void setup(void)
Definition: HcalBaseDQMonitor.cc:127
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
void HcalDetDiagLEDMonitor::setup ( )
virtual

Reimplemented from HcalBaseDQMonitor.

Member Data Documentation

HcalDetDiagLEDData* HcalDetDiagLEDMonitor::calib_data[5][5][72]
private

Definition at line 123 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), GetCalib(), LoadReference(), and SaveReference().

HcalDetDiagLEDData HcalDetDiagLEDMonitor::calib_data[5][5][72]
private

Definition at line 201 of file HcalDetDiagLEDMonitor.cc.

edm::InputTag HcalDetDiagLEDMonitor::calibDigiLabel_
private

Definition at line 137 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), and HcalDetDiagLEDMonitor().

EtaPhiHists HcalDetDiagLEDMonitor::ChannelsLEDEnergy
private

Definition at line 125 of file HcalDetDiagLEDMonitor.h.

Referenced by fill_energy(), get_energy(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelsLEDEnergy
private

Definition at line 203 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelsLEDEnergyRef
private

Definition at line 126 of file HcalDetDiagLEDMonitor.h.

Referenced by fill_energy(), get_energy(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelsLEDEnergyRef
private

Definition at line 204 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusLEDMean
private

Definition at line 130 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusLEDMean
private

Definition at line 208 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusLEDRMS
private

Definition at line 131 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusLEDRMS
private

Definition at line 209 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusMissingChannels
private

Definition at line 127 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusMissingChannels
private

Definition at line 205 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusTimeMean
private

Definition at line 132 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusTimeMean
private

Definition at line 210 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusTimeRMS
private

Definition at line 133 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusTimeRMS
private

Definition at line 211 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusUnstableChannels
private

Definition at line 128 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusUnstableChannels
private

Definition at line 206 of file HcalDetDiagLEDMonitor.cc.

EtaPhiHists HcalDetDiagLEDMonitor::ChannelStatusUnstableLEDsignal
private

Definition at line 129 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fill_channel_status(), and setup().

EtaPhiHists* HcalDetDiagLEDMonitor::ChannelStatusUnstableLEDsignal
private

Definition at line 207 of file HcalDetDiagLEDMonitor.cc.

int HcalDetDiagLEDMonitor::dataset_seq_number
private

Definition at line 81 of file HcalDetDiagLEDMonitor.h.

Referenced by HcalDetDiagLEDMonitor(), and SaveReference().

edm::InputTag HcalDetDiagLEDMonitor::digiLabel_
private

Definition at line 135 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), and HcalDetDiagLEDMonitor().

HcalElectronicsMap HcalDetDiagLEDMonitor::emap
private

Definition at line 44 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fillHistos(), SaveReference(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::Energy
private

Definition at line 94 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::Energy2Dhbhehf
private

Definition at line 100 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::Energy2Dho
private

Definition at line 101 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::EnergyCorr
private

Definition at line 106 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::EnergyHF
private

Definition at line 96 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::EnergyRMS
private

Definition at line 102 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::EnergyRMSHF
private

Definition at line 104 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

HcalLogicalMapGenerator * HcalDetDiagLEDMonitor::gen
private

Definition at line 43 of file HcalDetDiagLEDMonitor.h.

Referenced by setup().

HcalDetDiagLEDData* HcalDetDiagLEDMonitor::hb_data[85][72][4]
private

Definition at line 119 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), CheckStatus(), fillHistos(), LoadReference(), and SaveReference().

HcalDetDiagLEDData HcalDetDiagLEDMonitor::hb_data[85][72][4]
private

Definition at line 197 of file HcalDetDiagLEDMonitor.cc.

MonitorElement * HcalDetDiagLEDMonitor::HBMphi
private

Definition at line 108 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HBPphi
private

Definition at line 107 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

HcalDetDiagLEDData* HcalDetDiagLEDMonitor::he_data[85][72][4]
private

Definition at line 120 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), CheckStatus(), fillHistos(), LoadReference(), and SaveReference().

HcalDetDiagLEDData HcalDetDiagLEDMonitor::he_data[85][72][4]
private

Definition at line 198 of file HcalDetDiagLEDMonitor.cc.

MonitorElement * HcalDetDiagLEDMonitor::HEMphi
private

Definition at line 110 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HEPphi
private

Definition at line 109 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

HcalDetDiagLEDData* HcalDetDiagLEDMonitor::hf_data[85][72][4]
private

Definition at line 122 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), CheckStatus(), fillHistos(), LoadReference(), and SaveReference().

HcalDetDiagLEDData HcalDetDiagLEDMonitor::hf_data[85][72][4]
private

Definition at line 200 of file HcalDetDiagLEDMonitor.cc.

MonitorElement * HcalDetDiagLEDMonitor::HFMphi
private

Definition at line 112 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HFPphi
private

Definition at line 111 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HO0phi
private

Definition at line 113 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HO1Mphi
private

Definition at line 116 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HO1Pphi
private

Definition at line 114 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HO2Mphi
private

Definition at line 117 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::HO2Pphi
private

Definition at line 115 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

HcalDetDiagLEDData* HcalDetDiagLEDMonitor::ho_data[85][72][4]
private

Definition at line 121 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), CheckStatus(), fillHistos(), LoadReference(), and SaveReference().

HcalDetDiagLEDData HcalDetDiagLEDMonitor::ho_data[85][72][4]
private

Definition at line 199 of file HcalDetDiagLEDMonitor.cc.

int HcalDetDiagLEDMonitor::ievt_
private

Definition at line 78 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), CheckStatus(), endRun(), GetStatistics(), HcalDetDiagLEDMonitor(), and SaveReference().

bool HcalDetDiagLEDMonitor::IsReference
private

Definition at line 82 of file HcalDetDiagLEDMonitor.h.

Referenced by HcalDetDiagLEDMonitor(), and LoadReference().

std::map<unsigned int, int> HcalDetDiagLEDMonitor::KnownBadCells_
private

Definition at line 216 of file HcalDetDiagLEDMonitor.cc.

Referenced by beginRun(), CheckStatus(), and fillHistos().

double HcalDetDiagLEDMonitor::LEDMeanTreshold
private

Definition at line 84 of file HcalDetDiagLEDMonitor.h.

Referenced by HcalDetDiagLEDMonitor().

double HcalDetDiagLEDMonitor::LEDRmsTreshold
private

Definition at line 85 of file HcalDetDiagLEDMonitor.h.

Referenced by HcalDetDiagLEDMonitor().

HcalLogicalMap * HcalDetDiagLEDMonitor::lmap
private

Definition at line 45 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::meEVT_
private

Definition at line 92 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::meRUN_
private

Definition at line 170 of file HcalDetDiagLEDMonitor.cc.

Referenced by analyze(), and setup().

std::string HcalDetDiagLEDMonitor::OutputFilePath
private

Definition at line 90 of file HcalDetDiagLEDMonitor.h.

Referenced by HcalDetDiagLEDMonitor(), and SaveReference().

std::string HcalDetDiagLEDMonitor::ReferenceData
private

Definition at line 88 of file HcalDetDiagLEDMonitor.h.

Referenced by HcalDetDiagLEDMonitor(), and LoadReference().

std::string HcalDetDiagLEDMonitor::ReferenceRun
private

Definition at line 89 of file HcalDetDiagLEDMonitor.h.

Referenced by LoadReference(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::RefRun_
private

Definition at line 93 of file HcalDetDiagLEDMonitor.h.

Referenced by setup().

int HcalDetDiagLEDMonitor::run_number
private

Definition at line 80 of file HcalDetDiagLEDMonitor.h.

Referenced by analyze(), HcalDetDiagLEDMonitor(), and SaveReference().

MonitorElement * HcalDetDiagLEDMonitor::Time
private

Definition at line 95 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fillHistos(), SaveReference(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::Time2Dhbhehf
private

Definition at line 98 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::Time2Dho
private

Definition at line 99 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::TimeHF
private

Definition at line 97 of file HcalDetDiagLEDMonitor.h.

Referenced by CheckStatus(), fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::TimeRMS
private

Definition at line 103 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

MonitorElement * HcalDetDiagLEDMonitor::TimeRMSHF
private

Definition at line 105 of file HcalDetDiagLEDMonitor.h.

Referenced by fillHistos(), and setup().

edm::InputTag HcalDetDiagLEDMonitor::triggerLabel_
private

Definition at line 136 of file HcalDetDiagLEDMonitor.h.

bool HcalDetDiagLEDMonitor::UseDB
private

Definition at line 86 of file HcalDetDiagLEDMonitor.h.

std::string HcalDetDiagLEDMonitor::XmlFilePath
private

Definition at line 168 of file HcalDetDiagLEDMonitor.cc.

Referenced by HcalDetDiagLEDMonitor(), and SaveReference().

Generated for CMSSW Reference Manual by doxygen 1.8.5