#include <DPGAnalysis/SiStripTools/plugins/APVCyclePhaseProducerFromL1ABC.cc>

Inheritance diagram for APVCyclePhaseProducerFromL1ABC:

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

	~APVCyclePhaseProducerFromL1ABC ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

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

virtual	~ProducerBase () noexcept(false)

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

	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)

Private Member Functions
virtual void	beginJob () override

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

virtual void	endJob () override

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

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

Private Attributes
edm::EventNumber_t	_curroffevent

long long	_curroffset

const std::vector< std::string >	_defpartnames

const std::vector< int >	_defphases

const unsigned int	_firstgoodrun

TH1F **	_hbx

TH1F **	_hdbx

TH1F **	_hdorbit

edm::EDGetTokenT< L1AcceptBunchCrossingCollection >	_l1abccollectionToken

std::map< edm::EventNumber_t, long long >	_offsets

const int	_orbitoffsetSOR

const bool	_wantHistos

RunHistogramManager	m_rhm

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	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::EDProducer
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

Description: EDproducer for APVCyclePhaseCollection which uses the configuration file to assign a phase to the run

Implementation: <Notes on="" implementation>="">

Definition at line 54 of file APVCyclePhaseProducerFromL1ABC.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 98 of file APVCyclePhaseProducerFromL1ABC.cc.

References _hbx, _hdbx, _hdorbit, _wantHistos, m_rhm, and RunHistogramManager::makeTH1F().

                                                                                             :
   _l1abccollectionToken(mayConsume<L1AcceptBunchCrossingCollection>(iConfig.getParameter<edm::InputTag>("l1ABCCollection"))),
   _defpartnames(iConfig.getParameter<std::vector<std::string> >("defaultPartitionNames")),
   _defphases(iConfig.getParameter<std::vector<int> >("defaultPhases")),
   _orbitoffsetSOR(iConfig.getParameter<int>("StartOfRunOrbitOffset")),
   _wantHistos(iConfig.getUntrackedParameter<bool>("wantHistos",false)),
   m_rhm(consumesCollector()),
   _hbx(0),_hdbx(0),_hdorbit(0),_firstgoodrun(110878),
   _offsets(), _curroffset(0), _curroffevent(0)
 {
 
   produces<APVCyclePhaseCollection,edm::InEvent>();
 
    //now do what ever other initialization is needed
 
   if(_wantHistos) {
     _hbx = m_rhm.makeTH1F("l1abcbx","BX number from L1ABC",4096,-0.5,4095.5);
     _hdbx = m_rhm.makeTH1F("dbx","BX number difference",4096*2-1,-4095.5,4095.5);
     _hdorbit = m_rhm.makeTH1F("dorbit","Orbit Number difference",9999,-4999.5,4999.5);
   }
 
 
 }

APVCyclePhaseProducerFromL1ABC::~APVCyclePhaseProducerFromL1ABC ( )

Definition at line 123 of file APVCyclePhaseProducerFromL1ABC.cc.

 {
 
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void APVCyclePhaseProducerFromL1ABC::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 287 of file APVCyclePhaseProducerFromL1ABC.cc.

288 {

289 }

void APVCyclePhaseProducerFromL1ABC::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Definition at line 138 of file APVCyclePhaseProducerFromL1ABC.cc.

References _firstgoodrun, _hbx, _hdbx, _hdorbit, _offsets, _wantHistos, RunHistogramManager::beginRun(), m_rhm, and edm::RunBase::run().

 {
 
   // reset offset vector
 
   _offsets.clear();
   edm::LogInfo("AbsoluteBXOffsetReset") << "Absolute BX offset map reset";
 
 
   if(_wantHistos) {
 
     m_rhm.beginRun(iRun);
 
     if(_hbx && *_hbx) {
       (*_hbx)->GetXaxis()->SetTitle("BX");     (*_hbx)->GetYaxis()->SetTitle("Events");
     }
 
     if(_hdbx && *_hdbx) {
       (*_hdbx)->GetXaxis()->SetTitle("#DeltaBX");     (*_hdbx)->GetYaxis()->SetTitle("Events");
     }
 
     if(_hdorbit && *_hdorbit) {
       (*_hdorbit)->GetXaxis()->SetTitle("#Deltaorbit");     (*_hdorbit)->GetYaxis()->SetTitle("Events");
     }
 
   }
 
   if(iRun.run() < _firstgoodrun) {
     edm::LogInfo("UnreliableMissingL1AcceptBunchCrossingCollection") <<
       "In this run L1AcceptBunchCrossingCollection is missing or unreliable: default phases will be used";
   }
 
 }

void APVCyclePhaseProducerFromL1ABC::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDProducer.

Definition at line 293 of file APVCyclePhaseProducerFromL1ABC.cc.

References DEFINE_FWK_MODULE.

293 {

294 }

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

overrideprivatevirtual

Definition at line 174 of file APVCyclePhaseProducerFromL1ABC.cc.

References _offsets, and PFRecoTauDiscriminationByIsolation_cfi::offset.

 {
   // summary of absolute bx offset vector
 
   edm::LogInfo("L1AcceptBunchCrossingAbsoluteBXOffsetSummary") << "Absolute BX offset summary:";
   for(std::map<edm::EventNumber_t, long long>::const_iterator offset=_offsets.begin();offset!=_offsets.end();++offset) {
     edm::LogVerbatim("L1AcceptBunchCrossingAbsoluteBXOffsetSummary") << offset->first << " " << offset->second;
   }
 
 }

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

overrideprivatevirtual

Definition at line 187 of file APVCyclePhaseProducerFromL1ABC.cc.

References _curroffevent, _curroffset, _defpartnames, _defphases, _firstgoodrun, _hbx, _hdbx, _hdorbit, _l1abccollectionToken, _offsets, _orbitoffsetSOR, _wantHistos, edm::EventBase::bunchCrossing(), edm::EventID::event(), Exception, edm::Event::getByToken(), edm::EventBase::id(), eostools::move(), edm::EventBase::orbitNumber(), phases, edm::Event::put(), and edm::Event::run().

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

                                                                                      {
 
   using namespace edm;
 
   std::unique_ptr<APVCyclePhaseCollection> apvphases(new APVCyclePhaseCollection() );
 
 
   const std::vector<int>& phases = _defphases;
   const std::vector<std::string>& partnames = _defpartnames;
 
   // Look for the L1AcceptBunchCrossingCollection
 
   int phasechange = 0;
 
   if(iEvent.run() >= _firstgoodrun ) {
 
     Handle<L1AcceptBunchCrossingCollection > pIn;
     iEvent.getByToken(_l1abccollectionToken,pIn);
 
     // offset computation
 
     long long orbitoffset = _orbitoffsetSOR;
     int bxoffset = 0;
 
     for(L1AcceptBunchCrossingCollection::const_iterator l1abc=pIn->begin();l1abc!=pIn->end();++l1abc) {
       if(l1abc->l1AcceptOffset()==0) {
     if(l1abc->eventType()!=0) {
       orbitoffset = (long long)iEvent.orbitNumber() - (long long)l1abc->orbitNumber() ;
       bxoffset = iEvent.bunchCrossing() - l1abc->bunchCrossing();
 
       if(_wantHistos) {
         if(_hbx && *_hbx) (*_hbx)->Fill(l1abc->bunchCrossing());
         if(_hdbx && *_hdbx) (*_hdbx)->Fill(bxoffset);
         if(_hdorbit && *_hdorbit) (*_hdorbit)->Fill(orbitoffset);
       }
     }
     else {
       edm::LogWarning("L1AcceptBunchCrossingNoType") << "L1AcceptBunchCrossing with no type found: ";
       for(L1AcceptBunchCrossingCollection::const_iterator debu=pIn->begin();debu!=pIn->end();++debu) {
         edm::LogPrint("L1AcceptBunchCrossingNoType") << *debu;
       }
     }
       }
     }
 
     long long absbxoffset = orbitoffset*3564 + bxoffset;
 
     if(orbitoffset != _orbitoffsetSOR) phasechange = (orbitoffset*3564)%70;
 
      if(_offsets.size()==0) {
        _curroffset = absbxoffset;
        _curroffevent = iEvent.id().event();
        _offsets[iEvent.id().event()] = absbxoffset;
      }
      else {
        if(_curroffset != absbxoffset || iEvent.id().event() < _curroffevent ) {
 
      if( _curroffset != absbxoffset) {
        edm::LogInfo("L1AcceptBunchCrossingAbsoluteBXOffsetChanged") << "Absolute BX offset changed from "
                                     << _curroffset << " to "
                                     << absbxoffset << " at orbit "
                                     << iEvent.orbitNumber() << " and BX "
                                     << iEvent.bunchCrossing();
        for(L1AcceptBunchCrossingCollection::const_iterator l1abc=pIn->begin();l1abc!=pIn->end();++l1abc) {
          edm::LogVerbatim("L1AcceptBunchCrossingAbsoluteBXOffsetChanged") << *l1abc;
        }
      }
 
      _curroffset = absbxoffset;
      _curroffevent = iEvent.id().event();
      _offsets[iEvent.id().event()] = absbxoffset;
        }
      }
 
 
   }
 
 
 
   if(phases.size() < partnames.size() ) {
     // throw exception
     throw cms::Exception("InvalidAPVCyclePhases") << " Inconsistent phases/partitions vector sizes: "
                          << phases.size() << " "
                          << partnames.size();
   }
 
   for(unsigned int ipart=0;ipart<partnames.size();++ipart) {
     if(phases[ipart]>=0) {
       apvphases->get()[partnames[ipart]] = (phases[ipart]+phasechange)%70;
 
     }
   }
 
 
   iEvent.put(std::move(apvphases));
 
 }