Inheritance diagram for RawPCCProducer:

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

	~RawPCCProducer () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

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

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)

	~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

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

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

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

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

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

Private Attributes
bool	applyCorr_

std::vector< int >	clustersPerBXInput_

std::vector< int >	clustersPerBXOutput_

std::vector< float >	corrClustersPerBXOutput_

std::vector< float >	correctionScaleFactors_

std::ofstream	csvfile

std::string	csvOutLabel_

std::vector< float >	errorPerBX_

std::vector< int >	events_

std::vector< int >	goodMods_

std::vector< int >	modID_

std::vector< int >	modVeto_

std::unique_ptr< LumiInfo >	outputLumiInfo

std::string	outputProductName_

std::string	pccSource_

edm::EDGetTokenT< reco::PixelClusterCounts >	pccToken

std::string	prodInst_

bool	saveCSVFile_

float	statErrOnLumi_

std::string	takeAverageValue_

float	totalLumi_

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

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

class: RawPCCProducer.cc

description: Creates a LumiInfo object that will contain the luminosity per bunch crossing, along with the total luminosity and the statistical error.

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)

Definition at line 33 of file RawPCCProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 75 of file RawPCCProducer.cc.

References applyCorr_, clustersPerBXOutput_, corrClustersPerBXOutput_, correctionScaleFactors_, csvOutLabel_, edm::ParameterSet::getParameter(), modVeto_, LumiConstants::numBX, outputProductName_, pccSource_, pccToken, prodInst_, saveCSVFile_, AlCaHLTBitMon_QueryRunRegistry::string, and takeAverageValue_.

 {
     pccSource_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getParameter<std::string>("inputPccLabel");
     prodInst_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getParameter<std::string>("ProdInst");
     takeAverageValue_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getUntrackedParameter<std::string>("OutputValue",std::string("Totals")); 
     outputProductName_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getUntrackedParameter<std::string>("outputProductName","alcaLumi");
     modVeto_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getParameter<std::vector<int>>("modVeto");
     applyCorr_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getUntrackedParameter<bool>("ApplyCorrections",false);
     saveCSVFile_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getUntrackedParameter<bool>("saveCSVFile",false);
     csvOutLabel_ = iConfig.getParameter<edm::ParameterSet>("RawPCCProducerParameters").getUntrackedParameter<std::string>("label",std::string("rawPCC.csv"));
 
     edm::InputTag PCCInputTag_(pccSource_, prodInst_);
     clustersPerBXOutput_.resize(LumiConstants::numBX,0);//vector containing clusters per bxid 
     corrClustersPerBXOutput_.resize(LumiConstants::numBX,0);//vector containing clusters per bxid 
     pccToken=consumes<reco::PixelClusterCounts, edm::InLumi>(PCCInputTag_);
     produces<LumiInfo, edm::Transition::EndLuminosityBlock>(outputProductName_);
 
     if(!applyCorr_){
         correctionScaleFactors_.resize(LumiConstants::numBX,1.0);
     }
 }

RawPCCProducer::~RawPCCProducer ( )

override

Definition at line 98 of file RawPCCProducer.cc.

98 {

99 }

Member Function Documentation

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

finalprivate

Definition at line 108 of file RawPCCProducer.cc.

References csvfile, csvOutLabel_, edm::LuminosityBlockBase::luminosityBlock(), outputLumiInfo, edm::LuminosityBlockBase::run(), and saveCSVFile_.

                                                                                                        {
 
     outputLumiInfo = std::make_unique<LumiInfo>(); 
 
     if(saveCSVFile_){
         csvfile.open(csvOutLabel_, std::ios_base::app);
         csvfile<<std::to_string(lumiSeg.run())<<",";
         csvfile<<std::to_string(lumiSeg.luminosityBlock())<<",";
     }
 }

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

finalprivate

Definition at line 120 of file RawPCCProducer.cc.

