#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	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

	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

void	endLuminosityBlockProduce (edm::LuminosityBlock &, edm::EventSetup const &) override final

virtual void	produce (edm::Event &, edm::EventSetup const &) override final

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

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

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)

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

Static Public Member Functions inherited from edm::one::EDProducerBase
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_)

void	dqmPlots (const EcalPedestals &newpeds, DQMStore::IBooker &ibooker)

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::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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 146 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 160 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_, dqmPlots(), 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(int(i/100))+"/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(int(i/100))+"/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);
     }
 
 
 
     dqmPlots(pedestals, ibooker_);
 
     // 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::dqmPlots	(	const EcalPedestals &	newpeds,
		DQMStore::IBooker &	ibooker
	)

private

Definition at line 198 of file ECALpedestalPCLHarvester.cc.

References DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), EBDetId::detIdFromDenseIndex(), EEDetId::detIdFromDenseIndex(), dqmDir_, MonitorElement::Fill(), cond::hash, EBDetId::ieta(), EBDetId::iphi(), EEDetId::ix(), EEDetId::iy(), EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, RecoTauDiscriminantConfiguration::mean, RecoTauDiscriminantConfiguration::rms, DQMStore::IBooker::setCurrentFolder(), and EEDetId::zside().

Referenced by dqmEndJob().

                                                                                               {
 
      ibooker.cd();
      ibooker.setCurrentFolder(dqmDir_+"/Summary");
 
      MonitorElement * pmeb = ibooker.book2D("meaneb","Pedestal Means",360, 1., 361., 171, -85., 86.);
      MonitorElement * preb = ibooker.book2D("rmseb","Pedestal RMS",360, 1., 361., 171, -85., 86.);
 
      MonitorElement * pmeep = ibooker.book2D("meaneep","Pedestal Means",100,1,101,100,1,101);
      MonitorElement * preep = ibooker.book2D("rmseep","Pedestal RMS",100,1,101,100,1,101);
 
      MonitorElement * pmeem = ibooker.book2D("meaneem","Pedestal Means",100,1,101,100,1,101);
      MonitorElement * preem = ibooker.book2D("rmseem","Pedestal RMS",100,1,101,100,1,101);
 
      MonitorElement * pmebd = ibooker.book2D("meanebdiff","Pedestal Means Diff",360, 1., 361., 171, -85., 86.);
      MonitorElement * prebd = ibooker.book2D("rmsebdiff","Pedestal RMS Diff",360, 1., 361., 171, -85., 86.);
 
      MonitorElement * pmeepd = ibooker.book2D("meaneepdiff","Pedestal Means Diff",100,1,101,100,1,101);
      MonitorElement * preepd = ibooker.book2D("rmseepdiff","Pedestal RMS Diff",100,1,101,100,1,101);
 
      MonitorElement * pmeemd = ibooker.book2D("meaneemdiff","Pedestal Means Diff",100,1,101,100,1,101);
      MonitorElement * preemd = ibooker.book2D("rmseemdiff","Pedestal RMS Diff",100,1,101,100,1,101);
 
      
      for (int hash =0; hash<EBDetId::kSizeForDenseIndexing;++hash){
          
          EBDetId di= EBDetId::detIdFromDenseIndex(hash); 
          float mean= newpeds[di].mean_x12;
          float rms = newpeds[di].rms_x12;
          
          float cmean = (*currentPedestals_)[di].mean_x12;
          float crms  = (*currentPedestals_)[di].rms_x12;
                   
          pmeb->Fill(di.iphi(),di.ieta(),mean);
          preb->Fill(di.iphi(),di.ieta(),rms);
          pmebd->Fill(di.iphi(),di.ieta(),(mean-cmean)/cmean);
          prebd->Fill(di.iphi(),di.ieta(),(rms-crms)/crms);
 
      }
 
      
      for (int hash =0; hash<EEDetId::kSizeForDenseIndexing;++hash){
          
          EEDetId di= EEDetId::detIdFromDenseIndex(hash); 
          float mean= newpeds[di].mean_x12;
          float rms = newpeds[di].rms_x12;
          float cmean = (*currentPedestals_)[di].mean_x12;
          float crms  = (*currentPedestals_)[di].rms_x12;
 
          if (di.zside() >0){
              pmeep->Fill(di.ix(),di.iy(),mean);
              preep->Fill(di.ix(),di.iy(),rms);
              pmeepd->Fill(di.ix(),di.iy(),(mean-cmean)/cmean);
              preepd->Fill(di.ix(),di.iy(),(rms-crms)/crms);
          } else{
              pmeem->Fill(di.ix(),di.iy(),mean);
              preem->Fill(di.ix(),di.iy(),rms);
              pmeemd->Fill(di.ix(),di.iy(),(mean-cmean)/cmean);
              preemd->Fill(di.ix(),di.iy(),(rms-crms)/crms);
          }
          
      }
 
      
 
 }

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

virtual

Reimplemented from DQMEDHarvester.

Definition at line 133 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 125 of file ECALpedestalPCLHarvester.cc.

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

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