#include <DQM/SiStripMonitorHardware/plugins/SiStripCMMonitor.cc>

Inheritance diagram for SiStripCMMonitorPlugin:

Classes
struct	Statistics

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

	~SiStripCMMonitorPlugin () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

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

void	fillMaps (uint32_t aDetId, unsigned short aChInModule, std::pair< uint16_t, uint16_t > aMedians)

void	updateCabling (const edm::EventSetup &eventSetup)

Private Attributes
const SiStripFedCabling *	cabling_

uint32_t	cablingCacheId_

CMHistograms	cmHists_

std::map< unsigned int, Statistics >	CommonModes_

std::map< unsigned int, Statistics >	CommonModesAPV0minusAPV1_

bool	doTkHistoMap_

edm::EventNumber_t	evt_

std::vector< unsigned int >	fedIdVec_

bool	fillAllDetailedHistograms_

bool	fillWithEvtNum_

bool	fillWithLocalEvtNum_

std::string	folderName_

std::pair< uint16_t, uint16_t >	prevMedians_ [FEDNumbering::MAXSiStripFEDID+1][sistrip::FEDCH_PER_FED]

unsigned int	printDebug_

edm::InputTag	rawDataTag_

edm::EDGetTokenT< FEDRawDataCollection >	rawDataToken_

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

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Description: DQM source application to monitor common mode for SiStrip data

Definition at line 62 of file SiStripCMMonitor.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 123 of file SiStripCMMonitor.cc.

References cmHists_, CommonModes_, CommonModesAPV0minusAPV1_, doTkHistoMap_, evt_, sistrip::FEDCH_PER_FED, l1t::stage2::layer2::fedId, fillAllDetailedHistograms_, fillWithEvtNum_, folderName_, LogTrace, FEDNumbering::MAXSiStripFEDID, FEDNumbering::MINSiStripFEDID, prevMedians_, printDebug_, rawDataTag_, and rawDataToken_.

   : rawDataTag_(iConfig.getUntrackedParameter<edm::InputTag>("RawDataTag",edm::InputTag("source",""))),
     folderName_(iConfig.getUntrackedParameter<std::string>("HistogramFolderName","SiStrip/ReadoutView/CMMonitoring")),
     fedIdVec_(iConfig.getUntrackedParameter<std::vector<unsigned int> >("FedIdVec")),
     fillAllDetailedHistograms_(iConfig.getUntrackedParameter<bool>("FillAllDetailedHistograms",false)),
     fillWithEvtNum_(iConfig.getUntrackedParameter<bool>("FillWithEventNumber",false)),
     fillWithLocalEvtNum_(iConfig.getUntrackedParameter<bool>("FillWithLocalEventNumber",false)),
     printDebug_(iConfig.getUntrackedParameter<unsigned int>("PrintDebugMessages",1)),
     cablingCacheId_(0)
     
 {
   rawDataToken_ = consumes<FEDRawDataCollection>(rawDataTag_);
   //print config to debug log
   std::ostringstream debugStream;
   if (printDebug_>1) {
     debugStream << "[SiStripCMMonitorPlugin]Configuration for SiStripCMMonitorPlugin: " << std::endl
                 << "[SiStripCMMonitorPlugin]\tRawDataTag: " << rawDataTag_ << std::endl
                 << "[SiStripCMMonitorPlugin]\tHistogramFolderName: " << folderName_ << std::endl
                 << "[SiStripCMMonitorPlugin]\tFillAllDetailedHistograms? " << (fillAllDetailedHistograms_ ? "yes" : "no") << std::endl
         << "[SiStripCMMonitorPlugin]\tFillWithEventNumber?" << (fillWithEvtNum_ ? "yes" : "no") << std::endl
                 << "[SiStripCMMonitorPlugin]\tPrintDebugMessages? " << (printDebug_ ? "yes" : "no") << std::endl;
   }
     
  std::ostringstream* pDebugStream = (printDebug_>1 ? &debugStream : nullptr);
 
  cmHists_.initialise(iConfig,pDebugStream);
 
  doTkHistoMap_ = cmHists_.tkHistoMapEnabled();
 
  CommonModes_.clear();
  CommonModesAPV0minusAPV1_.clear();
 
  for (unsigned int fedId(FEDNumbering::MINSiStripFEDID); fedId <= FEDNumbering::MAXSiStripFEDID; fedId++){
    for (unsigned int iCh(0); iCh<sistrip::FEDCH_PER_FED; iCh++){
      prevMedians_[fedId][iCh] = std::pair<uint16_t,uint16_t>(0,0);
    }
  }
 
 
  if (printDebug_)
    LogTrace("SiStripMonitorHardware") << debugStream.str();
 
  evt_ = 0;
 
 }

