Inheritance diagram for CorrPCCProducer:

Public Member Functions
	CorrPCCProducer (const edm::ParameterSet &)

	~CorrPCCProducer () override

Public Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, const edm::EventSetup &iSetup) final

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

void	calculateCorrections (std::vector< float >, std::vector< float > &, float &)

void	dqmEndRun (edm::Run const &runSeg, const edm::EventSetup &iSetup) final

void	dqmEndRunProduce (const edm::Run &runSeg, const edm::EventSetup &iSetup)

void	endJob () final

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

void	estimateType1Frac (std::vector< float >, float &)

void	evaluateCorrectionResiduals (std::vector< float >)

float	getMaximum (std::vector< float >)

void	makeCorrectionTemplate ()

void	resetBlock ()

Private Attributes
unsigned int	approxLumiBlockSize_

std::vector< float >	correctionScaleFactors_

std::vector< float >	correctionTemplate_

TH1F *	corrlumiAvg_h

unsigned int	countLumi_

std::vector< float >	errOnLumiByBX_

std::vector< float >	events_

TFile *	histoFile

TList *	hlist

unsigned int	iBlock = 0

TH1F *	lumiAvg_h

std::map< std::pair< unsigned int, unsigned int >, unsigned int >	lumiInfoCounter

std::map< std::pair< unsigned int, unsigned int >, LumiInfo * >	lumiInfoMap

std::map< std::pair< unsigned int, unsigned int >, LumiInfo * >::iterator	lumiInfoMapIterator

std::map< unsigned int, LumiInfo * >	lumiInfoMapPerLS

edm::EDGetTokenT< LumiInfo >	lumiInfoToken

std::vector< unsigned int >	lumiSections

unsigned int	maxLS = 3500

float	mean_type1_residual

float	mean_type1_residual_unc

float	mean_type2_residual

float	mean_type2_residual_unc

unsigned int	minimumNumberOfEvents

unsigned int	nTrain

float	overallCorrection_

std::string	pccSrc_

float	pedestal

float	pedestal_unc

TGraphErrors *	pedestalGraph

edm::Service< cond::service::PoolDBOutputService >	poolDbService

std::string	prodInst_

std::vector< float >	rawlumiBX_

TH1F *	scaleFactorAvg_h

unsigned int	thisLS

std::vector< float >	totalLumiByBX_

std::vector< float >	totalLumiByBXAvg_

float	type1Frac

TGraphErrors *	type1FracGraph

TH1F *	type1FracHist

MonitorElement *	Type1FracMon

TGraphErrors *	type1resGraph

TH1F *	type1resHist

MonitorElement *	Type1ResMon

double	type2_a_

double	type2_b_

TGraphErrors *	type2resGraph

TH1F *	type2resHist

MonitorElement *	Type2ResMon

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Protected Attributes inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

class: CorrPCCProducer.cc

description: Computes the type 1 and type 2 corrections to the luminosity type 1 - first (spillover from previous BXs real clusters) type 2 - after (comes from real activation) pedestal - is a constant noise term for low lumi period

authors:Sam Higginbotham (shigg.nosp@m.inb@.nosp@m.cern..nosp@m.ch) and Chris Palmer (capal.nosp@m.mer@.nosp@m.cern..nosp@m.ch) , Jose Benitez (jose..nosp@m.beni.nosp@m.tez@c.nosp@m.ern..nosp@m.ch)

Definition at line 50 of file CorrPCCProducer.cc.

Constructor & Destructor Documentation

◆ CorrPCCProducer()

CorrPCCProducer::CorrPCCProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 133 of file CorrPCCProducer.cc.

References approxLumiBlockSize_, correctionScaleFactors_, correctionTemplate_, countLumi_, events_, edm::ParameterSet::getParameter(), histoFile, lumiInfoToken, makeCorrectionTemplate(), minimumNumberOfEvents, LumiConstants::numBX, pccSrc_, pedestalGraph, prodInst_, resetBlock(), totalLumiByBX_, totalLumiByBXAvg_, type1FracGraph, type1resGraph, type2_a_, type2_b_, and type2resGraph.

                                                                {
   pccSrc_ =
       iConfig.getParameter<edm::ParameterSet>("CorrPCCProducerParameters").getParameter<std::string>("inLumiObLabel");
   prodInst_ =
       iConfig.getParameter<edm::ParameterSet>("CorrPCCProducerParameters").getParameter<std::string>("ProdInst");
   approxLumiBlockSize_ =
       iConfig.getParameter<edm::ParameterSet>("CorrPCCProducerParameters").getParameter<int>("approxLumiBlockSize");
   type2_a_ = iConfig.getParameter<edm::ParameterSet>("CorrPCCProducerParameters").getParameter<double>("type2_a");
   type2_b_ = iConfig.getParameter<edm::ParameterSet>("CorrPCCProducerParameters").getParameter<double>("type2_b");
   countLumi_ = 0;
   minimumNumberOfEvents = 1000;
 
   totalLumiByBX_.resize(LumiConstants::numBX);
   totalLumiByBXAvg_.resize(LumiConstants::numBX);
   events_.resize(LumiConstants::numBX);
   correctionScaleFactors_.resize(LumiConstants::numBX);
   correctionTemplate_.resize(LumiConstants::numBX);
 
   resetBlock();
 
   makeCorrectionTemplate();
 
   edm::InputTag inputPCCTag_(pccSrc_, prodInst_);
 
   lumiInfoToken = consumes<LumiInfo, edm::InLumi>(inputPCCTag_);
 
   histoFile = new TFile("CorrectionHisto.root", "RECREATE");
 
   type1FracGraph = new TGraphErrors();
   type1resGraph = new TGraphErrors();
   type2resGraph = new TGraphErrors();
   type1FracGraph->SetName("Type1Fraction");
   type1resGraph->SetName("Type1Res");
   type2resGraph->SetName("Type2Res");
   type1FracGraph->GetYaxis()->SetTitle("Type 1 Fraction");
   type1resGraph->GetYaxis()->SetTitle("Type 1 Residual");
   type2resGraph->GetYaxis()->SetTitle("Type 2 Residual");
   type1FracGraph->GetXaxis()->SetTitle("Unique LS ID");
   type1resGraph->GetXaxis()->SetTitle("Unique LS ID");
   type2resGraph->GetXaxis()->SetTitle("Unique LS ID");
   type1FracGraph->SetMarkerStyle(8);
   type1resGraph->SetMarkerStyle(8);
   type2resGraph->SetMarkerStyle(8);
 
   pedestalGraph = new TGraphErrors();
   pedestalGraph->SetName("Pedestal");
   pedestalGraph->GetYaxis()->SetTitle("pedestal value (counts) per lumi section");
   pedestalGraph->GetXaxis()->SetTitle("Unique LS ID");
   pedestalGraph->SetMarkerStyle(8);
 }

