#include <LEDTask.h>

Inheritance diagram for LEDTask:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	endRun (edm::Run const &r, edm::EventSetup const &) override

	LEDTask (edm::ParameterSet const &)

	~LEDTask () override

Public Member Functions inherited from hcaldqm::DQTask
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

	DQTask (edm::ParameterSet const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

	~DQTask () override

Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from hcaldqm::DQModule
	DQModule (edm::ParameterSet const &)

virtual	~DQModule ()(false)

Protected Member Functions
virtual void	_dump ()

bool	_isApplicable (edm::Event const &) override

void	_process (edm::Event const &, edm::EventSetup const &) override

void	_resetMonitors (hcaldqm::UpdateFreq) override

Protected Member Functions inherited from hcaldqm::DQTask
virtual int	_getCalibType (edm::Event const &)

Protected Attributes
hcaldqm::Container2D	_cADCvsTS_SubdetPM

hcaldqm::ContainerSingle2D	_cLowSignal_CrateSlot

hcaldqm::Container2D	_cMissing_depth

hcaldqm::Container2D	_cMissing_FEDuTCA

hcaldqm::Container2D	_cMissing_FEDVME

hcaldqm::ContainerProf1D	_cShapeCut_FEDSlot

hcaldqm::ContainerProf2D	_cSignalMean_depth

hcaldqm::ContainerProf2D	_cSignalMean_FEDuTCA

hcaldqm::ContainerProf2D	_cSignalMean_FEDVME

hcaldqm::Container1D	_cSignalMean_Subdet

hcaldqm::ContainerProf2D	_cSignalRMS_depth

hcaldqm::ContainerProf2D	_cSignalRMS_FEDuTCA

hcaldqm::ContainerProf2D	_cSignalRMS_FEDVME

hcaldqm::Container1D	_cSignalRMS_Subdet

hcaldqm::Container1D	_cSumQ_SubdetPM

hcaldqm::ContainerProf2D	_cTDCTime_depth

hcaldqm::Container1D	_cTDCTime_SubdetPM

hcaldqm::ContainerProf2D	_cTimingMean_depth

hcaldqm::ContainerProf2D	_cTimingMean_FEDuTCA

hcaldqm::ContainerProf2D	_cTimingMean_FEDVME

hcaldqm::Container1D	_cTimingMean_Subdet

hcaldqm::ContainerProf2D	_cTimingRMS_depth

hcaldqm::ContainerProf2D	_cTimingRMS_FEDuTCA

hcaldqm::ContainerProf2D	_cTimingRMS_FEDVME

hcaldqm::Container1D	_cTimingRMS_Subdet

hcaldqm::electronicsmap::ElectronicsMap	_ehashmap

hcaldqm::filter::HashFilter	_filter_uTCA

hcaldqm::filter::HashFilter	_filter_VME

hcaldqm::Container2D	_LED_ADCvsBX_Subdet

hcaldqm::Container2D	_LED_ADCvsEvn_Subdet

std::map< HcalSubdetector, std::vector< HcalDetId > >	_ledCalibrationChannels

double	_lowHBHE

double	_lowHE

double	_lowHF

double	_lowHO

int	_nevents

edm::InputTag	_tagHBHE

edm::InputTag	_tagHE

edm::InputTag	_tagHF

edm::InputTag	_tagHO

edm::InputTag	_tagTrigger

edm::InputTag	_taguMN

edm::EDGetTokenT< HBHEDigiCollection >	_tokHBHE

edm::EDGetTokenT< QIE11DigiCollection >	_tokHE

edm::EDGetTokenT< QIE10DigiCollection >	_tokHF

edm::EDGetTokenT< HODigiCollection >	_tokHO

edm::EDGetTokenT< HcalTBTriggerData >	_tokTrigger

edm::EDGetTokenT< HcalUMNioDigi >	_tokuMN

hcaldqm::ContainerXXX< int >	_xEntries

hcaldqm::ContainerXXX< double >	_xSignalSum

hcaldqm::ContainerXXX< double >	_xSignalSum2

hcaldqm::ContainerXXX< double >	_xTimingSum

hcaldqm::ContainerXXX< double >	_xTimingSum2

Protected Attributes inherited from hcaldqm::DQTask
ContainerI	_cEvsPerLS

ContainerI	_cEvsTotal

ContainerS	_cProcessingTypeName

ContainerS	_cRunKeyName

ContainerI	_cRunKeyVal

edm::ESHandle< HcalDbService >	_dbService

HcalElectronicsMap const *	_emap

int	_procLSs

edm::InputTag	_tagRaw

edm::EDGetTokenT< FEDRawDataCollection >	_tokRaw

std::vector< uint32_t >	_vcdaqEids

ContainerXXX< uint32_t >	_xQuality

Protected Attributes inherited from hcaldqm::DQModule
int	_currentLS

int	_debug

int	_evsPerLS

int	_evsTotal

Logger	_logger

int	_maxLS

ModuleType	_mtype

std::string	_name

ProcessingType	_ptype

std::string	_runkeyName

int	_runkeyVal

std::string	_subsystem

Detailed Description

Definition at line 22 of file LEDTask.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file LEDTask.cc.

References _ledCalibrationChannels, _lowHBHE, _lowHE, _lowHF, _lowHO, _nevents, _tagHBHE, _tagHE, _tagHF, _tagHO, _tagTrigger, _taguMN, _tokHBHE, _tokHE, _tokHF, _tokHO, _tokTrigger, _tokuMN, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalOther, HcalOuter, and mps_fire::i.

                                          :
     DQTask(ps)
 {
     _nevents = ps.getUntrackedParameter<int>("nevents", 2000);
     //  tags
     _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
         edm::InputTag("hcalDigis"));
     _tagHE = ps.getUntrackedParameter<edm::InputTag>("tagHE",
         edm::InputTag("hcalDigis"));
     _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
         edm::InputTag("hcalDigis"));
     _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
         edm::InputTag("hcalDigis"));
     _tagTrigger = ps.getUntrackedParameter<edm::InputTag>("tagTrigger",
         edm::InputTag("tbunpacker"));
     _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN",
         edm::InputTag("hcalDigis"));
     _tokHBHE = consumes<HBHEDigiCollection>(_tagHBHE);
     _tokHE = consumes<QIE11DigiCollection>(_tagHE);
     _tokHO = consumes<HODigiCollection>(_tagHO);
     _tokHF = consumes<QIE10DigiCollection>(_tagHF);
     _tokTrigger = consumes<HcalTBTriggerData>(_tagTrigger);
     _tokuMN = consumes<HcalUMNioDigi>(_taguMN);
 
     //  constants
     _lowHBHE = ps.getUntrackedParameter<double>("lowHBHE",
         20);
     _lowHE = ps.getUntrackedParameter<double>("lowHE",
         20);
     _lowHO = ps.getUntrackedParameter<double>("lowHO",
         20);
     _lowHF = ps.getUntrackedParameter<double>("lowHF",
         20);
 
     // LED calibration channels
     std::vector<edm::ParameterSet> vLedCalibChannels = ps.getParameter<std::vector<edm::ParameterSet>>("ledCalibrationChannels");
     for (int i = 0; i <= 3; ++i) {
         HcalSubdetector this_subdet = HcalEmpty;
         switch (i) {
             case 0:
                 this_subdet = HcalBarrel;
                 break;
             case 1:
                 this_subdet = HcalEndcap;
                 break;
             case 2:
                 this_subdet = HcalOuter;
                 break;
             case 3:
                 this_subdet = HcalForward;
                 break;
             default:
                 this_subdet = HcalEmpty;
                 break;
         }
         std::vector<int32_t> subdet_calib_ietas = vLedCalibChannels[i].getUntrackedParameter<std::vector<int32_t>>("ieta");
         std::vector<int32_t> subdet_calib_iphis = vLedCalibChannels[i].getUntrackedParameter<std::vector<int32_t>>("iphi");
         std::vector<int32_t> subdet_calib_depths = vLedCalibChannels[i].getUntrackedParameter<std::vector<int32_t>>("depth");
         for (unsigned int ichannel = 0; ichannel < subdet_calib_ietas.size(); ++ichannel) {
             _ledCalibrationChannels[this_subdet].push_back(HcalDetId(HcalOther, subdet_calib_ietas[ichannel], subdet_calib_iphis[ichannel], subdet_calib_depths[ichannel]));
         }
     }
 
 }

LEDTask::~LEDTask ( )

inlineoverride

Definition at line 26 of file LEDTask.h.

References bookHistograms().

27 {}

Member Function Documentation

void LEDTask::_dump ( )

protectedvirtual

Definition at line 299 of file LEDTask.cc.

References _cMissing_depth, _cMissing_FEDuTCA, _cMissing_FEDVME, _cSignalMean_depth, _cSignalMean_FEDuTCA, _cSignalMean_FEDVME, _cSignalMean_Subdet, _cSignalRMS_depth, _cSignalRMS_FEDuTCA, _cSignalRMS_FEDVME, _cSignalRMS_Subdet, _cTimingMean_depth, _cTimingMean_FEDuTCA, _cTimingMean_FEDVME, _cTimingMean_Subdet, _cTimingRMS_depth, _cTimingRMS_FEDuTCA, _cTimingRMS_FEDVME, _cTimingRMS_Subdet, _ehashmap, hcaldqm::DQTask::_emap, hcaldqm::DQModule::_ptype, _xEntries, _xSignalSum, _xSignalSum2, _xTimingSum, _xTimingSum2, HcalElectronicsMap::allPrecisionId(), runTauDisplay::eid, hcaldqm::Container2D::fill(), hcaldqm::ContainerProf2D::fill(), hcaldqm::Container1D::fill(), hcaldqm::fLocal, hcaldqm::ContainerXXX< STDTYPE >::get(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), gen::n, hcaldqm::Container1D::reset(), and mathSSE::sqrt().

Referenced by _process(), and endRun().

 {
     _cSignalMean_Subdet.reset();
     _cSignalRMS_Subdet.reset();
     _cTimingMean_Subdet.reset();
     _cTimingRMS_Subdet.reset();
     _cSignalMean_depth.reset();
     _cSignalRMS_depth.reset();
     _cTimingMean_depth.reset();
     _cTimingRMS_depth.reset();
 
     if (_ptype == fLocal) { // hidefed2crate
         _cSignalMean_FEDVME.reset();
         _cSignalMean_FEDuTCA.reset();
         _cSignalRMS_FEDVME.reset();
         _cSignalRMS_FEDuTCA.reset();
         _cTimingMean_FEDVME.reset();
         _cTimingMean_FEDuTCA.reset();
         _cTimingRMS_FEDVME.reset();
         _cTimingRMS_FEDuTCA.reset();
     }
 
     std::vector<HcalGenericDetId> dids = _emap->allPrecisionId();
     for (std::vector<HcalGenericDetId>::const_iterator it=dids.begin();
         it!=dids.end(); ++it)
     {
         if (!it->isHcalDetId())
             continue;
         HcalDetId did = HcalDetId(it->rawId());
         HcalElectronicsId eid(_ehashmap.lookup(did));
         int n = _xEntries.get(did);
         double msig = _xSignalSum.get(did)/n; 
         double mtim = _xTimingSum.get(did)/n;
         double rsig = sqrt(_xSignalSum2.get(did)/n-msig*msig);
         double rtim = sqrt(_xTimingSum2.get(did)/n-mtim*mtim);
 
         //  channels missing or low signal
         if (n==0)
         {
             _cMissing_depth.fill(did);
             if (_ptype == fLocal) { // hidefed2crate
                 if (eid.isVMEid())
                     _cMissing_FEDVME.fill(eid);
                 else
                     _cMissing_FEDuTCA.fill(eid);
             }
             continue;
         }
         _cSignalMean_Subdet.fill(did, msig);
         _cSignalMean_depth.fill(did, msig);
         _cSignalRMS_Subdet.fill(did, rsig);
         _cSignalRMS_depth.fill(did, rsig);
         _cTimingMean_Subdet.fill(did, mtim);
         _cTimingMean_depth.fill(did, mtim);
         _cTimingRMS_Subdet.fill(did, rtim);
         _cTimingRMS_depth.fill(did, rtim);
         if (_ptype == fLocal) { // hidefed2crate
             if (eid.isVMEid())
             {
                 _cSignalMean_FEDVME.fill(eid, msig);
                 _cSignalRMS_FEDVME.fill(eid, rsig);
                 _cTimingMean_FEDVME.fill(eid, mtim);
                 _cTimingRMS_FEDVME.fill(eid, rtim);
             }
             else
             {
                 _cSignalMean_FEDuTCA.fill(eid, msig);
                 _cSignalRMS_FEDuTCA.fill(eid, rsig);
                 _cTimingMean_FEDuTCA.fill(eid, mtim);
                 _cTimingRMS_FEDuTCA.fill(eid, rtim);
             }
         }
     }
 }

bool LEDTask::_isApplicable ( edm::Event const & e )

overrideprotectedvirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 610 of file LEDTask.cc.

References hcaldqm::DQModule::_logger, hcaldqm::DQModule::_ptype, _tagTrigger, _tokTrigger, _tokuMN, DEFINE_FWK_MODULE, Logger::dqmthrow(), HcalUMNioDigi::eventType(), hcaldqm::constants::EVENTTYPE_LED, hcaldqm::fOnline, edm::Event::getByToken(), edm::InputTag::instance(), edm::InputTag::label(), and HcalTBTriggerData::wasLEDTrigger().

