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

Inheritance diagram for APVCyclePhaseMonitor:

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

	~APVCyclePhaseMonitor () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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

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)

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

std::map< std::string, TH1F * >	_hphases

std::map< std::string, TProfile * >	_hphasevsorbit

std::map< std::string, TH1F ** >	_hselectedphases

std::map< std::string, TH1F ** >	_hselectedphasessize

std::map< std::string, TH1F ** >	_hselectedphasesvector

std::map< std::string, TProfile ** >	_hselectedphasevectorvsorbit

std::map< std::string, TProfile ** >	_hselectedphasevsorbit

unsigned int	_nevents

std::vector< std::string >	_selectedparts

std::vector< std::string >	_selectedvectorparts

const unsigned int	m_LSfrac

const unsigned int	m_maxLS

RunHistogramManager	m_rhm

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

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

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: <one line="" class="" summary>="">

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

Definition at line 50 of file APVCyclePhaseMonitor.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 91 of file APVCyclePhaseMonitor.cc.

References _hselectedphases, _hselectedphasessize, _hselectedphasesvector, _hselectedphasevectorvsorbit, _hselectedphasevsorbit, _selectedparts, _selectedvectorparts, edm::ParameterSet::getParameter(), m_LSfrac, m_maxLS, m_rhm, RunHistogramManager::makeTH1F(), and RunHistogramManager::makeTProfile().

     : _apvphasecollectionToken(
           consumes<APVCyclePhaseCollection>(iConfig.getParameter<edm::InputTag>("apvCyclePhaseCollection"))),
       _selectedparts(
           iConfig.getUntrackedParameter<std::vector<std::string> >("selectedPartitions", std::vector<std::string>())),
       _selectedvectorparts(iConfig.getUntrackedParameter<std::vector<std::string> >("selectedVectorPartitions",
                                                                                     std::vector<std::string>())),
       m_maxLS(iConfig.getUntrackedParameter<unsigned int>("maxLSBeforeRebin", 125)),
       m_LSfrac(iConfig.getUntrackedParameter<unsigned int>("startingLSFraction", 16)),
       m_rhm(consumesCollector()),
       _hphases(),
       _hselectedphases(),
       _hselectedphasesvector(),
       _hselectedphasessize(),
       _hphasevsorbit(),
       _hselectedphasevsorbit(),
       _hselectedphasevectorvsorbit(),
       _nevents(0) {
   //now do what ever initialization is needed
 
   edm::LogInfo("UsedAPVCyclePhaseCollection")
       << " APVCyclePhaseCollection " << iConfig.getParameter<edm::InputTag>("apvCyclePhaseCollection") << " used";
 
   for (std::vector<std::string>::const_iterator part = _selectedparts.begin(); part != _selectedparts.end(); ++part) {
     char hname[300];
 
     sprintf(hname, "selected_phase_%s", part->c_str());
     edm::LogInfo("SelectedTH1FBeingBooked") << "TH1F " << hname << " being booked";
     _hselectedphases[*part] = m_rhm.makeTH1F(hname, hname, 70, -0.5, 69.5);
 
     sprintf(hname, "selected_phasevsorbit_%s", part->c_str());
     edm::LogInfo("SelectedTProfileBeingBooked") << "TProfile " << hname << " being booked";
     _hselectedphasevsorbit[*part] = m_rhm.makeTProfile(hname, hname, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
   }
 
   for (std::vector<std::string>::const_iterator part = _selectedvectorparts.begin(); part != _selectedvectorparts.end();
        ++part) {
     char hname[300];
 
     sprintf(hname, "selected_phase_vector_%s", part->c_str());
     edm::LogInfo("SelectedVectTH1FBeingBooked") << "TH1F " << hname << " being booked";
     _hselectedphasesvector[*part] = m_rhm.makeTH1F(hname, hname, 70, -0.5, 69.5);
 
     sprintf(hname, "selected_phase_size_%s", part->c_str());
     edm::LogInfo("SelectedVectSizeTH1FBeingBooked") << "TH1F " << hname << " being booked";
     _hselectedphasessize[*part] = m_rhm.makeTH1F(hname, hname, 10, -0.5, 9.5);
 
     sprintf(hname, "selected_phasevectorvsorbit_%s", part->c_str());
     edm::LogInfo("SelectedVectTProfileBeingBooked") << "TProfile " << hname << " being booked";
     _hselectedphasevectorvsorbit[*part] = m_rhm.makeTProfile(hname, hname, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
   }
 }

APVCyclePhaseMonitor::~APVCyclePhaseMonitor ( )

override

Definition at line 144 of file APVCyclePhaseMonitor.cc.

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

Member Function Documentation

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

overrideprivate

Definition at line 154 of file APVCyclePhaseMonitor.cc.

References _apvphasecollectionToken, _hphases, _hphasevsorbit, _hselectedphases, _hselectedphasessize, _hselectedphasesvector, _hselectedphasevectorvsorbit, _hselectedphasevsorbit, _nevents, _selectedparts, _selectedvectorparts, BTVHLTOfflineSource_cfi::dirname, HcalObjRepresent::Fill(), APVCyclePhaseCollection::get(), edm::Event::getByToken(), APVCyclePhaseCollection::getPhase(), APVCyclePhaseCollection::getPhases(), m_LSfrac, m_maxLS, TFileService::mkdir(), edm::EventBase::orbitNumber(), HLT_2018_cff::phase, phases, and edm::Event::run().

                                                                                       {
   using namespace edm;
 
   _nevents++;
 
   edm::Handle<APVCyclePhaseCollection> apvphases;
   iEvent.getByToken(_apvphasecollectionToken, apvphases);
 
   // improve the matchin between default and actual partitions
 
   edm::Service<TFileService> tfserv;
 
   for (std::map<std::string, int>::const_iterator phase = apvphases->get().begin(); phase != apvphases->get().end();
        ++phase) {
     if (_hphases.find(phase->first) == _hphases.end()) {
       char dirname[300];
       sprintf(dirname, "run_%d", iEvent.run());
       TFileDirectory subrun = tfserv->mkdir(dirname);
 
       char hname[300];
 
       sprintf(hname, "phase_%s", phase->first.c_str());
       edm::LogInfo("TH1FBeingBooked") << "TH1F " << hname << " being booked";
       _hphases[phase->first] = subrun.make<TH1F>(hname, hname, 70, -0.5, 69.5);
       _hphases[phase->first]->GetXaxis()->SetTitle("BX mod 70");
       _hphases[phase->first]->GetYaxis()->SetTitle("Events");
 
       sprintf(hname, "phasevsorbit_%s", phase->first.c_str());
       edm::LogInfo("TProfileBeingBooked") << "TProfile " << hname << " being booked";
       _hphasevsorbit[phase->first] = subrun.make<TProfile>(hname, hname, m_LSfrac * m_maxLS, 0, m_maxLS * 262144);
       _hphasevsorbit[phase->first]->SetCanExtend(TH1::kXaxis);
       _hphasevsorbit[phase->first]->GetXaxis()->SetTitle("time [orbit#]");
       _hphasevsorbit[phase->first]->GetYaxis()->SetTitle("Phase");
     }
     _hphases[phase->first]->Fill(phase->second);
     _hphasevsorbit[phase->first]->Fill(iEvent.orbitNumber(), phase->second);
   }
 
   // selected partitions
 
   for (std::vector<std::string>::const_iterator part = _selectedparts.begin(); part != _selectedparts.end(); ++part) {
     if (_hselectedphases.find(*part) != _hselectedphases.end() && _hselectedphases[*part] && *_hselectedphases[*part]) {
       (*_hselectedphases[*part])->Fill(apvphases->getPhase(*part));
     }
     if (_hselectedphasevsorbit.find(*part) != _hselectedphasevsorbit.end() && _hselectedphasevsorbit[*part] &&
         *_hselectedphasevsorbit[*part]) {
       (*_hselectedphasevsorbit[*part])->Fill(iEvent.orbitNumber(), apvphases->getPhase(*part));
     }
   }
 
   for (std::vector<std::string>::const_iterator part = _selectedvectorparts.begin(); part != _selectedvectorparts.end();
        ++part) {
     const std::vector<int> phases = apvphases->getPhases(*part);
 
     if (_hselectedphasessize.find(*part) != _hselectedphasessize.end() && _hselectedphasessize[*part] &&
         *_hselectedphasessize[*part]) {
       (*_hselectedphasessize[*part])->Fill(phases.size());
     }
 
     for (std::vector<int>::const_iterator phase = phases.begin(); phase != phases.end(); ++phase) {
       if (_hselectedphasesvector.find(*part) != _hselectedphasesvector.end() && _hselectedphasesvector[*part] &&
           *_hselectedphasesvector[*part]) {
         (*_hselectedphasesvector[*part])->Fill(*phase);
       }
       if (_hselectedphasevectorvsorbit.find(*part) != _hselectedphasevectorvsorbit.end() &&
           _hselectedphasevectorvsorbit[*part] && *_hselectedphasevectorvsorbit[*part]) {
         (*_hselectedphasevectorvsorbit[*part])->Fill(iEvent.orbitNumber(), *phase);
       }
     }
   }
 }

void APVCyclePhaseMonitor::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 278 of file APVCyclePhaseMonitor.cc.

278 {}

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

overrideprivate

Definition at line 226 of file APVCyclePhaseMonitor.cc.

References _hphases, _hphasevsorbit, _hselectedphases, _hselectedphasessize, _hselectedphasesvector, _hselectedphasevectorvsorbit, _hselectedphasevsorbit, RunHistogramManager::beginRun(), compare::hist, and m_rhm.

                                                                             {
   _hphases.clear();
   _hphasevsorbit.clear();
 
   m_rhm.beginRun(iRun);
 
   for (std::map<std::string, TH1F**>::const_iterator hist = _hselectedphases.begin(); hist != _hselectedphases.end();
        ++hist) {
     if (*(hist->second)) {
       (*(hist->second))->GetXaxis()->SetTitle("BX mod 70");
       (*(hist->second))->GetYaxis()->SetTitle("Events");
     }
   }
   for (std::map<std::string, TProfile**>::const_iterator prof = _hselectedphasevsorbit.begin();
        prof != _hselectedphasevsorbit.end();
        ++prof) {
     if (*(prof->second)) {
       (*(prof->second))->SetCanExtend(TH1::kXaxis);
       (*(prof->second))->GetXaxis()->SetTitle("time [orbit#]");
       (*(prof->second))->GetYaxis()->SetTitle("Phase");
     }
   }
   for (std::map<std::string, TH1F**>::const_iterator hist = _hselectedphasesvector.begin();
        hist != _hselectedphasesvector.end();
        ++hist) {
     if (*(hist->second)) {
       (*(hist->second))->GetXaxis()->SetTitle("BX mod 70");
       (*(hist->second))->GetYaxis()->SetTitle("Events");
     }
   }
   for (std::map<std::string, TH1F**>::const_iterator hist = _hselectedphasessize.begin();
        hist != _hselectedphasessize.end();
        ++hist) {
     if (*(hist->second)) {
       (*(hist->second))->GetXaxis()->SetTitle("Number of Phases");
       (*(hist->second))->GetYaxis()->SetTitle("Events");
     }
   }
   for (std::map<std::string, TProfile**>::const_iterator prof = _hselectedphasevectorvsorbit.begin();
        prof != _hselectedphasevectorvsorbit.end();
        ++prof) {
     if (*(prof->second)) {
       (*(prof->second))->SetCanExtend(TH1::kXaxis);
       (*(prof->second))->GetXaxis()->SetTitle("time [orbit#]");
       (*(prof->second))->GetYaxis()->SetTitle("Phase");
     }
   }
 }