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

Inheritance diagram for CommonModeAnalyzer:

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

	~CommonModeAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources, 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

void	updateDetCabling (const SiStripDetCablingRcd &iRcd)

Private Attributes
edm::EDGetTokenT < APVCyclePhaseCollection >	m_apvphasecollToken

std::vector< TH1D * >	m_cmdist

std::vector< TProfile ** >	m_cmvsbxrun

std::vector< TProfile * >	m_cmvsdbxincycle

std::vector< TProfile ** >	m_cmvsorbitrun

const SiStripDetCabling *	m_detCabling = nullptr
	The cabling object. More...

edm::ESGetToken < SiStripDetCabling, SiStripDetCablingRcd >	m_detCablingToken

edm::ESWatcher < SiStripDetCablingRcd >	m_detCablingWatcher

edm::EDGetTokenT< DetIdCollection >	m_digibadmodulecollectionToken

edm::EDGetTokenT < edm::DetSetVector < SiStripRawDigi > >	m_digicollectionToken

edm::EDGetTokenT < EventWithHistory >	m_historyProductToken

const bool	m_ignorebadfedmod

const bool	m_ignorenotconnected

std::vector< std::string >	m_labels

std::vector< TH1D * >	m_napvs

int	m_nevents

std::vector< TH1D * >	m_nmodules

const std::string	m_phasepart

RunHistogramManager	m_rhm

std::vector< DetIdSelector >	m_selections

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 65 of file CommonModeAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 117 of file CommonModeAnalyzer.cc.

References edm::ParameterSet::getParameter(), TFileService::kSharedResource, label, m_cmdist, m_cmvsbxrun, m_cmvsdbxincycle, m_cmvsorbitrun, m_napvs, m_nmodules, m_rhm, m_selections, TFileService::make(), RunHistogramManager::makeTProfile(), corrVsCorr::selection, and AlCaHLTBitMon_QueryRunRegistry::string.

     : m_digicollectionToken(
           consumes<edm::DetSetVector<SiStripRawDigi> >(iConfig.getParameter<edm::InputTag>("digiCollection"))),
       m_historyProductToken(consumes<EventWithHistory>(iConfig.getParameter<edm::InputTag>("historyProduct"))),
       m_apvphasecollToken(consumes<APVCyclePhaseCollection>(iConfig.getParameter<edm::InputTag>("apvPhaseCollection"))),
       m_digibadmodulecollectionToken(
           consumes<DetIdCollection>(iConfig.getParameter<edm::InputTag>("badModuleDigiCollection"))),
       m_phasepart(iConfig.getUntrackedParameter<std::string>("phasePartition", "None")),
       m_ignorebadfedmod(iConfig.getParameter<bool>("ignoreBadFEDMod")),
       m_ignorenotconnected(iConfig.getParameter<bool>("ignoreNotConnected")),
       m_nevents(0),
       m_selections(),
       m_labels(),
       m_cmdist(),
       m_nmodules(),
       m_napvs(),
       m_cmvsdbxincycle(),
       m_cmvsbxrun(),
       m_cmvsorbitrun(),
       m_rhm(consumesCollector()),
       m_detCablingWatcher(this, &CommonModeAnalyzer::updateDetCabling),
       m_detCablingToken(esConsumes()) {
   //now do what ever initialization is needed
   usesResource(TFileService::kSharedResource);
 
   edm::Service<TFileService> tfserv;
 
   std::vector<edm::ParameterSet> selconfigs = iConfig.getParameter<std::vector<edm::ParameterSet> >("selections");
 
   for (std::vector<edm::ParameterSet>::const_iterator selconfig = selconfigs.begin(); selconfig != selconfigs.end();
        ++selconfig) {
     std::string label = selconfig->getParameter<std::string>("label");
     DetIdSelector selection(*selconfig);
     m_selections.push_back(selection);
 
     {
       std::string hname = label + "_CommonMode";
       std::string htitle = label + " Common Mode";
       m_cmdist.push_back(tfserv->make<TH1D>(hname.c_str(), htitle.c_str(), 1024, -0.5, 1024 - 0.5));
       m_cmdist.back()->GetXaxis()->SetTitle("ADC");
     }
     {
       std::string hname = label + "_nmodules";
       std::string htitle = label + " Number of Modules with CM value";
       m_nmodules.push_back(tfserv->make<TH1D>(hname.c_str(), htitle.c_str(), 20000, -0.5, 20000 - 0.5));
       m_nmodules.back()->GetXaxis()->SetTitle("#modules");
     }
     {
       std::string hname = label + "_napvs";
       std::string htitle = label + " Number of APVs with CM value";
       m_napvs.push_back(tfserv->make<TH1D>(hname.c_str(), htitle.c_str(), 2000, -0.5, 80000 - 0.5));
       m_napvs.back()->GetXaxis()->SetTitle("#apvs");
     }
     {
       std::string hname = label + "_CMvsDBXinCycle";
       std::string htitle = label + " Common Mode vs DBX in Cycle";
       m_cmvsdbxincycle.push_back(tfserv->make<TProfile>(hname.c_str(), htitle.c_str(), 1000, -0.5, 1000 - 0.5));
       m_cmvsdbxincycle.back()->GetXaxis()->SetTitle("DBX in cycle");
       m_cmvsdbxincycle.back()->GetYaxis()->SetTitle("CM (ADC counts)");
     }
     {
       std::string hname = label + "_CMvsBX";
       std::string htitle = label + " Common Mode vs BX";
       m_cmvsbxrun.push_back(m_rhm.makeTProfile(hname.c_str(), htitle.c_str(), 3565, -0.5, 3565 - 0.5));
     }
     {
       std::string hname = label + "_CMvsOrbit";
       std::string htitle = label + " Common Mode vs Orbit";
       m_cmvsorbitrun.push_back(m_rhm.makeTProfile(hname.c_str(), htitle.c_str(), 4 * 500, 0, 500 * 262144));
     }
   }
 }