SiStripCMMonitorPlugin::~SiStripCMMonitorPlugin ( )

override

Definition at line 169 of file SiStripCMMonitor.cc.

170 {

171 }

Member Function Documentation

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

overrideprivate

Definition at line 199 of file SiStripCMMonitor.cc.

References edmScanValgrind::buffer, cabling_, cmHists_, sistrip::FEDChannel::cmMedian(), gather_cfg::cout, FEDRawData::data(), FedChannelConnection::detId(), edm::EventID::event(), evt_, sistrip::FEDCH_PER_FED, SiStripFedCabling::fedConnection(), FEDRawDataCollection::FEDData(), l1t::stage2::layer2::fedId, fillMaps(), fillWithEvtNum_, fillWithLocalEvtNum_, edm::EventSetup::get(), edm::Event::getByToken(), edm::EventBase::id(), sistrip::invalid32_, FedChannelConnection::isConnected(), sistrip::FEDChannel::length(), FEDNumbering::MAXSiStripFEDID, FEDNumbering::MINSiStripFEDID, FedChannelConnection::nApvPairs(), edm::EventBase::orbitNumber(), prevMedians_, printDebug_, edm::ESHandle< T >::product(), rawDataToken_, FEDRawData::size(), updateCabling(), and MuonErrorMatrixValues_cff::values.

 {
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
 
   //static bool firstEvent = true;
   //static bool isBeingFilled = false;
   //update cabling
   updateCabling(iSetup);
   
   //get raw data
   edm::Handle<FEDRawDataCollection> rawDataCollectionHandle;
   iEvent.getByToken(rawDataToken_,rawDataCollectionHandle);
   const FEDRawDataCollection& rawDataCollection = *rawDataCollectionHandle;
   
   //FED errors
   FEDErrors lFedErrors;
 
   //loop over siStrip FED IDs
   for (unsigned int fedId = FEDNumbering::MINSiStripFEDID; 
        fedId <= FEDNumbering::MAXSiStripFEDID; 
        fedId++) {//loop over FED IDs
     const FEDRawData& fedData = rawDataCollection.FEDData(fedId);
 
     //create an object to fill all errors
     lFedErrors.initialiseFED(fedId,cabling_,tTopo);
 
     //Do detailed check
     //first check if data exists
     bool lDataExist = lFedErrors.checkDataPresent(fedData);
     if (!lDataExist) {
       continue;
     }
 
     std::auto_ptr<const sistrip::FEDBuffer> buffer;
 
     if (!lFedErrors.fillFatalFEDErrors(fedData,0)) {
       continue;
     }
     else {
       //need to construct full object to go any further
       buffer.reset(new sistrip::FEDBuffer(fedData.data(),fedData.size(),true));
       bool channelLengthsOK = buffer->checkChannelLengthsMatchBufferLength();
       bool channelPacketCodesOK = buffer->checkChannelPacketCodes();
       bool feLengthsOK = buffer->checkFEUnitLengths();
       if ( !channelLengthsOK ||
        !channelPacketCodesOK ||
        !feLengthsOK ) {
     continue;
       }
     }
 
     std::ostringstream infoStream;
 
     
     if (printDebug_ > 1) {
       infoStream << " --- Processing FED #" << fedId << std::endl;
     }
 
 
     std::vector<CMHistograms::CMvalues> values;
 
     for (unsigned int iCh = 0; 
      iCh < sistrip::FEDCH_PER_FED; 
      iCh++) {//loop on channels
 
       const FedChannelConnection & lConnection = cabling_->fedConnection(fedId,iCh);
       bool connected = lConnection.isConnected();
 
       //std::cout << "FedID " << fedId << ", ch " << iCh << ", nAPVPairs " << lConnection.nApvPairs() << " apvPairNumber " << lConnection.apvPairNumber() << std::endl;
 
       if (!connected) {
     continue;
       }
 
       uint32_t lDetId = lConnection.detId();
       unsigned short nChInModule = lConnection.nApvPairs();
 
       if (!lDetId || lDetId == sistrip::invalid32_) continue;
 
       bool lFailUnpackerChannelCheck = !buffer->channelGood(iCh, true) && connected;
 
       if (lFailUnpackerChannelCheck) {
     continue;
       }
       
 
       //short lAPVPair = lConnection.apvPairNumber();
       //short lSubDet = DetId(lDetId).subdetId();
 
 //       if (firstEvent){
 //  infoStream << "Subdet " << lSubDet << ", " ;
 //  if (lSubDet == 3) {
 //    
 //    infoStream << "TIB layer " << tTopo->tibLayer(lDetId)  << ", fedID " << fedId << ", channel " << iCh << std::endl;
 //  }
 //  else if (lSubDet == 4) {
 //    
 //    infoStream << "TID side " << tTopo->tibSide(lDetId)  << " wheel " << tTopo->tibWheel(lDetId) << ", ring " << tTopo->tibRing(lDetId) << ", fedID " << fedId << ", channel " << iCh << std::endl;
 //  }
 //  else if (lSubDet == 5) {
 //    
 //    infoStream << "TOB side " << tTopo->tibRod(lDetId)[0]  << " layer " << tTopo->tibLayer(lDetId) << ", rod " << tTopo->tibRodNumber(lDetId) << ", fedID " << fedId << ", channel " << iCh << std::endl;
 //  }
 //  else if (lSubDet == 6) {
 //    
 //    infoStream << "TEC side " << tTopo->tibSide(lDetId)  << " wheel " << tTopo->tibWheel(lDetId) << ", petal " << tTopo->tibPetalNumber(lDetId) << ", ring " << tTopo->tibRing(lDetId) << ", fedID " << fedId << ", channel " << iCh << std::endl;
 //  }
 //  isBeingFilled=true;
 //       }
 
       std::ostringstream lMode;
       lMode << buffer->readoutMode();
       if (evt_ == 0 && printDebug_ > 1) 
     std::cout << "Readout mode: " << lMode.str() << std::endl;
       
 
       const sistrip::FEDChannel & lChannel = buffer->channel(iCh);
       std::pair<uint16_t,uint16_t> medians = std::pair<uint16_t,uint16_t>(0,0);
 
       if (lMode.str().find("Zero suppressed") != lMode.str().npos && lMode.str().find("lite") == lMode.str().npos) medians = std::pair<uint16_t,uint16_t>(lChannel.cmMedian(0),lChannel.cmMedian(1));
       
       CMHistograms::CMvalues lVal;
       lVal.ChannelID = iCh;
       lVal.Medians = std::pair<uint16_t,uint16_t>(medians.first,medians.second);
       lVal.PreviousMedians = prevMedians_[fedId][iCh];
 
 //       if (medians.second-medians.first > 26){
 //  std::ostringstream info;
 //  if (medians.second-medians.first > 44) info << " --- Second bump: event " << iEvent.id().event() << ", FED/Channel " << fedId << "/" << iCh << ", delta=" << medians.second-medians.first << std::endl;
 //  else info << " --- First bump: event " << iEvent.id().event() << ", FED/Channel " << fedId << "/" << iCh << ", delta=" << medians.second-medians.first << std::endl;
 //  edm::LogVerbatim("SiStripMonitorHardware") << info.str();
 //       }
 
       if (printDebug_ > 1) {
     if (lChannel.length() > 7) {
       infoStream << "Medians for channel #" << iCh << " (length " << lChannel.length() << "): " << medians.first << ", " << medians.second << std::endl;
     }
       }
 
       values.push_back(lVal);
 
       //if (iEvent.id().event() > 1000)
       fillMaps(lDetId,nChInModule,medians);
 
       prevMedians_[fedId][iCh] = std::pair<uint16_t,uint16_t>(medians.first,medians.second);
       
     }//loop on channels
     
     float lTime = 0;
     if (fillWithEvtNum_) {
       // casting from unsigned long long to a float here
       // doing it explicitely
       lTime = static_cast<float>(iEvent.id().event());
     } else {
       if (fillWithLocalEvtNum_) {
         // casting from unsigned long long to a float here
         // doing it explicitely
         lTime = static_cast<float>(evt_);
       } else {
         lTime = iEvent.orbitNumber()/11223.;
       }
     }
 
     cmHists_.fillHistograms(values,lTime,fedId);
 
     //if (printDebug_ > 0 && isBeingFilled && firstEvent) edm::LogVerbatim("SiStripMonitorHardware") << infoStream.str();
  
 
 
   }//loop on FEDs
 
 
   //if (isBeingFilled) 
   //firstEvent = false;
 
   evt_++;
 
 }//analyze method