◆ ~CorrPCCProducer()

CorrPCCProducer::~CorrPCCProducer ( )

override

Definition at line 185 of file CorrPCCProducer.cc.

185 {}

Member Function Documentation

◆ beginLuminosityBlock()

void CorrPCCProducer::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		const edm::EventSetup &	iSetup
	)

finalprivate

Definition at line 377 of file CorrPCCProducer.cc.

References countLumi_.

                                                                                                          {
   countLumi_++;
 }

◆ bookHistograms()

void CorrPCCProducer::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	context
	)

overrideprivatevirtual

Implements DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 670 of file CorrPCCProducer.cc.

References dqm::implementation::IBooker::book1D(), maxLS, dqm::implementation::NavigatorBase::setCurrentFolder(), Type1FracMon, Type1ResMon, and Type2ResMon.

                                                                                                                {
   ibooker.setCurrentFolder("AlCaReco/LumiPCC/");
   auto scope = DQMStore::IBooker::UseRunScope(ibooker);
   Type1FracMon = ibooker.book1D("type1Fraction", "Type1Fraction;Lumisection;Type 1 Fraction", maxLS, 0, maxLS);
   Type1ResMon = ibooker.book1D("type1Residual", "Type1Residual;Lumisection;Type 1 Residual", maxLS, 0, maxLS);
   Type2ResMon = ibooker.book1D("type2Residual", "Type2Residual;Lumisection;Type 2 Residual", maxLS, 0, maxLS);
 }

◆ calculateCorrections()

void CorrPCCProducer::calculateCorrections	(	std::vector< float >	uncorrected,
		std::vector< float > &	correctionScaleFactors_,
		float &	overallCorrection_
	)

private

Definition at line 296 of file CorrPCCProducer.cc.

References correctionScaleFactors_, correctionTemplate_, estimateType1Frac(), evaluateCorrectionResiduals(), getMaximum(), mps_fire::i, dqmiolumiharvest::j, SiStripPI::max, LumiConstants::numBX, overallCorrection_, pedestal, pedestal_unc, mathSSE::sqrt(), remoteMonitoring_LASER_era2018_cfg::threshold, dtNoiseDBValidation_cfg::trial, and type1Frac.

Referenced by dqmEndRunProduce().

                                                                       {
   type1Frac = 0;
 
   int nTrials = 4;
 
   for (int trial = 0; trial < nTrials; trial++) {
     estimateType1Frac(uncorrected, type1Frac);
     edm::LogInfo("INFO") << "type 1 fraction after iteration " << trial << " is  " << type1Frac;
   }
 
   //correction should never be negative
   type1Frac = std::max(0.0, (double)type1Frac);
 
   std::vector<float> corrected_tmp_;
   for (size_t i = 0; i < uncorrected.size(); i++) {
     corrected_tmp_.push_back(uncorrected.at(i));
   }
 
   //Apply all corrections
   for (size_t i = 0; i < LumiConstants::numBX - 1; i++) {
     // type 1 - first (spillover from previous BXs real clusters)
     float bin_i = corrected_tmp_.at(i);
     corrected_tmp_.at(i + 1) = corrected_tmp_.at(i + 1) - type1Frac * bin_i;
 
     // type 2 - after (comes from real activation)
     bin_i = corrected_tmp_.at(i);
     for (size_t j = i + 1; j < i + LumiConstants::numBX - 1; j++) {
       if (j < LumiConstants::numBX) {
         corrected_tmp_.at(j) = corrected_tmp_.at(j) - bin_i * correctionTemplate_.at(j - i);
       } else {
         corrected_tmp_.at(j - LumiConstants::numBX) =
             corrected_tmp_.at(j - LumiConstants::numBX) - bin_i * correctionTemplate_.at(j - i);
       }
     }
   }
 
   float lumiMax = getMaximum(corrected_tmp_);
   float threshold = lumiMax * 0.2;
 
   //here subtract the pedestal
   pedestal = 0.;
   pedestal_unc = 0.;
   int nped = 0;
   for (size_t i = 0; i < LumiConstants::numBX; i++) {
     if (corrected_tmp_.at(i) < threshold) {
       pedestal += corrected_tmp_.at(i);
       nped++;
     }
   }
   if (nped > 0) {
     pedestal_unc = sqrt(pedestal) / nped;
     pedestal = pedestal / nped;
   }
   for (size_t i = 0; i < LumiConstants::numBX; i++) {
     corrected_tmp_.at(i) = corrected_tmp_.at(i) - pedestal;
   }
 
   evaluateCorrectionResiduals(corrected_tmp_);
 
   float integral_uncorr_clusters = 0;
   float integral_corr_clusters = 0;
 
   //Calculate Per-BX correction factor and overall correction factor
   for (size_t ibx = 0; ibx < corrected_tmp_.size(); ibx++) {
     if (corrected_tmp_.at(ibx) > threshold) {
       integral_uncorr_clusters += uncorrected.at(ibx);
       integral_corr_clusters += corrected_tmp_.at(ibx);
     }
     if (corrected_tmp_.at(ibx) != 0.0 && uncorrected.at(ibx) != 0.0) {
       correctionScaleFactors_.at(ibx) = corrected_tmp_.at(ibx) / uncorrected.at(ibx);
     } else {
       correctionScaleFactors_.at(ibx) = 0.0;
     }
   }
 
   overallCorrection_ = integral_corr_clusters / integral_uncorr_clusters;
 }

◆ dqmEndRun()

void CorrPCCProducer::dqmEndRun	(	edm::Run const &	runSeg,
		const edm::EventSetup &	iSetup
	)

finalprivatevirtual

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 413 of file CorrPCCProducer.cc.

References dqmEndRunProduce().

                                                                                  {
   // TODO: why was this code not put here in the first place?
   dqmEndRunProduce(runSeg, iSetup);
 }

◆ dqmEndRunProduce()

void CorrPCCProducer::dqmEndRunProduce	(	const edm::Run &	runSeg,
		const edm::EventSetup &	iSetup
	)

private

Definition at line 419 of file CorrPCCProducer.cc.

References approxLumiBlockSize_, simKBmtfDigis_cfi::bx, calculateCorrections(), correctionScaleFactors_, corrlumiAvg_h, gather_cfg::cout, errOnLumiByBX_, events_, dqmdumpme::first, dqmMemoryStats::float, histoFile, mps_fire::i, iBlock, edm::RunBase::id(), edm::Service< T >::isAvailable(), lumiAvg_h, lumiInfoCounter, lumiInfoMap, lumiInfoMapIterator, lumiInfoMapPerLS, lumiSections, maxLS, mean_type1_residual, mean_type1_residual_unc, mean_type2_residual, mean_type2_residual_unc, LumiConstants::numBX, overallCorrection_, pedestal, pedestal_unc, pedestalGraph, poolDbService, rawlumiBX_, resetBlock(), edm::RunID::run(), edm::RunBase::run(), scaleFactorAvg_h, dqm::impl::MonitorElement::setBinContent(), dqm::impl::MonitorElement::setBinError(), LumiCorrections::setCorrectionsBX(), LumiCorrections::setOverallCorrection(), LumiCorrections::setType1Fraction(), LumiCorrections::setType1Residual(), LumiCorrections::setType2Residual(), jetUpdater_cfi::sort, thisLS, totalLumiByBX_, totalLumiByBXAvg_, type1Frac, type1FracGraph, type1FracHist, Type1FracMon, type1resGraph, type1resHist, Type1ResMon, type2resGraph, type2resHist, Type2ResMon, and cond::service::PoolDBOutputService::writeOneIOV().

Referenced by dqmEndRun().

                                                                                         {
   if (lumiSections.empty()) {
     return;
   }
 
   std::sort(lumiSections.begin(), lumiSections.end());
 
   edm::LogInfo("INFO") << "Number of Lumisections " << lumiSections.size() << " in run " << runSeg.run();
 
   //Determining integer number of blocks
   float nBlocks_f = float(lumiSections.size()) / approxLumiBlockSize_;
   unsigned int nBlocks = 1;
   if (nBlocks_f > 1) {
     if (nBlocks_f - lumiSections.size() / approxLumiBlockSize_ < 0.5) {
       nBlocks = lumiSections.size() / approxLumiBlockSize_;
     } else {
       nBlocks = lumiSections.size() / approxLumiBlockSize_ + 1;
     }
   }
 
   float nLSPerBlock = float(lumiSections.size()) / nBlocks;
 
   std::vector<std::pair<unsigned int, unsigned int>> lsKeys;
   lsKeys.clear();
 
   //Constructing nBlocks IOVs
   for (unsigned iKey = 0; iKey < nBlocks; iKey++) {
     lsKeys.push_back(std::make_pair(lumiSections[(unsigned int)(iKey * nLSPerBlock)],
                                     lumiSections[(unsigned int)((iKey + 1) * nLSPerBlock) - 1]));
   }
 
   lsKeys[0].first = 1;
 
   for (unsigned int lumiSection = 0; lumiSection < lumiSections.size(); lumiSection++) {
     thisLS = lumiSections[lumiSection];
 
     std::pair<unsigned int, unsigned int> lsKey;
 
     bool foundKey = false;
 
     for (unsigned iKey = 0; iKey < nBlocks; iKey++) {
       if ((thisLS >= lsKeys[iKey].first) && (thisLS <= lsKeys[iKey].second)) {
         lsKey = lsKeys[iKey];
         foundKey = true;
         break;
       }
     }
 
     if (!foundKey) {
       edm::LogInfo("WARNING") << "Didn't find key " << thisLS;
       continue;
     }
 
     if (lumiInfoMap.count(lsKey) == 0) {
       lumiInfoMap[lsKey] = new LumiInfo();
     }
 
     //Sum all lumi in IOV of lsKey
     totalLumiByBX_ = lumiInfoMap[lsKey]->getInstLumiAllBX();
     events_ = lumiInfoMap[lsKey]->getErrorLumiAllBX();
 
     rawlumiBX_ = lumiInfoMapPerLS[thisLS]->getInstLumiAllBX();
     errOnLumiByBX_ = lumiInfoMapPerLS[thisLS]->getErrorLumiAllBX();
 
     for (unsigned int bx = 0; bx < LumiConstants::numBX; bx++) {
       totalLumiByBX_[bx] += rawlumiBX_[bx];
       events_[bx] += errOnLumiByBX_[bx];
     }
     lumiInfoMap[lsKey]->setInstLumiAllBX(totalLumiByBX_);
     lumiInfoMap[lsKey]->setErrorLumiAllBX(events_);
     lumiInfoCounter[lsKey]++;
   }
 
   cond::Time_t thisIOV = 1;
 
   char* histname1 = new char[100];
   char* histname2 = new char[100];
   char* histname3 = new char[100];
   char* histTitle1 = new char[100];
   char* histTitle2 = new char[100];
   char* histTitle3 = new char[100];
   sprintf(histTitle1, "Type1Fraction_%d", runSeg.run());
   sprintf(histTitle2, "Type1Res_%d", runSeg.run());
   sprintf(histTitle3, "Type2Res_%d", runSeg.run());
   type1FracHist = new TH1F(histTitle1, histTitle1, 1000, -0.5, 0.5);
   type1resHist = new TH1F(histTitle2, histTitle2, 4000, -0.2, 0.2);
   type2resHist = new TH1F(histTitle3, histTitle3, 4000, -0.2, 0.2);
   delete[] histTitle1;
   delete[] histTitle2;
   delete[] histTitle3;
 
   for (lumiInfoMapIterator = lumiInfoMap.begin(); (lumiInfoMapIterator != lumiInfoMap.end()); ++lumiInfoMapIterator) {
     totalLumiByBX_ = lumiInfoMapIterator->second->getInstLumiAllBX();
     events_ = lumiInfoMapIterator->second->getErrorLumiAllBX();
 
     if (events_.empty()) {
       continue;
     }
 
     edm::LuminosityBlockID lu(runSeg.id().run(), edm::LuminosityBlockNumber_t(lumiInfoMapIterator->first.first));
     thisIOV = (cond::Time_t)(lu.value());
 
     sprintf(histname1,
             "CorrectedLumiAvg_%d_%d_%d_%d",
             runSeg.run(),
             iBlock,
             lumiInfoMapIterator->first.first,
             lumiInfoMapIterator->first.second);
     sprintf(histname2,
             "ScaleFactorsAvg_%d_%d_%d_%d",
             runSeg.run(),
             iBlock,
             lumiInfoMapIterator->first.first,
             lumiInfoMapIterator->first.second);
     sprintf(histname3,
             "RawLumiAvg_%d_%d_%d_%d",
             runSeg.run(),
             iBlock,
             lumiInfoMapIterator->first.first,
             lumiInfoMapIterator->first.second);
 
     corrlumiAvg_h = new TH1F(histname1, "", LumiConstants::numBX, 1, LumiConstants::numBX);
     scaleFactorAvg_h = new TH1F(histname2, "", LumiConstants::numBX, 1, LumiConstants::numBX);
     lumiAvg_h = new TH1F(histname3, "", LumiConstants::numBX, 1, LumiConstants::numBX);
 
     //Averaging by the number of events
     for (unsigned int i = 0; i < LumiConstants::numBX; i++) {
       if (events_.at(i) != 0) {
         totalLumiByBXAvg_[i] = totalLumiByBX_[i] / events_[i];
       } else {
         totalLumiByBXAvg_[i] = 0.0;
       }
     }
 
     calculateCorrections(totalLumiByBXAvg_, correctionScaleFactors_, overallCorrection_);
 
     for (unsigned int bx = 0; bx < LumiConstants::numBX; bx++) {
       corrlumiAvg_h->SetBinContent(bx, totalLumiByBXAvg_[bx] * correctionScaleFactors_[bx]);
       if (events_.at(bx) != 0) {
         corrlumiAvg_h->SetBinError(bx,
                                    totalLumiByBXAvg_[bx] * correctionScaleFactors_[bx] / TMath::Sqrt(events_.at(bx)));
       } else {
         corrlumiAvg_h->SetBinError(bx, 0.0);
       }
 
       scaleFactorAvg_h->SetBinContent(bx, correctionScaleFactors_[bx]);
       lumiAvg_h->SetBinContent(bx, totalLumiByBXAvg_[bx]);
     }
 
     //Writing the corrections to SQL lite file for db
     LumiCorrections pccCorrections;
     pccCorrections.setOverallCorrection(overallCorrection_);
     pccCorrections.setType1Fraction(type1Frac);
     pccCorrections.setType1Residual(mean_type1_residual);
     pccCorrections.setType2Residual(mean_type2_residual);
     pccCorrections.setCorrectionsBX(correctionScaleFactors_);
 
     if (poolDbService.isAvailable()) {
       poolDbService->writeOneIOV(pccCorrections, thisIOV, "LumiCorrectionsRcd");
     } else {
       throw std::runtime_error("PoolDBService required.");
     }
 
     histoFile->cd();
     corrlumiAvg_h->Write();
     scaleFactorAvg_h->Write();
     lumiAvg_h->Write();
 
     delete corrlumiAvg_h;
     delete scaleFactorAvg_h;
     delete lumiAvg_h;
 
     type1FracHist->Fill(type1Frac);
     type1resHist->Fill(mean_type1_residual);
     type2resHist->Fill(mean_type2_residual);
 
     for (unsigned int ils = lumiInfoMapIterator->first.first; ils < lumiInfoMapIterator->first.second + 1; ils++) {
       if (ils > maxLS) {
         std::cout << "ils out of maxLS range!!" << std::endl;
         break;
       }
       Type1FracMon->setBinContent(ils, type1Frac);
       Type1FracMon->setBinError(ils, mean_type1_residual_unc);
       Type1ResMon->setBinContent(ils, mean_type1_residual);
       Type1ResMon->setBinError(ils, mean_type1_residual_unc);
       Type2ResMon->setBinContent(ils, mean_type2_residual);
       Type2ResMon->setBinError(ils, mean_type2_residual_unc);
     }
 
     type1FracGraph->SetPoint(iBlock, thisIOV + approxLumiBlockSize_ / 2.0, type1Frac);
     type1resGraph->SetPoint(iBlock, thisIOV + approxLumiBlockSize_ / 2.0, mean_type1_residual);
     type2resGraph->SetPoint(iBlock, thisIOV + approxLumiBlockSize_ / 2.0, mean_type2_residual);
     type1FracGraph->SetPointError(iBlock, approxLumiBlockSize_ / 2.0, mean_type1_residual_unc);
     type1resGraph->SetPointError(iBlock, approxLumiBlockSize_ / 2.0, mean_type1_residual_unc);
     type2resGraph->SetPointError(iBlock, approxLumiBlockSize_ / 2.0, mean_type2_residual_unc);
     pedestalGraph->SetPoint(iBlock, thisIOV + approxLumiBlockSize_ / 2.0, pedestal);
     pedestalGraph->SetPointError(iBlock, approxLumiBlockSize_ / 2.0, pedestal_unc);
 
     edm::LogInfo("INFO")
         << "iBlock type1Frac mean_type1_residual mean_type2_residual mean_type1_residual_unc mean_type2_residual_unc "
         << iBlock << " " << type1Frac << " " << mean_type1_residual << " " << mean_type2_residual << " "
         << mean_type1_residual_unc << " " << mean_type2_residual_unc;
 
     type1Frac = 0.0;
     mean_type1_residual = 0.0;
     mean_type2_residual = 0.0;
     mean_type1_residual_unc = 0;
     mean_type2_residual_unc = 0;
     pedestal = 0.;
     pedestal_unc = 0.;
 
     iBlock++;
 
     resetBlock();
   }
   histoFile->cd();
   type1FracHist->Write();
   type1resHist->Write();
   type2resHist->Write();
 
   delete type1FracHist;
   delete type1resHist;
   delete type2resHist;
 
   delete[] histname1;
   delete[] histname2;
   delete[] histname3;
 
   for (lumiInfoMapIterator = lumiInfoMap.begin(); (lumiInfoMapIterator != lumiInfoMap.end()); ++lumiInfoMapIterator) {
     delete lumiInfoMapIterator->second;
   }
   for (unsigned int lumiSection = 0; lumiSection < lumiSections.size(); lumiSection++) {
     thisLS = lumiSections[lumiSection];
     delete lumiInfoMapPerLS[thisLS];
   }
   lumiInfoMap.clear();
   lumiInfoMapPerLS.clear();
   lumiSections.clear();
   lumiInfoCounter.clear();
 }

