#include <ECALpedestalPCLHarvester.h>

Inheritance diagram for ECALpedestalPCLHarvester:

Public Member Functions
	ECALpedestalPCLHarvester (const edm::ParameterSet &ps)

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

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

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

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

	DQMEDHarvester (void)

virtual void	dqmEndLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endJob () final

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

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::WatchLuminosityBlocks, edm::one::SharedResources >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	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 &	itemsToGetFromEvent () 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)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
bool	checkStatusCode (const DetId &id)

bool	checkVariation (const EcalPedestalsMap &oldPedestals, const EcalPedestalsMap &newPedestals)

virtual void	dqmEndJob (DQMStore::IBooker &ibooker_, DQMStore::IGetter &igetter_)

Private Attributes
const EcalChannelStatus *	channelStatus_

bool	checkAnomalies_

std::vector< int >	chStatusToExclude_

const EcalPedestals *	currentPedestals_

std::string	dqmDir_

int	minEntries_

double	nSigma_

double	thresholdAnomalies_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

Detailed Description

Description: Fill DQM histograms with pedestals. Intended to be used on laser data from the TestEnablesEcalHcal dataset

Definition at line 33 of file ECALpedestalPCLHarvester.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file ECALpedestalPCLHarvester.cc.

References checkAnomalies_, chStatusToExclude_, dqmDir_, edm::ParameterSet::getParameter(), minEntries_, nSigma_, AlCaHLTBitMon_QueryRunRegistry::string, and thresholdAnomalies_.

                                                                            :
     currentPedestals_(0),channelStatus_(0){
 
     chStatusToExclude_= StringToEnumValue<EcalChannelStatusCode::Code>(ps.getParameter<std::vector<std::string> >("ChannelStatusToExclude"));
     minEntries_=ps.getParameter<int>("MinEntries");
     checkAnomalies_  = ps.getParameter<bool>("checkAnomalies"); 
     nSigma_  = ps.getParameter<double>("nSigma");
     thresholdAnomalies_ = ps.getParameter<double>("thresholdAnomalies"); 
     dqmDir_         = ps.getParameter<std::string>("dqmDir");
 }

Member Function Documentation

bool ECALpedestalPCLHarvester::checkStatusCode ( const DetId & id )

private

Definition at line 143 of file ECALpedestalPCLHarvester.cc.

References channelStatus_, chStatusToExclude_, spr::find(), EcalCondObjectContainer< T >::find(), and EcalCondObjectContainer< T >::getMap().

Referenced by dqmEndJob().

                                                              {
 
     EcalChannelStatusMap::const_iterator dbstatusPtr;
     dbstatusPtr = channelStatus_->getMap().find(id.rawId());
     EcalChannelStatusCode::Code  dbstatus = dbstatusPtr->getStatusCode();
 
     std::vector<int>::const_iterator res =
         std::find( chStatusToExclude_.begin(), chStatusToExclude_.end(), dbstatus );
     if ( res != chStatusToExclude_.end() ) return false;
 
     return true;
 }

bool ECALpedestalPCLHarvester::checkVariation	(	const EcalPedestalsMap &	oldPedestals,
		const EcalPedestalsMap &	newPedestals
	)

private

Definition at line 157 of file ECALpedestalPCLHarvester.cc.

References funct::abs(), EBDetId::detIdFromDenseIndex(), EEDetId::detIdFromDenseIndex(), EcalCondObjectContainer< T >::find(), mps_fire::i, EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, EcalPedestal::mean_x12, nSigma_, EcalPedestal::rms_x12, and thresholdAnomalies_.

Referenced by dqmEndJob().

                                                                                     {
 
     uint32_t nAnomaliesEB =0;
     uint32_t nAnomaliesEE =0;
 
     for (uint16_t i =0; i< EBDetId::kSizeForDenseIndexing; ++i) {
        
         DetId id = EBDetId::detIdFromDenseIndex(i);
         const EcalPedestal& newped=* newPedestals.find(id.rawId());   
         const EcalPedestal& oldped=* oldPedestals.find(id.rawId());
 
         if  (std::abs(newped.mean_x12 -oldped.mean_x12)  >  nSigma_ * oldped.rms_x12) nAnomaliesEB++;
 
     }
 
     for (uint16_t i =0; i< EEDetId::kSizeForDenseIndexing; ++i) {
        
         DetId id = EEDetId::detIdFromDenseIndex(i);
         const EcalPedestal& newped=* newPedestals.find(id.rawId());   
         const EcalPedestal& oldped=* oldPedestals.find(id.rawId());
 
         if  (std::abs(newped.mean_x12 -oldped.mean_x12)  >  nSigma_ * oldped.rms_x12) nAnomaliesEE++;
 
     }
     
     if (nAnomaliesEB > thresholdAnomalies_ *  EBDetId::kSizeForDenseIndexing || 
         nAnomaliesEE > thresholdAnomalies_ *  EEDetId::kSizeForDenseIndexing)  
         return true;
 
     return false;
 
 
 }