CommonModeAnalyzer::~CommonModeAnalyzer ( )

override

Definition at line 190 of file CommonModeAnalyzer.cc.

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

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 200 of file CommonModeAnalyzer.cc.

References edm::EventBase::bunchCrossing(), edm::ESWatcher< T >::check(), HcalObjRepresent::Fill(), spr::find(), edm::Event::getByToken(), SiStripDetCabling::getConnections(), hcalSimParameters_cfi::he, APVCyclePhaseCollection::invalid, isSelected(), LogDebug, m_apvphasecollToken, m_cmdist, m_cmvsbxrun, m_cmvsdbxincycle, m_cmvsorbitrun, m_detCabling, m_detCablingWatcher, m_digibadmodulecollectionToken, m_digicollectionToken, m_historyProductToken, m_ignorebadfedmod, m_ignorenotconnected, m_napvs, m_nevents, m_nmodules, m_phasepart, m_selections, mod(), APVCyclePhaseCollection::multiphase, APVCyclePhaseCollection::nopartition, edm::EventBase::orbitNumber(), and digitizers_cfi::strip.

                                                                                     {
   using namespace edm;
 
   m_detCablingWatcher.check(iSetup);
 
   m_nevents++;
 
   edm::Handle<EventWithHistory> he;
   iEvent.getByToken(m_historyProductToken, he);
 
   edm::Handle<APVCyclePhaseCollection> apvphase;
   iEvent.getByToken(m_apvphasecollToken, apvphase);
 
   Handle<DetIdCollection> badmodules;
   iEvent.getByToken(m_digibadmodulecollectionToken, badmodules);
 
   int thephase = APVCyclePhaseCollection::invalid;
   if (apvphase.isValid() && !apvphase.failedToGet()) {
     thephase = apvphase->getPhase(m_phasepart);
   }
   bool isphaseok = (thephase != APVCyclePhaseCollection::invalid && thephase != APVCyclePhaseCollection::multiphase &&
                     thephase != APVCyclePhaseCollection::nopartition);
 
   Handle<edm::DetSetVector<SiStripRawDigi> > digis;
   iEvent.getByToken(m_digicollectionToken, digis);
 
   // loop on detector with digis
 
   std::vector<int> nmodules(m_selections.size(), 0);
   std::vector<int> napvs(m_selections.size(), 0);
 
   for (edm::DetSetVector<SiStripRawDigi>::const_iterator mod = digis->begin(); mod != digis->end(); mod++) {
     std::vector<const FedChannelConnection*> conns = m_detCabling->getConnections(mod->detId());
 
     if (!m_ignorebadfedmod || std::find(badmodules->begin(), badmodules->end(), mod->detId()) == badmodules->end()) {
       for (unsigned int isel = 0; isel < m_selections.size(); ++isel) {
         if (m_selections[isel].isSelected(mod->detId())) {
           unsigned int strip = 0;
           ++nmodules[isel];
           for (edm::DetSet<SiStripRawDigi>::const_iterator digi = mod->begin(); digi != mod->end(); digi++, strip++) {
             LogDebug("StripNumber") << "Strip number " << strip;
             if (!m_ignorenotconnected ||
                 ((conns.size() > strip / 2) && conns[strip / 2] && conns[strip / 2]->isConnected())) {
               ++napvs[isel];
               m_cmdist[isel]->Fill(digi->adc());
               if (isphaseok)
                 m_cmvsdbxincycle[isel]->Fill(he->deltaBXinCycle(thephase), digi->adc());
               if (m_cmvsbxrun[isel] && *(m_cmvsbxrun[isel]))
                 (*(m_cmvsbxrun[isel]))->Fill(iEvent.bunchCrossing(), digi->adc());
               if (m_cmvsorbitrun[isel] && *(m_cmvsorbitrun[isel]))
                 (*(m_cmvsorbitrun[isel]))->Fill(iEvent.orbitNumber(), digi->adc());
             } else if (digi->adc() > 0) {
               edm::LogWarning("NonZeroCMWhenDisconnected")
                   << " Non zero CM in " << mod->detId() << " APV " << strip << " with " << conns.size()
                   << " connections and connection pointer" << conns[strip / 2];
             }
           }
         }
       }
     }
   }
   for (unsigned int isel = 0; isel < m_selections.size(); ++isel) {
     m_nmodules[isel]->Fill(nmodules[isel]);
     m_napvs[isel]->Fill(napvs[isel]);
   }
 }

void CommonModeAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 289 of file CommonModeAnalyzer.cc.

289 {}

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

overrideprivate

Definition at line 267 of file CommonModeAnalyzer.cc.

References RunHistogramManager::beginRun(), m_cmvsbxrun, m_cmvsorbitrun, and m_rhm.

                                                                           {
   m_rhm.beginRun(iRun);
 
   for (std::vector<TProfile**>::const_iterator cmvsbx = m_cmvsbxrun.begin(); cmvsbx != m_cmvsbxrun.end(); ++cmvsbx) {
     if (*cmvsbx && *(*cmvsbx)) {
       (*(*cmvsbx))->GetXaxis()->SetTitle("BX");
       (*(*cmvsbx))->GetYaxis()->SetTitle("CM (ADC counts)");
     }
   }
   for (std::vector<TProfile**>::const_iterator cmvsorbit = m_cmvsorbitrun.begin(); cmvsorbit != m_cmvsorbitrun.end();
        ++cmvsorbit) {
     if (*cmvsorbit && *(*cmvsorbit)) {
       (*(*cmvsorbit))->GetXaxis()->SetTitle("orbit");
       (*(*cmvsorbit))->GetYaxis()->SetTitle("CM (ADC counts)");
       (*(*cmvsorbit))->SetCanExtend(TH1::kXaxis);
     }
   }
 }

void CommonModeAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 292 of file CommonModeAnalyzer.cc.