◆ endJob()

void CorrPCCProducer::endJob ( void )

finalprivatevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 678 of file CorrPCCProducer.cc.

References histoFile, pedestalGraph, type1FracGraph, type1resGraph, and type2resGraph.

Referenced by o2olib.O2ORunMgr::executeJob().

                              {
   histoFile->cd();
   type1FracGraph->Write();
   type1resGraph->Write();
   type2resGraph->Write();
   pedestalGraph->Write();
   histoFile->Write();
   histoFile->Close();
   delete type1FracGraph;
   delete type1resGraph;
   delete type2resGraph;
   delete pedestalGraph;
 }

◆ endLuminosityBlock()

void CorrPCCProducer::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		const edm::EventSetup &	iSetup
	)

finalprivate

Definition at line 382 of file CorrPCCProducer.cc.

References simKBmtfDigis_cfi::bx, errOnLumiByBX_, events_, edm::LuminosityBlock::getByToken(), LumiInfo::getErrorLumiAllBX(), LumiInfo::getInstLumiAllBX(), lumiInfoMapPerLS, lumiInfoToken, edm::LuminosityBlockBase::luminosityBlock(), lumiSections, minimumNumberOfEvents, LumiConstants::numBX, edm::Handle< T >::product(), thisLS, and totalLumiByBX_.

                                                                                                        {
   thisLS = lumiSeg.luminosityBlock();
 
   edm::Handle<LumiInfo> PCCHandle;
   lumiSeg.getByToken(lumiInfoToken, PCCHandle);
 
   const LumiInfo& inLumiOb = *(PCCHandle.product());
 
   errOnLumiByBX_ = inLumiOb.getErrorLumiAllBX();
 
   unsigned int totalEvents = 0;
   for (unsigned int bx = 0; bx < LumiConstants::numBX; bx++) {
     totalEvents += errOnLumiByBX_[bx];
   }
 
   if (totalEvents < minimumNumberOfEvents) {
     edm::LogInfo("INFO") << "number of events in this LS is too few " << totalEvents;
     //return;
   } else {
     edm::LogInfo("INFO") << "Skipping Lumisection " << thisLS;
   }
 
   lumiInfoMapPerLS[thisLS] = new LumiInfo();
   totalLumiByBX_ = inLumiOb.getInstLumiAllBX();
   events_ = inLumiOb.getErrorLumiAllBX();
   lumiInfoMapPerLS[thisLS]->setInstLumiAllBX(totalLumiByBX_);
   lumiInfoMapPerLS[thisLS]->setErrorLumiAllBX(events_);
   lumiSections.push_back(thisLS);
 }

◆ estimateType1Frac()

void CorrPCCProducer::estimateType1Frac	(	std::vector< float >	uncorrPCCPerBX,
		float &	type1Frac
	)

private

Definition at line 209 of file CorrPCCProducer.cc.

References correctionTemplate_, evaluateCorrectionResiduals(), mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, mean_type1_residual, LumiConstants::numBX, and type1Frac.

Referenced by calculateCorrections().

                                                                                          {
   std::vector<float> corrected_tmp_;
   for (size_t i = 0; i < uncorrPCCPerBX.size(); i++) {
     corrected_tmp_.push_back(uncorrPCCPerBX.at(i));
   }
 
   //Apply initial type 1 correction
   for (size_t k = 0; k < LumiConstants::numBX - 1; k++) {
     float bin_k = corrected_tmp_.at(k);
     corrected_tmp_.at(k + 1) = corrected_tmp_.at(k + 1) - type1Frac * bin_k;
   }
 
   //Apply type 2 correction
   for (size_t i = 0; i < LumiConstants::numBX - 1; i++) {
     for (size_t j = i + 1; j < i + LumiConstants::numBX - 1; j++) {
       float bin_i = corrected_tmp_.at(i);
       if (j < LumiConstants::numBX) {
         corrected_tmp_.at(j) = corrected_tmp_.at(j) - bin_i * correctionTemplate_.at(j - i);
       } else {
         corrected_tmp_.at(j - LumiConstants::numBX) =
             corrected_tmp_.at(j - LumiConstants::numBX) - bin_i * correctionTemplate_.at(j - i);
       }
     }
   }
 
   //Apply additional iteration for type 1 correction
   evaluateCorrectionResiduals(corrected_tmp_);
   type1Frac += mean_type1_residual;
 }

◆ evaluateCorrectionResiduals()

void CorrPCCProducer::evaluateCorrectionResiduals ( std::vector< float > corrected_tmp_ )

private

Definition at line 240 of file CorrPCCProducer.cc.

References getMaximum(), histoFile, mean_type1_residual, mean_type1_residual_unc, mean_type2_residual, mean_type2_residual_unc, nTrain, LumiConstants::numBX, and remoteMonitoring_LASER_era2018_cfg::threshold.

Referenced by calculateCorrections(), and estimateType1Frac().

                                                                                  {
   float lumiMax = getMaximum(corrected_tmp_);
   float threshold = lumiMax * 0.2;
 
   mean_type1_residual = 0;
   mean_type2_residual = 0;
   nTrain = 0;
   float lumi = 0;
   std::vector<float> afterGlow;
   TH1F type1("type1", "", 1000, -0.5, 0.5);
   TH1F type2("type2", "", 1000, -0.5, 0.5);
   for (size_t ibx = 2; ibx < LumiConstants::numBX - 5; ibx++) {
     lumi = corrected_tmp_.at(ibx);
     afterGlow.clear();
     afterGlow.push_back(corrected_tmp_.at(ibx + 1));
     afterGlow.push_back(corrected_tmp_.at(ibx + 2));
 
     //Where type 1 and type 2 residuals are computed
     if (lumi > threshold && afterGlow[0] < threshold && afterGlow[1] < threshold) {
       for (int index = 3; index < 6; index++) {
         float thisAfterGlow = corrected_tmp_.at(ibx + index);
         if (thisAfterGlow < threshold) {
           afterGlow.push_back(thisAfterGlow);
         } else {
           break;
         }
       }
       float thisType1 = 0;
       float thisType2 = 0;
       if (afterGlow.size() > 1) {
         int nAfter = 0;
         for (unsigned int index = 1; index < afterGlow.size(); index++) {
           thisType2 += afterGlow[index];
           type2.Fill(afterGlow[index] / lumi);
           nAfter++;
         }
         thisType2 /= nAfter;
       }
       thisType1 = (afterGlow[0] - thisType2) / lumi;
 
       type1.Fill(thisType1);
       nTrain += 1;
     }
   }
 
   mean_type1_residual = type1.GetMean();  //Calculate the mean value of the type 1 residual
   mean_type2_residual = type2.GetMean();  //Calculate the mean value of the type 2 residual
   mean_type1_residual_unc = type1.GetMeanError();
   mean_type2_residual_unc = type2.GetMeanError();
 
   histoFile->cd();
   type1.Write();
   type2.Write();
 }

◆ getMaximum()

float CorrPCCProducer::getMaximum ( std::vector< float > lumi_vector )

private

Definition at line 197 of file CorrPCCProducer.cc.

References mps_fire::i.

Referenced by calculateCorrections(), and evaluateCorrectionResiduals().

                                                               {
   float max_lumi = 0;
   for (size_t i = 0; i < lumi_vector.size(); i++) {
     if (lumi_vector.at(i) > max_lumi)
       max_lumi = lumi_vector.at(i);
   }
   return max_lumi;
 }

◆ makeCorrectionTemplate()

void CorrPCCProducer::makeCorrectionTemplate ( )

private

Definition at line 189 of file CorrPCCProducer.cc.

References simKBmtfDigis_cfi::bx, correctionTemplate_, JetChargeProducer_cfi::exp, LumiConstants::numBX, type2_a_, and type2_b_.

Referenced by CorrPCCProducer().

                                              {
   for (unsigned int bx = 1; bx < LumiConstants::numBX; bx++) {
     correctionTemplate_.at(bx) = type2_a_ * exp(-(float(bx) - 1) * type2_b_);
   }
 }

◆ resetBlock()

void CorrPCCProducer::resetBlock ( )

private

Definition at line 661 of file CorrPCCProducer.cc.

References simKBmtfDigis_cfi::bx, correctionScaleFactors_, events_, LumiConstants::numBX, totalLumiByBX_, and totalLumiByBXAvg_.

Referenced by CorrPCCProducer(), and dqmEndRunProduce().

                                  {
   for (unsigned int bx = 0; bx < LumiConstants::numBX; bx++) {
     totalLumiByBX_[bx] = 0;
     totalLumiByBXAvg_[bx] = 0;
     events_[bx] = 0;
     correctionScaleFactors_[bx] = 1.0;
   }
 }

Member Data Documentation

◆ approxLumiBlockSize_

unsigned int CorrPCCProducer::approxLumiBlockSize_

private

Definition at line 119 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), and dqmEndRunProduce().

◆ correctionScaleFactors_

std::vector<float> CorrPCCProducer::correctionScaleFactors_

private

Definition at line 81 of file CorrPCCProducer.cc.

Referenced by calculateCorrections(), CorrPCCProducer(), dqmEndRunProduce(), and resetBlock().

◆ correctionTemplate_

std::vector<float> CorrPCCProducer::correctionTemplate_

private

Definition at line 80 of file CorrPCCProducer.cc.

Referenced by calculateCorrections(), CorrPCCProducer(), estimateType1Frac(), and makeCorrectionTemplate().

◆ corrlumiAvg_h

TH1F* CorrPCCProducer::corrlumiAvg_h

private

Definition at line 94 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ countLumi_

unsigned int CorrPCCProducer::countLumi_

private

Definition at line 118 of file CorrPCCProducer.cc.

Referenced by beginLuminosityBlock(), and CorrPCCProducer().

◆ errOnLumiByBX_

std::vector<float> CorrPCCProducer::errOnLumiByBX_

private

Definition at line 76 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and endLuminosityBlock().

◆ events_

std::vector<float> CorrPCCProducer::events_

private

Definition at line 79 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), endLuminosityBlock(), and resetBlock().

◆ histoFile

TFile* CorrPCCProducer::histoFile

private

Definition at line 110 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), endJob(), and evaluateCorrectionResiduals().

◆ hlist

TList* CorrPCCProducer::hlist

private

Definition at line 109 of file CorrPCCProducer.cc.

◆ iBlock

unsigned int CorrPCCProducer::iBlock = 0

private

Definition at line 84 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ lumiAvg_h

TH1F* CorrPCCProducer::lumiAvg_h

private

Definition at line 96 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ lumiInfoCounter

std::map<std::pair<unsigned int, unsigned int>, unsigned int> CorrPCCProducer::lumiInfoCounter

private

Definition at line 92 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ lumiInfoMap

std::map<std::pair<unsigned int, unsigned int>, LumiInfo*> CorrPCCProducer::lumiInfoMap

private

Definition at line 91 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ lumiInfoMapIterator

std::map<std::pair<unsigned int, unsigned int>, LumiInfo*>::iterator CorrPCCProducer::lumiInfoMapIterator

private

Definition at line 89 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ lumiInfoMapPerLS

std::map<unsigned int, LumiInfo*> CorrPCCProducer::lumiInfoMapPerLS

private

Definition at line 87 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and endLuminosityBlock().

◆ lumiInfoToken

edm::EDGetTokenT<LumiInfo> CorrPCCProducer::lumiInfoToken

private

Definition at line 71 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), and endLuminosityBlock().

◆ lumiSections

std::vector<unsigned int> CorrPCCProducer::lumiSections

private

Definition at line 88 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and endLuminosityBlock().

◆ maxLS

unsigned int CorrPCCProducer::maxLS = 3500

private

Definition at line 101 of file CorrPCCProducer.cc.

Referenced by bookHistograms(), and dqmEndRunProduce().

◆ mean_type1_residual

float CorrPCCProducer::mean_type1_residual

private

Definition at line 113 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), estimateType1Frac(), and evaluateCorrectionResiduals().

◆ mean_type1_residual_unc

float CorrPCCProducer::mean_type1_residual_unc

private

Definition at line 115 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and evaluateCorrectionResiduals().

◆ mean_type2_residual

float CorrPCCProducer::mean_type2_residual

private

Definition at line 114 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and evaluateCorrectionResiduals().

◆ mean_type2_residual_unc

float CorrPCCProducer::mean_type2_residual_unc

private

Definition at line 116 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and evaluateCorrectionResiduals().

◆ minimumNumberOfEvents

unsigned int CorrPCCProducer::minimumNumberOfEvents

private

Definition at line 85 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), and endLuminosityBlock().

◆ nTrain

unsigned int CorrPCCProducer::nTrain

private

Definition at line 117 of file CorrPCCProducer.cc.

Referenced by evaluateCorrectionResiduals().

◆ overallCorrection_

float CorrPCCProducer::overallCorrection_

private

Definition at line 82 of file CorrPCCProducer.cc.

Referenced by calculateCorrections(), and dqmEndRunProduce().

◆ pccSrc_

std::string CorrPCCProducer::pccSrc_

private

Definition at line 72 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer().

◆ pedestal

float CorrPCCProducer::pedestal

private

Definition at line 125 of file CorrPCCProducer.cc.

Referenced by calculateCorrections(), and dqmEndRunProduce().

◆ pedestal_unc

float CorrPCCProducer::pedestal_unc

private

Definition at line 126 of file CorrPCCProducer.cc.

Referenced by calculateCorrections(), and dqmEndRunProduce().

◆ pedestalGraph

TGraphErrors* CorrPCCProducer::pedestalGraph

private

Definition at line 127 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), and endJob().

◆ poolDbService

edm::Service<cond::service::PoolDBOutputService> CorrPCCProducer::poolDbService

private

Definition at line 129 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ prodInst_

std::string CorrPCCProducer::prodInst_

private

Definition at line 73 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer().

◆ rawlumiBX_

std::vector<float> CorrPCCProducer::rawlumiBX_

private

Definition at line 75 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ scaleFactorAvg_h

TH1F* CorrPCCProducer::scaleFactorAvg_h

private

Definition at line 95 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ thisLS

unsigned int CorrPCCProducer::thisLS

private

Definition at line 120 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce(), and endLuminosityBlock().

◆ totalLumiByBX_

std::vector<float> CorrPCCProducer::totalLumiByBX_

private

Definition at line 77 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), endLuminosityBlock(), and resetBlock().

◆ totalLumiByBXAvg_

std::vector<float> CorrPCCProducer::totalLumiByBXAvg_

private

Definition at line 78 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), and resetBlock().

◆ type1Frac

float CorrPCCProducer::type1Frac

private

Definition at line 112 of file CorrPCCProducer.cc.

Referenced by calculateCorrections(), dqmEndRunProduce(), and estimateType1Frac().

◆ type1FracGraph

TGraphErrors* CorrPCCProducer::type1FracGraph

private

Definition at line 106 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), and endJob().

◆ type1FracHist

TH1F* CorrPCCProducer::type1FracHist

private

Definition at line 97 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ Type1FracMon

MonitorElement* CorrPCCProducer::Type1FracMon

private

Definition at line 102 of file CorrPCCProducer.cc.

Referenced by bookHistograms(), and dqmEndRunProduce().

◆ type1resGraph

TGraphErrors* CorrPCCProducer::type1resGraph

private

Definition at line 107 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), and endJob().

◆ type1resHist

TH1F* CorrPCCProducer::type1resHist

private

Definition at line 98 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ Type1ResMon

MonitorElement* CorrPCCProducer::Type1ResMon

private

Definition at line 103 of file CorrPCCProducer.cc.

Referenced by bookHistograms(), and dqmEndRunProduce().

◆ type2_a_

double CorrPCCProducer::type2_a_

private

Definition at line 122 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), and makeCorrectionTemplate().

◆ type2_b_

double CorrPCCProducer::type2_b_

private

Definition at line 123 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), and makeCorrectionTemplate().

◆ type2resGraph

TGraphErrors* CorrPCCProducer::type2resGraph

private

Definition at line 108 of file CorrPCCProducer.cc.

Referenced by CorrPCCProducer(), dqmEndRunProduce(), and endJob().

◆ type2resHist

TH1F* CorrPCCProducer::type2resHist

private

Definition at line 99 of file CorrPCCProducer.cc.

Referenced by dqmEndRunProduce().

◆ Type2ResMon

MonitorElement* CorrPCCProducer::Type2ResMon

private

Definition at line 104 of file CorrPCCProducer.cc.

Referenced by bookHistograms(), and dqmEndRunProduce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ CorrPCCProducer()

◆ ~CorrPCCProducer()

Member Function Documentation

◆ beginLuminosityBlock()

◆ bookHistograms()

◆ calculateCorrections()

◆ dqmEndRun()

◆ dqmEndRunProduce()

◆ endJob()

◆ endLuminosityBlock()

◆ estimateType1Frac()

◆ evaluateCorrectionResiduals()

◆ getMaximum()

◆ makeCorrectionTemplate()

◆ resetBlock()

Member Data Documentation

◆ approxLumiBlockSize_

◆ correctionScaleFactors_

◆ correctionTemplate_

◆ corrlumiAvg_h

◆ countLumi_

◆ errOnLumiByBX_

◆ events_

◆ histoFile

◆ hlist

◆ iBlock

◆ lumiAvg_h

◆ lumiInfoCounter

◆ lumiInfoMap

◆ lumiInfoMapIterator

◆ lumiInfoMapPerLS

◆ lumiInfoToken

◆ lumiSections

◆ maxLS

◆ mean_type1_residual

◆ mean_type1_residual_unc

◆ mean_type2_residual

◆ mean_type2_residual_unc

◆ minimumNumberOfEvents

◆ nTrain

◆ overallCorrection_

◆ pccSrc_

◆ pedestal

◆ pedestal_unc

◆ pedestalGraph

◆ poolDbService

◆ prodInst_

◆ rawlumiBX_

◆ scaleFactorAvg_h

◆ thisLS

◆ totalLumiByBX_

◆ totalLumiByBXAvg_

◆ type1Frac

◆ type1FracGraph

◆ type1FracHist

◆ Type1FracMon

◆ type1resGraph

◆ type1resHist

◆ Type1ResMon

◆ type2_a_

◆ type2_b_

◆ type2resGraph

◆ type2resHist

◆ Type2ResMon