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

Inheritance diagram for EventTimeDistribution:

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

	~EventTimeDistribution () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () 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	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

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

void	endJob () override

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

Private Attributes
edm::EDGetTokenT < APVCyclePhaseCollection >	_apvphasecollToken

TH1F **	_bx

TH1F **	_bxincycle

TH2F **	_bxincyclevsbx

TH1F **	_dbx

TH2F **	_dbxvsbx

TH2F **	_dbxvsbxincycle

edm::EDGetTokenT < EventWithHistory >	_historyProductToken

unsigned int	_nevents

TH1F **	_orbit

TH2F **	_orbitvsbxincycle

const std::string	_phasepart

RunHistogramManager	_rhm

const bool	_wantbxincyclevsbx

const bool	_wantdbxvsbx

const bool	_wantdbxvsbxincycle

const bool	_wantorbitvsbxincycle

std::vector< TH1F ** >	m_dbxhistos

std::vector< std::pair < unsigned int, unsigned int > >	m_dbxindices

TH1F **	m_ewhdepth

const bool	m_ewhdepthHisto

const unsigned int	m_LSfrac

const unsigned int	m_maxLS

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
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)

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)

Detailed Description

Description: <one line="" class="" summary>="">

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

Definition at line 52 of file EventTimeDistribution.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 104 of file EventTimeDistribution.cc.

References _bx, _bxincycle, _bxincyclevsbx, _dbx, _dbxvsbx, _dbxvsbxincycle, _orbit, _orbitvsbxincycle, _rhm, _wantbxincyclevsbx, _wantdbxvsbx, _wantdbxvsbxincycle, _wantorbitvsbxincycle, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), LogDebug, m_dbxhistos, m_dbxindices, m_ewhdepth, m_ewhdepthHisto, m_LSfrac, m_maxLS, RunHistogramManager::makeTH1F(), RunHistogramManager::makeTH2F(), and submitPVValidationJobs::params.

     : _historyProductToken(consumes<EventWithHistory>(iConfig.getParameter<edm::InputTag>("historyProduct"))),
       _apvphasecollToken(consumes<APVCyclePhaseCollection>(iConfig.getParameter<edm::InputTag>("apvPhaseCollection"))),
       _phasepart(iConfig.getUntrackedParameter<std::string>("phasePartition", "None")),
       _wantdbxvsbxincycle(iConfig.getUntrackedParameter<bool>("wantDBXvsBXincycle", false)),
       _wantdbxvsbx(iConfig.getUntrackedParameter<bool>("wantDBXvsBX", false)),
       _wantbxincyclevsbx(iConfig.getUntrackedParameter<bool>("wantBXincyclevsBX", false)),
       _wantorbitvsbxincycle(iConfig.getUntrackedParameter<bool>("wantOrbitvsBXincycle", false)),
       _nevents(0),
       m_maxLS(iConfig.getUntrackedParameter<unsigned int>("maxLSBeforeRebin", 100)),
       m_LSfrac(iConfig.getUntrackedParameter<unsigned int>("startingLSFraction", 4)),
       m_ewhdepthHisto(iConfig.getUntrackedParameter<bool>("wantEWHDepthHisto", false)),
       _rhm(consumesCollector()),
       _dbxvsbxincycle(nullptr),
       _dbxvsbx(nullptr),
       _bxincyclevsbx(nullptr),
       _orbitvsbxincycle(nullptr),
       m_ewhdepth(nullptr) {
   //now do what ever initialization is needed
 
   std::vector<edm::ParameterSet> dbxhistoparams(iConfig.getUntrackedParameter<std::vector<edm::ParameterSet> >(
       "dbxHistosParams", std::vector<edm::ParameterSet>()));
 
   for (std::vector<edm::ParameterSet>::const_iterator params = dbxhistoparams.begin(); params != dbxhistoparams.end();
        ++params) {
     m_dbxindices.push_back(std::pair<unsigned int, unsigned int>(params->getParameter<unsigned int>("firstEvent"),
                                                                  params->getParameter<unsigned int>("secondEvent")));
     char hname[300];
     sprintf(hname,
             "dbx_%d_%d",
             params->getParameter<unsigned int>("firstEvent"),
             params->getParameter<unsigned int>("secondEvent"));
     char htitle[300];
     sprintf(htitle,
             "dbx(%d,%d)",
             params->getParameter<unsigned int>("firstEvent"),
             params->getParameter<unsigned int>("secondEvent"));
 
     m_dbxhistos.push_back(_rhm.makeTH1F(hname,
                                         htitle,
                                         params->getParameter<int>("nbins"),
                                         params->getParameter<double>("min"),
                                         params->getParameter<double>("max")));
     LogDebug("DBXHistoPreBooking") << "Booked DBX histo named " << hname << " untitled " << htitle;
   }
 
   _dbx = _rhm.makeTH1F("dbx", "dbx", 1000, -0.5, 999.5);
   _bx = _rhm.makeTH1F("bx", "BX number", 3564, -0.5, 3563.5);
   _bxincycle = _rhm.makeTH1F("bxcycle", "bxcycle", 70, -0.5, 69.5);
   _orbit = _rhm.makeTH1F("orbit", "orbit", m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
   if (_wantdbxvsbxincycle)
     _dbxvsbxincycle = _rhm.makeTH2F("dbxvsbxincycle", "dbxvsbxincycle", 70, -0.5, 69.5, 1000, -0.5, 999.5);
   if (_wantdbxvsbx)
     _dbxvsbx = _rhm.makeTH2F("dbxvsbx", "dbxvsbx", 3564, -0.5, 3563.5, 1000, -0.5, 999.5);
   if (_wantbxincyclevsbx)
     _bxincyclevsbx = _rhm.makeTH2F("bxincyclevsbx", "bxincyclevsbx", 3564, -0.5, 3563.5, 70, -0.5, 69.5);
   if (_wantorbitvsbxincycle)
     _orbitvsbxincycle =
         _rhm.makeTH2F("orbitvsbxincycle", "orbitvsbxincycle", 70, -0.5, 69.5, m_maxLS, 0, m_maxLS * 262144);
   if (m_ewhdepthHisto)
     m_ewhdepth = _rhm.makeTH1F("ewhdepth", "EventWithHistory Depth", 11, -0.5, 10.5);
 
   edm::LogInfo("UsedAPVCyclePhaseCollection")
       << " APVCyclePhaseCollection " << iConfig.getParameter<edm::InputTag>("apvPhaseCollection") << " used";
 }

EventTimeDistribution::~EventTimeDistribution ( )

override

Definition at line 170 of file EventTimeDistribution.cc.

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

Member Function Documentation

void EventTimeDistribution::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 180 of file EventTimeDistribution.cc.

References _apvphasecollToken, _bxincyclevsbx, _dbxvsbx, _dbxvsbxincycle, _historyProductToken, _nevents, _orbitvsbxincycle, _phasepart, edm::EventBase::bunchCrossing(), HcalObjRepresent::Fill(), edm::Event::getByToken(), hcalSimParameters_cfi::he, dqmdumpme::indices, APVCyclePhaseCollection::invalid, LogDebug, m_dbxhistos, m_dbxindices, m_ewhdepth, APVCyclePhaseCollection::multiphase, APVCyclePhaseCollection::nopartition, and edm::EventBase::orbitNumber().

                                                                                        {
   using namespace edm;
 
   _nevents++;
 
   edm::Handle<EventWithHistory> he;
   iEvent.getByToken(_historyProductToken, he);
 
   // improve the matchin between default and actual partitions
 
   (*_dbx)->Fill(he->deltaBX());
   std::vector<std::pair<unsigned int, unsigned int> >::const_iterator indices = m_dbxindices.begin();
   for (std::vector<TH1F**>::const_iterator dbxhist = m_dbxhistos.begin(); dbxhist != m_dbxhistos.end();
        ++dbxhist, ++indices) {
     (*(*dbxhist))->Fill(he->deltaBX(indices->first, indices->second));
   }
 
   (*_bx)->Fill(iEvent.bunchCrossing() % 3564);
   (*_orbit)->Fill(iEvent.orbitNumber());
   if (_dbxvsbx && *_dbxvsbx)
     (*_dbxvsbx)->Fill(iEvent.bunchCrossing() % 3564, he->deltaBX());
   if (m_ewhdepth && *m_ewhdepth)
     (*m_ewhdepth)->Fill(he->depth());
 
   edm::Handle<APVCyclePhaseCollection> apvphase;
   iEvent.getByToken(_apvphasecollToken, apvphase);
 
   long long tbx = he->absoluteBX();
   if (apvphase.isValid() && !apvphase.failedToGet()) {
     const int thephase = apvphase->getPhase(_phasepart);
     if (thephase != APVCyclePhaseCollection::invalid && thephase != APVCyclePhaseCollection::multiphase &&
         thephase != APVCyclePhaseCollection::nopartition) {
       tbx -= thephase;
       (*_bxincycle)->Fill(tbx % 70);
       if (_dbxvsbxincycle && *_dbxvsbxincycle)
         (*_dbxvsbxincycle)->Fill(tbx % 70, he->deltaBX());
       if (_bxincyclevsbx && *_bxincyclevsbx)
         (*_bxincyclevsbx)->Fill(iEvent.bunchCrossing() % 3564, tbx % 70);
       if (_orbitvsbxincycle && *_orbitvsbxincycle)
         (*_orbitvsbxincycle)->Fill(tbx % 70, iEvent.orbitNumber());
 
     } else {
       LogDebug("InvalidPhase") << "Invalid APVCyclePhase value : " << _phasepart << " " << thephase;
     }
   }
 }

void EventTimeDistribution::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 287 of file EventTimeDistribution.cc.

287 {}

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

overrideprivate

Definition at line 227 of file EventTimeDistribution.cc.

References _bx, _bxincycle, _bxincyclevsbx, _dbx, _dbxvsbx, _dbxvsbxincycle, _orbit, _orbitvsbxincycle, _rhm, RunHistogramManager::beginRun(), LogDebug, m_dbxhistos, and m_ewhdepth.

                                                                              {
   _rhm.beginRun(iRun);
 
   if (*_dbx) {
     (*_dbx)->GetXaxis()->SetTitle("#DeltaBX");
   }
 
   LogDebug("NomberOfHistos") << m_dbxhistos.size();
   for (std::vector<TH1F**>::const_iterator dbxhist = m_dbxhistos.begin(); dbxhist != m_dbxhistos.end(); ++dbxhist) {
     LogDebug("HistoPointer") << *dbxhist;
     if (*(*dbxhist)) {
       (*(*dbxhist))->GetXaxis()->SetTitle("#DeltaBX");
     }
   }
   LogDebug("LabelDone") << "all labels set";
 
   if (*_bx) {
     (*_bx)->GetXaxis()->SetTitle("BX");
   }
 
   if (*_bxincycle) {
     (*_bxincycle)->GetXaxis()->SetTitle("Event BX mod(70)");
   }
 
   if (*_orbit) {
     (*_orbit)->SetCanExtend(TH1::kXaxis);
     (*_orbit)->GetXaxis()->SetTitle("time [Orb#]");
   }
 
   LogDebug("StdPlotsDone") << "all labels in std plots set";
 
   if (_dbxvsbxincycle && *_dbxvsbxincycle) {
     (*_dbxvsbxincycle)->GetXaxis()->SetTitle("Event BX mod(70)");
     (*_dbxvsbxincycle)->GetYaxis()->SetTitle("#DeltaBX");
   }
 
   if (_dbxvsbx && *_dbxvsbx) {
     (*_dbxvsbx)->GetXaxis()->SetTitle("BX");
     (*_dbxvsbx)->GetYaxis()->SetTitle("#DeltaBX");
   }
 
   if (_bxincyclevsbx && *_bxincyclevsbx) {
     (*_bxincyclevsbx)->GetXaxis()->SetTitle("BX");
     (*_bxincyclevsbx)->GetYaxis()->SetTitle("Event BX mod(70)");
   }
 
   if (_orbitvsbxincycle && *_orbitvsbxincycle) {
     (*_orbitvsbxincycle)->SetCanExtend(TH1::kYaxis);
     (*_orbitvsbxincycle)->GetXaxis()->SetTitle("Event BX mod(70)");
     (*_orbitvsbxincycle)->GetYaxis()->SetTitle("time [Orb#]");
   }
 
   if (m_ewhdepth && *m_ewhdepth) {
     (*m_ewhdepth)->GetXaxis()->SetTitle("Depth");
   }
 }