References m_nevents.

292 { edm::LogInfo("EndOfJob") << m_nevents << " analyzed events"; }

CommonModeAnalyzer::m_nevents

int m_nevents

Definition: CommonModeAnalyzer.cc:88

edm::LogInfo

Log< level::Info, false > LogInfo

Definition: MessageLogger.h:125

void CommonModeAnalyzer::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 286 of file CommonModeAnalyzer.cc.

286 {}

void CommonModeAnalyzer::updateDetCabling ( const SiStripDetCablingRcd & iRcd )

private

Definition at line 294 of file CommonModeAnalyzer.cc.

References edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), m_detCabling, and m_detCablingToken.

                                                                           {
   m_detCabling = &iRcd.get(m_detCablingToken);
 }

Member Data Documentation

edm::EDGetTokenT<APVCyclePhaseCollection> CommonModeAnalyzer::m_apvphasecollToken

private

Definition at line 83 of file CommonModeAnalyzer.cc.

Referenced by analyze().

std::vector<TH1D*> CommonModeAnalyzer::m_cmdist

private

Definition at line 92 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

std::vector<TProfile**> CommonModeAnalyzer::m_cmvsbxrun

private

Definition at line 96 of file CommonModeAnalyzer.cc.

Referenced by analyze(), beginRun(), and CommonModeAnalyzer().

std::vector<TProfile*> CommonModeAnalyzer::m_cmvsdbxincycle

private

Definition at line 95 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

std::vector<TProfile**> CommonModeAnalyzer::m_cmvsorbitrun

private

Definition at line 97 of file CommonModeAnalyzer.cc.

Referenced by analyze(), beginRun(), and CommonModeAnalyzer().

const SiStripDetCabling* CommonModeAnalyzer::m_detCabling = nullptr

private

The cabling object.

Definition at line 103 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and updateDetCabling().

edm::ESGetToken<SiStripDetCabling, SiStripDetCablingRcd> CommonModeAnalyzer::m_detCablingToken

private

Definition at line 102 of file CommonModeAnalyzer.cc.

Referenced by updateDetCabling().

edm::ESWatcher<SiStripDetCablingRcd> CommonModeAnalyzer::m_detCablingWatcher

private

Definition at line 101 of file CommonModeAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<DetIdCollection> CommonModeAnalyzer::m_digibadmodulecollectionToken

private

Definition at line 84 of file CommonModeAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<edm::DetSetVector<SiStripRawDigi> > CommonModeAnalyzer::m_digicollectionToken

private

Definition at line 81 of file CommonModeAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<EventWithHistory> CommonModeAnalyzer::m_historyProductToken

private

Definition at line 82 of file CommonModeAnalyzer.cc.

Referenced by analyze().

const bool CommonModeAnalyzer::m_ignorebadfedmod

private

Definition at line 86 of file CommonModeAnalyzer.cc.

Referenced by analyze().

const bool CommonModeAnalyzer::m_ignorenotconnected

private

Definition at line 87 of file CommonModeAnalyzer.cc.

Referenced by analyze().

std::vector<std::string> CommonModeAnalyzer::m_labels

private

Definition at line 91 of file CommonModeAnalyzer.cc.

std::vector<TH1D*> CommonModeAnalyzer::m_napvs

private

Definition at line 94 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

int CommonModeAnalyzer::m_nevents

private

Definition at line 88 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and endJob().

std::vector<TH1D*> CommonModeAnalyzer::m_nmodules

private

Definition at line 93 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

const std::string CommonModeAnalyzer::m_phasepart

private

Definition at line 85 of file CommonModeAnalyzer.cc.

Referenced by analyze().

RunHistogramManager CommonModeAnalyzer::m_rhm

private

Definition at line 99 of file CommonModeAnalyzer.cc.

Referenced by beginRun(), and CommonModeAnalyzer().

std::vector<DetIdSelector> CommonModeAnalyzer::m_selections

private

Definition at line 90 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation