#include <LaserTask.h>

Inheritance diagram for LaserTask:

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

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

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

	LaserTask (edm::ParameterSet const &)

virtual	~LaserTask ()

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

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

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

	DQTask (edm::ParameterSet const &)

virtual	~DQTask ()

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

virtual	~DQModule ()(false)

Protected Member Functions
virtual void	_dump ()

virtual bool	_isApplicable (edm::Event const &)

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

virtual void	_resetMonitors (hcaldqm::UpdateFreq)

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

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
hcaldqm::Container1D	_cADC_SubdetPM

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::ContainerProf1D	_cSignalvsBX_SubdetPM

hcaldqm::ContainerProf1D	_cSignalvsEvent_SubdetPM

hcaldqm::ContainerProf1D	_cSignalvsLS_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::ContainerProf1D	_cTimingvsBX_SubdetPM

hcaldqm::ContainerProf1D	_cTimingvsEvent_SubdetPM

hcaldqm::ContainerProf1D	_cTimingvsLS_SubdetPM

hcaldqm::electronicsmap::ElectronicsMap	_ehashmap

hcaldqm::filter::HashFilter	_filter_uTCA

hcaldqm::filter::HashFilter	_filter_VME

uint32_t	_laserType

double	_lowHBHE

double	_lowHEP17

double	_lowHF

double	_lowHO

int	_nevents

edm::InputTag	_tagHBHE

edm::InputTag	_tagHEP17

edm::InputTag	_tagHF

edm::InputTag	_tagHO

edm::InputTag	_taguMN

edm::EDGetTokenT< HBHEDigiCollection >	_tokHBHE

edm::EDGetTokenT< QIE11DigiCollection >	_tokHEP17

edm::EDGetTokenT< QIE10DigiCollection >	_tokHF

edm::EDGetTokenT< HODigiCollection >	_tokHO

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

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)

Detailed Description

Definition at line 21 of file LaserTask.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file LaserTask.cc.

References _laserType, _lowHBHE, _lowHEP17, _lowHF, _lowHO, _nevents, _tagHBHE, _tagHEP17, _tagHF, _tagHO, _taguMN, _tokHBHE, _tokHEP17, _tokHF, _tokHO, _tokuMN, and edm::ParameterSet::getUntrackedParameter().

                                              :
     DQTask(ps)
 {
     _nevents = ps.getUntrackedParameter<int>("nevents", 2000);
 
     //  tags
     _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
         edm::InputTag("hcalDigis"));
     _tagHEP17 = ps.getUntrackedParameter<edm::InputTag>("tagHEP17",
         edm::InputTag("hcalDigis"));
     _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
         edm::InputTag("hcalDigis"));
     _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
         edm::InputTag("hcalDigis"));
     _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN",
         edm::InputTag("hcalDigis"));
     _tokHBHE = consumes<HBHEDigiCollection>(_tagHBHE);
     _tokHEP17 = consumes<QIE11DigiCollection>(_tagHEP17);
     _tokHO = consumes<HODigiCollection>(_tagHO);
     _tokHF = consumes<QIE10DigiCollection>(_tagHF);
     _tokuMN = consumes<HcalUMNioDigi>(_taguMN);
 
     //  constants
     _lowHBHE = ps.getUntrackedParameter<double>("lowHBHE",
         20);
     _lowHEP17 = ps.getUntrackedParameter<double>("lowHEP17",
         20);
     _lowHO = ps.getUntrackedParameter<double>("lowHO",
         20);
     _lowHF = ps.getUntrackedParameter<double>("lowHF",
         20);
     _laserType = (uint32_t)ps.getUntrackedParameter<uint32_t>("laserType");
 }

virtual LaserTask::~LaserTask ( )

inlinevirtual

Definition at line 25 of file LaserTask.h.

References bookHistograms().

26 {}

Member Function Documentation

void LaserTask::_dump ( )

protectedvirtual

Definition at line 265 of file LaserTask.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::Container1D::fill(), hcaldqm::fOffline, hcaldqm::ContainerXXX< STDTYPE >::get(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), gen::n, hcaldqm::Container1D::reset(), and mathSSE::sqrt().

Referenced by endLuminosityBlock(), 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 != fOffline) { // 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(*it));
         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 != fOffline) { // 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 != fOffline) { // 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 LaserTask::_isApplicable ( edm::Event const & e )

protectedvirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 557 of file LaserTask.cc.

References _laserType, hcaldqm::DQModule::_ptype, _tokuMN, DEFINE_FWK_MODULE, HcalUMNioDigi::eventType(), hcaldqm::constants::EVENTTYPE_LASER, hcaldqm::fOnline, edm::Event::getByToken(), and HcalUMNioDigi::valueUserWord().

Referenced by endRun().

 {
     if (_ptype!=fOnline)
         return true;
     else 
     {
         //  fOnline mode
         edm::Handle<HcalUMNioDigi> cumn;
         if (!e.getByToken(_tokuMN, cumn))
             return false;
         
         //  event type check first
         uint8_t eventType = cumn->eventType();
         if (eventType!=constants::EVENTTYPE_LASER)
             return false;
 
         //  check if this analysis task is of the right laser type
         uint32_t laserType = cumn->valueUserWord(0);
         if (laserType==_laserType) return true;
     }
 
     return false;
 }

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

protectedvirtual

Implements hcaldqm::DQTask.

Definition at line 340 of file LaserTask.cc.

References _cADC_SubdetPM, _cShapeCut_FEDSlot, _cSignalvsBX_SubdetPM, _cSignalvsEvent_SubdetPM, _cSignalvsLS_SubdetPM, _cTimingvsBX_SubdetPM, _cTimingvsEvent_SubdetPM, _cTimingvsLS_SubdetPM, hcaldqm::DQModule::_currentLS, hcaldqm::DQTask::_dbService, _ehashmap, hcaldqm::DQModule::_logger, _lowHBHE, _lowHEP17, _lowHF, _lowHO, hcaldqm::DQModule::_ptype, _tagHBHE, _tagHEP17, _tagHF, _tagHO, _tokHBHE, _tokHEP17, _tokHF, _tokHO, _xEntries, _xSignalSum, _xSignalSum2, _xTimingSum, _xTimingSum2, ecalMGPA::adc(), hcaldqm::utilities::aveTS(), edm::SortedCollection< T, SORT >::begin(), edm::DataFrameContainer::begin(), edm::EventBase::bunchCrossing(), QIE11DataFrame::detid(), QIE10DataFrame::detid(), Logger::dqmthrow(), runTauDisplay::eid, edm::SortedCollection< T, SORT >::end(), edm::DataFrameContainer::end(), edm::EventID::event(), edm::Event::eventAuxiliary(), hcaldqm::Container1D::fill(), hcaldqm::fLocal, hcaldqm::fOffline, hcaldqm::ContainerXXX< STDTYPE >::get(), edm::Event::getByToken(), HcalEndcap, HcalForward, mps_fire::i, edm::EventAuxiliary::id(), edm::InputTag::instance(), edm::InputTag::label(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), QIE11DataFrame::samples(), QIE10DataFrame::samples(), HBHEDataFrame::size(), HODataFrame::size(), HcalDetId::subdet(), and hcaldqm::utilities::sumQ().

Referenced by endRun().

 {
     edm::Handle<HBHEDigiCollection>     chbhe;
     edm::Handle<QIE11DigiCollection>        chep17;
     edm::Handle<HODigiCollection>       cho;
     edm::Handle<QIE10DigiCollection>        chf;
 
     if (!e.getByToken(_tokHBHE, chbhe))
         _logger.dqmthrow("Collection HBHEDigiCollection isn't available "
             + _tagHBHE.label() + " " + _tagHBHE.instance());
     if (!e.getByToken(_tokHEP17, chep17))
         _logger.dqmthrow("Collection QIE11DigiCollection isn't available "
             + _tagHEP17.label() + " " + _tagHEP17.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());
 
 //  int currentEvent = e.eventAuxiliary().id().event();
     int bx = e.bunchCrossing();
     
     for (HBHEDigiCollection::const_iterator it=chbhe->begin();
         it!=chbhe->end(); ++it)
     {
         const HBHEDataFrame digi = (const HBHEDataFrame)(*it);
         double sumQ = hcaldqm::utilities::sumQ<HBHEDataFrame>(digi, 2.5, 0, 
             digi.size()-1);
         if (sumQ<_lowHBHE)
             continue;
         HcalDetId did = digi.id();
         HcalElectronicsId eid = digi.elecId();
 
         double aveTS = hcaldqm::utilities::aveTS<HBHEDataFrame>(digi, 2.5, 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)++;
 
         for (int i=0; i<digi.size(); i++)
         {
             if (_ptype != fOffline) { // hidefed2crate
                 _cShapeCut_FEDSlot.fill(eid, i, 
                     digi.sample(i).nominal_fC()-2.5);
             }
             _cADC_SubdetPM.fill(did, digi.sample(i).adc());
         }
 
         //  select based on local global
         if (_ptype==fLocal)
         {
             int currentEvent = e.eventAuxiliary().id().event();
             _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
             _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
         }
         else
         {
             _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
             _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
             _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
             _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
         }
     }
     for (QIE11DigiCollection::const_iterator it=chep17->begin(); it!=chep17->end();
         ++it)
     {
         const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
         HcalDetId const& did = digi.detid();
         if (did.subdet() != HcalEndcap) {
             continue;
         }
         uint32_t rawid = _ehashmap.lookup(did);
         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<_lowHEP17)
             continue;
 
 
         //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)++;
 
         for (int i=0; i<digi.samples(); i++)
         {
             _cShapeCut_FEDSlot.fill(eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, digi, i));
             _cADC_SubdetPM.fill(did, digi[i].adc());
         }
 
         //  select based on local global
         if (_ptype==fLocal)
         {
             int currentEvent = e.eventAuxiliary().id().event();
             _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
             _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
         }
         else
         {
             _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
             _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
             _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
             _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
         }
     }
     for (HODigiCollection::const_iterator it=cho->begin();
         it!=cho->end(); ++it)
     {
         const HODataFrame digi = (const HODataFrame)(*it);
         double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(digi, 8.5, 0, 
             digi.size()-1);
         if (sumQ<_lowHO)
             continue;
         HcalDetId did = digi.id();
         HcalElectronicsId eid = digi.elecId();
 
         double aveTS = hcaldqm::utilities::aveTS<HODataFrame>(digi, 8.5, 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)++;
 
         for (int i=0; i<digi.size(); i++)
         {
             if (_ptype != fOffline) { // hidefed2crate
                 _cShapeCut_FEDSlot.fill(eid, i, 
                     digi.sample(i).nominal_fC()-8.5);
             }
             _cADC_SubdetPM.fill(did, digi.sample(i).adc());
         }
 
         //  select based on local global
         if (_ptype==fLocal)
         {
             int currentEvent = e.eventAuxiliary().id().event();
             _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
             _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
         }
         else
         {
             _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
             _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
             _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
             _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
         }
     }
     for (QIE10DigiCollection::const_iterator it=chf->begin();
         it!=chf->end(); ++it)
     {
         const QIE10DataFrame digi = (const QIE10DataFrame)(*it);
         HcalDetId did = digi.detid();
         if (did.subdet() != HcalForward) {
             continue;
         }
         HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
 
         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);
         //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
         if (sumQ<_lowHF)
             continue;
 
         //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)++;
 
         for (int i=0; i<digi.samples(); i++)
         {
             if (_ptype != fOffline) { // hidefed2crate
                 _cShapeCut_FEDSlot.fill(eid, (int)i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i));
             }
             _cADC_SubdetPM.fill(did, digi[i].adc());
         }
 
         //  select based on local global
         if (_ptype==fLocal)
         {
             int currentEvent = e.eventAuxiliary().id().event();
             _cTimingvsEvent_SubdetPM.fill(did, currentEvent, aveTS);
             _cSignalvsEvent_SubdetPM.fill(did, currentEvent, sumQ);
         }
         else
         {
             _cTimingvsLS_SubdetPM.fill(did, _currentLS, aveTS);
             _cSignalvsLS_SubdetPM.fill(did, _currentLS, sumQ);
             _cTimingvsBX_SubdetPM.fill(did, bx, aveTS);
             _cSignalvsBX_SubdetPM.fill(did, bx, sumQ);
         }
     }
 }

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

protectedvirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 260 of file LaserTask.cc.

Referenced by endRun().

 {
     DQTask::_resetMonitors(uf);
 }

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

virtual

Reimplemented from hcaldqm::DQTask.

Definition at line 41 of file LaserTask.cc.

Referenced by ~LaserTask().

 {
     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_1000),
         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);
 
     _cADC_SubdetPM.initialize(_name, "ADC",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_128),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
 
     if (_ptype != fOffline) { // 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_3000),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_3000),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_3000),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_3000),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_3000),0);
     }
     _cTimingvsEvent_SubdetPM.initialize(_name, "TimingvsEvent",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::EventNumber(_nevents),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
     _cSignalvsEvent_SubdetPM.initialize(_name, "SignalvsEvent",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::EventNumber(_nevents),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),0);
     _cTimingvsLS_SubdetPM.initialize(_name, "TimingvsLS",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::LumiSection(_maxLS),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
     _cSignalvsLS_SubdetPM.initialize(_name, "SignalvsLS",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::LumiSection(_maxLS),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),0);
     _cTimingvsBX_SubdetPM.initialize(_name, "TimingvsBX",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),0);
     _cSignalvsBX_SubdetPM.initialize(_name, "SignalvsBX",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),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_3000),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_1000),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);
     }
 
     
     //  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);
 
     if (_ptype==fLocal)
     {
         _cTimingvsEvent_SubdetPM.book(ib, _emap, _subsystem);
         _cSignalvsEvent_SubdetPM.book(ib, _emap, _subsystem);
     }
     else
     {   
         _cTimingvsLS_SubdetPM.book(ib, _emap, _subsystem);
         _cSignalvsLS_SubdetPM.book(ib, _emap, _subsystem);
         _cTimingvsBX_SubdetPM.book(ib, _emap, _subsystem);
         _cSignalvsBX_SubdetPM.book(ib, _emap, _subsystem);
     }
 
     if (_ptype != fOffline) { // 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);
     }
     _cADC_SubdetPM.book(ib, _emap, _subsystem);
 
     _cMissing_depth.book(ib, _emap,_subsystem);
     if (_ptype != fOffline) { // hidefed2crate
         _cShapeCut_FEDSlot.book(ib, _emap, _subsystem);
         _cMissing_FEDVME.book(ib, _emap, _filter_uTCA,_subsystem);
         _cMissing_FEDuTCA.book(ib, _emap, _filter_VME,_subsystem);
     }
 
     _xSignalSum.book(_emap);
     _xSignalSum2.book(_emap);
     _xEntries.book(_emap);
     _xTimingSum.book(_emap);
     _xTimingSum2.book(_emap);
 
     _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
 }

void LaserTask::endLuminosityBlock	(	edm::LuminosityBlock const &	lb,
		edm::EventSetup const &	es
	)

virtual

Reimplemented from hcaldqm::DQTask.

Definition at line 547 of file LaserTask.cc.

References _dump(), hcaldqm::DQModule::_ptype, and hcaldqm::fLocal.

Referenced by endRun().

 {
     if (_ptype==fLocal)
         return;
     this->_dump();
 
     DQTask::endLuminosityBlock(lb, es);
 }

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

inlinevirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 30 of file LaserTask.h.

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

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