void ECALpedestalPCLHarvester::dqmEndJob	(	DQMStore::IBooker &	ibooker_,
		DQMStore::IGetter &	igetter_
	)

privatevirtual

Implements DQMEDHarvester.

Definition at line 31 of file ECALpedestalPCLHarvester.cc.

References checkAnomalies_, checkStatusCode(), checkVariation(), currentPedestals_, cond::service::PoolDBOutputService::currentTime(), EBDetId::detIdFromDenseIndex(), EEDetId::detIdFromDenseIndex(), dqmDir_, EcalCondObjectContainer< T >::find(), DQMStore::IGetter::get(), MonitorElement::getEntries(), MonitorElement::getMean(), MonitorElement::getRMS(), mps_fire::i, edm::Service< T >::isAvailable(), EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, RecoTauDiscriminantConfiguration::mean, EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, minEntries_, RecoTauDiscriminantConfiguration::rms, EcalPedestal::rms_x1, EcalPedestal::rms_x12, EcalPedestal::rms_x6, EcalCondObjectContainer< T >::setValue(), AlCaHLTBitMon_QueryRunRegistry::string, and cond::service::PoolDBOutputService::writeOne().

                                                                                              {
 
 
     // calculate pedestals and fill db record
     EcalPedestals pedestals;
 
 
     for (uint16_t i =0; i< EBDetId::kSizeForDenseIndexing; ++i) {
         std::string hname = dqmDir_+"/eb_" + std::to_string(i);
         MonitorElement* ch= igetter_.get(hname);
         double mean = ch->getMean();
         double rms  = ch->getRMS();
 
         DetId id = EBDetId::detIdFromDenseIndex(i);
         EcalPedestal ped;
         EcalPedestal oldped=* currentPedestals_->find(id.rawId());
 
         ped.mean_x12=mean;
         ped.rms_x12=rms;
 
         // if bad channel or low stat skip
         if(ch->getEntries()< minEntries_ || !checkStatusCode(id)){
 
             ped.mean_x12=oldped.mean_x12;
             ped.rms_x12=oldped.rms_x12;
 
         }
 
         ped.mean_x6=oldped.mean_x6;
         ped.rms_x6=oldped.rms_x6;
         ped.mean_x1=oldped.mean_x1;
         ped.rms_x1=oldped.rms_x1;
 
         pedestals.setValue(id.rawId(),ped);
     }
 
 
     for (uint16_t i =0; i< EEDetId::kSizeForDenseIndexing; ++i) {
 
         std::string hname = dqmDir_+"/ee_" + std::to_string(i);
      
         MonitorElement* ch= igetter_.get(hname);
         double mean = ch->getMean();
         double rms  = ch->getRMS();
 
         DetId id = EEDetId::detIdFromDenseIndex(i);
         EcalPedestal ped;
         EcalPedestal oldped= *currentPedestals_->find(id.rawId());
 
         ped.mean_x12=mean;
         ped.rms_x12=rms;
 
         // if bad channel or low stat skip
         if(ch->getEntries()< minEntries_ || !checkStatusCode(id)){
             ped.mean_x12=oldped.mean_x12;
             ped.rms_x12=oldped.rms_x12;
         }
 
         ped.mean_x6=oldped.mean_x6;
         ped.rms_x6=oldped.rms_x6;
         ped.mean_x1=oldped.mean_x1;
         ped.rms_x1=oldped.rms_x1;
 
         pedestals.setValue(id.rawId(),ped);
     }
 
 
     // check if there are large variations wrt exisiting pedstals
 
     if (checkAnomalies_){
         if (checkVariation(*currentPedestals_, pedestals)) {
             edm::LogError("Large Variations found wrt to old pedestals, no file created");
             return;
         }
     }
 
     // write out pedestal record
     edm::Service<cond::service::PoolDBOutputService> poolDbService;
 
     if( poolDbService.isAvailable() )
         poolDbService->writeOne( &pedestals, poolDbService->currentTime(),
                                  "EcalPedestalsRcd"  );
     else
         throw std::runtime_error("PoolDBService required.");
 }

void ECALpedestalPCLHarvester::endRun	(	edm::Run const &	run,
		edm::EventSetup const &	isetup
	)

virtual

Reimplemented from DQMEDHarvester.

Definition at line 130 of file ECALpedestalPCLHarvester.cc.

References channelStatus_, currentPedestals_, edm::EventSetup::get(), and edm::ESHandle< T >::product().

                                                                                     {
 
 
     edm::ESHandle<EcalChannelStatus> chStatus;
     isetup.get<EcalChannelStatusRcd>().get(chStatus);
     channelStatus_=chStatus.product();
 
     edm::ESHandle<EcalPedestals> peds;
     isetup.get<EcalPedestalsRcd>().get(peds);
     currentPedestals_=peds.product();
 
 }

void ECALpedestalPCLHarvester::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 122 of file ECALpedestalPCLHarvester.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                      {
 
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }