CommonModeAnalyzer Class Reference

#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::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
void	updateDetCabling (const edm::EventSetup &setup)

Private Attributes
edm::EDGetTokenT< APVCyclePhaseCollection >	m_apvphasecollToken
uint32_t	m_cacheIdDet
	DB cache ID used to establish if the cabling has changed during the run. More...
std::vector< TH1D * >	m_cmdist
std::vector< TProfile ** >	m_cmvsbxrun
std::vector< TProfile * >	m_cmvsdbxincycle
std::vector< TProfile ** >	m_cmvsorbitrun
const SiStripDetCabling *	m_detCabling
	The cabling object. More...
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::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 66 of file CommonModeAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 120 of file CommonModeAnalyzer.cc.

References edm::ParameterSet::getParameter(), 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_cacheIdDet(0),
  m_detCabling(nullptr)
{
   //now do what ever initialization is needed

 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 191 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  )

overrideprivate

Definition at line 209 of file CommonModeAnalyzer.cc.

References edm::EDCollection< T >::begin(), edm::EventBase::bunchCrossing(), EventWithHistory::deltaBXinCycle(), edm::EDCollection< T >::end(), edm::HandleBase::failedToGet(), HcalObjRepresent::Fill(), spr::find(), edm::Event::getByToken(), SiStripDetCabling::getConnections(), APVCyclePhaseCollection::getPhase(), hcalSimParameters_cfi::he, APVCyclePhaseCollection::invalid, edm::HandleBase::isValid(), LogDebug, m_apvphasecollToken, m_cmdist, m_cmvsbxrun, m_cmvsdbxincycle, m_cmvsorbitrun, m_detCabling, 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(), digitizers_cfi::strip, and updateDetCabling().

{
   using namespace edm;

   updateDetCabling( 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::EDAnalyzer.

Definition at line 307 of file CommonModeAnalyzer.cc.

{

}

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

overrideprivate

Definition at line 279 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::EDAnalyzer.

Definition at line 314 of file CommonModeAnalyzer.cc.

References m_nevents.

                           {

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

}

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

overrideprivate

Definition at line 300 of file CommonModeAnalyzer.cc.

{
}

void CommonModeAnalyzer::updateDetCabling ( const edm::EventSetup &  setup )

private

Definition at line 319 of file CommonModeAnalyzer.cc.

References EnergyCorrector::c, DEFINE_FWK_MODULE, edm::EventSetup::get(), m_cacheIdDet, m_detCabling, and edm::ESHandle< T >::product().

Referenced by analyze().

{
  uint32_t cache_id = setup.get<SiStripDetCablingRcd>().cacheIdentifier();//.get( cabling_ );

  if ( m_cacheIdDet != cache_id ) { // If the cache ID has changed since the last update...
      // Update the cabling object
    edm::ESHandle<SiStripDetCabling> c;
    setup.get<SiStripDetCablingRcd>().get( c );
    m_detCabling = c.product();
    m_cacheIdDet = cache_id;
  } // end of new cache ID check
}

Member Data Documentation

edm::EDGetTokenT<APVCyclePhaseCollection> CommonModeAnalyzer::m_apvphasecollToken

private

Definition at line 85 of file CommonModeAnalyzer.cc.

Referenced by analyze().

uint32_t CommonModeAnalyzer::m_cacheIdDet

private

DB cache ID used to establish if the cabling has changed during the run.

Definition at line 103 of file CommonModeAnalyzer.cc.

Referenced by updateDetCabling().

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

private

Definition at line 94 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

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

private

Definition at line 98 of file CommonModeAnalyzer.cc.

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

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

private

Definition at line 97 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

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

private

Definition at line 99 of file CommonModeAnalyzer.cc.

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

const SiStripDetCabling* CommonModeAnalyzer::m_detCabling

private

The cabling object.

Definition at line 104 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and updateDetCabling().

edm::EDGetTokenT<DetIdCollection> CommonModeAnalyzer::m_digibadmodulecollectionToken

private

Definition at line 86 of file CommonModeAnalyzer.cc.

Referenced by analyze().

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

private

Definition at line 83 of file CommonModeAnalyzer.cc.

Referenced by analyze().

edm::EDGetTokenT<EventWithHistory> CommonModeAnalyzer::m_historyProductToken

private

Definition at line 84 of file CommonModeAnalyzer.cc.

Referenced by analyze().

const bool CommonModeAnalyzer::m_ignorebadfedmod

private

Definition at line 88 of file CommonModeAnalyzer.cc.

Referenced by analyze().

const bool CommonModeAnalyzer::m_ignorenotconnected

private

Definition at line 89 of file CommonModeAnalyzer.cc.

Referenced by analyze().

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

private

Definition at line 93 of file CommonModeAnalyzer.cc.

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

private

Definition at line 96 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

int CommonModeAnalyzer::m_nevents

private

Definition at line 90 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and endJob().

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

private

Definition at line 95 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().

const std::string CommonModeAnalyzer::m_phasepart

private

Definition at line 87 of file CommonModeAnalyzer.cc.

Referenced by analyze().

RunHistogramManager CommonModeAnalyzer::m_rhm

private

Definition at line 101 of file CommonModeAnalyzer.cc.

Referenced by beginRun(), and CommonModeAnalyzer().

std::vector<DetIdSelector> CommonModeAnalyzer::m_selections

private

Definition at line 92 of file CommonModeAnalyzer.cc.

Referenced by analyze(), and CommonModeAnalyzer().