Referenced by endRun().

 {
     if (_ptype!=fOnline)
     {
         //  local
         edm::Handle<HcalTBTriggerData> ctrigger;
         if (!e.getByToken(_tokTrigger, ctrigger))
             _logger.dqmthrow("Collection HcalTBTriggerData isn't available "
                 + _tagTrigger.label() + " " + _tagTrigger.instance());
         return ctrigger->wasLEDTrigger();
     } else {
         //  fOnline mode
         edm::Handle<HcalUMNioDigi> cumn;
         if (!e.getByToken(_tokuMN, cumn)) {
             return false;
         }
         
         return (cumn->eventType() == constants::EVENTTYPE_LED);
     }
 
     return false;
 }

void LEDTask::_process	(	edm::Event const &	e,
		edm::EventSetup const &	es
	)

overrideprotectedvirtual

Implements hcaldqm::DQTask.

Definition at line 374 of file LEDTask.cc.

References _cADCvsTS_SubdetPM, _cLowSignal_CrateSlot, _cShapeCut_FEDSlot, _cSumQ_SubdetPM, _cTDCTime_depth, _cTDCTime_SubdetPM, hcaldqm::DQTask::_dbService, _dump(), _ehashmap, hcaldqm::DQModule::_evsTotal, _LED_ADCvsBX_Subdet, _LED_ADCvsEvn_Subdet, _ledCalibrationChannels, hcaldqm::DQModule::_logger, _lowHBHE, _lowHE, _lowHF, _lowHO, hcaldqm::DQModule::_ptype, _tagHBHE, _tagHE, _tagHF, _tagHO, _tokHBHE, _tokHE, _tokHF, _tokHO, _xEntries, _xSignalSum, _xSignalSum2, _xTimingSum, _xTimingSum2, ecalMGPA::adc(), hcaldqm::utilities::aveTS(), begin, edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), edm::EventBase::bunchCrossing(), hcaldqm::constants::CALIBEVENTS_MIN, HcalDetId::depth(), QIE11DataFrame::detid(), QIE10DataFrame::detid(), Logger::dqmthrow(), runTauDisplay::eid, HBHEDataFrame::elecId(), HODataFrame::elecId(), end, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), edm::EventID::event(), edm::Event::eventAuxiliary(), hcaldqm::Container2D::fill(), hcaldqm::ContainerSingle2D::fill(), hcaldqm::ContainerProf2D::fill(), hcaldqm::Container1D::fill(), spr::find(), hcaldqm::fLocal, hcaldqm::fOnline, hcaldqm::ContainerXXX< STDTYPE >::get(), edm::Event::getByToken(), HcalCalibration, HcalEndcap, HcalForward, HcalOther, mps_fire::i, HODataFrame::id(), HBHEDataFrame::id(), edm::EventAuxiliary::id(), HcalDetId::ieta(), edm::InputTag::instance(), HcalDetId::iphi(), edm::InputTag::label(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), QIE11DataFrame::samples(), QIE10DataFrame::samples(), HODataFrame::size(), HBHEDataFrame::size(), AlCaHLTBitMon_QueryRunRegistry::string, HcalDetId::subdet(), hcaldqm::utilities::sumQ(), ntuplemaker::time, and Logger::warn().

Referenced by endRun().

 {
     edm::Handle<HBHEDigiCollection>     chbhe;
     edm::Handle<HODigiCollection>       cho;
     edm::Handle<QIE10DigiCollection>        chf;
     edm::Handle<QIE11DigiCollection>        che;
 
     if (!e.getByToken(_tokHBHE, chbhe))
         _logger.dqmthrow("Collection HBHEDigiCollection isn't available "
             + _tagHBHE.label() + " " + _tagHBHE.instance());
     if (!e.getByToken(_tokHO, cho))
         _logger.dqmthrow("Collection HODigiCollection isn't available "
             + _tagHO.label() + " " + _tagHO.instance());
     if (!e.getByToken(_tokHF, chf))
         _logger.dqmthrow("Collection QIE10DigiCollection isn't available "
             + _tagHF.label() + " " + _tagHF.instance());
     if (!e.getByToken(_tokHE, che))
         _logger.dqmthrow("Collection QIE11DigiCollection isn't available "
             + _tagHE.label() + " " + _tagHE.instance());
 
 //  int currentEvent = e.eventAuxiliary().id().event();
 
     for (HBHEDigiCollection::const_iterator it=chbhe->begin();
         it!=chbhe->end(); ++it)
     {
         const HBHEDataFrame digi = (const HBHEDataFrame)(*it);
         HcalDetId did = digi.id();
         HcalElectronicsId eid = digi.elecId();
 
         // Get total charge and apply charge cut
         CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HBHEDataFrame>(_dbService, did, digi);
         //double sumQ = hcaldqm::utilities::sumQ<HBHEDataFrame>(digi, 2.5, 0, digi.size()-1);
         double sumQ = hcaldqm::utilities::sumQDB<HBHEDataFrame>(_dbService, digi_fC, did, digi, 0, digi.size()-1);
         if (sumQ >= _lowHBHE) {
             //double aveTS = hcaldqm::utilities::aveTS<HBHEDataFrame>(digi, 2.5, 0,digi.size()-1);
             double aveTS = hcaldqm::utilities::aveTSDB<HBHEDataFrame>(_dbService, digi_fC, did, digi, 0, digi.size()-1);
 
             _xSignalSum.get(did)+=sumQ;
             _xSignalSum2.get(did)+=sumQ*sumQ;
             _xTimingSum.get(did)+=aveTS;
             _xTimingSum2.get(did)+=aveTS*aveTS;
             _xEntries.get(did)++;
 
             if (_ptype == fLocal) { // hidefed2crate
                 for (int i=0; i<digi.size(); i++) {
                     //_cShapeCut_FEDSlot.fill(eid, i, digi.sample(i).nominal_fC()-2.5);
                     _cShapeCut_FEDSlot.fill(eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<HBHEDataFrame>(_dbService, digi_fC, did, digi, i));
                 }
             }
 
             if (_ptype == fOnline) {
                 for (int iTS = 0; iTS < digi.size(); ++iTS) {
                     _cADCvsTS_SubdetPM.fill(did, iTS, digi.sample(iTS).adc());                  
                 }
                 _cSumQ_SubdetPM.fill(did, sumQ);
             }
         }
     }
 
     for (QIE11DigiCollection::const_iterator it=che->begin(); it!=che->end();
         ++it)
     {
         const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
         HcalDetId const& did = digi.detid();
         if (did.subdet() != HcalEndcap) {
             // LED monitoring from calibration channels
             if (did.subdet() == HcalOther) {
                 HcalOtherDetId hodid(digi.detid());
                 if (hodid.subdet() == HcalCalibration) {
                     if (std::find(_ledCalibrationChannels[HcalEndcap].begin(), _ledCalibrationChannels[HcalEndcap].end(), did) != _ledCalibrationChannels[HcalEndcap].end()) {
                         for (int i=0; i<digi.samples(); i++) {
                             if (_ptype == fOnline) {
                                 _LED_ADCvsBX_Subdet.fill(HcalDetId(HcalEndcap, 16, 1, 1), e.bunchCrossing(), digi[i].adc());
                             } else if (_ptype == fLocal) {
                                 _LED_ADCvsEvn_Subdet.fill(HcalDetId(HcalEndcap, 16, 1, 1), e.eventAuxiliary().id().event(), digi[i].adc());
                             }
                         }
                     }
                 }
             }
 
             continue;
         }
         uint32_t rawid = _ehashmap.lookup(did);
         if (!rawid) {
           std::string unknown_id_string="Detid "+std::to_string(int(did))+", ieta "+std::to_string(did.ieta());
           unknown_id_string+=", iphi "+std::to_string(did.iphi())+", depth "+std::to_string(did.depth());
           unknown_id_string+=", is not in emap. Skipping.";
           _logger.warn(unknown_id_string);
           continue;
         }
         HcalElectronicsId const& eid(rawid);
 
         CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, digi);
         //double sumQ = hcaldqm::utilities::sumQ_v10<QIE11DataFrame>(digi, 2.5, 0, digi.samples()-1);
         double sumQ = hcaldqm::utilities::sumQDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
         if (sumQ >= _lowHE) {
             //double aveTS = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(digi, 2.5, 0,digi.samples()-1);
             double aveTS = hcaldqm::utilities::aveTSDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.size()-1);
 
             _xSignalSum.get(did)+=sumQ;
             _xSignalSum2.get(did)+=sumQ*sumQ;
             _xTimingSum.get(did)+=aveTS;
             _xTimingSum2.get(did)+=aveTS*aveTS;
             _xEntries.get(did)++;
 
             if (_ptype == fLocal) { // hidefed2crate
                 for (int i=0; i<digi.samples(); i++) {
                     //_cShapeCut_FEDSlot.fill(eid, i, digi.sample(i).nominal_fC()-2.5);
                     _cShapeCut_FEDSlot.fill(eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, digi, i));
                 }
             }
             if (_ptype == fOnline) {
                 for (int iTS = 0; iTS < digi.samples(); ++iTS) {
                     _cADCvsTS_SubdetPM.fill(did, iTS, digi[iTS].adc());                 
                     if (digi[iTS].tdc() <50) {
                         double time = iTS*25. + (digi[iTS].tdc() / 2.);
                         _cTDCTime_SubdetPM.fill(did, time);
                         _cTDCTime_depth.fill(did, time);
                     }
                 }
                 _cSumQ_SubdetPM.fill(did, sumQ);
 
                 // Low signal in SOI
                 short soi = -1;
                 for (int i=0; i<digi.samples(); i++) {
                     if (digi[i].soi()) {
                         soi = i;
                         break;
                     }
                 }
                 if (digi[soi].adc() < 30) {
                     _cLowSignal_CrateSlot.fill(eid);
                 }   
             }
         }           
     }
     for (HODigiCollection::const_iterator it=cho->begin();
         it!=cho->end(); ++it)
     {
         const HODataFrame digi = (const HODataFrame)(*it);
         HcalDetId did = digi.id();
         HcalElectronicsId eid = digi.elecId();
         //double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(digi, 8.5, 0, digi.size()-1);
         CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HODataFrame>(_dbService, did, digi);
         double sumQ = hcaldqm::utilities::sumQDB<HODataFrame>(_dbService, digi_fC, did, digi, 0, digi.size()-1);
         if (sumQ >= _lowHO) {
             //double aveTS = hcaldqm::utilities::aveTS<HODataFrame>(digi, 8.5, 0, digi.size()-1);
             double aveTS = hcaldqm::utilities::aveTSDB<HODataFrame>(_dbService, digi_fC, did, digi, 0, digi.size()-1);
 
             _xSignalSum.get(did)+=sumQ;
             _xSignalSum2.get(did)+=sumQ*sumQ;
             _xTimingSum.get(did)+=aveTS;
             _xTimingSum2.get(did)+=aveTS*aveTS;
             _xEntries.get(did)++;
 
             if (_ptype == fLocal) { // hidefed2crate
                 for (int i=0; i<digi.size(); i++) {
                     //_cShapeCut_FEDSlot.fill(eid, i, digi.sample(i).nominal_fC()-8.5);
                     _cShapeCut_FEDSlot.fill(eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<HODataFrame>(_dbService, digi_fC, did, digi, i));
                 }
             }
             if (_ptype == fOnline) {
                 for (int iTS = 0; iTS < digi.size(); ++iTS) {
                     _cADCvsTS_SubdetPM.fill(did, iTS, digi.sample(iTS).adc());                  
                 }
                 _cSumQ_SubdetPM.fill(did, sumQ);
             }
         }
     }
 
     for (QIE10DigiCollection::const_iterator it=chf->begin();
         it!=chf->end(); ++it)
     {
         const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
         HcalDetId did = digi.detid();
         if (did.subdet() != HcalForward) {
             // LED monitoring from calibration channels
             if (did.subdet() == HcalOther) {
                 HcalOtherDetId hodid(digi.detid());
                 if (hodid.subdet() == HcalCalibration) {
                     if (std::find(_ledCalibrationChannels[HcalForward].begin(), _ledCalibrationChannels[HcalForward].end(), did) != _ledCalibrationChannels[HcalForward].end()) {
                         for (int i=0; i<digi.samples(); i++) {
                             if (_ptype == fOnline) {
                                 _LED_ADCvsBX_Subdet.fill(HcalDetId(HcalForward, 29, 1, 1), e.bunchCrossing(), digi[i].adc());
                             } else if (_ptype == fLocal) {
                                 _LED_ADCvsEvn_Subdet.fill(HcalDetId(HcalForward, 29, 1, 1), e.eventAuxiliary().id().event(), digi[i].adc());
                             }
                         }
                     }
                 }
             }
             continue;
         }
         HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
         //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
         CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, digi);
         double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples()-1);
         if (sumQ >= _lowHF) {
             //double aveTS = hcaldqm::utilities::aveTS_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
             double aveTS = hcaldqm::utilities::aveTSDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.size()-1);
 
             _xSignalSum.get(did)+=sumQ;
             _xSignalSum2.get(did)+=sumQ*sumQ;
             _xTimingSum.get(did)+=aveTS;
             _xTimingSum2.get(did)+=aveTS*aveTS;
             _xEntries.get(did)++;
 
             if (_ptype == fLocal) { // hidefed2crate
                 for (int i = 0; i < digi.samples(); ++i) {
                     // Note: this used to be digi.sample(i).nominal_fC() - 2.5, but this branch doesn't exist in QIE10DataFrame.
                     // Instead, use lookup table.
                     //_cShapeCut_FEDSlot.fill(eid, i, constants::adc2fC[digi[i].adc()]);
                     _cShapeCut_FEDSlot.fill(eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i));
                 }
             }
             if (_ptype == fOnline) {
                 for (int iTS = 0; iTS < digi.samples(); ++iTS) {
                     _cADCvsTS_SubdetPM.fill(did, iTS, digi[iTS].adc());                 
                     if (digi[iTS].le_tdc() <50) {
                         double time = iTS*25. + (digi[iTS].le_tdc() / 2.);
                         _cTDCTime_SubdetPM.fill(did, time);
                         _cTDCTime_depth.fill(did, time);
                     }
                 }
                 _cSumQ_SubdetPM.fill(did, sumQ);
             }
         }           
     }
 
     if (_ptype==fOnline && _evsTotal>0 &&
         _evsTotal%constants::CALIBEVENTS_MIN==0)
         this->_dump();
 }

void LEDTask::_resetMonitors ( hcaldqm::UpdateFreq uf )

overrideprotectedvirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 294 of file LEDTask.cc.

Referenced by endRun().

 {
     DQTask::_resetMonitors(uf);
 }

void LEDTask::bookHistograms	(	DQMStore::IBooker &	ib,
		edm::Run const &	r,
		edm::EventSetup const &	es
	)

override

Definition at line 72 of file LEDTask.cc.

Referenced by ~LEDTask().

 {
     if (_ptype==fLocal)
         if (r.runAuxiliary().run()==1)
             return;
 
     DQTask::bookHistograms(ib, r, es);
     
     edm::ESHandle<HcalDbService> dbService;
     es.get<HcalDbRecord>().get(dbService);
     _emap = dbService->getHcalMapping();
 
     std::vector<uint32_t> vhashVME;
     std::vector<uint32_t> vhashuTCA;
     std::vector<uint32_t> vhashC36;
     vhashVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
         constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
     vhashuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
         FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
     _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
         vhashVME);
     _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
         vhashuTCA);
 
     //  INITIALIZE
     _cSignalMean_Subdet.initialize(_name, "SignalMean",
         hcaldqm::hashfunctions::fSubdet, 
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
     _cSignalRMS_Subdet.initialize(_name, "SignalRMS",
         hcaldqm::hashfunctions::fSubdet, 
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
     _cTimingMean_Subdet.initialize(_name, "TimingMean",
         hcaldqm::hashfunctions::fSubdet, 
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
     _cTimingRMS_Subdet.initialize(_name, "TimingRMS",
         hcaldqm::hashfunctions::fSubdet, 
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200), 
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
 
     if (_ptype == fLocal) { // hidefed2crate
         _cSignalMean_FEDVME.initialize(_name, "SignalMean",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),0);
         _cSignalMean_FEDuTCA.initialize(_name, "SignalMean",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),0);
         _cSignalRMS_FEDVME.initialize(_name, "SignalRMS",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_10000, true),0);
         _cSignalRMS_FEDuTCA.initialize(_name, "SignalRMS",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_10000, true),0);
         _cTimingMean_FEDVME.initialize(_name, "TimingMean",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
         _cTimingMean_FEDuTCA.initialize(_name, "TimingMean",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
         _cTimingRMS_FEDVME.initialize(_name, "TimingRMS",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
         _cTimingRMS_FEDuTCA.initialize(_name, "TimingRMS",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
 
         _cShapeCut_FEDSlot.initialize(_name, "Shape", 
             hcaldqm::hashfunctions::fFEDSlot,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000),0);
     }
 
     _cSignalMean_depth.initialize(_name, "SignalMean",
         hcaldqm::hashfunctions::fdepth, 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta), 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),0);
     _cSignalRMS_depth.initialize(_name, "SignalRMS",
         hcaldqm::hashfunctions::fdepth, 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta), 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_10000, true),0);
     _cTimingMean_depth.initialize(_name, "TimingMean",
         hcaldqm::hashfunctions::fdepth, 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta), 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
     _cTimingRMS_depth.initialize(_name, "TimingRMS",
         hcaldqm::hashfunctions::fdepth,
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta), 
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
 
     _cMissing_depth.initialize(_name, "Missing",
         hcaldqm::hashfunctions::fdepth,
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
     if (_ptype != fOffline) { // hidefed2crate
         _cMissing_FEDVME.initialize(_name, "Missing",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cMissing_FEDuTCA.initialize(_name, "Missing",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
     }
 
     // Plots for LED in global
     if (_ptype == fOnline) {
         _cADCvsTS_SubdetPM.initialize(_name, "ADCvsTS",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
         _cLowSignal_CrateSlot.initialize(_name, "LowSignal", 
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fCrateuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
         _cSumQ_SubdetPM.initialize(_name, "SumQ", 
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
         _cTDCTime_SubdetPM.initialize(_name, "TDCTime", 
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
         _cTDCTime_depth.initialize(_name, "TDCTime", 
             hcaldqm::hashfunctions::fdepth,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),          
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),0);
         _LED_ADCvsBX_Subdet.initialize(_name, "LED_ADCvsBX", 
             hcaldqm::hashfunctions::fSubdet, 
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX_36),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_256_4),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);     
     } else if (_ptype == fLocal) {
         _LED_ADCvsEvn_Subdet.initialize(_name, "LED_ADCvsEvn", 
             hcaldqm::hashfunctions::fSubdet, 
             new hcaldqm::quantity::EventNumber(_nevents),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_256_4),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
     }
     
     //  initialize compact containers
     _xSignalSum.initialize(hcaldqm::hashfunctions::fDChannel);
     _xSignalSum2.initialize(hcaldqm::hashfunctions::fDChannel);
     _xTimingSum.initialize(hcaldqm::hashfunctions::fDChannel);
     _xTimingSum2.initialize(hcaldqm::hashfunctions::fDChannel);
     _xEntries.initialize(hcaldqm::hashfunctions::fDChannel);
 
     //  BOOK
     _cSignalMean_Subdet.book(ib, _emap, _subsystem);
     _cSignalRMS_Subdet.book(ib, _emap, _subsystem);
     _cTimingMean_Subdet.book(ib, _emap, _subsystem);
     _cTimingRMS_Subdet.book(ib, _emap, _subsystem);
 
     _cSignalMean_depth.book(ib, _emap, _subsystem);
     _cSignalRMS_depth.book(ib, _emap, _subsystem);
     _cTimingMean_depth.book(ib, _emap, _subsystem);
     _cTimingRMS_depth.book(ib, _emap, _subsystem);
 
     _cMissing_depth.book(ib, _emap, _subsystem);
     if (_ptype == fLocal) { // hidefed2crate
         _cSignalMean_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cSignalMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cSignalRMS_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cSignalRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cTimingMean_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cTimingMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cTimingRMS_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cTimingRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cShapeCut_FEDSlot.book(ib, _emap, _subsystem);
         _cMissing_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cMissing_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
     }
     if (_ptype == fOnline) {
         _cADCvsTS_SubdetPM.book(ib, _emap, _subsystem);
         _cLowSignal_CrateSlot.book(ib, _subsystem);
         _cSumQ_SubdetPM.book(ib, _emap, _subsystem);
         _cTDCTime_SubdetPM.book(ib, _emap, _subsystem);
         _cTDCTime_depth.book(ib, _emap, _subsystem);
     }
 
     if (_ptype == fOnline) {
         _LED_ADCvsBX_Subdet.book(ib, _emap, _subsystem);
     } else if (_ptype == fLocal) {
         _LED_ADCvsEvn_Subdet.book(ib, _emap, _subsystem);
     }
 
     _xSignalSum.book(_emap);
     _xSignalSum2.book(_emap);
     _xEntries.book(_emap);
     _xTimingSum.book(_emap);
     _xTimingSum2.book(_emap);
 
     _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
 }

void LEDTask::endRun	(	edm::Run const &	r,
		edm::EventSetup const &
	)

inlineoverride

Definition at line 31 of file LEDTask.h.

References _dump(), _isApplicable(), _process(), hcaldqm::DQModule::_ptype, _resetMonitors(), hcaldqm::fLocal, edm::RunAuxiliary::run(), and edm::Run::runAuxiliary().

         {
             if (_ptype==hcaldqm::fLocal)
                 if (r.runAuxiliary().run()==1)
                     return;
             this->_dump();
         }