void SiStripCMMonitorPlugin::bookHistograms	(	DQMStore::IBooker &	,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 179 of file SiStripCMMonitor.cc.

References cmHists_, fillAllDetailedHistograms_, folderName_, edm::EventSetup::get(), edm::ESHandle< T >::product(), and DQMStore::IBooker::setCurrentFolder().

 {
   ibooker.setCurrentFolder(folderName_);
 
   edm::ESHandle<TkDetMap> tkDetMapHandle;
   eSetup.get<TrackerTopologyRcd>().get(tkDetMapHandle);
   const TkDetMap* tkDetMap = tkDetMapHandle.product();
 
   cmHists_.bookTopLevelHistograms(ibooker, tkDetMap);
 
   if (fillAllDetailedHistograms_) cmHists_.bookAllFEDHistograms(ibooker);
 }

void SiStripCMMonitorPlugin::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

overrideprivate

Definition at line 192 of file SiStripCMMonitor.cc.

 {
 
 }

void SiStripCMMonitorPlugin::fillMaps	(	uint32_t	aDetId,
		unsigned short	aChInModule,
		std::pair< uint16_t, uint16_t >	aMedians
	)

private

Definition at line 455 of file SiStripCMMonitor.cc.

References CommonModes_, CommonModesAPV0minusAPV1_, SiStripCMMonitorPlugin::Statistics::Counter, DEFINE_FWK_MODULE, doTkHistoMap_, SiStripCMMonitorPlugin::Statistics::Mean, SiStripCMMonitorPlugin::Statistics::Rms, and edm::second().

Referenced by analyze().

 {
 
   if (doTkHistoMap_){//if TkHistMap is enabled
     std::pair<std::map<unsigned int,Statistics>::iterator,bool> alreadyThere[2];
 
     Statistics lStat;
     lStat.Mean = (aMedians.first+aMedians.second)*1./(2*aChInModule);
     lStat.Rms = (aMedians.first+aMedians.second)*(aMedians.first+aMedians.second)*1./(4*aChInModule);
     lStat.Counter = 1./aChInModule;
 
     alreadyThere[0] = CommonModes_.insert(std::pair<unsigned int,Statistics>(aDetId,lStat));
     if (!alreadyThere[0].second) {
       ((alreadyThere[0].first)->second).Mean += (aMedians.first+aMedians.second)*1./(2*aChInModule);
       ((alreadyThere[0].first)->second).Rms += (aMedians.first+aMedians.second)*(aMedians.first+aMedians.second)*1./(4*aChInModule);
       ((alreadyThere[0].first)->second).Counter += 1./aChInModule;
     }
 
     lStat.Mean = (aMedians.first-aMedians.second)*1./aChInModule;
     lStat.Rms = (aMedians.first-aMedians.second)*(aMedians.first-aMedians.second)*1./aChInModule;
     lStat.Counter = 1./aChInModule;
 
     alreadyThere[1] = CommonModesAPV0minusAPV1_.insert(std::pair<unsigned int,Statistics>(aDetId,lStat));
     if (!alreadyThere[1].second) {
       ((alreadyThere[1].first)->second).Mean += (aMedians.first-aMedians.second)*1./aChInModule;
       ((alreadyThere[1].first)->second).Rms += (aMedians.first-aMedians.second)*(aMedians.first-aMedians.second)*1./aChInModule;
       ((alreadyThere[1].first)->second).Counter += 1./aChInModule;
     }
 
   }
 
 }

void SiStripCMMonitorPlugin::updateCabling ( const edm::EventSetup & eventSetup )

private

Definition at line 443 of file SiStripCMMonitor.cc.

References cabling_, cablingCacheId_, edm::EventSetup::get(), and edm::ESHandle< T >::product().

Referenced by analyze().

 {
   uint32_t currentCacheId = eventSetup.get<SiStripFedCablingRcd>().cacheIdentifier();
   if (cablingCacheId_ != currentCacheId) {
     edm::ESHandle<SiStripFedCabling> cablingHandle;
     eventSetup.get<SiStripFedCablingRcd>().get(cablingHandle);
     cabling_ = cablingHandle.product();
     cablingCacheId_ = currentCacheId;
   }
 }