References applyCorr_, clustersPerBXInput_, clustersPerBXOutput_, corrClustersPerBXOutput_, correctionScaleFactors_, csvfile, errorPerBX_, events_, lumiContext::fill, spr::find(), objects.autophobj::float, edm::EventSetup::get(), edm::LuminosityBlock::getByToken(), LumiCorrections::getCorrectionsBX(), goodMods_, mps_fire::i, createfilelist::int, modID_, modVeto_, LumiConstants::numBX, outputLumiInfo, pccToken, edm::Handle< T >::product(), edm::ESHandle< T >::product(), reco::PixelClusterCounts::readCounts(), reco::PixelClusterCounts::readEvents(), reco::PixelClusterCounts::readModID(), saveCSVFile_, statErrOnLumi_, takeAverageValue_, and totalLumi_.

                                                                                                      {
     totalLumi_=0.0;
     statErrOnLumi_=0.0;
 
     std::fill(clustersPerBXOutput_.begin(),  clustersPerBXOutput_.end(),0);
     std::fill(corrClustersPerBXOutput_.begin(),     corrClustersPerBXOutput_.end(),   0);
 
     goodMods_.clear();
 
     edm::Handle<reco::PixelClusterCounts> pccHandle; 
     lumiSeg.getByToken(pccToken,pccHandle);
 
     const reco::PixelClusterCounts& inputPcc = *(pccHandle.product()); 
 
     modID_ = inputPcc.readModID();
     events_= inputPcc.readEvents();
     clustersPerBXInput_ = inputPcc.readCounts();
 
     //making list of modules to sum over
     for (unsigned int i=0;i<modID_.size();i++){
         if (std::find(modVeto_.begin(),modVeto_.end(), modID_.at(i)) == modVeto_.end()){
             goodMods_.push_back(i);
         }
     }
 
     //summing over good modules only 
     for (int bx=0;bx<int(LumiConstants::numBX);bx++){
         for (unsigned int i=0;i<goodMods_.size();i++){
             clustersPerBXOutput_.at(bx)+=clustersPerBXInput_.at(goodMods_.at(i)*int(LumiConstants::numBX)+bx);
 
         }
     }
 
     if(applyCorr_){
         edm::ESHandle< LumiCorrections > corrHandle;
         iSetup.get<LumiCorrectionsRcd>().get(corrHandle);
         const LumiCorrections *pccCorrections = corrHandle.product();
         correctionScaleFactors_ = pccCorrections->getCorrectionsBX();
     }
 
     for (unsigned int i=0;i<clustersPerBXOutput_.size();i++){
         if (events_.at(i)!=0){
             corrClustersPerBXOutput_[i]=clustersPerBXOutput_[i]*correctionScaleFactors_[i];
         }else{
             corrClustersPerBXOutput_[i]=0.0;
         }
         totalLumi_+=corrClustersPerBXOutput_[i];
         statErrOnLumi_+=float(events_[i]);
     }
 
     errorPerBX_.clear();
     errorPerBX_.assign(events_.begin(),events_.end()); 
 
     if(takeAverageValue_=="Average"){
         unsigned int NActiveBX=0;
         for (int bx=0;bx<int(LumiConstants::numBX);bx++){
             if(events_[bx]>0){
                 NActiveBX++;  
                 // Counting where events are will only work 
                 // for ZeroBias or AlwaysTrue triggers.
                 // Random triggers will get all BXs.
                 corrClustersPerBXOutput_[bx]/=float(events_[bx]);
                 errorPerBX_[bx]=1/TMath::Sqrt(float(events_[bx]));
             }
         }
         if (statErrOnLumi_!=0) {
             totalLumi_=totalLumi_/statErrOnLumi_*float(NActiveBX);
             statErrOnLumi_=1/TMath::Sqrt(statErrOnLumi_)*totalLumi_;
         }
     }
 
 
     outputLumiInfo->setTotalInstLumi(totalLumi_);
     outputLumiInfo->setTotalInstStatError(statErrOnLumi_);
 
     outputLumiInfo->setErrorLumiAllBX(errorPerBX_);
     outputLumiInfo->setInstLumiAllBX(corrClustersPerBXOutput_);
 
     if(saveCSVFile_){
         csvfile<<std::to_string(totalLumi_);
         
         if(totalLumi_>0){
             for(unsigned int bx=0;bx<LumiConstants::numBX;bx++){
                 csvfile<<","<<std::to_string(corrClustersPerBXOutput_[bx]);
             }
             csvfile<<std::endl;   
         } else if (totalLumi_<0) {
             edm::LogInfo("WARNING")<<"WHY IS LUMI NEGATIVE?!?!?!? "<<totalLumi_;
         }
 
         csvfile.close();
     }
 }

void RawPCCProducer::endLuminosityBlockProduce	(	edm::LuminosityBlock &	lumiSeg,
		const edm::EventSetup &	iSetup
	)

finalprivate

Definition at line 215 of file RawPCCProducer.cc.

References DEFINE_FWK_MODULE, eostools::move(), outputLumiInfo, outputProductName_, edm::LuminosityBlock::put(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                       {
     lumiSeg.put(std::move(outputLumiInfo), std::string(outputProductName_)); 
 
 }

void RawPCCProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

finalprivate

Definition at line 102 of file RawPCCProducer.cc.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                          {
 
 }

Member Data Documentation

bool RawPCCProducer::applyCorr_

private

Definition at line 66 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

std::vector<int> RawPCCProducer::clustersPerBXInput_

private

Definition at line 53 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

std::vector<int> RawPCCProducer::clustersPerBXOutput_

private

Definition at line 56 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

std::vector<float> RawPCCProducer::corrClustersPerBXOutput_

private

Definition at line 57 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

std::vector<float> RawPCCProducer::correctionScaleFactors_

private

Definition at line 67 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

std::ofstream RawPCCProducer::csvfile

private

Definition at line 71 of file RawPCCProducer.cc.

Referenced by beginLuminosityBlock(), and endLuminosityBlock().

std::string RawPCCProducer::csvOutLabel_

private

Definition at line 63 of file RawPCCProducer.cc.

Referenced by beginLuminosityBlock(), and RawPCCProducer().

std::vector<float> RawPCCProducer::errorPerBX_

private

Definition at line 58 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

std::vector<int> RawPCCProducer::events_

private

Definition at line 55 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

std::vector<int> RawPCCProducer::goodMods_

private

Definition at line 59 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

std::vector<int> RawPCCProducer::modID_

private

Definition at line 54 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

std::vector<int> RawPCCProducer::modVeto_

private

Definition at line 51 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

std::unique_ptr<LumiInfo> RawPCCProducer::outputLumiInfo

private

Definition at line 69 of file RawPCCProducer.cc.

Referenced by beginLuminosityBlock(), endLuminosityBlock(), and endLuminosityBlockProduce().

std::string RawPCCProducer::outputProductName_

private

Definition at line 52 of file RawPCCProducer.cc.

Referenced by endLuminosityBlockProduce(), and RawPCCProducer().

std::string RawPCCProducer::pccSource_

private

Definition at line 47 of file RawPCCProducer.cc.

Referenced by RawPCCProducer().

edm::EDGetTokenT<reco::PixelClusterCounts> RawPCCProducer::pccToken

private

Definition at line 46 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

std::string RawPCCProducer::prodInst_

private

Definition at line 48 of file RawPCCProducer.cc.

Referenced by RawPCCProducer().

bool RawPCCProducer::saveCSVFile_

private

Definition at line 64 of file RawPCCProducer.cc.

Referenced by beginLuminosityBlock(), endLuminosityBlock(), and RawPCCProducer().

float RawPCCProducer::statErrOnLumi_

private

Definition at line 61 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

std::string RawPCCProducer::takeAverageValue_

private

Definition at line 49 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock(), and RawPCCProducer().

float RawPCCProducer::totalLumi_

private

Definition at line 60 of file RawPCCProducer.cc.

Referenced by endLuminosityBlock().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation