#include <AlcaBeamSpotProducer.h>

Inheritance diagram for AlcaBeamSpotProducer:

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

	~AlcaBeamSpotProducer ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks >
	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)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

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

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

Private Attributes
int	beginLumiOfBSFit_

int	countEvt_

int	countLumi_

int	endLumiOfBSFit_

int	fitNLumi_

int	ftmprun

int	ftmprun0

int	ftotalevents

double	inputBeamWidth_

std::time_t	refBStime [2]

int	resetFitNLumi_

bool	runallfitters_

bool	runbeamwidthfit_

BeamFitter *	theBeamFitter

bool	write2DB_

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

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: AlcaBeamSpotProducer.h package: Calibration/TkAlCaRecoProducers

author: Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov)

Definition at line 26 of file AlcaBeamSpotProducer.h.

Constructor & Destructor Documentation

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

explicit

class: AlcaBeamSpotProducer.cc package: RecoVertex/BeamSpotProducer

author: Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov) Geng-Yuan Jeng, UC Riverside (Geng-.nosp@m.Yuan.nosp@m..Jeng.nosp@m.@cer.nosp@m.n.ch)

Definition at line 35 of file AlcaBeamSpotProducer.cc.

References beginLumiOfBSFit_, edm::EDConsumerBase::consumesCollector(), countLumi_, endLumiOfBSFit_, fitNLumi_, ftmprun, ftmprun0, ftotalevents, edm::ParameterSet::getParameter(), BeamFitter::resetCutFlow(), resetFitNLumi_, BeamFitter::resetLSRange(), BeamFitter::resetPVFitter(), BeamFitter::resetRefTime(), BeamFitter::resetTrkVector(), runallfitters_, runbeamwidthfit_, theBeamFitter, and write2DB_.

                                                                         {
   // get parameter
   write2DB_        = iConfig.getParameter<edm::ParameterSet>("AlcaBeamSpotProducerParameters").getParameter<bool>("WriteToDB");
   runallfitters_   = iConfig.getParameter<edm::ParameterSet>("AlcaBeamSpotProducerParameters").getParameter<bool>("RunAllFitters");
   fitNLumi_        = iConfig.getParameter<edm::ParameterSet>("AlcaBeamSpotProducerParameters").getUntrackedParameter<int>("fitEveryNLumi",-1);
   resetFitNLumi_   = iConfig.getParameter<edm::ParameterSet>("AlcaBeamSpotProducerParameters").getUntrackedParameter<int>("resetEveryNLumi",-1);
   runbeamwidthfit_ = iConfig.getParameter<edm::ParameterSet>("AlcaBeamSpotProducerParameters").getParameter<bool>("RunBeamWidthFit");
   
   theBeamFitter = new BeamFitter(iConfig, consumesCollector());
   theBeamFitter->resetTrkVector();
   theBeamFitter->resetLSRange();
   theBeamFitter->resetCutFlow();
   theBeamFitter->resetRefTime();
   theBeamFitter->resetPVFitter();
 
   ftotalevents = 0;
   ftmprun0 = ftmprun = -1;
   countLumi_ = 0;
   beginLumiOfBSFit_ = endLumiOfBSFit_ = -1;
   
   produces<reco::BeamSpot, edm::InLumi>("alcaBeamSpot");
 }

AlcaBeamSpotProducer::~AlcaBeamSpotProducer ( )

Definition at line 59 of file AlcaBeamSpotProducer.cc.

References theBeamFitter.

                                            {
   delete theBeamFitter;
 }

Member Function Documentation

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

finaloverrideprivatevirtual

Definition at line 71 of file AlcaBeamSpotProducer.cc.

References beginLumiOfBSFit_, edm::LuminosityBlockBase::beginTime(), countLumi_, ftmprun, ftmprun0, edm::LuminosityBlockBase::luminosityBlock(), refBStime, resetFitNLumi_, edm::LuminosityBlockBase::run(), and edm::Timestamp::value().

                                                                                                              {
   const edm::TimeValue_t fbegintimestamp = lumiSeg.beginTime().value();
   const std::time_t ftmptime = fbegintimestamp >> 32;
 
   if ( countLumi_ == 0 || (resetFitNLumi_ > 0 && countLumi_%resetFitNLumi_ == 0) ) {
     ftmprun0 = lumiSeg.run();
     ftmprun = ftmprun0;
     beginLumiOfBSFit_ = lumiSeg.luminosityBlock();
     refBStime[0] = ftmptime;
   }
     
   countLumi_++;
   
 }

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

