#include <RecHitTask.h>

Inheritance diagram for RecHitTask:

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

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

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

	RecHitTask (edm::ParameterSet const &)

	~RecHitTask () override

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

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

	DQTask (edm::ParameterSet const &)

	~DQTask () override

Public Member Functions inherited from DQMEDAnalyzer
void	beginRun (edm::Run const &, edm::EventSetup const &) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

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

	~EDAnalyzerBase () override

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 Types
enum	RecoFlag { fUni =0, fTCDS =1, fUnknownIds = 2, nRecoFlag =3 }

Protected Member Functions
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 &)

virtual bool	_isApplicable (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	_cDAAsymmetry_cut_SubdetPM

hcaldqm::ContainerProf2D	_cDAAsymmetryMean_cut_depth

hcaldqm::Container2D	_cDAAsymmetryVsCharge_SubdetPM

hcaldqm::ContainerProf2D	_cEnergy_depth

hcaldqm::Container1D	_cEnergy_Subdet

hcaldqm::ContainerProf1D	_cEnergyvsBX_SubdetPM

hcaldqm::ContainerProf1D	_cEnergyvsieta_Subdet

hcaldqm::ContainerProf1D	_cEnergyvsiphi_SubdetPM

hcaldqm::ContainerProf1D	_cEnergyvsLS_SubdetPM

hcaldqm::Container2D	_cOccupancy_depth

hcaldqm::Container2D	_cOccupancy_ElectronicsuTCA

hcaldqm::Container2D	_cOccupancy_ElectronicsVME

hcaldqm::Container2D	_cOccupancy_FEDuTCA

hcaldqm::Container2D	_cOccupancy_FEDVME

hcaldqm::Container2D	_cOccupancyCut_depth

hcaldqm::Container2D	_cOccupancyCut_ElectronicsuTCA

hcaldqm::Container2D	_cOccupancyCut_ElectronicsVME

hcaldqm::Container2D	_cOccupancyCut_FEDuTCA

hcaldqm::Container2D	_cOccupancyCut_FEDVME

hcaldqm::ContainerProf1D	_cOccupancyCutvsBX_Subdet

hcaldqm::Container1D	_cOccupancyCutvsieta_Subdet

hcaldqm::Container1D	_cOccupancyCutvsiphi_SubdetPM

hcaldqm::Container2D	_cOccupancyCutvsiphivsLS_SubdetPM

hcaldqm::ContainerProf1D	_cOccupancyCutvsLS_Subdet

hcaldqm::Container1D	_cOccupancyvsieta_Subdet

hcaldqm::Container1D	_cOccupancyvsiphi_SubdetPM

hcaldqm::ContainerProf1D	_cOccupancyvsLS_Subdet

hcaldqm::ContainerSingle2D	_cSummaryvsLS

hcaldqm::Container2D	_cSummaryvsLS_FED

hcaldqm::ContainerProf2D	_cTimingCut_depth

hcaldqm::ContainerProf2D	_cTimingCut_ElectronicsuTCA

hcaldqm::ContainerProf2D	_cTimingCut_ElectronicsVME

hcaldqm::ContainerProf2D	_cTimingCut_FEDuTCA

hcaldqm::ContainerProf2D	_cTimingCut_FEDVME

hcaldqm::Container1D	_cTimingCut_HBHEPartition

hcaldqm::Container1D	_cTimingCut_SubdetPM

hcaldqm::ContainerProf1D	_cTimingCutvsBX_SubdetPM

hcaldqm::ContainerProf1D	_cTimingCutvsieta_Subdet

hcaldqm::ContainerProf1D	_cTimingCutvsiphi_SubdetPM

hcaldqm::ContainerProf1D	_cTimingCutvsLS_FED

hcaldqm::ContainerProf1D	_cTimingCutvsLS_SubdetPM

hcaldqm::Container2D	_cTimingvsEnergy_SubdetPM

double	_cutE_HBHE

double	_cutE_HF

double	_cutE_HO

hcaldqm::electronicsmap::ElectronicsMap	_ehashmap

hcaldqm::filter::HashFilter	_filter_FEDsuTCA

hcaldqm::filter::HashFilter	_filter_FEDsVME

hcaldqm::filter::HashFilter	_filter_HF

hcaldqm::filter::HashFilter	_filter_uTCA

hcaldqm::filter::HashFilter	_filter_VME

std::vector< HcalGenericDetId >	_gids

bool	_hfPreRecHitsAvailable

edm::InputTag	_tagHBHE

edm::InputTag	_tagHF

edm::InputTag	_tagHO

edm::InputTag	_tagPreHF

double	_thresh_unihf

edm::EDGetTokenT< HBHERecHitCollection >	_tokHBHE

edm::EDGetTokenT< HFRecHitCollection >	_tokHF

edm::EDGetTokenT< HORecHitCollection >	_tokHO

edm::EDGetTokenT< HFPreRecHitCollection >	_tokPreHF

bool	_unknownIdsPresent

std::vector< hcaldqm::flag::Flag >	_vflags

std::vector< uint32_t >	_vhashFEDs

hcaldqm::ContainerXXX< uint32_t >	_xUni

hcaldqm::ContainerXXX< uint32_t >	_xUniHF

MonitorElement *	meUnknownIds1LS

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

module: RecHitTask.h Author: VK Description: HCAL RECO Data Tier Evaluation

Online: Offline:

Definition at line 26 of file RecHitTask.h.

Member Enumeration Documentation

enum RecHitTask::RecoFlag

protected

Enumerator
fUni
fTCDS
fUnknownIds
nRecoFlag

Definition at line 61 of file RecHitTask.h.

         {
             fUni=0,
             fTCDS=1,
             fUnknownIds = 2,
             nRecoFlag=3
         };

Constructor & Destructor Documentation

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

Definition at line 8 of file RecHitTask.cc.

References _cutE_HBHE, _cutE_HF, _cutE_HO, _hfPreRecHitsAvailable, _tagHBHE, _tagHF, _tagHO, _tagPreHF, _thresh_unihf, _tokHBHE, _tokHF, _tokHO, _tokPreHF, _vflags, nanoDQM_cfi::Flag, fTCDS, fUni, fUnknownIds, edm::ParameterSet::getUntrackedParameter(), and nRecoFlag.

                                                :
     DQTask(ps)
 {
     _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
         edm::InputTag("hbhereco"));
     _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
         edm::InputTag("horeco"));
     _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
         edm::InputTag("hfreco"));
     _tagPreHF = ps.getUntrackedParameter<edm::InputTag>("tagPreHF",
         edm::InputTag(""));
     _hfPreRecHitsAvailable = ps.getUntrackedParameter<bool>("hfPreRecHitsAvailable", false);
 
     _tokHBHE = consumes<HBHERecHitCollection>(_tagHBHE);
     _tokHO = consumes<HORecHitCollection>(_tagHO);
     _tokHF = consumes<HFRecHitCollection>(_tagHF);
     _tokPreHF = consumes<HFPreRecHitCollection>(_tagPreHF);
 
     _cutE_HBHE = ps.getUntrackedParameter<double>("cutE_HBHE", 5);
     _cutE_HO = ps.getUntrackedParameter<double>("cutE_HO", 5);
     _cutE_HF = ps.getUntrackedParameter<double>("cutE_HF", 5);
     _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
 
     //  order must be the same as in RecoFlag enum
     _vflags.resize(nRecoFlag);
     _vflags[fUni]=flag::Flag("UniSlotHF");
     _vflags[fTCDS]=flag::Flag("TCDS");
     _vflags[fUnknownIds] = flag::Flag("UnknownIds");
 }

RecHitTask::~RecHitTask ( )

inlineoverride

Definition at line 30 of file RecHitTask.h.

References _process(), _resetMonitors(), beginLuminosityBlock(), bookHistograms(), and endLuminosityBlock().

30 {}

Member Function Documentation

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

overrideprotectedvirtual

Implements hcaldqm::DQTask.

Definition at line 415 of file RecHitTask.cc.

References _cDAAsymmetry_cut_SubdetPM, _cDAAsymmetryMean_cut_depth, _cDAAsymmetryVsCharge_SubdetPM, _cEnergy_depth, _cEnergy_Subdet, _cEnergyvsBX_SubdetPM, _cEnergyvsieta_Subdet, _cEnergyvsiphi_SubdetPM, _cEnergyvsLS_SubdetPM, _cOccupancy_depth, _cOccupancy_ElectronicsuTCA, _cOccupancy_ElectronicsVME, _cOccupancy_FEDuTCA, _cOccupancy_FEDVME, _cOccupancyCut_depth, _cOccupancyCut_ElectronicsuTCA, _cOccupancyCut_ElectronicsVME, _cOccupancyCut_FEDuTCA, _cOccupancyCut_FEDVME, _cOccupancyCutvsBX_Subdet, _cOccupancyCutvsieta_Subdet, _cOccupancyCutvsiphi_SubdetPM, _cOccupancyCutvsiphivsLS_SubdetPM, _cOccupancyCutvsLS_Subdet, _cOccupancyvsieta_Subdet, _cOccupancyvsiphi_SubdetPM, _cOccupancyvsLS_Subdet, _cTimingCut_depth, _cTimingCut_ElectronicsuTCA, _cTimingCut_ElectronicsVME, _cTimingCut_FEDuTCA, _cTimingCut_FEDVME, _cTimingCut_HBHEPartition, _cTimingCut_SubdetPM, _cTimingCutvsBX_SubdetPM, _cTimingCutvsieta_Subdet, _cTimingCutvsiphi_SubdetPM, _cTimingCutvsLS_FED, _cTimingCutvsLS_SubdetPM, _cTimingvsEnergy_SubdetPM, hcaldqm::DQModule::_currentLS, _cutE_HBHE, _cutE_HF, _cutE_HO, _ehashmap, _filter_HF, _hfPreRecHitsAvailable, hcaldqm::DQModule::_logger, hcaldqm::DQModule::_ptype, _tagHBHE, _tagHF, _tagHO, _tagPreHF, _tokHBHE, _tokHF, _tokHO, _tokPreHF, _unknownIdsPresent, _xUniHF, edm::SortedCollection< T, SORT >::begin(), edm::EventBase::bunchCrossing(), Logger::dqmthrow(), runTauDisplay::eid, edm::SortedCollection< T, SORT >::end(), hcaldqm::Container2D::fill(), hcaldqm::Container1D::fill(), MonitorElement::Fill(), hcaldqm::filter::HashFilter::filter(), hcaldqm::fOffline, hcaldqm::fOnline, hcaldqm::ContainerXXX< STDTYPE >::get(), edm::Event::getByToken(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HcalDetId::ieta(), edm::InputTag::instance(), HcalElectronicsId::isUTCAid(), HcalElectronicsId::isVMEid(), edm::InputTag::label(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), meUnknownIds1LS, DetId::rawId(), and HcalDetId::subdet().

Referenced by ~RecHitTask().

 {
     edm::Handle<HBHERecHitCollection> chbhe;
     edm::Handle<HORecHitCollection> cho;
     edm::Handle<HFRecHitCollection> chf;
 
     if (!(e.getByToken(_tokHBHE, chbhe)))
         _logger.dqmthrow("Collection HBHERecHitCollection not available "
             + _tagHBHE.label() + " " + _tagHBHE.instance());
     if (!(e.getByToken(_tokHO, cho)))
         _logger.dqmthrow("Collection HORecHitCollection not available "
             + _tagHO.label() + " " + _tagHO.instance());
     if (!(e.getByToken(_tokHF, chf)))
         _logger.dqmthrow("Collection HFRecHitCollection not available "
             + _tagHF.label() + " " + _tagHF.instance());
 
     edm::Handle<HFPreRecHitCollection> cprehf;
     if (_hfPreRecHitsAvailable) {
         if (!(e.getByToken(_tokPreHF, cprehf)))
             _logger.dqmthrow("Collection HFPreRecHitCollection not available "
                 + _tagPreHF.label() + " " + _tagPreHF.instance());
     }
 
     //  extract some info per event
     int bx = e.bunchCrossing();
 
     //  To fill histograms outside of the loop, you need to determine if there were
     //  any valid det ids first
     uint32_t rawidValid = 0;
     uint32_t rawidHBValid = 0;
     uint32_t rawidHEValid = 0;
 
     double ehbm = 0; double ehbp = 0;
     double ehem = 0; double ehep = 0;
     int nChsHB = 0; int nChsHE = 0;
     int nChsHBCut = 0; int nChsHECut = 0;
     for (HBHERecHitCollection::const_iterator it=chbhe->begin();
         it!=chbhe->end(); ++it)
     {
         double energy = it->energy();
         double timing = it->time();
 
         //  Explicit check on the DetIds present in the Collection
         HcalDetId did = it->id();
         uint32_t rawid = _ehashmap.lookup(did);
         /*
          * Needs to be removed as DetIds that belong to the HEP17 after combination
          * are not present in the emap
          * Removed until further notice!
          *
          */
         //if (rawid==0)
         //{meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
         
         HcalElectronicsId const& eid(rawid);
         rawidValid = did.rawId();
         if (did.subdet()==HcalBarrel)
             rawidHBValid = did.rawId();
         else if (did.subdet()==HcalEndcap)
             rawidHEValid = did.rawId();
 
         _cEnergy_Subdet.fill(did, energy);
         _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
         _cOccupancy_depth.fill(did);
         did.subdet()==HcalBarrel?did.ieta()>0?ehbp+=energy:ehbm+=energy:
             did.ieta()>0?ehep+=energy:ehem+=energy;
 
         //  ONLINE ONLY!
         if (_ptype==fOnline)
         {
             _cOccupancyvsiphi_SubdetPM.fill(did);
             _cOccupancyvsieta_Subdet.fill(did);
         }
         //  ^^^ONLINE ONLY!
         //  
         if (_ptype != fOffline) { // hidefed2crate
             //  Also, for these electronics plots, require that the channel was found in the emap.
             if (rawid != 0) {
                 if (eid.isVMEid())
                 {
                     _cOccupancy_FEDVME.fill(eid);
                     _cOccupancy_ElectronicsVME.fill(eid);
                 } else if (eid.isUTCAid()) {
                     _cOccupancy_FEDuTCA.fill(eid);
                     _cOccupancy_ElectronicsuTCA.fill(eid);
                 }
             }
         }
 
         if (energy>_cutE_HBHE)
         {
             //  ONLINE ONLY!
             if (_ptype==fOnline)
             {
                 _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
                 _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
                 _cEnergyvsieta_Subdet.fill(did, energy);
                 _cEnergyvsiphi_SubdetPM.fill(did, energy);
                 _cTimingCutvsieta_Subdet.fill(did, timing);
                 _cTimingCutvsiphi_SubdetPM.fill(did, timing);
                 _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
                 _cOccupancyCutvsiphi_SubdetPM.fill(did);
                 _cOccupancyCutvsieta_Subdet.fill(did);
                 _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
             }
             //  ^^^ONLINE ONLY!
             _cEnergy_depth.fill(did, energy);
             _cTimingCut_SubdetPM.fill(did, timing);
             _cTimingCut_HBHEPartition.fill(did, timing);
             _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
 
             //  ONLINE 
             if (_ptype==fOnline) {
                 if (rawid != 0) {
                     _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
                 }
                 _cTimingCut_depth.fill(did, timing);
                 //  ^^^ONLINE
             } else {
                 if (_ptype != fOffline) { // hidefed2crate
                     if (rawid != 0) {
                         _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
                     }
                 }
                 _cTimingCut_depth.fill(did, timing);
             }
             _cOccupancyCut_depth.fill(did);
             if (_ptype != fOffline) { // hidefed2crate
                 if (rawid != 0) {
                     if (eid.isVMEid()){ 
                         //  ONLINE 
                         if (_ptype==fOnline)
                         {
                             _cTimingCut_FEDVME.fill(eid, timing);
                             _cTimingCut_ElectronicsVME.fill(eid, timing);
                         } // ^^^ ONLINE
                         else
                         {
                             _cTimingCut_FEDVME.fill(eid, timing);
                             _cTimingCut_ElectronicsVME.fill(eid, timing);
                         }
                         //  ^^^ONLINE
 
                         _cOccupancyCut_FEDVME.fill(eid);
                         _cOccupancyCut_ElectronicsVME.fill(eid);
                     } else if (eid.isUTCAid()) {
                         if (_ptype==fOnline)
                         {
                             //  time constraints are explicit!
                             _cTimingCut_FEDuTCA.fill(eid, timing);
                             _cTimingCut_ElectronicsuTCA.fill(eid, timing);
                         }
                         else
                         {
                             _cTimingCut_FEDuTCA.fill(eid, timing);
                             _cTimingCut_ElectronicsuTCA.fill(eid, timing);
                         }
                         _cOccupancyCut_FEDuTCA.fill(eid);
                         _cOccupancyCut_ElectronicsuTCA.fill(eid);
                     }
                 }
             }
             did.subdet()==HcalBarrel?nChsHBCut++:nChsHECut++;
         }
         did.subdet()==HcalBarrel?nChsHB++:nChsHE++;
     }
 
     if (rawidHBValid!=0 && rawidHEValid!=0)
     {
         _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, 
             nChsHB);
         _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS,
             nChsHE);
     
         //  ONLINE ONLY!
         if (_ptype==fOnline)
         {
             _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid),
                 bx, nChsHBCut);
             _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid),
                 bx, nChsHECut);
             _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid), 
                 _currentLS, nChsHBCut);
             _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid), 
                 _currentLS, nChsHECut);
         }
         //  ^^^ONLINE ONLY!
     }
 
     //  reset
     rawidValid = 0;
 
     int nChsHO = 0; int nChsHOCut = 0;
     double ehop = 0; double ehom = 0;
     for (HORecHitCollection::const_iterator it=cho->begin();
         it!=cho->end(); ++it)
     {
         double energy = it->energy();
         double timing = it->time();
         
         //  Explicit check on the DetIds present in the Collection
         HcalDetId did = it->id();
         uint32_t rawid = _ehashmap.lookup(did);
         if (rawid==0)
         {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
         HcalElectronicsId const& eid(rawid);
         if (did.subdet()==HcalOuter)
             rawidValid = did.rawId();
 
         _cEnergy_Subdet.fill(did, energy);
         _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
         _cOccupancy_depth.fill(did);
         did.ieta()>0?ehop+=energy:ehom+=energy;
 
         //  IMPORTANT: ONLINE ONLY!
         if (_ptype==fOnline)
         {
             _cOccupancyvsiphi_SubdetPM.fill(did);
             _cOccupancyvsieta_Subdet.fill(did);
         }
         //  ONLINE ONLY!
 
         if (_ptype != fOffline) { // hidefed2crate
             if (eid.isVMEid())
             {
                 _cOccupancy_FEDVME.fill(eid);
                 _cOccupancy_ElectronicsVME.fill(eid);
             }
             else
             {
                 _cOccupancy_FEDuTCA.fill(eid);
                 _cOccupancy_ElectronicsuTCA.fill(eid);
             }
         }
 
         if (energy>_cutE_HO)
         {
             //  ONLINE ONLY!
             if (_ptype==fOnline)
             {
                 _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
                 _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
                 _cEnergyvsieta_Subdet.fill(did, energy);
                 _cEnergyvsiphi_SubdetPM.fill(did, energy);
                 _cTimingCutvsieta_Subdet.fill(did, timing);
                 _cTimingCutvsiphi_SubdetPM.fill(did, timing);
                 _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
                 _cOccupancyCutvsiphi_SubdetPM.fill(did);
                 _cOccupancyCutvsieta_Subdet.fill(did);
                 _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
             }
             //  ^^^ONLINE ONLY!
             
             _cEnergy_depth.fill(did, energy);
             _cTimingCut_SubdetPM.fill(did, timing);
             _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
             if (_ptype != fOffline) { // hidefed2crate
                 _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
             }
             _cOccupancyCut_depth.fill(did);
             _cTimingCut_depth.fill(did, timing);
             if (_ptype != fOffline) { // hidefed2crate
                 if (eid.isVMEid())
                 {
                     _cTimingCut_FEDVME.fill(eid, timing);
                     _cTimingCut_ElectronicsVME.fill(eid, timing);
                     _cOccupancyCut_FEDVME.fill(eid);
                     _cOccupancyCut_ElectronicsVME.fill(eid);
                 }
                 else
                 {
                     _cTimingCut_FEDuTCA.fill(eid, timing);
                     _cTimingCut_ElectronicsuTCA.fill(eid, timing);
                     _cOccupancyCut_FEDuTCA.fill(eid);
                     _cOccupancyCut_ElectronicsuTCA.fill(eid);
                 }
             }
             nChsHOCut++;
         }
         nChsHO++;
     }
 
     if (rawidValid!=0)
     {
         _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
             nChsHO);
         //  ONLINE ONLY!
         if (_ptype==fOnline)
         {
             _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid),
                 bx, nChsHOCut);
             _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), 
                 _currentLS, nChsHOCut);
         }
         //  ^^^ONLINE ONLY!
     }
 
     //reset 
     rawidValid = 0;
 
     int nChsHF = 0; int nChsHFCut = 0;
     double ehfp = 0; double ehfm = 0;
     for (HFRecHitCollection::const_iterator it=chf->begin();
         it!=chf->end(); ++it)
     {
         double energy = it->energy();
         double timing = it->time();
 
         //  Explicit check on the DetIds present in the Collection
         HcalDetId did = it->id();
         uint32_t rawid = _ehashmap.lookup(did);
         if (rawid==0)
         {meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}
         HcalElectronicsId const& eid(rawid);
         if (did.subdet()==HcalForward)
             rawidValid = did.rawId();
 
         _cEnergy_Subdet.fill(did, energy);
         _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
         _cOccupancy_depth.fill(did);
         did.ieta()>0?ehfp+=energy:ehfm+=energy;
         
         //  IMPORTANT:
         //  only for Online Processing
         //
         if (_ptype==fOnline)
         {
             _cOccupancyvsiphi_SubdetPM.fill(did);
             _cOccupancyvsieta_Subdet.fill(did);
         }
         //  ONLINE ONLY!
 
         if (_ptype != fOffline) { // hidefed2crate
             if (eid.isVMEid())
             {
                 _cOccupancy_FEDVME.fill(eid);
                 _cOccupancy_ElectronicsVME.fill(eid);
             }
             else
             {
                 _cOccupancy_FEDuTCA.fill(eid);
                 _cOccupancy_ElectronicsuTCA.fill(eid);
             }
         }
 
         if (energy>_cutE_HF)
         {
             //  ONLINE ONLY!
             if (_ptype==fOnline)
             {
                 _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
                 _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
                 _cEnergyvsieta_Subdet.fill(did, energy);
                 _cEnergyvsiphi_SubdetPM.fill(did, energy);
                 _cTimingCutvsieta_Subdet.fill(did, timing);
                 _cTimingCutvsiphi_SubdetPM.fill(did, timing);
                 _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
                 _cOccupancyCutvsiphi_SubdetPM.fill(did);
                 _cOccupancyCutvsieta_Subdet.fill(did);
                 _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
                 if (_ptype != fOffline) { // hidefed2crate
                     _xUniHF.get(eid)++;
                 }
             }
             //  ^^^ONLINE ONLY!
             _cEnergy_depth.fill(did, energy);
             _cTimingCut_SubdetPM.fill(did, timing);
             _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
             if (_ptype != fOffline) { // hidefed2crate
                 _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
             }
             _cOccupancyCut_depth.fill(did);
             _cTimingCut_depth.fill(did, timing);
             if (_ptype != fOffline) { // hidefed2crate
                 if (eid.isVMEid())
                 {
                     _cTimingCut_FEDVME.fill(eid, timing);
                     _cTimingCut_ElectronicsVME.fill(eid, timing);
                     _cOccupancyCut_FEDVME.fill(eid);
                     _cOccupancyCut_ElectronicsVME.fill(eid);
                 }
                 else
                 {
                     _cTimingCut_FEDuTCA.fill(eid, timing);
                     _cTimingCut_ElectronicsuTCA.fill(eid, timing);
                     _cOccupancyCut_FEDuTCA.fill(eid);
                     _cOccupancyCut_ElectronicsuTCA.fill(eid);
                 }
             }
             nChsHFCut++;
         }
         nChsHF++;
     }
 
     if (rawidValid!=0)
     {
         _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS,
             nChsHF);
         //  ONLINE ONLY!
         if (_ptype==fOnline)
         {
             _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid),
                 bx, nChsHFCut);
             _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), 
                 _currentLS, nChsHFCut);
         }
         //  ^^^ONLINE ONLY!
     }
 
     // Loop over HFPreRecHits to get charge and charge asymmetry
     if (_hfPreRecHitsAvailable) {
         for (HFPreRecHitCollection::const_iterator it=cprehf->begin();
             it!=cprehf->end(); ++it)
         {
             HcalDetId did = it->id();
             if (_filter_HF.filter(did)) {
                 continue;
             }
             std::pair<float, bool> chargeAsymmetry = it->chargeAsymmetry(0.);
             std::pair<float, bool> chargeAsymmetryCut = it->chargeAsymmetry(20.);
 
             if (chargeAsymmetry.second) {
                 _cDAAsymmetryVsCharge_SubdetPM.fill(did, chargeAsymmetry.first, it->charge());
             }
             if (chargeAsymmetryCut.second) {
                 _cDAAsymmetryMean_cut_depth.fill(did, chargeAsymmetryCut.first);
                 _cDAAsymmetry_cut_SubdetPM.fill(did, chargeAsymmetryCut.first);
             }
         }
     }
 }

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

overrideprotectedvirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 401 of file RecHitTask.cc.

References _unknownIdsPresent, and hcaldqm::f1LS.

Referenced by ~RecHitTask().

 {
     switch(uf)
     {
         case hcaldqm::f1LS:
             _unknownIdsPresent = false;
             break;
         default:
             break;
     }
 
     DQTask::_resetMonitors(uf);
 }

void RecHitTask::beginLuminosityBlock	(	edm::LuminosityBlock const &	lb,
		edm::EventSetup const &	es
	)

overridevirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 848 of file RecHitTask.cc.

Referenced by ~RecHitTask().

 {
     DQTask::beginLuminosityBlock(lb, es);
 }

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

overridevirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 38 of file RecHitTask.cc.

Referenced by ~RecHitTask().

 {
     DQTask::bookHistograms(ib,r,es);
 
     //  GET WHAT YOU NEED
     edm::ESHandle<HcalDbService> dbs;
     es.get<HcalDbRecord>().get(dbs);
     _emap = dbs->getHcalMapping();
 
     std::vector<uint32_t> vVME;
     std::vector<uint32_t> vuTCA;
     vVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN, 
         constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
     vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
         FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
     _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
         vVME);
     _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics,
         vuTCA);
     std::vector<uint32_t> vhashHF; 
     vhashHF.push_back(hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29,1,1)));
     _filter_HF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet, vhashHF);
 
     //  INITIALIZE FIRST
     //  Energy
     _cEnergy_Subdet.initialize(_name, "Energy", hcaldqm::hashfunctions::fSubdet,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
     _cEnergy_depth.initialize(_name, "Energy", hcaldqm::hashfunctions::fdepth,
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),0);
 
     //  Timing
     _cTimingCut_SubdetPM.initialize(_name, "TimingCut", 
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
     _cTimingvsEnergy_SubdetPM.initialize(_name, "TimingvsEnergy",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
 
     _cTimingCut_HBHEPartition.initialize(_name, "TimingCut",
         hcaldqm::hashfunctions::fHBHEPartition,
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
     _cTimingCut_depth.initialize(_name, "TimingCut",
         hcaldqm::hashfunctions::fdepth,
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
         new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
     _cTimingCutvsLS_SubdetPM.initialize(_name, "TimingCutvsLS",
         hcaldqm::hashfunctions::fSubdetPM,
         new hcaldqm::quantity::LumiSection(_maxLS),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
 
     //  Occupancy
     _cOccupancy_depth.initialize(_name, "Occupancy",
         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);
 
     _cOccupancyvsLS_Subdet.initialize(_name, "OccupancyvsLS",
         hcaldqm::hashfunctions::fSubdet,
         new hcaldqm::quantity::LumiSection(_maxLS),
         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
 
     _cOccupancyCut_depth.initialize(_name, "OccupancyCut",
         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 (_hfPreRecHitsAvailable) {
         _cDAAsymmetryVsCharge_SubdetPM.initialize(_name, "ChargeVsAsymmetry", 
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry), 
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cDAAsymmetryMean_cut_depth.initialize(_name, "AsymmetryMean",
             hcaldqm::hashfunctions::fdepth,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),0);
         _cDAAsymmetry_cut_SubdetPM.initialize(_name, "Asymmetry",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN), 0);
     }
 
     //  INITIALIZE HISTOGRAMS to be used only in Online
     if (_ptype==fOnline)
     {
         _cEnergyvsieta_Subdet.initialize(_name, "Energyvsieta",
             hcaldqm::hashfunctions::fSubdet,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
         _cEnergyvsiphi_SubdetPM.initialize(_name, "Energyvsiphi",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
         _cEnergyvsLS_SubdetPM.initialize(_name, "EnergyvsLS",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::LumiSection(_maxLS),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
         _cEnergyvsBX_SubdetPM.initialize(_name, "EnergyvsBX",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),0);
         _cTimingCutvsieta_Subdet.initialize(_name, "TimingCutvsieta",
             hcaldqm::hashfunctions::fSubdet,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cTimingCutvsiphi_SubdetPM.initialize(_name, "TimingCutvsiphi",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cTimingCutvsBX_SubdetPM.initialize(_name, "TimingCutvsBX",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cOccupancyvsiphi_SubdetPM.initialize(_name, "Occupancyvsiphi",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyvsieta_Subdet.initialize(_name, "Occupancyvsieta",
             hcaldqm::hashfunctions::fSubdet,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyCutvsiphi_SubdetPM.initialize(_name, "OccupancyCutvsiphi",
             hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyCutvsieta_Subdet.initialize(_name, "OccupancyCutvsieta",
             hcaldqm::hashfunctions::fSubdet,
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyCutvsBX_Subdet.initialize(_name, "OccupancyCutvsBX",
             hcaldqm::hashfunctions::fSubdet,
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
             "OccupancyCutvsiphivsLS", hcaldqm::hashfunctions::fSubdetPM,
             new hcaldqm::quantity::LumiSection(_maxLS),
             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyCutvsLS_Subdet.initialize(_name, "OccupancyCutvsLS",
             hcaldqm::hashfunctions::fSubdet,
             new hcaldqm::quantity::LumiSection(_maxLS),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),0);
     }
 
     // FED-based plots
     if (_ptype != fOffline) { // hidefed2crate
         std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
         std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
         std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);
     
         //  push the rawIds of each fed into the vector
         for (std::vector<int>::const_iterator it=vFEDsVME.begin();
             it!=vFEDsVME.end(); ++it)
             _vhashFEDs.push_back(HcalElectronicsId(
                 FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it)-FED_VME_MIN).rawId());
         for (std::vector<int>::const_iterator it=vFEDsuTCA.begin(); 
             it!=vFEDsuTCA.end(); ++it)
         {
             std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
             _vhashFEDs.push_back(HcalElectronicsId(
                 cspair.first, cspair.second, FIBER_uTCA_MIN1,
                 FIBERCH_MIN, false).rawId());
         }
     
         _cTimingCut_FEDVME.initialize(_name, "TimingCut",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberVMEFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cTimingCut_FEDuTCA.initialize(_name, "TimingCut",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cTimingCutvsLS_FED.initialize(_name, "TimingCutvsLS",
             hcaldqm::hashfunctions::fFED,
             new hcaldqm::quantity::LumiSection(_maxLS),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cTimingCut_ElectronicsVME.initialize(_name, "TimingCut",
             hcaldqm::hashfunctions::fElectronics,
             new hcaldqm::quantity::FEDQuantity(vFEDsVME),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
         _cTimingCut_ElectronicsuTCA.initialize(_name, "TimingCut",
             hcaldqm::hashfunctions::fElectronics,
             new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),0);
     
         _cOccupancy_FEDVME.initialize(_name, "Occupancy",
             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);
         _cOccupancy_FEDuTCA.initialize(_name, "Occupancy",
             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);
         _cOccupancy_ElectronicsVME.initialize(_name, "Occupancy",
             hcaldqm::hashfunctions::fElectronics,
             new hcaldqm::quantity::FEDQuantity(vFEDsVME),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancy_ElectronicsuTCA.initialize(_name, "Occupancy",
             hcaldqm::hashfunctions::fElectronics,
             new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
     
         _cOccupancyCut_FEDVME.initialize(_name, "OccupancyCut",
             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);
         _cOccupancyCut_FEDuTCA.initialize(_name, "OccupancyCut",
             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);
         _cOccupancyCut_ElectronicsVME.initialize(_name, "OccupancyCut",
             hcaldqm::hashfunctions::fElectronics,
             new hcaldqm::quantity::FEDQuantity(vFEDsVME),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSpigot),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         _cOccupancyCut_ElectronicsuTCA.initialize(_name, "OccupancyCut",
             hcaldqm::hashfunctions::fElectronics,
             new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
             new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
         if (_ptype == fOnline) {
             _cSummaryvsLS_FED.initialize(_name, "SummaryvsLS",
                 hcaldqm::hashfunctions::fFED,
                 new hcaldqm::quantity::LumiSection(_maxLS),
                 new hcaldqm::quantity::FlagQuantity(_vflags),
                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
             _cSummaryvsLS.initialize(_name, "SummaryvsLS",
                 new hcaldqm::quantity::LumiSection(_maxLS),
                 new hcaldqm::quantity::FEDQuantity(vFEDs),
                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),0);
             
             _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
             _xUni.initialize(hcaldqm::hashfunctions::fFED);
         }
     }
 
     //  BOOK HISTOGRAMS
     char cutstr[200];
     sprintf(cutstr, "_EHBHE%dHO%dHF%d", int(_cutE_HBHE),
         int(_cutE_HO), int(_cutE_HF));
     char cutstr2[200];
     sprintf(cutstr2, "_EHF%d", int(_cutE_HF));
 
     //  Energy
     _cEnergy_Subdet.book(ib, _emap, _subsystem);
     _cEnergy_depth.book(ib, _emap, _subsystem);
 
     //  Timing
     _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
     _cTimingvsEnergy_SubdetPM.book(ib, _emap, _subsystem);
     if (_ptype != fOffline) { // hidefed2crate
         _cTimingCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cTimingCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
     }
     _cTimingCut_HBHEPartition.book(ib, _emap, _subsystem);
     _cTimingCut_depth.book(ib, _emap, _subsystem);
     _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
     if (_ptype != fOffline) { // hidefed2crate
         _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
     }
 
     //  Occupancy
     _cOccupancy_depth.book(ib, _emap, _subsystem);
     if (_ptype != fOffline) { // hidefed2crate
         _cOccupancy_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cOccupancy_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
     }
     _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
     _cOccupancyCut_depth.book(ib, _emap, _subsystem);
     if (_ptype != fOffline) { // hidefed2crate
         _cOccupancyCut_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
         _cOccupancyCut_ElectronicsVME.book(ib, _emap, _filter_uTCA, _subsystem);
         _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
     }
 
     if (_hfPreRecHitsAvailable) {
         _cDAAsymmetryVsCharge_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
         _cDAAsymmetryMean_cut_depth.book(ib, _emap, _filter_HF, _subsystem);
         _cDAAsymmetry_cut_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
     }
 
     //  BOOK HISTOGRAMS to be used only in Online
     if (_ptype==fOnline)
     {
         _cEnergyvsLS_SubdetPM.book(ib, _emap, _subsystem);
         _cEnergyvsieta_Subdet.book(ib, _emap, _subsystem);
         _cEnergyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
         _cEnergyvsBX_SubdetPM.book(ib, _emap, _subsystem);
         _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
         _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
         _cTimingCutvsBX_SubdetPM.book(ib, _emap, _subsystem);
         _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
         _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
         _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
         _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
         _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
         _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
         _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
         _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
         _cSummaryvsLS.book(ib, _subsystem);
 
         std::vector<uint32_t> vhashFEDHF;
         vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN,
             FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
         vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN,
             FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
         vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN,
             FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
         vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN+6,
             FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
         vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN+6,
             FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
         vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN+6,
             FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
         HashFilter filter_FEDHF;
         filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED,
             vhashFEDHF);
 
         _gids = _emap->allPrecisionId();
         if (_ptype != fOffline) { // hidefed2crate
             _xUniHF.book(_emap, filter_FEDHF);
             _xUni.book(_emap);
         }
     }
 
     //  initialize hash map
     _ehashmap.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap);
 
     //  book some mes...
     ib.setCurrentFolder(_subsystem+"/"+_name);
     meUnknownIds1LS = ib.book1D("UnknownIds", "UnknownIds",
         1, 0, 1);
     _unknownIdsPresent = false;
     meUnknownIds1LS->setLumiFlag();
 }

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

overridevirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 854 of file RecHitTask.cc.

References _cSummaryvsLS, _cSummaryvsLS_FED, hcaldqm::DQModule::_currentLS, hcaldqm::DQModule::_ptype, hcaldqm::DQModule::_runkeyVal, hcaldqm::flag::Flag::_state, _thresh_unihf, _unknownIdsPresent, hcaldqm::DQTask::_vcdaqEids, _vflags, _vhashFEDs, _xUni, _xUniHF, hcaldqm::ContainerXXX< STDTYPE >::begin(), DEFINE_FWK_MODULE, runTauDisplay::eid, hcaldqm::ContainerXXX< STDTYPE >::end(), hcaldqm::flag::fBAD, hcaldqm::flag::fGOOD, spr::find(), hcaldqm::flag::fNCDAQ, hcaldqm::fOffline, hcaldqm::fOnline, hcaldqm::flag::fPROBLEMATIC, fUni, fUnknownIds, hcaldqm::ContainerXXX< STDTYPE >::get(), hcaldqm::utilities::isFEDHF(), hcaldqm::ContainerXXX< STDTYPE >::reset(), hcaldqm::flag::Flag::reset(), hcaldqm::Container2D::setBinContent(), hcaldqm::ContainerSingle2D::setBinContent(), globals_cff::x1, and globals_cff::x2.

Referenced by ~RecHitTask().

 {
     if (_ptype!=fOnline)
         return;
 
     //
     //  GENERATE STATUS ONLY FOR ONLINE
     //
 //  for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
 //      it!=gids.end(); ++it)
 //  {}
 
     for (uintCompactMap::const_iterator it=_xUniHF.begin();
         it!=_xUniHF.end(); ++it)
     {
         uint32_t hash1 = it->first;
         HcalElectronicsId eid1(hash1);
         double x1 = it->second;
 
         for (uintCompactMap::const_iterator jt=_xUniHF.begin();
             jt!=_xUniHF.end(); ++jt)
         {
             if (jt==it)
                 continue;
             double x2 = jt->second;
             if (x2==0)
                 continue;
             if (x1/x2<_thresh_unihf)
                 _xUni.get(eid1)++;
         }
     }
 
     if (_ptype != fOffline) { // hidefed2crate
         for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
         it!=_vhashFEDs.end(); ++it)
         {
             flag::Flag fSum("RECO");
             HcalElectronicsId eid = HcalElectronicsId(*it);
 
             std::vector<uint32_t>::const_iterator cit=std::find(
                 _vcdaqEids.begin(), _vcdaqEids.end(), *it);
             if (cit==_vcdaqEids.end())
             {
                 //  not @cDAQ
                 for (uint32_t iflag=0; iflag<_vflags.size(); iflag++)
                     _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
                         int(flag::fNCDAQ));
                 _cSummaryvsLS.setBinContent(eid, _currentLS, int(flag::fNCDAQ));
                 continue;
             }
 
             //  FED is @cDAQ
             if (hcaldqm::utilities::isFEDHF(eid) && (_runkeyVal==0 || _runkeyVal==4))
             {
                 if (_xUni.get(eid)>0)
                     _vflags[fUni]._state = flag::fPROBLEMATIC;
                 else
                     _vflags[fUni]._state = flag::fGOOD;
             }
 
             if (_unknownIdsPresent)
                 _vflags[fUnknownIds]._state = flag::fBAD;
             else 
                 _vflags[fUnknownIds]._state = flag::fGOOD;
 
             int iflag=0;
             for (std::vector<flag::Flag>::iterator ft=_vflags.begin();
                 ft!=_vflags.end(); ++ft)
             {
                 _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag),
                     int(ft->_state));
                 fSum+=(*ft);
                 iflag++;
 
                 //  reset after using
                 ft->reset();
             }
             _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
         }
         _xUniHF.reset(); _xUni.reset();
     }
 
     //  in the end always do the DQTask::endLumi
     DQTask::endLuminosityBlock(lb, es);
 }