#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< one::DQMLuminosityBlockElements >
	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	_lowHF

double	_lowHO

int	_nevents

edm::InputTag	_tagHO

edm::InputTag	_tagQIE10

edm::InputTag	_tagQIE11

edm::InputTag	_tagTrigger

edm::InputTag	_taguMN

edm::EDGetTokenT< HODigiCollection >	_tokHO

edm::EDGetTokenT< QIE10DigiCollection >	_tokQIE10

edm::EDGetTokenT< QIE11DigiCollection >	_tokQIE11

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 8 of file LEDTask.cc.

References _ledCalibrationChannels, _lowHBHE, _lowHF, _lowHO, _nevents, _tagHO, _tagQIE10, _tagQIE11, _tagTrigger, _taguMN, _tokHO, _tokQIE10, _tokQIE11, _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
   _tagQIE11 = ps.getUntrackedParameter<edm::InputTag>("tagHE",
     edm::InputTag("hcalDigis"));
   _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
     edm::InputTag("hcalDigis"));
   _tagQIE10 = 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"));
   _tokQIE11 = consumes<QIE11DigiCollection>(_tagQIE11);
   _tokHO = consumes<HODigiCollection>(_tagHO);
   _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
   _tokTrigger = consumes<HcalTBTriggerData>(_tagTrigger);
   _tokuMN = consumes<HcalUMNioDigi>(_taguMN);
 
   //    constants
   _lowHBHE = ps.getUntrackedParameter<double>("lowHBHE", 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 25 of file LEDTask.h.

References bookHistograms().

25 {}

Member Function Documentation

void LEDTask::_dump ( )

protectedvirtual

Definition at line 340 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 621 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 408 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, _lowHF, _lowHO, hcaldqm::DQModule::_ptype, _tagHO, _tagQIE10, _tagQIE11, _tokHO, _tokQIE10, _tokQIE11, _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, 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(), HcalBarrel, HcalCalibration, HcalEndcap, HcalForward, HcalOther, mps_fire::i, HODataFrame::id(), edm::EventAuxiliary::id(), HcalDetId::ieta(), edm::InputTag::instance(), HcalDetId::iphi(), edm::InputTag::label(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), QIE11DataFrame::samples(), QIE10DataFrame::samples(), HODataFrame::size(), AlCaHLTBitMon_QueryRunRegistry::string, HcalDetId::subdet(), hcaldqm::utilities::sumQ(), protons_cff::time, and Logger::warn().

Referenced by endRun().

                                                                    {
   edm::Handle<HODigiCollection>   c_ho;
   edm::Handle<QIE10DigiCollection>    c_QIE10;
   edm::Handle<QIE11DigiCollection>    c_QIE11;
 
   if (!e.getByToken(_tokHO, c_ho))
     _logger.dqmthrow("Collection HODigiCollection isn't available "
       + _tagHO.label() + " " + _tagHO.instance());
   if (!e.getByToken(_tokQIE10, c_QIE10))
     _logger.dqmthrow("Collection QIE10DigiCollection isn't available "
       + _tagQIE10.label() + " " + _tagQIE10.instance());
   if (!e.getByToken(_tokQIE11, c_QIE11))
     _logger.dqmthrow("Collection QIE11DigiCollection isn't available "
       + _tagQIE11.label() + " " + _tagQIE11.instance());
 
   //    int currentEvent = e.eventAuxiliary().id().event();
 
   for (QIE11DigiCollection::const_iterator it = c_QIE11->begin(); it != c_QIE11->end(); ++it) {
     const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
     HcalDetId const& did = digi.detid();
     if ((did.subdet() != HcalBarrel) && (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());
               }
             }
           } else if (std::find(_ledCalibrationChannels[HcalBarrel].begin(), _ledCalibrationChannels[HcalBarrel].end(), did) !=
               _ledCalibrationChannels[HcalBarrel].end()) {
             for (int i = 0; i < digi.samples(); i++) {
               if (_ptype == fOnline) {
                 _LED_ADCvsBX_Subdet.fill(HcalDetId(HcalBarrel, 1, 1, 1), e.bunchCrossing(), digi[i].adc());
               } else if (_ptype == fLocal) {
                 _LED_ADCvsEvn_Subdet.fill(
                     HcalDetId(HcalBarrel, 1, 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 >= _lowHBHE) {
       //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 || _ptype == fLocal) {
         for (int iTS = 0; iTS < digi.samples(); ++iTS) {
           _cADCvsTS_SubdetPM.fill(did, iTS, digi[iTS].adc()); 
         }       
       }
       if (_ptype == fOnline) {
         for (int iTS = 0; iTS < digi.samples(); ++iTS) {
           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 = c_ho->begin(); it != c_ho->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 || _ptype == fLocal) {
         for (int iTS = 0; iTS < digi.size(); ++iTS) {
           _cADCvsTS_SubdetPM.fill(did, iTS, digi.sample(iTS).adc());
         }
       }
       if (_ptype == fOnline) {
         _cSumQ_SubdetPM.fill(did, sumQ);
       }
     }
   }
 
   for (QIE10DigiCollection::const_iterator it = c_QIE10->begin(); it != c_QIE10->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 || _ptype == fLocal) {
         for (int iTS = 0; iTS < digi.samples(); ++iTS) {
           _cADCvsTS_SubdetPM.fill(did, iTS, digi[iTS].adc());
         }
       }
       if (_ptype == fOnline) {
         for (int iTS = 0; iTS < digi.samples(); ++iTS) {
           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 338 of file LEDTask.cc.

Referenced by endRun().

338 { DQTask::_resetMonitors(uf); }

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

override

Definition at line 65 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 || _ptype == fLocal) {
     _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);   
   }
   if (_ptype == fOnline) {
     _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 || _ptype == fLocal) {
     _cADCvsTS_SubdetPM.book(ib, _emap, _subsystem);   
   }
   if (_ptype == fOnline) {
     _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 28 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();
   }