finaloverrideprivatevirtual

Definition at line 87 of file AlcaBeamSpotProducer.cc.

87 {

88 }

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

finaloverrideprivatevirtual

Definition at line 91 of file AlcaBeamSpotProducer.cc.

References beginLumiOfBSFit_, countLumi_, endLumiOfBSFit_, edm::LuminosityBlockBase::endTime(), reco::BeamSpot::Fake, fitNLumi_, ftmprun, ftmprun0, BeamFitter::getBeamSpot(), BeamFitter::getBSvector(), BeamFitter::getFitLSRange(), edm::LuminosityBlockBase::luminosityBlock(), edm::LuminosityBlock::put(), refBStime, BeamFitter::resetCutFlow(), resetFitNLumi_, BeamFitter::resetLSRange(), BeamFitter::resetPVFitter(), BeamFitter::resetRefTime(), BeamFitter::resetTrkVector(), query::result, BeamFitter::runPVandTrkFitter(), BeamFitter::setFitLSRange(), BeamFitter::setRefTime(), BeamFitter::setRun(), reco::BeamSpot::setType(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamFitter, and edm::Timestamp::value().

                                                                                                             {
   const edm::TimeValue_t fendtimestamp = lumiSeg.endTime().value();
   const std::time_t fendtime = fendtimestamp >> 32;
   refBStime[1] = fendtime;
     
   endLumiOfBSFit_ = lumiSeg.luminosityBlock();
     
   if ( fitNLumi_ == -1 && resetFitNLumi_ == -1 ) return;
     
   if (fitNLumi_ > 0 && countLumi_%fitNLumi_!=0) return;
 
   theBeamFitter->setFitLSRange(beginLumiOfBSFit_,endLumiOfBSFit_);
   theBeamFitter->setRefTime(refBStime[0],refBStime[1]);
   theBeamFitter->setRun(ftmprun0);
     
   std::pair<int,int> LSRange = theBeamFitter->getFitLSRange();
 
   reco::BeamSpot bs;
   if (theBeamFitter->runPVandTrkFitter()){
     bs = theBeamFitter->getBeamSpot();
     edm::LogInfo("AlcaBeamSpotProducer")
         << "\n RESULTS OF DEFAULT FIT " << std::endl
         << " for runs: " << ftmprun0 << " - " << ftmprun << std::endl
         << " for lumi blocks : " << LSRange.first << " - " << LSRange.second << std::endl
         << " lumi counter # " << countLumi_ << std::endl
         << bs << std::endl
         << "fit done. \n" << std::endl; 
   }
   else { // Fill in empty beam spot if beamfit fails
     bs.setType(reco::BeamSpot::Fake);
     edm::LogInfo("AlcaBeamSpotProducer")
         << "\n Empty Beam spot fit" << std::endl
         << " for runs: " << ftmprun0 << " - " << ftmprun << std::endl
         << " for lumi blocks : " << LSRange.first << " - " << LSRange.second << std::endl
         << " lumi counter # " << countLumi_ << std::endl
         << bs << std::endl
         << "fit failed \n" << std::endl;
   }
 
   std::auto_ptr<reco::BeamSpot> result(new reco::BeamSpot);
   *result = bs;
   lumiSeg.put(result, std::string("alcaBeamSpot"));
     
   if (resetFitNLumi_ > 0 && countLumi_%resetFitNLumi_ == 0) {
     std::vector<BSTrkParameters> theBSvector = theBeamFitter->getBSvector();
     edm::LogInfo("AlcaBeamSpotProducer")
         << "Total number of tracks accumulated = " << theBSvector.size() << std::endl
         << "Reset track collection for beam fit" <<std::endl;
     theBeamFitter->resetTrkVector();
     theBeamFitter->resetLSRange();
     theBeamFitter->resetCutFlow();
     theBeamFitter->resetRefTime();
     theBeamFitter->resetPVFitter();
     countLumi_=0;
   }
 }

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

finaloverrideprivatevirtual

Implements edm::one::EDProducerBase.

Definition at line 64 of file AlcaBeamSpotProducer.cc.

References ftmprun, ftotalevents, edm::EventBase::id(), BeamFitter::readEvent(), edm::EventID::run(), and theBeamFitter.

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

                                                                                {
   ftotalevents++;
   theBeamFitter->readEvent(iEvent);
   ftmprun = iEvent.id().run();
 }