#include <ECALpedestalPCLHarvester.h>

Inheritance diagram for ECALpedestalPCLHarvester:

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

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

Public Member Functions inherited from DQMEDHarvester
void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final

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

	DQMEDHarvester (void)

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

void	endJob () final

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

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

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

	~DQMEDHarvester () override=default

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

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies 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)

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

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

bool	isGood (const DetId &id)

Private Attributes
const EcalChannelStatus *	channelStatus_

bool	checkAnomalies_

std::vector< int >	chStatusToExclude_

const EcalPedestals *	currentPedestals_

std::string	dqmDir_

int	entriesEB_ [EBDetId::kSizeForDenseIndexing]

int	entriesEE_ [EEDetId::kSizeForDenseIndexing]

const EcalPedestals *	g6g1Pedestals_

std::string	labelG6G1_

int	minEntries_

double	nSigma_

float	threshDiffEB_

float	threshDiffEE_

double	thresholdAnomalies_

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

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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 35 of file ECALpedestalPCLHarvester.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file ECALpedestalPCLHarvester.cc.

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

                                                                            :
     currentPedestals_(nullptr),channelStatus_(nullptr){
 
     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");
     labelG6G1_      = ps.getParameter<std::string>("labelG6G1");
     threshDiffEB_   = ps.getParameter<double>("threshDiffEB");
     threshDiffEE_   = ps.getParameter<double>("threshDiffEE");
 }

Member Function Documentation

bool ECALpedestalPCLHarvester::checkStatusCode ( const DetId & id )

private

Definition at line 161 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 186 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_
	)

overrideprivatevirtual

Implements DQMEDHarvester.

Definition at line 32 of file ECALpedestalPCLHarvester.cc.

References funct::abs(), checkAnomalies_, checkStatusCode(), checkVariation(), currentPedestals_, cond::service::PoolDBOutputService::currentTime(), EBDetId::detIdFromDenseIndex(), EEDetId::detIdFromDenseIndex(), diffTreeTool::diff, dqmDir_, dqmPlots(), entriesEB_, entriesEE_, EcalCondObjectContainer< T >::find(), g6g1Pedestals_, DQMStore::IGetter::get(), MonitorElement::getEntries(), MonitorElement::getMean(), MonitorElement::getRMS(), mps_fire::i, edm::Service< T >::isAvailable(), EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, SiStripPI::mean, EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, minEntries_, SiStripPI::rms, EcalPedestal::rms_x1, EcalPedestal::rms_x12, EcalPedestal::rms_x6, EcalCondObjectContainer< T >::setValue(), AlCaHLTBitMon_QueryRunRegistry::string, threshDiffEB_, threshDiffEE_, 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();
         entriesEB_[i] = ch->getEntries();
 
         DetId id = EBDetId::detIdFromDenseIndex(i);
         EcalPedestal ped;
         EcalPedestal oldped =* currentPedestals_->find(id.rawId());
         EcalPedestal g6g1ped=* g6g1Pedestals_->find(id.rawId());
 
         ped.mean_x12=mean;
         ped.rms_x12=rms;
 
         float diff = std::abs(mean-oldped.mean_x12);
 
         // if bad channel or low stat skip or the difference is too large wrt to previous record
         if(ch->getEntries()< minEntries_ || !checkStatusCode(id) || diff>threshDiffEB_){
 
             ped.mean_x12=oldped.mean_x12;
             ped.rms_x12=oldped.rms_x12;
 
         }
         
         // copy g6 and g1 from the corressponding record
         ped.mean_x6= g6g1ped.mean_x6;
         ped.rms_x6 = g6g1ped.rms_x6;
         ped.mean_x1= g6g1ped.mean_x1;
         ped.rms_x1 = g6g1ped.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();
         entriesEE_[i] = ch->getEntries(); 
  
         DetId id = EEDetId::detIdFromDenseIndex(i);
         EcalPedestal ped;
         EcalPedestal oldped = *currentPedestals_->find(id.rawId());
         EcalPedestal g6g1ped= *g6g1Pedestals_->find(id.rawId());
                 
         ped.mean_x12=mean;
         ped.rms_x12=rms;
 
         float diff = std::abs(mean-oldped.mean_x12);
 
         // if bad channel or low stat skip or the difference is too large wrt to previous record
         if(ch->getEntries()< minEntries_ || !checkStatusCode(id)|| diff>threshDiffEE_){
             ped.mean_x12=oldped.mean_x12;
             ped.rms_x12=oldped.rms_x12;
         }
 
         // copy g6 and g1 pedestals from corresponding record
         ped.mean_x6= g6g1ped.mean_x6;
         ped.rms_x6 = g6g1ped.rms_x6;
         ped.mean_x1= g6g1ped.mean_x1;
         ped.rms_x1 = g6g1ped.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 224 of file ECALpedestalPCLHarvester.cc.

References funct::abs(), DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), EBDetId::detIdFromDenseIndex(), EEDetId::detIdFromDenseIndex(), dqmDir_, entriesEB_, entriesEE_, MonitorElement::Fill(), cond::hash, EBDetId::ieta(), EBDetId::iphi(), isGood(), EEDetId::ix(), EEDetId::iy(), EBDetId::kSizeForDenseIndexing, EEDetId::kSizeForDenseIndexing, SiStripPI::mean, SiStripPI::rms, DQMStore::IBooker::setCurrentFolder(), and EEDetId::zside().

Referenced by dqmEndJob().

                                                                                               {
 
      ibooker.cd();
      ibooker.setCurrentFolder(dqmDir_+"/Summary");
 
      MonitorElement * pmeb = ibooker.book2D("meaneb","Pedestal Means EB",360, 1., 361., 171, -85., 86.);
      MonitorElement * preb = ibooker.book2D("rmseb","Pedestal RMS EB ",360, 1., 361., 171, -85., 86.);
 
      MonitorElement * pmeep = ibooker.book2D("meaneep","Pedestal Means EEP",100,1,101,100,1,101);
      MonitorElement * preep = ibooker.book2D("rmseep","Pedestal RMS EEP",100,1,101,100,1,101);
 
      MonitorElement * pmeem = ibooker.book2D("meaneem","Pedestal Means EEM",100,1,101,100,1,101);
      MonitorElement * preem = ibooker.book2D("rmseem","Pedestal RMS EEM",100,1,101,100,1,101);
 
      MonitorElement * pmebd = ibooker.book2D("meanebdiff","Abs Rel Pedestal Means Diff EB",360,1., 361., 171, -85., 86.);
      MonitorElement * prebd = ibooker.book2D("rmsebdiff","Abs Rel Pedestal RMS Diff E ",360, 1., 361., 171, -85., 86.);
 
      MonitorElement * pmeepd = ibooker.book2D("meaneepdiff","Abs Rel Pedestal Means Diff EEP",100,1,101,100,1,101);
      MonitorElement * preepd = ibooker.book2D("rmseepdiff","Abs Rel Pedestal RMS Diff EEP",100,1,101,100,1,101);
 
      MonitorElement * pmeemd = ibooker.book2D("meaneemdiff","Abs Rel Pedestal Means Diff EEM",100,1,101,100,1,101);
      MonitorElement * preemd = ibooker.book2D("rmseemdiff","Abs RelPedestal RMS Diff EEM",100,1,101,100,1,101);
 
      MonitorElement * poeb = ibooker.book2D("occeb","Occupancy EB",360, 1., 361., 171, -85., 86.);
      MonitorElement * poeep = ibooker.book2D("occeep","Occupancy EEP",100,1,101,100,1,101);
      MonitorElement * poeem = ibooker.book2D("occeem","Occupancy EEM",100,1,101,100,1,101);
 
 
      MonitorElement * hdiffeb = ibooker.book1D("diffeb","Pedestal Differences EB",100,-2.5,2.5);
      MonitorElement * hdiffee = ibooker.book1D("diffee","Pedestal Differences EE",100,-2.5,2.5);
 
      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;
                   
          if (!isGood(di) ) continue;   // only good channels are plotted
  
          pmeb->Fill(di.iphi(),di.ieta(),mean);
          preb->Fill(di.iphi(),di.ieta(),rms);
          if (cmean) pmebd->Fill(di.iphi(),di.ieta(),std::abs(mean-cmean)/cmean);
          if (crms) prebd->Fill(di.iphi(),di.ieta(),std::abs(rms-crms)/crms);
          poeb->Fill(di.iphi(),di.ieta(),entriesEB_[hash]);
          hdiffeb->Fill(mean-cmean);
      }
 
      
      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 (!isGood(di) ) continue;   // only good channels are plotted
          
          if (di.zside() >0){
              pmeep->Fill(di.ix(),di.iy(),mean);
              preep->Fill(di.ix(),di.iy(),rms);
              poeep->Fill(di.ix(),di.iy(),entriesEE_[hash]);
              if (cmean) pmeepd->Fill(di.ix(),di.iy(),std::abs(mean-cmean)/cmean);
              if (crms)  preepd->Fill(di.ix(),di.iy(),std::abs(rms-crms)/crms);
          } else{
              pmeem->Fill(di.ix(),di.iy(),mean);
              preem->Fill(di.ix(),di.iy(),rms);
              if (cmean) pmeemd->Fill(di.ix(),di.iy(),std::abs(mean-cmean)/cmean);
              poeem->Fill(di.ix(),di.iy(),entriesEE_[hash]); 
              if (crms)preemd->Fill(di.ix(),di.iy(),std::abs(rms-crms)/crms);
          }
          hdiffee->Fill(mean-cmean);
          
      }
 
      
 
 }

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

override

Definition at line 144 of file ECALpedestalPCLHarvester.cc.

References channelStatus_, currentPedestals_, g6g1Pedestals_, edm::EventSetup::get(), labelG6G1_, 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();
 
     edm::ESHandle<EcalPedestals> g6g1peds;
     isetup.get<EcalPedestalsRcd>().get(labelG6G1_,g6g1peds);
     g6g1Pedestals_=peds.product();
 
 }

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

static

Definition at line 136 of file ECALpedestalPCLHarvester.cc.

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

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

bool ECALpedestalPCLHarvester::isGood ( const DetId & id )

private

Definition at line 174 of file ECALpedestalPCLHarvester.cc.

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

Referenced by dqmPlots().

                                                     {
     
     EcalChannelStatusMap::const_iterator dbstatusPtr;
     dbstatusPtr = channelStatus_->getMap().find(id.rawId());
     if (dbstatusPtr == channelStatus_->getMap().end())
         edm::LogError("Invalid DetId supplied");
     EcalChannelStatusCode::Code  dbstatus = dbstatusPtr->getStatusCode();
     if (dbstatus ==0 ) return true;
     return false;
 }