HcalDetDiagLaserMonitor Class Reference

#include <DQM/HcalDetDiagLaserMonitor/src/HcalDetDiagLaserMonitor.cc>

Inheritance diagram for HcalDetDiagLaserMonitor:

Public Member Functions
	HcalDetDiagLaserMonitor (const edm::ParameterSet &)
	~HcalDetDiagLaserMonitor ()
Public Member Functions inherited from HcalBaseDQMonitor
	HcalBaseDQMonitor (const edm::ParameterSet &ps)
	HcalBaseDQMonitor ()
virtual	~HcalBaseDQMonitor ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) 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	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override
void	bookHistograms (DQMStore::IBooker &ib, const edm::Run &run, const edm::EventSetup &c) override
void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override
void	endRun (const edm::Run &run, const edm::EventSetup &c) override
void	fillHistos (int sd)
void	fillProblems (int sd)
bool	get_ave_rbx (int sd, int side, int rbx, float *ave, float *rms)
bool	get_ave_subdet (int sd, float *ave_t, float *ave_e, float *ave_t_r, float *ave_e_r)
HcalDetDiagLaserData *	GetCalib (std::string sd, int eta, int phi)
void	LoadDataset ()
void	LoadReference ()
void	SaveRaddamData ()
void	SaveReference ()

Private Attributes
std::string	baseFolder_
HcalDetDiagLaserData	calib_data [5][5][72]
bool	createHTMLonly
int	dataset_seq_number
std::string	DatasetName
const HcalElectronicsMap *	emap
MonitorElement *	Energy2Dhbhehf
MonitorElement *	Energy2Dho
HcalDetDiagLaserData	hb_data [85][72][4]
MonitorElement *	hb_time_rbx
MonitorElement *	hbheEnergy
MonitorElement *	hbheEnergyRMS
bool	HBHEseq_
MonitorElement *	hbheTime
MonitorElement *	hbheTimeRMS
HcalDetDiagLaserData	he_data [85][72][4]
MonitorElement *	he_time_rbx
HcalDetDiagLaserData	hf_data [85][72][4]
MonitorElement *	hf_time_rbx
MonitorElement *	hfEnergy
MonitorElement *	hfEnergyRMS
bool	HFseq_
MonitorElement *	hfTime
MonitorElement *	hfTimeRMS
HcalDetDiagLaserData	ho_data [85][72][4]
MonitorElement *	ho_time_rbx
MonitorElement *	hoEnergy
MonitorElement *	hoEnergyRMS
bool	HOseq_
MonitorElement *	hoTime
MonitorElement *	hoTimeRMS
MonitorElement *	htmlFolder
std::string	htmlOutputPath
int	ievt_
int	ievt_hbhe_
int	ievt_hf_
int	ievt_ho_
edm::InputTag	inputLabelDigi_
bool	IsReference
std::map< unsigned int, int >	KnownBadCells_
double	LaserEnergyThreshold
double	LaserTimingThreshold
int	lastHBHEorbit_
int	lastHForbit_
int	lastHOorbit_
bool	LocalRun
MonitorElement *	meEVT_
MonitorElement *	meRUN_
int	nChecksHBHE_
int	nChecksHF_
int	nChecksHO_
int	nHB
int	nHBHEchecks
int	nHE
int	nHF
int	nHFchecks
int	nHO
int	nHOchecks
bool	Online_
std::string	OutputFilePath
bool	Overwrite
std::string	prefixME_
EtaPhiHists *	ProblemCellsByDepth_energy
EtaPhiHists *	ProblemCellsByDepth_energy_val
EtaPhiHists *	ProblemCellsByDepth_timing
EtaPhiHists *	ProblemCellsByDepth_timing_val
std::vector< std::string >	problemnames_
MonitorElement *	Raddam [56]
HcalRaddamData	Raddam_data [56]
bool	RaddamRun
double	RaddamThreshold1
double	RaddamThreshold2
MonitorElement *	refEnergy2Dhbhehf
MonitorElement *	refEnergy2Dho
std::string	ReferenceData
std::string	ReferenceRun
MonitorElement *	RefRun_
MonitorElement *	refTime2Dhbhehf
MonitorElement *	refTime2Dho
int	run_number
TH1F *	S1 [56]
TH1F *	S2 [56]
MonitorElement *	Time2Dhbhehf
MonitorElement *	Time2Dho
edm::EDGetTokenT < HcalCalibDigiCollection >	tok_calib_
edm::EDGetTokenT < HBHEDigiCollection >	tok_hbhe_
edm::EDGetTokenT < HFDigiCollection >	tok_hf_
edm::EDGetTokenT < HODigiCollection >	tok_ho_
edm::EDGetTokenT < FEDRawDataCollection >	tok_raw_
edm::EDGetTokenT < HcalTBTriggerData >	tok_tb_
std::string	XmlFilePath

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T...>	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from HcalBaseDQMonitor
virtual void	cleanup (void)
virtual void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &c)
void	getLogicalMap (const edm::EventSetup &c)
bool	IsAllowedCalibType ()
bool	LumiInOrder (int lumisec)
virtual void	reset (void)
virtual void	setup (DQMStore::IBooker &)
void	SetupEtaPhiHists (DQMStore::IBooker &ib, 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_
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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 249 of file HcalDetDiagLaserMonitor.cc.

Constructor & Destructor Documentation

HcalDetDiagLaserMonitor::HcalDetDiagLaserMonitor ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 386 of file HcalDetDiagLaserMonitor.cc.

References createHTMLonly, dataset_seq_number, DatasetName, HcalBaseDQMonitor::debug_, emap, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), htmlOutputPath, ievt_, inputLabelDigi_, HLT_25ns10e33_v2_cff::InputTag, IsReference, LaserEnergyThreshold, LaserTimingThreshold, LocalRun, nHB, nHBHEchecks, nHE, nHF, nHFchecks, nHO, nHOchecks, Online_, OutputFilePath, Overwrite, prefixME_, ProblemCellsByDepth_energy, ProblemCellsByDepth_energy_val, ProblemCellsByDepth_timing, ProblemCellsByDepth_timing_val, RaddamRun, RaddamThreshold1, RaddamThreshold2, ReferenceData, run_number, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, tok_calib_, tok_hbhe_, tok_hf_, tok_ho_, tok_raw_, tok_tb_, and XmlFilePath.

                                                                              :
 HcalBaseDQMonitor(iConfig)
 {
 

  ievt_=-1;
  emap=0;
  dataset_seq_number=1;
  run_number=-1;
  IsReference=false;
  LocalRun=RaddamRun=false;
  createHTMLonly=false;
  nHB=nHE=nHO=nHF=0;
  nHBHEchecks=nHOchecks=nHFchecks=0;

  inputLabelDigi_  = iConfig.getUntrackedParameter<edm::InputTag>("digiLabel",edm::InputTag("hcalDigis"));

  ReferenceData    = iConfig.getUntrackedParameter<std::string>("LaserReferenceData" ,"");
  OutputFilePath   = iConfig.getUntrackedParameter<std::string>("OutputFilePath", "");
  DatasetName      = iConfig.getUntrackedParameter<std::string>("LaserDatasetName", "");
  htmlOutputPath   = iConfig.getUntrackedParameter<std::string>("htmlOutputPath", "");
  XmlFilePath      = iConfig.getUntrackedParameter<std::string>("XmlFilePath", "");
  Online_          = iConfig.getUntrackedParameter<bool>  ("online",false);
  Overwrite        = iConfig.getUntrackedParameter<bool>  ("Overwrite",true);
  prefixME_        = iConfig.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
  if (prefixME_.size()>0 && prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
  prefixME_.append("/");
  subdir_          = iConfig.getUntrackedParameter<std::string>("TaskFolder","DetDiagPedestalMonitor_Hcal/");
  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
    subdir_.append("/");
  subdir_=prefixME_+subdir_;
  debug_           = iConfig.getUntrackedParameter<int>("debug",0);

  LaserTimingThreshold = iConfig.getUntrackedParameter<double>("LaserTimingThreshold",0.2);
  LaserEnergyThreshold = iConfig.getUntrackedParameter<double>("LaserEnergyThreshold",0.1);
  RaddamThreshold1     = iConfig.getUntrackedParameter<double>("RaddamThreshold1",10.0);
  RaddamThreshold2     = iConfig.getUntrackedParameter<double>("RaddamThreshold2",0.95);

  // register for data access
   tok_tb_ = consumes<HcalTBTriggerData>(iConfig.getParameter<edm::InputTag>("hcalTBTriggerDataTag"));
  tok_raw_  = consumes<FEDRawDataCollection>(iConfig.getUntrackedParameter<edm::InputTag>("RawDataLabel",edm::InputTag("source")));
  tok_calib_  = consumes<HcalCalibDigiCollection>(iConfig.getUntrackedParameter<edm::InputTag>("calibDigiLabel",edm::InputTag("hcalDigis")));
  tok_hbhe_ = consumes<HBHEDigiCollection>(inputLabelDigi_);
  tok_ho_  = consumes<HODigiCollection>(inputLabelDigi_);
  tok_hf_ = consumes<HFDigiCollection>(inputLabelDigi_);


  ProblemCellsByDepth_timing=0;
  ProblemCellsByDepth_energy=0;
  ProblemCellsByDepth_timing_val=0;
  ProblemCellsByDepth_energy_val=0;

}

HcalDetDiagLaserMonitor::~HcalDetDiagLaserMonitor

(

)

Definition at line 575 of file HcalDetDiagLaserMonitor.cc.

References ProblemCellsByDepth_energy, ProblemCellsByDepth_energy_val, ProblemCellsByDepth_timing, and ProblemCellsByDepth_timing_val.

                                                 {

  if ( ProblemCellsByDepth_timing ) delete ProblemCellsByDepth_timing;
  if ( ProblemCellsByDepth_energy ) delete ProblemCellsByDepth_energy;
  if ( ProblemCellsByDepth_timing_val ) delete ProblemCellsByDepth_timing_val;
  if ( ProblemCellsByDepth_energy_val ) delete ProblemCellsByDepth_energy_val;

}

Member Function Documentation

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

overrideprivatevirtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 585 of file HcalDetDiagLaserMonitor.cc.

References hcaldqm::constants::adc2fC, HcalDetDiagLaserData::add_statistics(), HcalBaseDQMonitor::analyze(), calib, calib_data, createHTMLonly, HcalRaddamData::CUT1EVNT, HcalRaddamData::CUT2EVNT, FEDRawData::data(), data, HLT_25ns10e33_v2_cff::depth, alignCSCRings::e, eta, MonitorElement::Fill(), HcalObjRepresent::Fill(), fillHistos(), fillProblems(), edm::Event::getByToken(), hb_data, HBHEseq_, hc_HBHEHPD, hc_HFPMT, hc_HOHPD, hc_RADDAM, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, he_data, hf_data, HFseq_, ho_data, HOseq_, i, edm::EventBase::id(), ievt_, ievt_hbhe_, ievt_hf_, ievt_ho_, edm::HandleBase::isValid(), j, relval_2017::k, lastHBHEorbit_, lastHForbit_, lastHOorbit_, LocalRun, FEDNumbering::MAXHCALFEDID, meEVT_, meRUN_, FEDNumbering::MINHCALFEDID, N, nChecksHBHE_, nChecksHF_, nChecksHO_, Online_, edm::EventBase::orbitNumber(), phi, Raddam, Raddam_data, RaddamRun, RaddamThreshold1, RaddamThreshold2, lumiPlot::rawdata, HcalBaseDQMonitor::reset(), edm::EventID::run(), run_number, HcalRaddamData::s1_adc, S2, HcalRaddamData::s2_adc, SaveReference(), FEDRawData::size(), tok_calib_, tok_hbhe_, tok_hf_, tok_ho_, tok_raw_, tok_tb_, HcalRaddamData::TOTEVNT, and relativeConstraints::value.

                                                                                        {
 if(createHTMLonly) return;
  HcalBaseDQMonitor::analyze(iEvent,iSetup); // base class increments ievt_, etc. counters
 
int  eta,phi,depth,nTS;
   if(ievt_==-1){ 
       ievt_=0;HBHEseq_=HOseq_=HFseq_=false; lastHBHEorbit_=lastHOorbit_=lastHForbit_=-1;nChecksHBHE_=nChecksHO_=nChecksHF_=0; 
       ievt_hbhe_=0,ievt_ho_=0,ievt_hf_=0;
   }

   bool LaserEvent=false;
   bool LaserRaddam=false;
   int orbit=iEvent.orbitNumber();
   meRUN_->Fill(iEvent.id().run());
   // for local runs 
   edm::Handle<HcalTBTriggerData> trigger_data;
   iEvent.getByToken(tok_tb_, trigger_data);
   if(trigger_data.isValid()){
       if(trigger_data->wasLaserTrigger()) LaserEvent=true;
       LocalRun=true;
   }
   if(!LocalRun && Online_){
      if(HBHEseq_ && (orbit-lastHBHEorbit_)>(11223*10) && ievt_hbhe_>40){
         HBHEseq_=false;
         fillHistos(HcalBarrel);
         fillProblems(HcalBarrel);
         fillProblems(HcalEndcap);
         nChecksHBHE_++;
         ievt_hbhe_=0;
         for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) hb_data[i][j][k].reset();
         for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) he_data[i][j][k].reset();
      }
      if(HOseq_ && (orbit-lastHOorbit_)>(11223*10) && ievt_ho_>40){
         HOseq_=false;
         fillHistos(HcalOuter);
         fillProblems(HcalOuter);
         nChecksHO_++; 
         ievt_ho_=0;
         for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) ho_data[i][j][k].reset();
      }
      if(HFseq_ && (orbit-lastHForbit_)>(11223*10) && ievt_hf_>40){
         HFseq_=false;
         fillHistos(HcalForward);
         fillProblems(HcalForward);
         nChecksHF_++; 
         ievt_hf_=0;
         if(nChecksHF_==1 || (nChecksHF_>1 && ((nChecksHF_-1)%12)==0)){
             SaveReference();
         }
         for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) hf_data[i][j][k].reset();
      }
   }

   // Abort Gap laser 
   if(LocalRun==false || LaserEvent==false){
     edm::Handle<FEDRawDataCollection> rawdata;
     iEvent.getByToken(tok_raw_ ,rawdata);
       // edm::Handle<FEDRawDataCollection> rawdata;
       // iEvent.getByType(rawdata);
       //checking FEDs for calibration information
       for (int i=FEDNumbering::MINHCALFEDID;i<=FEDNumbering::MAXHCALFEDID; i++) {
          const FEDRawData& fedData = rawdata->FEDData(i) ;
          if ( fedData.size() < 24 ) continue ;
          int value = ((const HcalDCCHeader*)(fedData.data()))->getCalibType() ;
      if(value==hc_HBHEHPD){ HBHEseq_=true; HOseq_=HFseq_=false; lastHBHEorbit_=orbit; ievt_hbhe_++; }
      if(value==hc_HOHPD){   HOseq_=true; HBHEseq_=HFseq_=false; lastHOorbit_=orbit;   ievt_ho_++;   }
      if(value==hc_HFPMT){   HFseq_=true; HBHEseq_=HOseq_=false; lastHForbit_=orbit;   ievt_hf_++;   }
          
          if(value==hc_HBHEHPD || value==hc_HOHPD || value==hc_HFPMT){ LaserEvent=true; break;}
      if(value==hc_RADDAM){ LaserEvent=true; LaserRaddam=true; break;} 
       }
   }   
   if(!LaserEvent) return;
   edm::Handle<HBHEDigiCollection> hbhe; 
   iEvent.getByToken(tok_hbhe_,hbhe);
   edm::Handle<HODigiCollection> ho; 
   iEvent.getByToken(tok_ho_,ho);
   edm::Handle<HFDigiCollection> hf;
   iEvent.getByToken(tok_hf_,hf);
   edm::Handle<HcalCalibDigiCollection> calib;
   iEvent.getByToken(tok_calib_, calib); 

   if(LocalRun && LaserEvent){
      int N=0; 
      if(hf.isValid()){
         for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
             eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth();
             float e=0;
         for(int i=0;i<digi->size();i++) e+=adc2fC[digi->sample(i).adc()&0xff]-2.5;
             if(e>40){ N++;}
         }
      }
      if(N>50 && N<57){ RaddamRun=true; /*LaserRaddam=true;*/}
   }
   if(RaddamRun){
      if(hf.isValid()){
         for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
             eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
             int N;
         for(N=0;N<56;N++)if(eta==RADDAM_CH[N].eta && phi==RADDAM_CH[N].phi && depth==RADDAM_CH[N].depth) break;
         if(N==56) continue; 
             float max1=0,max2=0;
             int   nmax1=0,nmax2=0;
         for(int i=0;i<nTS;i++){
                  if(max1<adc2fC[digi->sample(i).adc()&0xff]){ nmax1=i; max1=adc2fC[digi->sample(i).adc()&0xff]; }
             }
             Raddam_data[N].TOTEVNT++;
         for(int i=0;i<nTS;i++){
                  if(i==nmax1) continue;
                  if(max2<adc2fC[digi->sample(i).adc()&0xff]){ nmax2=i; max2=adc2fC[digi->sample(i).adc()&0xff]; }
             }
             if(nmax1>nmax2){
                int tmp1=nmax2;
                nmax2=nmax1;nmax1=tmp1;
             }
         if(nmax1==0 || nmax2==(nTS-1)) continue;
             if(nmax2!=(nmax1+1)) continue;
     
             if(max1<RaddamThreshold1 || max2<RaddamThreshold1) continue;
             Raddam_data[N].CUT1EVNT++;
             max1-=2.5; max2-=2.5;
             float S2=max1+max2;
             float S4=S2+adc2fC[digi->sample(nmax1-1).adc()&0xff]+adc2fC[digi->sample(nmax2+1).adc()&0xff]-5.0;
             if((S2/S4)<RaddamThreshold2) continue;
             Raddam_data[N].CUT2EVNT++;
             Raddam_data[N].s1_adc[digi->sample(nmax1).adc()&0xff]++;
             Raddam_data[N].s2_adc[digi->sample(nmax2).adc()&0xff]++;
         }
      }
   }

   meEVT_->Fill(++ievt_);
   run_number=iEvent.id().run();
   double data[20];
   if(!LaserRaddam){
      if(hbhe.isValid()){
         for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
             eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
             double ped=(adc2fC[digi->sample(0).adc()&0xff]+adc2fC[digi->sample(1).adc()&0xff])/2.0;
         if(digi->id().subdet()==HcalBarrel){
        for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
        hb_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
         }   
             if(digi->id().subdet()==HcalEndcap){
        for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
        he_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
         }
         }
      }
      if(ho.isValid()){
         for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
             eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
             double ped=(adc2fC[digi->sample(0).adc()&0xff]+adc2fC[digi->sample(1).adc()&0xff])/2.0;
         if((eta>=11 && eta<=15 && phi>=59 && phi<=70) || (eta>=5 && eta<=10 && phi>=47 && phi<=58)){
            for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
         }else{
            for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
         }
             ho_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
         }
      }
      if(hf.isValid()){
         for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
             eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
             double ped=adc2fC[digi->sample(0).adc()&0xff];
         for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
         hf_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
         }   
      }
      if(calib.isValid())for(HcalCalibDigiCollection::const_iterator digi=calib->begin();digi!=calib->end();digi++){
         if(digi->id().cboxChannel()!=0 || digi->id().hcalSubdet()==0) continue; 
         nTS=digi->size();
         float e=0;
         for(int i=0;i<nTS;i++){ data[i]=adc2fC[digi->sample(i).adc()&0xff]; e+=data[i];} 
         if(e<15000) calib_data[digi->id().hcalSubdet()][digi->id().ieta()+2][digi->id().iphi()-1].add_statistics(data,nTS);
      }
   }else{ //Raddam
      if(hf.isValid()){
         for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
             eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
         int N;
         for(N=0;N<56;N++)if(eta==RADDAM_CH[N].eta && phi==RADDAM_CH[N].phi && depth==RADDAM_CH[N].depth) break;
         if(N==56) continue;      
         for(int i=0;i<nTS;i++) Raddam[N]->Fill(i,adc2fC[digi->sample(i).adc()&0xff]-2.5);
         
         }   
      }
   }
}

void HcalDetDiagLaserMonitor::beginLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg, const edm::EventSetup &  c  )

overrideprivatevirtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 2061 of file HcalDetDiagLaserMonitor.cc.

{}

void HcalDetDiagLaserMonitor::bookHistograms ( DQMStore::IBooker &  ib, const edm::Run &  run, const edm::EventSetup &  c  )

overrideprivatevirtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 441 of file HcalDetDiagLaserMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookInt(), DQMStore::IBooker::bookString(), createHTMLonly, dataset_seq_number, DatasetName, EtaPhiHists::depth, emap, Energy2Dhbhehf, Energy2Dho, eta, MonitorElement::Fill(), edm::EventSetup::get(), HcalCondObjectContainer< Item >::getAllChannels(), reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), hb_time_rbx, hbheEnergy, hbheEnergyRMS, hbheTime, hbheTimeRMS, DetId::Hcal, HcalChannelStatus::HcalCellMask, HcalChannelStatus::HcalCellOff, he_time_rbx, hf_time_rbx, hfEnergy, hfEnergyRMS, hfTime, hfTimeRMS, ho_time_rbx, hoEnergy, hoEnergyRMS, hoTime, hoTimeRMS, htmlFolder, htmlOutputPath, i, KnownBadCells_, LoadDataset(), LoadReference(), hlt_dqm_clientPB-live_cfg::me, meEVT_, meRUN_, mergeVDriftHistosByStation::name, nHB, nHE, nHF, nHO, AlCaHLTBitMon_ParallelJobs::p, phi, prefixME_, ProblemCellsByDepth_energy, ProblemCellsByDepth_energy_val, ProblemCellsByDepth_timing, ProblemCellsByDepth_timing_val, problemnames_, edm::ESHandle< class >::product(), Raddam, refEnergy2Dhbhehf, refEnergy2Dho, ReferenceRun, RefRun_, refTime2Dhbhehf, refTime2Dho, run_number, MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), EtaPhiHists::setup(), HcalBaseDQMonitor::setup(), mps_update::status, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, Time2Dhbhehf, and Time2Dho.

                                                                                                          {
  edm::ESHandle<HcalChannelQuality> p;
  c.get<HcalChannelQualityRcd>().get("withTopo",p);
  const HcalChannelQuality* chanquality= p.product();
  std::vector<DetId> mydetids = chanquality->getAllChannels();
  KnownBadCells_.clear();

  for (std::vector<DetId>::const_iterator i = mydetids.begin();i!=mydetids.end();++i){
     if (i->det()!=DetId::Hcal) continue; // not an hcal cell
     HcalDetId id=HcalDetId(*i);
     int status=(chanquality->getValues(id))->getValue();
     if((status & HcalChannelStatus::HcalCellOff) || (status & HcalChannelStatus::HcalCellMask)){
     KnownBadCells_[id.rawId()]=status;
     }
  } 

  edm::ESHandle<HcalDbService> conditions_;
  c.get<HcalDbRecord>().get(conditions_);
  emap=conditions_->getHcalMapping();
  
  HcalBaseDQMonitor::setup(ib);
  std::string name;
 
  ib.setCurrentFolder(subdir_);   
  meEVT_ = ib.bookInt("HcalDetDiagLaserMonitor Event Number");
  meRUN_ = ib.bookInt("HcalDetDiagLaserMonitor Run Number");

  ReferenceRun="UNKNOWN";
  LoadReference();
  LoadDataset();
  if(DatasetName.size()>0 && createHTMLonly){
     char str[200]; sprintf(str,"%sHcalDetDiagLaserData_run%i_%i/",htmlOutputPath.c_str(),run_number,dataset_seq_number);
     htmlFolder=ib.bookString("HcalDetDiagLaserMonitor HTML folder",str);
     MonitorElement *me;
     ib.setCurrentFolder(prefixME_+"HcalInfo");
     me=ib.bookInt("HBpresent");
     if(nHB>0) me->Fill(1);
     me=ib.bookInt("HEpresent");
     if(nHE>0) me->Fill(1);
     me=ib.bookInt("HOpresent");
     if(nHO>0) me->Fill(1);
     me=ib.bookInt("HFpresent");
     if(nHF>0) me->Fill(1);
  }
  ProblemCellsByDepth_timing = new EtaPhiHists();
  ProblemCellsByDepth_timing->setup(ib," Problem Bad Laser Timing");
  for(unsigned int i=0;i<ProblemCellsByDepth_timing->depth.size();i++)
          problemnames_.push_back(ProblemCellsByDepth_timing->depth[i]->getName());
  ProblemCellsByDepth_energy = new EtaPhiHists();
  ProblemCellsByDepth_energy->setup(ib," Problem Bad Laser Energy");
  for(unsigned int i=0;i<ProblemCellsByDepth_energy->depth.size();i++)
          problemnames_.push_back(ProblemCellsByDepth_energy->depth[i]->getName());

  ib.setCurrentFolder(subdir_+"Summary Plots");
     
  name="HBHE Laser Energy Distribution";                hbheEnergy        = ib.book1D(name,name,200,0,3000);
  name="HBHE Laser Timing Distribution";                hbheTime          = ib.book1D(name,name,200,0,10);
  name="HBHE Laser Energy RMS_div_Energy Distribution"; hbheEnergyRMS     = ib.book1D(name,name,200,0,0.5);
  name="HBHE Laser Timing RMS Distribution";            hbheTimeRMS       = ib.book1D(name,name,200,0,1);
  name="HO Laser Energy Distribution";                  hoEnergy          = ib.book1D(name,name,200,0,3000);
  name="HO Laser Timing Distribution";                  hoTime            = ib.book1D(name,name,200,0,10);
  name="HO Laser Energy RMS_div_Energy Distribution";   hoEnergyRMS       = ib.book1D(name,name,200,0,0.5);
  name="HO Laser Timing RMS Distribution";              hoTimeRMS         = ib.book1D(name,name,200,0,1);
  name="HF Laser Energy Distribution";                  hfEnergy          = ib.book1D(name,name,200,0,3000);
  name="HF Laser Timing Distribution";                  hfTime            = ib.book1D(name,name,200,0,10);
  name="HF Laser Energy RMS_div_Energy Distribution";   hfEnergyRMS       = ib.book1D(name,name,200,0,0.7);
  name="HF Laser Timing RMS Distribution";              hfTimeRMS         = ib.book1D(name,name,200,0,1);
     
  name="Laser Timing HBHEHF";                           Time2Dhbhehf      = ib.book2D(name,name,87,-43,43,74,0,73);
  name="Laser Timing HO";                               Time2Dho          = ib.book2D(name,name,33,-16,16,74,0,73);
  name="Laser Energy HBHEHF";                           Energy2Dhbhehf    = ib.book2D(name,name,87,-43,43,74,0,73);
  name="Laser Energy HO";                               Energy2Dho        = ib.book2D(name,name,33,-16,16,74,0,73);
  name="HBHEHF Laser (Timing-Ref)+1";                   refTime2Dhbhehf   = ib.book2D(name,name,87,-43,43,74,0,73);
  name="HO Laser (Timing-Ref)+1";                       refTime2Dho       = ib.book2D(name,name,33,-16,16,74,0,73);
  name="HBHEHF Laser Energy_div_Ref";                   refEnergy2Dhbhehf = ib.book2D(name,name,87,-43,43,74,0,73);
  name="HO Laser Energy_div_Ref";                       refEnergy2Dho     = ib.book2D(name,name,33,-16,16,74,0,73);
     
  name="HB RBX average Time-Ref";                       hb_time_rbx       = ib.book1D(name,name,36,0.5,36.5);
  name="HE RBX average Time-Ref";                       he_time_rbx       = ib.book1D(name,name,36,0.5,36.5);
  name="HO RBX average Time-Ref";                       ho_time_rbx       = ib.book1D(name,name,36,0.5,36.5);
  name="HF RoBox average Time-Ref";                     hf_time_rbx       = ib.book1D(name,name,24,0.5,24.5);
  
  char str[200];
  for(int i=1;i<=18;i++){ sprintf(str,"HBM%02i",i);     hb_time_rbx->setBinLabel(i,str);    }
  for(int i=1;i<=18;i++){ sprintf(str,"HBP%02i",i);     hb_time_rbx->setBinLabel(i+18,str); }
  for(int i=1;i<=18;i++){ sprintf(str,"HEM%02i",i);     he_time_rbx->setBinLabel(i,str);    }
  for(int i=1;i<=18;i++){ sprintf(str,"HEP%02i",i);     he_time_rbx->setBinLabel(i+18,str); }
  for(int i=1;i<=12;i++){ sprintf(str,"HFM%02i",i);     hf_time_rbx->setBinLabel(i,str);    }
  for(int i=1;i<=12;i++){ sprintf(str,"HFP%02i",i);     hf_time_rbx->setBinLabel(i+12,str); }
  for(int i=1;i<=6;i++){  sprintf(str,"HO2M%02i",i*2);  ho_time_rbx->setBinLabel(i,str);    }
  for(int i=1;i<=6;i++){  sprintf(str,"HO1M%02i",i*2);  ho_time_rbx->setBinLabel(i+6,str);  }
  for(int i=1;i<=12;i++){ sprintf(str,"HO0%02i",i);     ho_time_rbx->setBinLabel(i+12,str); }
  for(int i=1;i<=6;i++){  sprintf(str,"HO1P%02i",i*2);  ho_time_rbx->setBinLabel(i+24,str); }
  for(int i=1;i<=6;i++){  sprintf(str,"HO2P%02i",i*2);  ho_time_rbx->setBinLabel(i+30,str); }

  Time2Dhbhehf->setAxisTitle("i#eta",1);
  Time2Dhbhehf->setAxisTitle("i#phi",2);
  Time2Dho->setAxisTitle("i#eta",1);
  Time2Dho->setAxisTitle("i#phi",2);
  Energy2Dhbhehf->setAxisTitle("i#eta",1);
  Energy2Dhbhehf->setAxisTitle("i#phi",2);
  Energy2Dho->setAxisTitle("i#eta",1);
  Energy2Dho->setAxisTitle("i#phi",2);
  refTime2Dhbhehf->setAxisTitle("i#eta",1);
  refTime2Dhbhehf->setAxisTitle("i#phi",2);
  refTime2Dho->setAxisTitle("i#eta",1);
  refTime2Dho->setAxisTitle("i#phi",2);
  refEnergy2Dhbhehf->setAxisTitle("i#eta",1);
  refEnergy2Dhbhehf->setAxisTitle("i#phi",2);
  refEnergy2Dho->setAxisTitle("i#eta",1);
  refEnergy2Dho->setAxisTitle("i#phi",2);

  refTime2Dhbhehf->setAxisRange(0,2,3);
  refTime2Dho->setAxisRange(0,2,3);
  refEnergy2Dhbhehf->setAxisRange(0.5,1.5,3);
  refEnergy2Dho->setAxisRange(0.5,1.5,3);

  ib.setCurrentFolder(subdir_);
  RefRun_= ib.bookString("HcalDetDiagLaserMonitor Reference Run",ReferenceRun);

  ib.setCurrentFolder(subdir_+"Raddam Plots");
  for(int i=0;i<56;i++){
     sprintf(str,"RADDAM (%i %i)",RADDAM_CH[i].eta,RADDAM_CH[i].phi);
     Raddam[i] = ib.book1D(str,str,10,-0.5,9.5); 
  }

  ib.setCurrentFolder(subdir_+"Plots for client");
  ProblemCellsByDepth_timing_val = new EtaPhiHists();
  ProblemCellsByDepth_timing_val->setup(ib," Laser Timing difference");
  ProblemCellsByDepth_energy_val = new EtaPhiHists();
  ProblemCellsByDepth_energy_val->setup(ib," Laser Energy difference");
}

void HcalDetDiagLaserMonitor::endLuminosityBlock ( const edm::LuminosityBlock &  lumiSeg, const edm::EventSetup &  c  )

overrideprivatevirtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 2062 of file HcalDetDiagLaserMonitor.cc.

{}

void HcalDetDiagLaserMonitor::endRun ( const edm::Run &  run, const edm::EventSetup &  c  )

overrideprivatevirtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 2045 of file HcalDetDiagLaserMonitor.cc.

References createHTMLonly, fillHistos(), fillProblems(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, ievt_, LocalRun, Online_, RaddamRun, SaveRaddamData(), and SaveReference().

                                                                              {
    if(RaddamRun){
       SaveRaddamData();
    }
    if((LocalRun || !Online_ || createHTMLonly) && ievt_>10){
       fillHistos(HcalBarrel);
       fillHistos(HcalOuter);
       fillHistos(HcalForward);
       fillProblems(HcalBarrel);
       fillProblems(HcalEndcap);
       fillProblems(HcalOuter); 
       fillProblems(HcalForward); 
       if(!RaddamRun)SaveReference();
    }
}

void HcalDetDiagLaserMonitor::fillHistos ( int  sd )

private

Definition at line 1145 of file HcalDetDiagLaserMonitor.cc.

References funct::abs(), HLT_25ns10e33_v2_cff::depth, Energy2Dhbhehf, Energy2Dho, eta, MonitorElement::Fill(), get_ave_rbx(), HcalDetDiagLaserData::get_average_amp(), HcalDetDiagLaserData::get_average_time(), hb_data, hb_time_rbx, hbheEnergy, hbheEnergyRMS, hbheTime, hbheTimeRMS, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, he_data, he_time_rbx, hf_data, hf_time_rbx, hfEnergy, hfEnergyRMS, hfTime, hfTimeRMS, ho_data, ho_time_rbx, hoEnergy, hoEnergyRMS, hoTime, hoTimeRMS, i, isSiPM(), bookConverter::max, min(), phi, refEnergy2Dhbhehf, refEnergy2Dho, refTime2Dhbhehf, refTime2Dho, MonitorElement::Reset(), RMS, plotscripts::rms(), MonitorElement::setAxisRange(), MonitorElement::setBinContent(), MonitorElement::setBinError(), Time2Dhbhehf, and Time2Dho.

Referenced by analyze(), and endRun().

                                              {
   if(sd==HcalBarrel || sd==HcalEndcap){
      hbheEnergy->Reset();
      hbheTime->Reset();
      hbheEnergyRMS->Reset();
      hbheTimeRMS->Reset();
      hb_time_rbx->Reset();
      he_time_rbx->Reset();
   }
   if(sd==HcalOuter){
      hoEnergy->Reset();
      hoTime->Reset();
      hoEnergyRMS->Reset();
      hoTimeRMS->Reset();
      Time2Dho->Reset();
      Energy2Dho->Reset();
      refTime2Dho->Reset(); 
      refEnergy2Dho->Reset();
      ho_time_rbx->Reset();
   }
   if(sd==HcalForward){
      hfEnergy->Reset();
      hfTime->Reset();
      hfEnergyRMS->Reset();
      hfTimeRMS->Reset();
      hf_time_rbx->Reset();
   }
   if(sd==HcalBarrel || sd==HcalEndcap){
     // HB histograms
     for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){ 
        double T=0,nT=0,E=0,nE=0;
        for(int depth=1;depth<=2;depth++){
           if(hb_data[eta+42][phi-1][depth-1].get_statistics()>10){
          double ave=0;
          double rms=0;
          double time=0; 
          double time_rms=0;
          hb_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
          hb_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
          hbheEnergy->Fill(ave);
          if(ave>0)hbheEnergyRMS->Fill(rms/ave);
          hbheTime->Fill(time);
          hbheTimeRMS->Fill(time_rms);
          T+=time; nT++; E+=ave; nE++;
           }
        } 
        if(nT>0){Time2Dhbhehf->setBinContent(eta+44,phi+1,T/nT);Energy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE); }
     } 
     // HE histograms
     for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
        double T=0,nT=0,E=0,nE=0;
        for(int depth=1;depth<=3;depth++){
           if(he_data[eta+42][phi-1][depth-1].get_statistics()>10){
              double ave=0; double rms=0; double time=0; double time_rms=0;
          he_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
          he_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
          hbheEnergy->Fill(ave);
          if(ave>0)hbheEnergyRMS->Fill(rms/ave);
          hbheTime->Fill(time);
          hbheTimeRMS->Fill(time_rms);
          T+=time; nT++; E+=ave; nE++;
           }
        }
        if(nT>0 && abs(eta)>16 ){Time2Dhbhehf->setBinContent(eta+44,phi+1,T/nT); Energy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE); }     
        if(nT>0 && abs(eta)>20 ){Time2Dhbhehf->setBinContent(eta+44,phi+2,T/nT); Energy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);}  
     } 
   }
   if(sd==HcalForward){
     // HF histograms
     for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
        double T=0,nT=0,E=0,nE=0;
        for(int depth=1;depth<=2;depth++){
           if(hf_data[eta+42][phi-1][depth-1].get_statistics()>10){
              double ave=0; double rms=0; double time=0; double time_rms=0;
          hf_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
          hf_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
          hfEnergy->Fill(ave);
          if(ave>0)hfEnergyRMS->Fill(rms/ave);
          hfTime->Fill(time);
          T+=time; nT++; E+=ave; nE++;
          hfTimeRMS->Fill(time_rms);
           }
        }   
        if(nT>0 && abs(eta)>29 ){ Time2Dhbhehf->setBinContent(eta+44,phi+1,T/nT);   Time2Dhbhehf->setBinContent(eta+44,phi+2,T/nT);}     
        if(nT>0 && abs(eta)>29 ){ Energy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE); Energy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);}   
     }
   } 
   if(sd==HcalOuter){
     // HO histograms
     for(int eta=-10;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
        if(eta>10 && !isSiPM(eta,phi,4)) continue;
        double T=0,nT=0,E=0,nE=0;
        for(int depth=4;depth<=4;depth++){
           if(ho_data[eta+42][phi-1][depth-1].get_statistics()>10){
              double ave=0; double rms=0; double time=0; double time_rms=0;
          ho_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
          ho_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
          hoEnergy->Fill(ave);
          if(ave>0)hoEnergyRMS->Fill(rms/ave);
          hoTime->Fill(time);
          T+=time; nT++; E+=ave; nE++;
          hoTimeRMS->Fill(time_rms);
           }
        }
        if(nT>0){ Time2Dho->Fill(eta,phi,T/nT); Energy2Dho->Fill(eta,phi+1,E/nE) ;}
     } 
   }

   // compare with reference...
   if(sd==HcalBarrel || sd==HcalEndcap){
     for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){ 
        double T=0,nT=0,E=0,nE=0;
        for(int depth=1;depth<=2;depth++){
           if(hb_data[eta+42][phi-1][depth-1].get_statistics()>10){
         double val=0,rms=0,time=0,time_rms=0;
         double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
         if(!hb_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
             if(!hb_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
         if(!hb_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
         E+=VAL/val; nE++;
         T+=TIME-time; nT++;
       }  
        }
        if(nE>0) refEnergy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE);  
        if(nT>0){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01;  refTime2Dhbhehf->setBinContent(eta+44,phi+1,TTT); } 
     } 
     for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
        double T=0,nT=0,E=0,nE=0;
        for(int depth=1;depth<=3;depth++){
           if(he_data[eta+42][phi-1][depth-1].get_statistics()>10){
         double val=0,rms=0,time=0,time_rms=0;
         double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
         if(!he_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
             if(!he_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
         if(!he_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
         E+=VAL/val; nE++;
         T+=TIME-time; nT++;
       }  
        }
        if(nE>0 && abs(eta)>16) refEnergy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE);  
        if(nT>0 && abs(eta)>16){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01;  refTime2Dhbhehf->setBinContent(eta+44,phi+1,TTT); } 
        if(nE>0 && abs(eta)>20) refEnergy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);  
        if(nT>0 && abs(eta)>20){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01;  refTime2Dhbhehf->setBinContent(eta+44,phi+2,TTT); }  
     } 
   }
   if(sd==HcalForward){
     for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
        double T=0,nT=0,E=0,nE=0;
        for(int depth=1;depth<=2;depth++){
           if(hf_data[eta+42][phi-1][depth-1].get_statistics()>10){
         double val=0,rms=0,time=0,time_rms=0;
         double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
         if(!hf_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
             if(!hf_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
         if(!hf_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
         E+=VAL/val; nE++;
         T+=TIME-time; nT++;
       }  
        }
        if(nE>0 && abs(eta)>29) refEnergy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE);  
        if(nT>0 && abs(eta)>29){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+1,TTT); }
        if(nE>0 && abs(eta)>29) refEnergy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);  
        if(nT>0 && abs(eta)>29){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+2,TTT); } 
     }
   } 
   if(sd==HcalOuter){
     for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
        if(eta>10 && !isSiPM(eta,phi,4)) continue;
        double T=0,nT=0,E=0,nE=0;
        for(int depth=4;depth<=4;depth++){
           if(ho_data[eta+42][phi-1][depth-1].get_statistics()>10){
         double val=0,rms=0,time=0,time_rms=0;
         double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
         if(!ho_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
             if(!ho_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
         if(!ho_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
         E+=VAL/val; nE++;
         T+=TIME-time; nT++;
       }  
         }
        if(nE>0) refEnergy2Dho->Fill(eta,phi,E/nE);  
        if(nT>0){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dho->Fill(eta,phi,TTT);}  
     } 
   }
  float min,max;
  if(sd==HcalBarrel || sd==HcalEndcap){
    min=100;max=-100;
    for(int i=1;i<=18;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalBarrel,-1,i,&ave,&rms)){
        hb_time_rbx->setBinContent(i,ave);
        hb_time_rbx->setBinError(i,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=18;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalBarrel,1,i,&ave,&rms)){
        hb_time_rbx->setBinContent(i+18,ave);
        hb_time_rbx->setBinError(i+18,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    if(max>-100)hb_time_rbx->setAxisRange(min-1,max+1,2);
    min=100;max=-100;
    for(int i=1;i<=18;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalEndcap,-1,i,&ave,&rms)){
        he_time_rbx->setBinContent(i,ave);
        he_time_rbx->setBinError(i,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=18;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalEndcap,1,i,&ave,&rms)){
        he_time_rbx->setBinContent(i+18,ave);
        he_time_rbx->setBinError(i+18,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    if(max>-100)he_time_rbx->setAxisRange(min-1,max+1,2);
  }
  if(sd==HcalOuter){
    min=100;max=-100;
    for(int i=1;i<=6;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalOuter,-2,i,&ave,&rms)){
        ho_time_rbx->setBinContent(i,ave);
        ho_time_rbx->setBinError(i,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=6;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalOuter,-1,i,&ave,&rms)){
        ho_time_rbx->setBinContent(i+6,ave);
        ho_time_rbx->setBinError(i+6,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=12;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalOuter,0,i,&ave,&rms)){
        ho_time_rbx->setBinContent(i+12,ave);
        ho_time_rbx->setBinError(i+12,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=6;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalOuter,1,i,&ave,&rms)){
        ho_time_rbx->setBinContent(i+24,ave);
        ho_time_rbx->setBinError(i+24,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=6;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalOuter,2,i,&ave,&rms)){
        ho_time_rbx->setBinContent(i+30,ave);
        ho_time_rbx->setBinError(i+30,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    if(max>-100)ho_time_rbx->setAxisRange(min-1,max+1,2);
  }
  if(sd==HcalForward){
    min=100;max=-100;
    for(int i=1;i<=12;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalForward,-1,i,&ave,&rms)){
        hf_time_rbx->setBinContent(i,ave);
        hf_time_rbx->setBinError(i,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    for(int i=1;i<=12;i++){
      float ave=-10,rms=-10;
      if(get_ave_rbx(HcalForward,1,i,&ave,&rms)){
        hf_time_rbx->setBinContent(i+12,ave);
        hf_time_rbx->setBinError(i+12,rms);
        if(ave<min) min=ave;
        if(ave>max) max=ave;
      }
    }
    if(max>-100)hf_time_rbx->setAxisRange(min-1,max+1,2);
  }
} 

void HcalDetDiagLaserMonitor::fillProblems ( int  sd )

private

Definition at line 835 of file HcalDetDiagLaserMonitor.cc.

References HcalElectronicsMap::allElectronicsIdPrecision(), CalcEtaBin(), HcalDetId::depth(), EtaPhiHists::depth, HLT_25ns10e33_v2_cff::depth, DetId::det(), alignCSCRings::e, emap, eta, get_ave_subdet(), 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(), HcalDetId::iphi(), KnownBadCells_, LaserEnergyThreshold, LaserTimingThreshold, HcalElectronicsMap::lookup(), HcalDetDiagLaserData::nBadEnergy, HcalDetDiagLaserData::nBadTime, HcalDetDiagLaserData::nChecks, phi, ProblemCellsByDepth_energy, ProblemCellsByDepth_energy_val, ProblemCellsByDepth_timing, ProblemCellsByDepth_timing_val, DetId::rawId(), RMS, plotscripts::rms(), and DetId::subdetId().

Referenced by analyze(), and endRun().

                                                {
float ave_t,ave_e,ave_t_r,ave_e_r;
   if(!get_ave_subdet(sd,&ave_t,&ave_e,&ave_t_r,&ave_e_r)) return;

   for(int i=0;i<4;i++){
      ProblemCellsByDepth_energy->depth[i]->Reset();
      ProblemCellsByDepth_timing->depth[i]->Reset();
   }

   std::vector <HcalElectronicsId> AllElIds = emap->allElectronicsIdPrecision();
   for (std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
     DetId detid=emap->lookup(*eid);
     if (detid.det()!=DetId::Hcal) continue;
     HcalGenericDetId gid(emap->lookup(*eid));
     if(!(!(gid.null()) && 
            (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap  ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
     int eta=0,phi=0,depth=0;
     HcalDetId hid(detid);
     if(KnownBadCells_.find(hid.rawId())==KnownBadCells_.end()) continue;
     eta=hid.ieta();
     phi=hid.iphi();
     depth=hid.depth();

     int e=CalcEtaBin(sd,eta,depth)+1;
     if(detid.subdetId()==HcalBarrel && sd==HcalBarrel){
    double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
    if(!hb_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
        if(!hb_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
    if(!hb_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
        hb_data[eta+42][phi-1][depth-1].nChecks++;
        float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
        if(diff_t>LaserTimingThreshold){
             hb_data[eta+42][phi-1][depth-1].nBadTime++; 
             ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
        }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0); 
        if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
          float diff_e=((VAL/ave_e))/(val/ave_e_r);
          if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
               hb_data[eta+42][phi-1][depth-1].nBadEnergy++;
               ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
          }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
        }
     }
     if(detid.subdetId()==HcalEndcap && sd==HcalEndcap){
    double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
    if(!he_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
        if(!he_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
    if(!he_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
        he_data[eta+42][phi-1][depth-1].nChecks++;
        float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
        if(diff_t>LaserTimingThreshold){
          he_data[eta+42][phi-1][depth-1].nBadTime++;
          ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
        }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0); 
        if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
          float diff_e=((VAL/ave_e))/(val/ave_e_r);
          if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
            he_data[eta+42][phi-1][depth-1].nBadEnergy++;
            ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
          }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
        }
     }
     if(detid.subdetId()==HcalOuter && sd==HcalOuter){
    double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
    if(!ho_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
        if(!ho_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
    if(!ho_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
        ho_data[eta+42][phi-1][depth-1].nChecks++;
        float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
        if(diff_t>LaserTimingThreshold){
           ho_data[eta+42][phi-1][depth-1].nBadTime++;
           ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
        }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0); 
        if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
          float diff_e=((VAL/ave_e))/(val/ave_e_r);
          if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
            ho_data[eta+42][phi-1][depth-1].nBadEnergy++;
            ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
          }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
        }
     }
     if(detid.subdetId()==HcalForward && sd==HcalForward){
    double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
    if(!hf_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
        if(!hf_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
    if(!hf_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
        hf_data[eta+42][phi-1][depth-1].nChecks++;
        float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
        if(diff_t>LaserTimingThreshold){
           hf_data[eta+42][phi-1][depth-1].nBadTime++;
           ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
        }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0); 
        if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
          float diff_e=((VAL/ave_e))/(val/ave_e_r);
          if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
            hf_data[eta+42][phi-1][depth-1].nBadEnergy++;
            ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
          }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
        }
     }
   }
   for (std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
     DetId detid=emap->lookup(*eid);
     if (detid.det()!=DetId::Hcal) continue;
     HcalGenericDetId gid(emap->lookup(*eid));
     if(!(!(gid.null()) && 
            (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap  ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
     int eta=0,phi=0,depth=0;
     HcalDetId hid(detid);
     eta=hid.ieta();
     phi=hid.iphi();
     depth=hid.depth();
   
     if(detid.subdetId()==HcalBarrel){
        if(hb_data[eta+42][phi-1][depth-1].nBadTime>0){
           int e=CalcEtaBin(HcalBarrel,eta,depth)+1; 
           double val=hb_data[eta+42][phi-1][depth-1].nBadTime/hb_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
        } 
        if(hb_data[eta+42][phi-1][depth-1].nBadEnergy>0){
           int e=CalcEtaBin(HcalBarrel,eta,depth)+1; 
           double val=hb_data[eta+42][phi-1][depth-1].nBadEnergy/hb_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
        } 
     }
     if(detid.subdetId()==HcalEndcap){
        if(he_data[eta+42][phi-1][depth-1].nBadTime>0){
           int e=CalcEtaBin(HcalEndcap,eta,depth)+1; 
           double val=he_data[eta+42][phi-1][depth-1].nBadTime/he_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
        } 
        if(he_data[eta+42][phi-1][depth-1].nBadEnergy>0){
           int e=CalcEtaBin(HcalEndcap,eta,depth)+1; 
           double val=he_data[eta+42][phi-1][depth-1].nBadEnergy/he_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
        } 
     }
     if(detid.subdetId()==HcalOuter){
        if(ho_data[eta+42][phi-1][depth-1].nBadTime>0){
           int e=CalcEtaBin(HcalOuter,eta,depth)+1; 
           double val=ho_data[eta+42][phi-1][depth-1].nBadTime/ho_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
        } 
        if(ho_data[eta+42][phi-1][depth-1].nBadEnergy>0){
           int e=CalcEtaBin(HcalOuter,eta,depth)+1; 
           double val=ho_data[eta+42][phi-1][depth-1].nBadEnergy/ho_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val); 
        } 
     }
     if(detid.subdetId()==HcalForward){
        if(hf_data[eta+42][phi-1][depth-1].nBadTime>0){
           int e=CalcEtaBin(HcalForward,eta,depth)+1; 
           double val=hf_data[eta+42][phi-1][depth-1].nBadTime/hf_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
        } 
        if(hf_data[eta+42][phi-1][depth-1].nBadEnergy>0){
           int e=CalcEtaBin(HcalForward,eta,depth)+1; 
           double val=hf_data[eta+42][phi-1][depth-1].nBadEnergy/hf_data[eta+42][phi-1][depth-1].nChecks;
           ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
        } 
     }
   }
}

bool HcalDetDiagLaserMonitor::get_ave_rbx ( int  sd, int  side, int  rbx, float *  ave, float *  rms  )

private

Definition at line 1005 of file HcalDetDiagLaserMonitor.cc.

References hb_data, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, he_data, hf_data, ho_data, i, j, relval_2017::k, gen::n, plotscripts::rms(), and mathSSE::sqrt().

Referenced by fillHistos().

                                                                                      {
   double xt=0,xxt=0; 
   int eta_min=0,eta_max=0,n=0;
   if(sd==HcalBarrel){
      if(side>0){eta_min=1; eta_max=29;}
      if(side<0){eta_min=-29; eta_max=-1;}
      if(rbx==1){
         for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++)for(int k=1;k<=3;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!hb_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!hb_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
     }
         for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=2;j++)for(int k=1;k<=3;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!hb_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!hb_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
     }
      
      }else{
         for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*4-1);j<=((rbx-1)*4+2);j++)for(int k=1;k<=3;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!hb_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!hb_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
     }
      }
   }
   if(sd==HcalEndcap){
      if(side>0){eta_min=1; eta_max=29;}
      if(side<0){eta_min=-29; eta_max=-1;}
      if(rbx==1){
         for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++)for(int k=1;k<=3;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!he_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!he_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
     }
         for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=2;j++)for(int k=1;k<=3;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!he_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!he_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
     }
      
      }else{
         for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*4-1);j<=((rbx-1)*4+2);j++)for(int k=1;k<=3;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!he_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!he_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
     }
      }
   }  
   if(sd==HcalForward){
      if(side>0){eta_min=29; eta_max=40;}
      if(side<0){eta_min=-40; eta_max=-29;}
      for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*6+1);j<=((rbx-1)*6+6);j++)for(int k=1;k<=2;k++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!hf_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!hf_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
   }   
   if(sd==HcalOuter){
      if(side==0){
        eta_min=-4,eta_max=4;
        if(rbx==1){
          for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
          for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=4;j++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
    
    }else{
          for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*6-1);j<=((rbx-1)*6+4);j++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
    }
      }
      if(side==-1){ eta_min=-10,eta_max=-5;} 
      if(side==-2){ eta_min=-15,eta_max=-11;} 
      if(side==1) { eta_min=5,  eta_max=10;} 
      if(side==2) { eta_min=11, eta_max=15;}
      if(side!=0){
        if(rbx==1){
          for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
          for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=10;j++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
    
    }else{
          for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*12-1);j<=((rbx-1)*12+10);j++){
       double val,rms,time,time_rms;
       double TIME,TIME_RMS;
       if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
       if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
           xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
      }
    }      
      } 
   }
   if(n<10) return false;
   *ave=xt/n;
   *rms=sqrt(xxt/n-(xt*xt)/(n*n));
   return true; 
}

bool HcalDetDiagLaserMonitor::get_ave_subdet ( int  sd, float *  ave_t, float *  ave_e, float *  ave_t_r, float *  ave_e_r  )

private

Definition at line 774 of file HcalDetDiagLaserMonitor.cc.

References HLT_25ns10e33_v2_cff::depth, eta, HcalDetDiagLaserData::get_average_amp(), HcalDetDiagLaserData::get_average_time(), hb_data, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, he_data, hf_data, ho_data, phi, and plotscripts::rms().

Referenced by fillProblems().

                                                                                                          {
double T=0,nT=0,E=0,nE=0,Tr=0,nTr=0,Er=0,nEr=0;
   if(sd==HcalBarrel) for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){ 
      for(int depth=1;depth<=2;depth++){
         if(hb_data[eta+42][phi-1][depth-1].get_statistics()>10){
        double ave=0,rms=0,time=0,time_rms=0;
        hb_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
        hb_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
            T+=time; nT++; E+=ave; nE++;
            if(hb_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
              Tr+=time; nTr++; Er+=ave; nEr++;}
         }
      } 
   } 
   // HE histograms
   if(sd==HcalEndcap) for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
      for(int depth=1;depth<=3;depth++){
         if(he_data[eta+42][phi-1][depth-1].get_statistics()>10){
            double ave=0; double rms=0; double time=0; double time_rms=0;
        he_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
        he_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
            T+=time; nT++; E+=ave; nE++;
            if(he_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
              Tr+=time; nTr++; Er+=ave; nEr++;}
        }
      }
   } 
   // HF histograms
   if(sd==HcalForward) for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
      for(int depth=1;depth<=2;depth++){
         if(hf_data[eta+42][phi-1][depth-1].get_statistics()>10){
            double ave=0; double rms=0; double time=0; double time_rms=0;
        hf_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
        hf_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
        T+=time; nT++; E+=ave; nE++;
            if(hf_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
              Tr+=time; nTr++; Er+=ave; nEr++;}
         }
      } 
   } 
   // HO histograms
   if(sd==HcalOuter) for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
      for(int depth=4;depth<=4;depth++){
         if(ho_data[eta+42][phi-1][depth-1].get_statistics()>10){
            double ave=0; double rms=0; double time=0; double time_rms=0;
        ho_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
        ho_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
        T+=time; nT++; E+=ave; nE++;
            if(ho_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
              Tr+=time; nTr++; Er+=ave; nEr++;}
         }
      }
   } 
   if(nT<200 || nE<200 || nTr<200 || nEr<200) return false;
   *ave_t=T/nT;
   *ave_e=E/nE;
   *ave_t_r=Tr/nTr;
   *ave_e_r=Er/nEr;
   return true;
}

HcalDetDiagLaserData* HcalDetDiagLaserMonitor::GetCalib ( std::string  sd, int  eta, int  phi  )

inlineprivate

Definition at line 255 of file HcalDetDiagLaserMonitor.cc.

References funct::abs(), calib_data, ETA, PHI, and SD.

                                                                  {
        int SD=0,ETA=0,PHI=0;
        if(sd.compare("HB")==0) SD=1; 
        if(sd.compare("HE")==0) SD=2; 
        if(sd.compare("HO")==0) SD=3; 
        if(sd.compare("HF")==0) SD=4; 
        if(SD==1 || SD==2){
          if(eta>0) ETA=1; else ETA=-1;
          if(phi==71 ||phi==72 || phi==1 || phi==2) PHI=71; else PHI=((phi-3)/4)*4+3;
        }else if(SD==3){
          if(abs(eta)<=4){
      ETA=0;
      if(phi==71 ||phi==72 || phi==1 || phi==2 || phi==3 || phi==4) PHI=71; else PHI=((phi-5)/6)*6+5;
        }else{
      if(abs(eta)>4  && abs(eta)<=10)  ETA=1;
      if(abs(eta)>10 && abs(eta)<=15)  ETA=2;
      if(eta<0) ETA=-ETA;
      if(phi==71 ||phi==72 || (phi>=1 && phi<=10)) PHI=71; else PHI=((phi-11)/12)*12+11;
        }
        }else if(SD==4){
          if(eta>0) ETA=1; else ETA=-1;
          if(phi>=1  && phi<=18) PHI=1;
          if(phi>=19 && phi<=36) PHI=19;
          if(phi>=37 && phi<=54) PHI=37;
          if(phi>=55 && phi<=72) PHI=55;
      }
      return &calib_data[SD][ETA+2][PHI-1];
      };   

void HcalDetDiagLaserMonitor::LoadDataset ( )

private

Definition at line 1818 of file HcalDetDiagLaserMonitor.cc.

References calib_data, createHTMLonly, dataset_seq_number, DatasetName, g4TestOverlap_cfg::Depth, reco::tau::disc::Eta(), patZpeak::events, f, MonitorElement::Fill(), hb_data, he_data, hf_data, ho_data, ievt_, meEVT_, meRUN_, nHB, nHE, nHF, nHO, colinearityKinematic::Phi, plotscripts::rms(), DTTTrigCorrFirst::run, run_number, HcalDetDiagLaserData::set_data(), HcalDetDiagLaserData::set_statistics(), HcalDetDiagLaserData::set_statistics1(), and lumiQTWidget::t.

Referenced by bookHistograms().

                                         {
double amp,rms,time,time_rms;
int Eta,Phi,Depth,Statistic;
char subdet[10];
TFile *f;
      if(DatasetName.size()==0) return;
      createHTMLonly=true;
      if(gSystem->AccessPathName(DatasetName.c_str())) return;
      f = new TFile(DatasetName.c_str(),"READ");
 
      if(!f->IsOpen()) return ;

      TTree*  t=0;
      t=(TTree*)f->Get("HCAL Laser data");
      if(!t) return;
      t->SetBranchAddress("Subdet",   subdet);
      t->SetBranchAddress("eta",      &Eta);
      t->SetBranchAddress("phi",      &Phi);
      t->SetBranchAddress("depth",    &Depth);
      t->SetBranchAddress("amp",      &amp);
      t->SetBranchAddress("rms",      &rms);
      t->SetBranchAddress("time",     &time);
      t->SetBranchAddress("time_rms", &time_rms);
      t->SetBranchAddress("statistic",&Statistic);
      for(int ievt=0;ievt<t->GetEntries();ievt++){
         t->GetEntry(ievt);
     if(strcmp(subdet,"HB")==0){ nHB++;
            hb_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
            hb_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
            hb_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
     }
         if(strcmp(subdet,"HE")==0){ nHE++; 
            he_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
            he_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
            he_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
     }
         if(strcmp(subdet,"HO")==0){ nHO++;
            ho_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
            ho_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
            ho_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
     }
         if(strcmp(subdet,"HF")==0){ nHF++;
            hf_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
            hf_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
            hf_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
        }
        if(strcmp(subdet,"CALIB_HB")==0){
            calib_data[1][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
            calib_data[1][Eta+2][Phi-1].set_statistics(Statistic);
            calib_data[1][Eta+2][Phi-1].set_statistics1(Statistic);
        }
        if(strcmp(subdet,"CALIB_HE")==0){
            calib_data[2][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
            calib_data[2][Eta+2][Phi-1].set_statistics(Statistic);
            calib_data[2][Eta+2][Phi-1].set_statistics1(Statistic);
        }
        if(strcmp(subdet,"CALIB_HO")==0){
            calib_data[3][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
            calib_data[3][Eta+2][Phi-1].set_statistics(Statistic);
            calib_data[3][Eta+2][Phi-1].set_statistics1(Statistic);
        }
        if(strcmp(subdet,"CALIB_HF")==0){
            calib_data[4][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
            calib_data[4][Eta+2][Phi-1].set_statistics(Statistic);
            calib_data[4][Eta+2][Phi-1].set_statistics1(Statistic); 
        }
      }
      TObjString *STR1=(TObjString *)f->Get("run number");
      if(STR1){ int run; sscanf(STR1->String(),"%i",&run); meRUN_->Fill(run); run_number=run;}

      TObjString *STR2=(TObjString *)f->Get("Total events processed");
      if(STR2){ int events; sscanf(STR2->String(),"%i",&events); meEVT_->Fill(events); ievt_=events;}

      TObjString *STR3=(TObjString *)f->Get("Dataset number");
      if(STR3){ int ds; sscanf(STR3->String(),"%i",&ds); dataset_seq_number=ds;}
      f->Close(); 
      f->Delete();
} 

void HcalDetDiagLaserMonitor::LoadReference ( )

private

Definition at line 1780 of file HcalDetDiagLaserMonitor.cc.

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

Referenced by bookHistograms().

                                           {
double amp,rms,time,time_rms;
int Eta,Phi,Depth;
char subdet[10];
TFile *f;
 if(gSystem->AccessPathName(ReferenceData.c_str())) return;
 f = new TFile(ReferenceData.c_str(),"READ");
 
 if(!f->IsOpen()) return ;
 TObjString *STR=(TObjString *)f->Get("run number");
 
 if(STR){ std::string Ref(STR->String()); ReferenceRun=Ref;}
 
      TTree*  t=(TTree*)f->Get("HCAL Laser data");
      if(!t) return;
      t->SetBranchAddress("Subdet",   subdet);
      t->SetBranchAddress("eta",      &Eta);
      t->SetBranchAddress("phi",      &Phi);
      t->SetBranchAddress("depth",    &Depth);
      t->SetBranchAddress("amp",      &amp);
      t->SetBranchAddress("rms",      &rms);
      t->SetBranchAddress("time",     &time);
      t->SetBranchAddress("time_rms", &time_rms);
      for(int ievt=0;ievt<t->GetEntries();ievt++){
         t->GetEntry(ievt);
     if(strcmp(subdet,"HB")==0) hb_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
     if(strcmp(subdet,"HE")==0) he_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
     if(strcmp(subdet,"HO")==0) ho_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
     if(strcmp(subdet,"HF")==0) hf_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
         if(strcmp(subdet,"CALIB_HB")==0) calib_data[1][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
         if(strcmp(subdet,"CALIB_HE")==0) calib_data[2][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
         if(strcmp(subdet,"CALIB_HO")==0) calib_data[3][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
         if(strcmp(subdet,"CALIB_HF")==0) calib_data[4][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
     }
      f->Close();
      f->Delete();
      IsReference=true;
} 

void HcalDetDiagLaserMonitor::SaveRaddamData ( )

private

Definition at line 1897 of file HcalDetDiagLaserMonitor.cc.

References hcaldqm::constants::adc2fC, HcalRaddamData::CUT1EVNT, HcalRaddamData::CUT2EVNT, dataset_seq_number, Raddam_ch::depth, HLT_25ns10e33_v2_cff::depth, Raddam_ch::eta, eta, i, j, N, gen::n, OutputFilePath, Overwrite, phi, Raddam_ch::phi, Raddam_data, run_number, S1, HcalRaddamData::s1_adc, HcalRaddamData::S1BINWIDTH, HcalRaddamData::S1CHI2, HcalRaddamData::S1FITMEAN, HcalRaddamData::S1FITMEANER, HcalRaddamData::S1FITSIGMA, HcalRaddamData::S1MEAN, HcalRaddamData::S1NDF, HcalRaddamData::S1RMS, S2, HcalRaddamData::s2_adc, HcalRaddamData::S2BINWIDTH, HcalRaddamData::S2CHI2, HcalRaddamData::S2FITMEAN, HcalRaddamData::S2FITMEANER, HcalRaddamData::S2FITSIGMA, HcalRaddamData::S2MEAN, HcalRaddamData::S2NDF, HcalRaddamData::S2RMS, AlCaHLTBitMon_QueryRunRegistry::string, combinedConstraintHelpers::sum2(), lumiQTWidget::t, interactiveExample::theFile, HcalRaddamData::TOTEVNT, xmlFile(), and XmlFilePath.

Referenced by endRun().

                                            {
float adc_range[20]={14,28,40,52,67,132,202,262,322,397,722,1072,1372,1672,2047,3672,5422,6922,8422,10297};
int   adc_bins[20]={1,2,3,4,5,5,10,15,20,25,25,50,75,100,125,125,250,375,500,625};
char str[100];
      TF1 *fitFunc = new TF1("fitFunc","gaus");
      if(fitFunc==0) return;
      for(int i=0;i<56;i++){
          float sum1=0,sum2=0,n=0;
          S1[i]=S2[i]=0;
          for(int j=0;j<128;j++){
            sum1+=(adc2fC[j]-2.5)*Raddam_data[i].s1_adc[j];
            sum2+=(adc2fC[j]-2.5)*Raddam_data[i].s2_adc[j];
            n+=Raddam_data[i].s1_adc[j];
          }
          if(n<100) continue;
          sum1=sum1/n;
          sum2=sum2/n;
          int N=0;
          int Ws1=1,Ws2=1;
          for(N=1;N<19;N++) if(sum1>adc_range[N-1] && sum1<adc_range[N]) break;
          Ws1=adc_bins[N+1];
          for(N=1;N<19;N++) if(sum2>adc_range[N-1] && sum2<adc_range[N]) break;
          Ws2=adc_bins[N+1];
          sprintf(str,"Raddam(%i,%i,%i) S1",RADDAM_CH[i].eta,RADDAM_CH[i].phi,RADDAM_CH[i].depth);
          S1[i]=new TH1F(str,str,10000/Ws1,0,10000);
          sprintf(str,"Raddam(%i,%i,%i) S2",RADDAM_CH[i].eta,RADDAM_CH[i].phi,RADDAM_CH[i].depth);
          S2[i]=new TH1F(str,str,10000/Ws1,0,10000);
          for(int j=0;j<128;j++){
            S1[i]->Fill(adc2fC[j]-2.5,Raddam_data[i].s1_adc[j]);
            S2[i]->Fill(adc2fC[j]-2.5,Raddam_data[i].s2_adc[j]); 
          }
          double parm[3];
          S1[i]->Fit(fitFunc);
          S1[i]->GetFunction("fitFunc")->GetParameters(parm);
          Raddam_data[i].S1MEAN=S1[i]->GetMean();
          Raddam_data[i].S1RMS=S1[i]->GetRMS();
          Raddam_data[i].S1FITMEAN=parm[1];
          Raddam_data[i].S1FITMEANER=S1[i]->GetFunction("fitFunc")->GetParError(1);
          Raddam_data[i].S1FITSIGMA=parm[2];
          Raddam_data[i].S1CHI2=S1[i]->GetFunction("fitFunc")->GetChisquare();
          Raddam_data[i].S1NDF=S1[i]->GetFunction("fitFunc")->GetNDF();
          Raddam_data[i].S1BINWIDTH=Ws1;
          S2[i]->Fit(fitFunc);
          S2[i]->GetFunction("fitFunc")->GetParameters(parm);
          Raddam_data[i].S2MEAN=S2[i]->GetMean();
          Raddam_data[i].S2RMS=S2[i]->GetRMS();
          Raddam_data[i].S2FITMEAN=parm[1];
          Raddam_data[i].S2FITMEANER=S2[i]->GetFunction("fitFunc")->GetParError(1);
          Raddam_data[i].S2FITSIGMA=parm[2];
          Raddam_data[i].S2CHI2=S2[i]->GetFunction("fitFunc")->GetChisquare();
          Raddam_data[i].S2NDF=S2[i]->GetFunction("fitFunc")->GetNDF();
          Raddam_data[i].S2BINWIDTH=Ws2;
      }
      if(XmlFilePath.size()>0){
          char TIME[40];
          Long_t t; t=time(0); strftime(TIME,30,"%F %T",localtime(&t));      
          //create XML file
          if(!Overwrite){
             sprintf(str,"HcalDetDiagRaddam_%i_%i.xml",run_number,dataset_seq_number);
          }else{
             sprintf(str,"HcalDetDiagRaddam.xml");
          }
          std::string xmlName=str;
          std::ofstream xmlFile;
          xmlFile.open(xmlName.c_str());
          xmlFile<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?>\n";
          xmlFile<<"<ROOT>\n";
          xmlFile<<"  <HEADER>\n";
          xmlFile<<"    <TYPE>\n";
          xmlFile<<"      <EXTENSION_TABLE_NAME>HCAL_RADDAM</EXTENSION_TABLE_NAME>\n";
          xmlFile<<"      <NAME>HCAL Raddam</NAME>\n";
          xmlFile<<"    </TYPE>\n";
          xmlFile<<"    <!-- run details -->\n";
          xmlFile<<"    <RUN>\n";
          xmlFile<<"      <RUN_TYPE>TEST LOCAL-RUN</RUN_TYPE>\n";
          xmlFile<<"      <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
          xmlFile<<"      <RUN_BEGIN_TIMESTAMP>"<<TIME<<"</RUN_BEGIN_TIMESTAMP>\n";
          xmlFile<<"      <COMMENT_DESCRIPTION>hcal raddam data</COMMENT_DESCRIPTION>\n";
          xmlFile<<"      <LOCATION>P5</LOCATION>\n";
          xmlFile<<"      <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
          xmlFile<<"    </RUN>\n";
          xmlFile<<"  </HEADER>\n";
          xmlFile<<"  <DATA_SET>\n";
          xmlFile<<"     <COMMENT_DESCRIPTION>Test Raddam data</COMMENT_DESCRIPTION>\n";
          xmlFile<<"     <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
          xmlFile<<"     <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
          xmlFile<<"     <VERSION>Test_Version_1</VERSION>\n";
 
          for(int i=0;i<56;i++){
             xmlFile<<"     <DATA>\n";
             xmlFile<<"        <SUBDET>HF</SUBDET>\n";
             xmlFile<<"        <IETA>"<<RADDAM_CH[i].eta<<"</IETA>\n";
             xmlFile<<"        <IPHI>"<<RADDAM_CH[i].phi<<"</IPHI>\n";
             xmlFile<<"        <DEPTH>"<<RADDAM_CH[i].depth<<"</DEPTH>\n";

             xmlFile<<"        <TOTEVNT>"<<Raddam_data[i].TOTEVNT<<"</TOTEVNT>\n";
             xmlFile<<"        <CUT1EVNT>"<<Raddam_data[i].CUT1EVNT<<"</CUT1EVNT>\n";
             xmlFile<<"        <CUT2EVNT>"<<Raddam_data[i].CUT2EVNT<<"</CUT2EVNT>\n";

             xmlFile<<"        <S1MEAN>"<<Raddam_data[i].S1MEAN <<"</S1MEAN>\n";
             xmlFile<<"        <S1RMS>"<<Raddam_data[i].S1RMS <<"</S1RMS>\n";
             xmlFile<<"        <S1FITMEAN>"<<Raddam_data[i].S1FITMEAN <<"</S1FITMEAN>\n";
             xmlFile<<"        <S1FITMEANER>"<<Raddam_data[i].S1FITMEANER <<"</S1FITMEANER>\n";
             xmlFile<<"        <S1FITSIGMA>"<<Raddam_data[i].S1FITSIGMA <<"</S1FITSIGMA>\n";
             xmlFile<<"        <S1CHI2>"<<Raddam_data[i].S1CHI2 <<"</S1CHI2>\n";
             xmlFile<<"        <S1NDF>"<<Raddam_data[i].S1NDF <<"</S1NDF>\n";
             xmlFile<<"        <S1BINWIDTH>"<<Raddam_data[i].S1BINWIDTH <<"</S1BINWIDTH>\n";

             xmlFile<<"        <S2MEAN>"<<Raddam_data[i].S2MEAN <<"</S2MEAN>\n";
             xmlFile<<"        <S2RMS>"<<Raddam_data[i].S2RMS <<"</S2RMS>\n";
             xmlFile<<"        <S2FITMEAN>"<<Raddam_data[i].S2FITMEAN <<"</S2FITMEAN>\n";
             xmlFile<<"        <S2FITMEANER>"<<Raddam_data[i].S2FITMEANER <<"</S2FITMEANER>\n";
             xmlFile<<"        <S2FITSIGMA>"<<Raddam_data[i].S2FITSIGMA <<"</S2FITSIGMA>\n";
             xmlFile<<"        <S2CHI2>"<<Raddam_data[i].S2CHI2 <<"</S2CHI2>\n";
             xmlFile<<"        <S2NDF>"<<Raddam_data[i].S2NDF <<"</S2NDF>\n";
             xmlFile<<"        <S2BINWIDTH>"<<Raddam_data[i].S2BINWIDTH <<"</S2BINWIDTH>\n";
             xmlFile<<"    </DATA>\n";
          }
          xmlFile<<"  </DATA_SET>\n";
          xmlFile<<"</ROOT>\n";

          xmlFile.close();
          sprintf(str,"zip %s.zip %s",xmlName.c_str(),xmlName.c_str());
          system(str);
          sprintf(str,"rm -f %s",xmlName.c_str());
          system(str);
          sprintf(str,"mv -f %s.zip %s",xmlName.c_str(),XmlFilePath.c_str());
          system(str);
      }
      if(OutputFilePath.size()>0){
         if(!Overwrite){
            sprintf(str,"%sHcalDetDiagRaddamData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
         }else{
            sprintf(str,"%sHcalDetDiagRaddamData.root",OutputFilePath.c_str());
         }
         TFile *theFile = new TFile(str, "RECREATE");
         if(!theFile->IsOpen()) return;
         theFile->cd();
         for(int i=0;i<56;i++){
            if(S1[i]!=0)S1[i]->Write();
            if(S2[i]!=0)S2[i]->Write();
         }
         theFile->Write();
         theFile->Close();
     theFile->Delete();
      } 
}

void HcalDetDiagLaserMonitor::SaveReference ( )

private

Definition at line 1448 of file HcalDetDiagLaserMonitor.cc.

References HcalElectronicsMap::allElectronicsIdPrecision(), calib_data, dataset_seq_number, g4TestOverlap_cfg::Depth, HcalDetId::depth(), HLT_25ns10e33_v2_cff::depth, DetId::det(), alignCSCRings::e, emap, eta, reco::tau::disc::Eta(), patZpeak::events, HcalDetDiagLaserData::get_average_amp1(), HcalDetDiagLaserData::get_average_time1(), HcalDetDiagLaserData::get_statistics1(), HcalDetDiagLaserData::get_status(), 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(), j, relval_2017::k, HcalElectronicsMap::lookup(), OutputFilePath, Overwrite, phi, colinearityKinematic::Phi, plotscripts::rms(), DTTTrigCorrFirst::run, run_number, sd, AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), lumiQTWidget::t, interactiveExample::theFile, MainPageGenerator::tree, xmlFile(), and XmlFilePath.

Referenced by analyze(), and endRun().

                                           {
double amp,rms,Time,time_rms;
int    Eta,Phi,Depth,Statistic,Status=0;
char   Subdet[10],str[500];
   if(OutputFilePath.size()>0){
       if(!Overwrite){
          sprintf(str,"%sHcalDetDiagLaserData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
       }else{
          sprintf(str,"%sHcalDetDiagLaserData.root",OutputFilePath.c_str());
       }
       TFile *theFile = new TFile(str, "RECREATE");
       if(!theFile->IsOpen()) return;
       theFile->cd();
       sprintf(str,"%d",run_number); TObjString run(str);    run.Write("run number");
       sprintf(str,"%d",ievt_);      TObjString events(str); events.Write("Total events processed");
       sprintf(str,"%d",dataset_seq_number);      TObjString dsnum(str);  dsnum.Write("Dataset number");
       Long_t t; t=time(0); strftime(str,30,"%F %T",localtime(&t)); TObjString tm(str);  tm.Write("Dataset creation time");

       TTree *tree   =new TTree("HCAL Laser data","HCAL Laser data");
       if(tree==0)   return;
       tree->Branch("Subdet",   &Subdet,         "Subdet/C");
       tree->Branch("eta",      &Eta,            "Eta/I");
       tree->Branch("phi",      &Phi,            "Phi/I");
       tree->Branch("depth",    &Depth,          "Depth/I");
       tree->Branch("statistic",&Statistic,      "Statistic/I");
       tree->Branch("status",   &Status,         "Status/I");
       tree->Branch("amp",      &amp,            "amp/D");
       tree->Branch("rms",      &rms,            "rms/D");
       tree->Branch("time",     &Time,           "time/D");
       tree->Branch("time_rms", &time_rms,       "time_rms/D");
       sprintf(Subdet,"HB");
       for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
          if((Statistic=hb_data[eta+42][phi-1][depth-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=depth;
         Status=hb_data[eta+42][phi-1][depth-1].get_status();
         hb_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
         hb_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
          }
       } 
       sprintf(Subdet,"HE");
       for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=3;depth++){
         if((Statistic=he_data[eta+42][phi-1][depth-1].get_statistics1())>10){
            Eta=eta; Phi=phi; Depth=depth;
        Status=he_data[eta+42][phi-1][depth-1].get_status();
        he_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
        he_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
        tree->Fill();
         }
       } 
       sprintf(Subdet,"HO");
       for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++) for(int depth=4;depth<=4;depth++){
         if((Statistic=ho_data[eta+42][phi-1][depth-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=depth;
         Status=ho_data[eta+42][phi-1][depth-1].get_status();
         ho_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
         ho_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
         }
       } 
       sprintf(Subdet,"HF");
       for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
         if((Statistic=hf_data[eta+42][phi-1][depth-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=depth;
         Status=hf_data[eta+42][phi-1][depth-1].get_status();
         hf_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
         hf_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
         }
       }
       sprintf(Subdet,"CALIB_HB");
       for(int eta=-1;eta<=1;eta++) for(int phi=1;phi<=72;phi++){
          if((Statistic=calib_data[1][eta+2][phi-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=0;
         Status=calib_data[1][eta+2][phi-1].get_status();
         calib_data[1][eta+2][phi-1].get_average_amp1(&amp,&rms);
         calib_data[1][eta+2][phi-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
          }
       } 
       sprintf(Subdet,"CALIB_HE");
       for(int eta=-1;eta<=1;eta++) for(int phi=1;phi<=72;phi++){
          if((Statistic=calib_data[2][eta+2][phi-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=0;
         Status=calib_data[2][eta+2][phi-1].get_status();
         calib_data[2][eta+2][phi-1].get_average_amp1(&amp,&rms);
         calib_data[2][eta+2][phi-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
          }
       } 
       sprintf(Subdet,"CALIB_HO");
       for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
          if((Statistic=calib_data[3][eta+2][phi-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=0;
         Status=calib_data[3][eta+2][phi-1].get_status();
         calib_data[3][eta+2][phi-1].get_average_amp1(&amp,&rms);
         calib_data[3][eta+2][phi-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
          }
       } 
       sprintf(Subdet,"CALIB_HF");
       for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
          if((Statistic=calib_data[4][eta+2][phi-1].get_statistics1())>10){
             Eta=eta; Phi=phi; Depth=0;
         Status=calib_data[4][eta+2][phi-1].get_status();
         calib_data[4][eta+2][phi-1].get_average_amp1(&amp,&rms);
         calib_data[4][eta+2][phi-1].get_average_time1(&Time,&time_rms);
         tree->Fill();
          }
       } 
       theFile->Write();
       theFile->Close();
       theFile->Delete();
   }
   if(XmlFilePath.size()>0){
      char TIME[40];
      Long_t t; t=time(0); strftime(TIME,30,"%F %T",localtime(&t));
      //create XML file
      if(!Overwrite){
         sprintf(str,"HcalDetDiagLaser_%i_%i.xml",run_number,dataset_seq_number);
      }else{
         sprintf(str,"HcalDetDiagLaser.xml");
      }
      std::string xmlName=str;
      std::ofstream xmlFile;
      xmlFile.open(xmlName.c_str());

      xmlFile<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"yes\"?>\n";
      xmlFile<<"<ROOT>\n";
      xmlFile<<"  <HEADER>\n";
      xmlFile<<"    <HINTS mode='only-det-root'/>\n";
      xmlFile<<"    <TYPE>\n";
      xmlFile<<"      <EXTENSION_TABLE_NAME>HCAL_DETMON_LED_LASER_V1</EXTENSION_TABLE_NAME>\n";
      xmlFile<<"      <NAME>HCAL Laser HBHE HPD [abort gap global]</NAME>\n";
      xmlFile<<"    </TYPE>\n";
      xmlFile<<"    <!-- run details -->\n";
      xmlFile<<"    <RUN>\n";
      xmlFile<<"      <RUN_TYPE>GLOBAL-RUN</RUN_TYPE>\n";
      xmlFile<<"      <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
      xmlFile<<"      <RUN_BEGIN_TIMESTAMP>2009-01-01 00:00:00</RUN_BEGIN_TIMESTAMP>\n";
      xmlFile<<"      <COMMENT_DESCRIPTION>hcal laser data</COMMENT_DESCRIPTION>\n";
      xmlFile<<"      <LOCATION>P5</LOCATION>\n";
      xmlFile<<"      <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
      xmlFile<<"    </RUN>\n";
      xmlFile<<"  </HEADER>\n";
      xmlFile<<"  <DATA_SET>\n";
      xmlFile<<"     <!-- optional dataset metadata -->\n\n";
      xmlFile<<"     <SET_NUMBER>"<<dataset_seq_number<<"</SET_NUMBER>\n";
      xmlFile<<"     <SET_BEGIN_TIMESTAMP>2009-01-01 00:00:00</SET_BEGIN_TIMESTAMP>\n";
      xmlFile<<"     <SET_END_TIMESTAMP>2009-01-01 00:00:00</SET_END_TIMESTAMP>\n";
      xmlFile<<"     <NUMBER_OF_EVENTS_IN_SET>"<<ievt_<<"</NUMBER_OF_EVENTS_IN_SET>\n";
      xmlFile<<"     <COMMENT_DESCRIPTION>Automatic DQM output</COMMENT_DESCRIPTION>\n";
      xmlFile<<"     <DATA_FILE_NAME>"<< xmlName <<"</DATA_FILE_NAME>\n";
      xmlFile<<"     <IMAGE_FILE_NAME>data plot url or file path</IMAGE_FILE_NAME>\n";
      xmlFile<<"     <!-- who and when created this dataset-->\n\n";
      xmlFile<<"     <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
      xmlFile<<"     <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
      xmlFile<<"     <!-- version (string) and subversion (number) -->\n";
      xmlFile<<"     <!-- fields are used to read data back from the database -->\n\n";
      xmlFile<<"     <VERSION>"<<run_number<<dataset_seq_number<<"</VERSION>\n";
      xmlFile<<"     <SUBVERSION>1</SUBVERSION>\n";
      xmlFile<<"     <!--  Assign predefined dataset attributes -->\n\n";
      xmlFile<<"     <PREDEFINED_ATTRIBUTES>\n";
      xmlFile<<"        <ATTRIBUTE>\n";
      xmlFile<<"           <NAME>HCAL Dataset Status</NAME>\n";
      xmlFile<<"           <VALUE>VALID</VALUE>\n";
      xmlFile<<"        </ATTRIBUTE>\n";
      xmlFile<<"     </PREDEFINED_ATTRIBUTES>\n";
      xmlFile<<"     <!-- multiple data block records -->\n\n";

      std::vector <HcalElectronicsId> AllElIds = emap->allElectronicsIdPrecision();
      for(std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
         DetId detid=emap->lookup(*eid);
     if (detid.det()!=DetId::Hcal) continue;
         HcalGenericDetId gid(emap->lookup(*eid));
         if(!(!(gid.null()) && 
            (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap  ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
             gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
         int eta,phi,depth; 
         std::string subdet="";
     HcalDetId hid(detid);
     eta=hid.ieta();
     phi=hid.iphi();
     depth=hid.depth(); 
         
         double e=0,e_rms=0,t=0,t_rms=0;
         if(detid.subdetId()==HcalBarrel){
             subdet="HB";
             Statistic=hb_data[eta+42][phi-1][depth-1].get_statistics1();
         Status   =hb_data[eta+42][phi-1][depth-1].get_status();
         hb_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
         hb_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
         }else if(detid.subdetId()==HcalEndcap){
             subdet="HE";
             Statistic=he_data[eta+42][phi-1][depth-1].get_statistics1();
         Status   =he_data[eta+42][phi-1][depth-1].get_status();
         he_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
         he_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
     }else if(detid.subdetId()==HcalForward){
             subdet="HF";
             Statistic=hf_data[eta+42][phi-1][depth-1].get_statistics1();
         Status   =hf_data[eta+42][phi-1][depth-1].get_status();
         hf_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
         hf_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
     }else if(detid.subdetId()==HcalOuter){
             subdet="HO";
             Statistic=ho_data[eta+42][phi-1][depth-1].get_statistics1();
         Status   =ho_data[eta+42][phi-1][depth-1].get_status();
         ho_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
         ho_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
         }else continue;
         xmlFile<<"       <DATA>\n";
         xmlFile<<"          <NUMBER_OF_EVENTS_USED>"<<Statistic<<"</NUMBER_OF_EVENTS_USED>\n";
         xmlFile<<"          <SIGNAL_MEAN>"<<e<<"</SIGNAL_MEAN>\n";
         xmlFile<<"          <SIGNAL_RMS>"<<e_rms<<"</SIGNAL_RMS>\n";
         xmlFile<<"          <TIME_MEAN>"<<t<<"</TIME_MEAN>\n";
         xmlFile<<"          <TIME_RMS>"<<t_rms<<"</TIME_RMS>\n";
         xmlFile<<"          <CHANNEL_STATUS_WORD>"<<Status<<"</CHANNEL_STATUS_WORD>\n";
         xmlFile<<"          <CHANNEL_OBJECTNAME>HcalDetId</CHANNEL_OBJECTNAME>\n";
         xmlFile<<"             <SUBDET>"<<subdet<<"</SUBDET>\n";
         xmlFile<<"             <IETA>"<<eta<<"</IETA>\n";
         xmlFile<<"             <IPHI>"<<phi<<"</IPHI>\n";
         xmlFile<<"             <DEPTH>"<<depth<<"</DEPTH>\n";
         xmlFile<<"             <TYPE>0</TYPE>\n";
         xmlFile<<"       </DATA>\n";
      }
      xmlFile<<"  </DATA_SET>\n";
      xmlFile<<"</ROOT>\n";
      xmlFile.close();


      //create CALIB XML file 
      sprintf(str,"HcalDetDiagLaserCalib_%i_%i.xml",run_number,dataset_seq_number);
      std::string xmlNameCalib=str;
      std::ofstream xmlFileCalib;
      xmlFileCalib.open(xmlNameCalib.c_str());

      xmlFileCalib<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"yes\"?>\n";
      xmlFileCalib<<"<ROOT>\n";
      xmlFileCalib<<"  <HEADER>\n";
      xmlFileCalib<<"    <HINTS mode='only-det-root'/>\n";
      xmlFileCalib<<"    <TYPE>\n";
      xmlFileCalib<<"      <EXTENSION_TABLE_NAME>HCAL_DETMON_LED_LASER_V1</EXTENSION_TABLE_NAME>\n";
      xmlFileCalib<<"      <NAME>HCAL Laser CALIB [abort gap global]</NAME>\n";
      xmlFileCalib<<"    </TYPE>\n";
      xmlFileCalib<<"    <!-- run details -->\n";
      xmlFileCalib<<"    <RUN>\n";
      xmlFileCalib<<"      <RUN_TYPE>Global-RUN</RUN_TYPE>\n";
      xmlFileCalib<<"      <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
      xmlFileCalib<<"      <RUN_BEGIN_TIMESTAMP>2009-01-01 00:00:00</RUN_BEGIN_TIMESTAMP>\n";
      xmlFileCalib<<"      <COMMENT_DESCRIPTION>hcal Laser CALIB data</COMMENT_DESCRIPTION>\n";
      xmlFileCalib<<"      <LOCATION>P5</LOCATION>\n";
      xmlFileCalib<<"      <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
      xmlFileCalib<<"    </RUN>\n";
      xmlFileCalib<<"  </HEADER>\n";
      xmlFileCalib<<"  <DATA_SET>\n";
      xmlFileCalib<<"     <!-- optional dataset metadata -->\n\n";
      xmlFileCalib<<"     <SET_NUMBER>"<<dataset_seq_number<<"</SET_NUMBER>\n";
      xmlFileCalib<<"     <SET_BEGIN_TIMESTAMP>2009-01-01 00:00:00</SET_BEGIN_TIMESTAMP>\n";
      xmlFileCalib<<"     <SET_END_TIMESTAMP>2009-01-01 00:00:00</SET_END_TIMESTAMP>\n";
      xmlFileCalib<<"     <NUMBER_OF_EVENTS_IN_SET>"<<ievt_<<"</NUMBER_OF_EVENTS_IN_SET>\n";
      xmlFileCalib<<"     <COMMENT_DESCRIPTION>Automatic DQM output</COMMENT_DESCRIPTION>\n";
      xmlFileCalib<<"     <DATA_FILE_NAME>"<< xmlNameCalib <<"</DATA_FILE_NAME>\n";
      xmlFileCalib<<"     <IMAGE_FILE_NAME>data plot url or file path</IMAGE_FILE_NAME>\n";
      xmlFileCalib<<"     <!-- who and when created this dataset-->\n\n";
      xmlFileCalib<<"     <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
      xmlFileCalib<<"     <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
      xmlFileCalib<<"     <!-- version (string) and subversion (number) -->\n";
      xmlFileCalib<<"     <!-- fields are used to read data back from the database -->\n\n";
      xmlFileCalib<<"     <VERSION>"<<run_number<<dataset_seq_number<<"</VERSION>\n";
      xmlFileCalib<<"     <SUBVERSION>1</SUBVERSION>\n";
      xmlFileCalib<<"     <!--  Assign predefined dataset attributes -->\n\n";
      xmlFileCalib<<"     <PREDEFINED_ATTRIBUTES>\n";
      xmlFileCalib<<"        <ATTRIBUTE>\n";
      xmlFileCalib<<"           <NAME>HCAL Dataset Status</NAME>\n";
      xmlFileCalib<<"           <VALUE>VALID</VALUE>\n";
      xmlFileCalib<<"        </ATTRIBUTE>\n";
      xmlFileCalib<<"     </PREDEFINED_ATTRIBUTES>\n";
      xmlFileCalib<<"     <!-- multiple data block records -->\n\n";

      for(int sd=1;sd<=4;sd++) for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
         std::string subdet="";
         if(sd==1) subdet="HB";
         if(sd==2) subdet="HE";
         if(sd==3) subdet="HO";
         if(sd==4) subdet="HF";
         if((calib_data[sd][eta+2][phi-1].get_statistics())>100){
             double e=0,e_rms=0,t=0,t_rms=0;
         Status=calib_data[sd][eta+2][phi-1].get_status();
             Statistic=calib_data[sd][eta+2][phi-1].get_statistics1(); 
         calib_data[sd][eta+2][phi-1].get_average_amp1(&e,&e_rms);
         calib_data[sd][eta+2][phi-1].get_average_time1(&t,&t_rms);
             xmlFileCalib<<"       <DATA>\n";
             xmlFileCalib<<"          <NUMBER_OF_EVENTS_USED>"<<Statistic<<"</NUMBER_OF_EVENTS_USED>\n";
             xmlFileCalib<<"          <SIGNAL_MEAN>"<<e<<"</SIGNAL_MEAN>\n";
             xmlFileCalib<<"          <SIGNAL_RMS>"<<e_rms<<"</SIGNAL_RMS>\n";
             xmlFileCalib<<"          <TIME_MEAN>"<<t<<"</TIME_MEAN>\n";
             xmlFileCalib<<"          <TIME_RMS>"<<t_rms<<"</TIME_RMS>\n";
             xmlFileCalib<<"          <CHANNEL_STATUS_WORD>"<<Status<<"</CHANNEL_STATUS_WORD>\n";
             xmlFileCalib<<"          <CHANNEL_OBJECTNAME>HcalDetId</CHANNEL_OBJECTNAME>\n";
             xmlFileCalib<<"             <SUBDET>"<<subdet<<"</SUBDET>\n";
             xmlFileCalib<<"             <IETA>"<<eta<<"</IETA>\n";
             xmlFileCalib<<"             <IPHI>"<<phi<<"</IPHI>\n";
             xmlFileCalib<<"             <DEPTH>"<<0<<"</DEPTH>\n";
             xmlFileCalib<<"             <TYPE>0</TYPE>\n";
             xmlFileCalib<<"       </DATA>\n";  
         }
      }
      xmlFileCalib<<"  </DATA_SET>\n";
      xmlFileCalib<<"</ROOT>\n";
      xmlFileCalib.close();

      sprintf(str,"zip %s.zip %s %s",xmlName.c_str(),xmlName.c_str(),xmlNameCalib.c_str());
      system(str);
      sprintf(str,"rm -f %s %s",xmlName.c_str(),xmlNameCalib.c_str());
      system(str);
      sprintf(str,"mv -f %s.zip %s",xmlName.c_str(),XmlFilePath.c_str());
      system(str);

   }
   for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++)   hb_data[i][j][k].reset1();
   for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++)   he_data[i][j][k].reset1();
   for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++)   ho_data[i][j][k].reset1();
   for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++)   hf_data[i][j][k].reset1();
   for(int i=0;i<5;i++)for(int j=0;j<5;j++)for(int k=0;k<72;k++)calib_data[i][j][k].reset1();
   ievt_=0;
   dataset_seq_number++;
}

Member Data Documentation

std::string HcalDetDiagLaserMonitor::baseFolder_

private

Definition at line 321 of file HcalDetDiagLaserMonitor.cc.

HcalDetDiagLaserData HcalDetDiagLaserMonitor::calib_data[5][5][72]

private

Definition at line 374 of file HcalDetDiagLaserMonitor.cc.

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

bool HcalDetDiagLaserMonitor::createHTMLonly

private

Definition at line 329 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), bookHistograms(), endRun(), HcalDetDiagLaserMonitor(), and LoadDataset().

int HcalDetDiagLaserMonitor::dataset_seq_number

private

Definition at line 313 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), HcalDetDiagLaserMonitor(), LoadDataset(), SaveRaddamData(), and SaveReference().

std::string HcalDetDiagLaserMonitor::DatasetName

private

Definition at line 327 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), HcalDetDiagLaserMonitor(), and LoadDataset().

const HcalElectronicsMap* HcalDetDiagLaserMonitor::emap

private

Definition at line 289 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), fillProblems(), HcalDetDiagLaserMonitor(), and SaveReference().

MonitorElement* HcalDetDiagLaserMonitor::Energy2Dhbhehf

private

Definition at line 348 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::Energy2Dho

private

Definition at line 349 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

HcalDetDiagLaserData HcalDetDiagLaserMonitor::hb_data[85][72][4]

private

Definition at line 370 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), fillHistos(), fillProblems(), get_ave_rbx(), get_ave_subdet(), LoadDataset(), LoadReference(), and SaveReference().

MonitorElement* HcalDetDiagLaserMonitor::hb_time_rbx

private

Definition at line 355 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hbheEnergy

private

Definition at line 333 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hbheEnergyRMS

private

Definition at line 335 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

bool HcalDetDiagLaserMonitor::HBHEseq_

private

Definition at line 381 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

MonitorElement* HcalDetDiagLaserMonitor::hbheTime

private

Definition at line 334 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hbheTimeRMS

private

Definition at line 336 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

HcalDetDiagLaserData HcalDetDiagLaserMonitor::he_data[85][72][4]

private

Definition at line 371 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), fillHistos(), fillProblems(), get_ave_rbx(), get_ave_subdet(), LoadDataset(), LoadReference(), and SaveReference().

MonitorElement* HcalDetDiagLaserMonitor::he_time_rbx

private

Definition at line 356 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

HcalDetDiagLaserData HcalDetDiagLaserMonitor::hf_data[85][72][4]

private

Definition at line 373 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), fillHistos(), fillProblems(), get_ave_rbx(), get_ave_subdet(), LoadDataset(), LoadReference(), and SaveReference().

MonitorElement* HcalDetDiagLaserMonitor::hf_time_rbx

private

Definition at line 358 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hfEnergy

private

Definition at line 341 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hfEnergyRMS

private

Definition at line 343 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

bool HcalDetDiagLaserMonitor::HFseq_

private

Definition at line 381 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

MonitorElement* HcalDetDiagLaserMonitor::hfTime

private

Definition at line 342 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hfTimeRMS

private

Definition at line 344 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

HcalDetDiagLaserData HcalDetDiagLaserMonitor::ho_data[85][72][4]

private

Definition at line 372 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), fillHistos(), fillProblems(), get_ave_rbx(), get_ave_subdet(), LoadDataset(), LoadReference(), and SaveReference().

MonitorElement* HcalDetDiagLaserMonitor::ho_time_rbx

private

Definition at line 357 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hoEnergy

private

Definition at line 337 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), fillHistos(), and Jet.Jet::hoEnergyFraction().

MonitorElement* HcalDetDiagLaserMonitor::hoEnergyRMS

private

Definition at line 339 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

bool HcalDetDiagLaserMonitor::HOseq_

private

Definition at line 381 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

MonitorElement* HcalDetDiagLaserMonitor::hoTime

private

Definition at line 338 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::hoTimeRMS

private

Definition at line 340 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement * HcalDetDiagLaserMonitor::htmlFolder

private

Definition at line 331 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms().

std::string HcalDetDiagLaserMonitor::htmlOutputPath

private

Definition at line 328 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and HcalDetDiagLaserMonitor().

int HcalDetDiagLaserMonitor::ievt_

private

Definition at line 311 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), endRun(), HcalDetDiagLaserMonitor(), LoadDataset(), and SaveReference().

int HcalDetDiagLaserMonitor::ievt_hbhe_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::ievt_hf_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::ievt_ho_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

edm::InputTag HcalDetDiagLaserMonitor::inputLabelDigi_

private

Definition at line 290 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor().

bool HcalDetDiagLaserMonitor::IsReference

private

Definition at line 314 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor(), and LoadReference().

std::map<unsigned int, int> HcalDetDiagLaserMonitor::KnownBadCells_

private

Definition at line 376 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillProblems().

double HcalDetDiagLaserMonitor::LaserEnergyThreshold

private

Definition at line 309 of file HcalDetDiagLaserMonitor.cc.

Referenced by fillProblems(), and HcalDetDiagLaserMonitor().

double HcalDetDiagLaserMonitor::LaserTimingThreshold

private

Definition at line 309 of file HcalDetDiagLaserMonitor.cc.

Referenced by fillProblems(), and HcalDetDiagLaserMonitor().

int HcalDetDiagLaserMonitor::lastHBHEorbit_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::lastHForbit_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::lastHOorbit_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

bool HcalDetDiagLaserMonitor::LocalRun

private

Definition at line 315 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), endRun(), and HcalDetDiagLaserMonitor().

MonitorElement* HcalDetDiagLaserMonitor::meEVT_

private

Definition at line 331 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), bookHistograms(), and LoadDataset().

MonitorElement * HcalDetDiagLaserMonitor::meRUN_

private

Definition at line 331 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), bookHistograms(), and LoadDataset().

int HcalDetDiagLaserMonitor::nChecksHBHE_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::nChecksHF_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::nChecksHO_

private

Definition at line 382 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze().

int HcalDetDiagLaserMonitor::nHB

private

Definition at line 316 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), HcalDetDiagLaserMonitor(), and LoadDataset().

int HcalDetDiagLaserMonitor::nHBHEchecks

private

Definition at line 308 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor().

int HcalDetDiagLaserMonitor::nHE

private

Definition at line 316 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), HcalDetDiagLaserMonitor(), and LoadDataset().

int HcalDetDiagLaserMonitor::nHF

private

Definition at line 316 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), HcalDetDiagLaserMonitor(), and LoadDataset().

int HcalDetDiagLaserMonitor::nHFchecks

private

Definition at line 308 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor().

int HcalDetDiagLaserMonitor::nHO

private

Definition at line 316 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), HcalDetDiagLaserMonitor(), and LoadDataset().

int HcalDetDiagLaserMonitor::nHOchecks

private

Definition at line 308 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor().

bool HcalDetDiagLaserMonitor::Online_

private

Definition at line 323 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), endRun(), and HcalDetDiagLaserMonitor().

std::string HcalDetDiagLaserMonitor::OutputFilePath

private

Definition at line 319 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor(), SaveRaddamData(), and SaveReference().

bool HcalDetDiagLaserMonitor::Overwrite

private

Definition at line 324 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor(), SaveRaddamData(), and SaveReference().

std::string HcalDetDiagLaserMonitor::prefixME_

private

Definition at line 322 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and HcalDetDiagLaserMonitor().

EtaPhiHists* HcalDetDiagLaserMonitor::ProblemCellsByDepth_energy

private

Definition at line 364 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), fillProblems(), HcalDetDiagLaserMonitor(), and ~HcalDetDiagLaserMonitor().

EtaPhiHists* HcalDetDiagLaserMonitor::ProblemCellsByDepth_energy_val

private

Definition at line 368 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), fillProblems(), HcalDetDiagLaserMonitor(), and ~HcalDetDiagLaserMonitor().

EtaPhiHists* HcalDetDiagLaserMonitor::ProblemCellsByDepth_timing

private

Definition at line 363 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), fillProblems(), HcalDetDiagLaserMonitor(), and ~HcalDetDiagLaserMonitor().

EtaPhiHists* HcalDetDiagLaserMonitor::ProblemCellsByDepth_timing_val

private

Definition at line 367 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), fillProblems(), HcalDetDiagLaserMonitor(), and ~HcalDetDiagLaserMonitor().

std::vector<std::string> HcalDetDiagLaserMonitor::problemnames_

private

Definition at line 365 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms().

MonitorElement* HcalDetDiagLaserMonitor::Raddam[56]

private

Definition at line 360 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and bookHistograms().

HcalRaddamData HcalDetDiagLaserMonitor::Raddam_data[56]

private

Definition at line 378 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and SaveRaddamData().

bool HcalDetDiagLaserMonitor::RaddamRun

private

Definition at line 315 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), endRun(), and HcalDetDiagLaserMonitor().

double HcalDetDiagLaserMonitor::RaddamThreshold1

private

Definition at line 309 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

double HcalDetDiagLaserMonitor::RaddamThreshold2

private

Definition at line 309 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

MonitorElement* HcalDetDiagLaserMonitor::refEnergy2Dhbhehf

private

Definition at line 352 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::refEnergy2Dho

private

Definition at line 353 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

std::string HcalDetDiagLaserMonitor::ReferenceData

private

Definition at line 317 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor(), and LoadReference().

std::string HcalDetDiagLaserMonitor::ReferenceRun

private

Definition at line 318 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and LoadReference().

MonitorElement* HcalDetDiagLaserMonitor::RefRun_

private

Definition at line 332 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms().

MonitorElement* HcalDetDiagLaserMonitor::refTime2Dhbhehf

private

Definition at line 350 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::refTime2Dho

private

Definition at line 351 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

int HcalDetDiagLaserMonitor::run_number

private

Definition at line 312 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), bookHistograms(), HcalDetDiagLaserMonitor(), LoadDataset(), SaveRaddamData(), and SaveReference().

TH1F* HcalDetDiagLaserMonitor::S1[56]

private

Definition at line 361 of file HcalDetDiagLaserMonitor.cc.

Referenced by SaveRaddamData().

TH1F * HcalDetDiagLaserMonitor::S2[56]

private

Definition at line 361 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and SaveRaddamData().

MonitorElement* HcalDetDiagLaserMonitor::Time2Dhbhehf

private

Definition at line 346 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

MonitorElement* HcalDetDiagLaserMonitor::Time2Dho

private

Definition at line 347 of file HcalDetDiagLaserMonitor.cc.

Referenced by bookHistograms(), and fillHistos().

edm::EDGetTokenT<HcalCalibDigiCollection> HcalDetDiagLaserMonitor::tok_calib_

private

Definition at line 293 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

edm::EDGetTokenT<HBHEDigiCollection> HcalDetDiagLaserMonitor::tok_hbhe_

private

Definition at line 295 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

edm::EDGetTokenT<HFDigiCollection> HcalDetDiagLaserMonitor::tok_hf_

private

Definition at line 297 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

edm::EDGetTokenT<HODigiCollection> HcalDetDiagLaserMonitor::tok_ho_

private

Definition at line 296 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

edm::EDGetTokenT<FEDRawDataCollection> HcalDetDiagLaserMonitor::tok_raw_

private

Definition at line 292 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

edm::EDGetTokenT<HcalTBTriggerData> HcalDetDiagLaserMonitor::tok_tb_

private

Definition at line 294 of file HcalDetDiagLaserMonitor.cc.

Referenced by analyze(), and HcalDetDiagLaserMonitor().

std::string HcalDetDiagLaserMonitor::XmlFilePath

private

Definition at line 320 of file HcalDetDiagLaserMonitor.cc.

Referenced by HcalDetDiagLaserMonitor(), SaveRaddamData(), and SaveReference().