Inheritance diagram for SiStripSpyMonitorModule:

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

	~SiStripSpyMonitorModule ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob () override

virtual void	endJob () override

edm::DetSetVector < SiStripRawDigi > ::detset::const_iterator	findTwoConsecutive (const edm::DetSetVector< SiStripRawDigi >::detset &channelDigis, const uint16_t threshold, uint16_t &aCounter)

bool	hasNegativePedSubtr (const edm::DetSetVector< SiStripRawDigi >::detset &channelDigis, uint16_t aPair)

bool	identifyTickmarks (const edm::DetSetVector< SiStripRawDigi >::detset &channelDigis, const uint16_t threshold)

Private Attributes
DQMStore *	dqm_

std::string	dqmStoreFileName_

unsigned int	evt_

bool	fillAllDetailedHistograms_

bool	fillWithEvtNum_

bool	fillWithLocalEvtNum_

uint16_t	firstHeaderBit_

uint16_t	firstTrailerBit_

std::string	folderName_

sistrip::SpyUtilities::FrameQuality	frameQuality_

SPYHistograms	histManager_

uint32_t	maxDigiRange_

uint32_t	maxDigitalHigh_

uint32_t	maxDigitalLow_

uint32_t	minDigiRange_

uint32_t	minDigitalHigh_

uint32_t	minDigitalLow_

std::ofstream	outfile_ [20]

std::map< std::string, unsigned int >	outfileMap_

std::vector< std::string >	outfileNames_

edm::InputTag	spyAPVeTag_

edm::EDGetTokenT< std::vector < uint32_t > >	spyAPVeToken_

edm::InputTag	spyL1Tag_

edm::EDGetTokenT< std::vector < uint32_t > >	spyL1Token_

edm::InputTag	spyPedSubtrDigisTag_

edm::EDGetTokenT < edm::DetSetVector < SiStripRawDigi > >	spyPedSubtrDigisToken_

edm::InputTag	spyScopeRawDigisTag_

edm::EDGetTokenT < edm::DetSetVector < SiStripRawDigi > >	spyScopeRawDigisToken_

edm::InputTag	spyTotCountTag_

edm::EDGetTokenT< std::vector < uint32_t > >	spyTotCountToken_

sistrip::SpyUtilities	utility_

bool	writeCabling_

bool	writeDQMStore_

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

Detailed Description

Definition at line 50 of file SiStripSpyMonitorModule.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 140 of file SiStripSpyMonitorModule.cc.

References evt_, frameQuality_, edm::ParameterSet::getUntrackedParameter(), histManager_, i, SPYHistograms::initialise(), sistrip::SpyUtilities::FrameQuality::maxDigiRange, sistrip::SpyUtilities::FrameQuality::maxTickHeight, sistrip::SpyUtilities::FrameQuality::maxZeroLight, sistrip::SpyUtilities::FrameQuality::minDigiRange, sistrip::SpyUtilities::FrameQuality::minTickHeight, sistrip::SpyUtilities::FrameQuality::minZeroLight, dbtoconf::out, outfile_, outfileMap_, outfileNames_, spyAPVeTag_, spyAPVeToken_, spyL1Tag_, spyL1Token_, spyPedSubtrDigisTag_, spyPedSubtrDigisToken_, spyScopeRawDigisTag_, spyScopeRawDigisToken_, spyTotCountTag_, and spyTotCountToken_.

   : spyScopeRawDigisTag_(iConfig.getUntrackedParameter<edm::InputTag>("SpyScopeRawDigisTag",edm::InputTag("SiStripSpyUnpacker","ScopeRawDigis"))),
     spyPedSubtrDigisTag_(iConfig.getUntrackedParameter<edm::InputTag>("SpyPedSubtrDigisTag",edm::InputTag("SiStripFEDEmulator",""))),
     spyL1Tag_(iConfig.getUntrackedParameter<edm::InputTag>("SpyL1Tag",edm::InputTag("SiStripSpyDigiConverter","L1ACount"))),
     spyTotCountTag_(iConfig.getUntrackedParameter<edm::InputTag>("SpyTotalEventCountTag",edm::InputTag("SiStripSpyDigiConverter","TotalEventCount"))),
     spyAPVeTag_(iConfig.getUntrackedParameter<edm::InputTag>("SpyAPVeTag",edm::InputTag("SiStripSpyDigiConverter","APVAddress"))),
     dqm_(0),
     folderName_(iConfig.getUntrackedParameter<std::string>("HistogramFolderName","SiStrip/ReadoutView/SpyMonitoringSummary")),
     fillAllDetailedHistograms_(iConfig.getUntrackedParameter<bool>("FillAllDetailedHistograms",false)),
     fillWithEvtNum_(iConfig.getUntrackedParameter<bool>("FillWithEventNumber",false)),
     fillWithLocalEvtNum_(iConfig.getUntrackedParameter<bool>("FillWithLocalEventNumber",false)),
     writeDQMStore_(iConfig.getUntrackedParameter<bool>("WriteDQMStore",false)),
     dqmStoreFileName_(iConfig.getUntrackedParameter<std::string>("DQMStoreFileName","DQMStore.root")),
     firstHeaderBit_(0),
     firstTrailerBit_(0),
     outfileNames_(iConfig.getUntrackedParameter<std::vector<std::string> >("OutputErrors")),
     writeCabling_(iConfig.getUntrackedParameter<bool>("WriteCabling",false))
 {
 
   spyScopeRawDigisToken_ = consumes<edm::DetSetVector<SiStripRawDigi> >(spyScopeRawDigisTag_);
   spyPedSubtrDigisToken_ = consumes<edm::DetSetVector<SiStripRawDigi> >(spyPedSubtrDigisTag_);
 
   spyL1Token_       = consumes<std::vector<uint32_t> >(spyL1Tag_);
   spyTotCountToken_ = consumes<std::vector<uint32_t> >(spyTotCountTag_);
   spyAPVeToken_     = consumes<std::vector<uint32_t> >(spyAPVeTag_);
 
   evt_ = 0;
   std::ostringstream pDebugStream;
   histManager_.initialise(iConfig,&pDebugStream);
   const unsigned int nFiles = outfileNames_.size();
 
   for (unsigned int i(0); i<nFiles; i++){
     std::ostringstream lName;
     lName << outfileNames_.at(i) <<  ".out";
     if (i<20) outfile_[i].open(lName.str().c_str(),std::ios::out);
     outfileMap_[outfileNames_.at(i)] = i;
   }
 
   frameQuality_.minDigiRange = static_cast<uint16_t>(iConfig.getUntrackedParameter<uint32_t>("MinDigiRange",100));
   frameQuality_.maxDigiRange = static_cast<uint16_t>(iConfig.getUntrackedParameter<uint32_t>("MaxDigiRange",1024));
   frameQuality_.minZeroLight = static_cast<uint16_t>(iConfig.getUntrackedParameter<uint32_t>("MinZeroLight",0));
   frameQuality_.maxZeroLight = static_cast<uint16_t>(iConfig.getUntrackedParameter<uint32_t>("MaxZeroLight",1024));
   frameQuality_.minTickHeight = static_cast<uint16_t>(iConfig.getUntrackedParameter<uint32_t>("MinTickHeight",0));
   frameQuality_.maxTickHeight = static_cast<uint16_t>(iConfig.getUntrackedParameter<uint32_t>("MaxTickHeight",1024));
   
 }

SiStripSpyMonitorModule::~SiStripSpyMonitorModule ( )

Definition at line 187 of file SiStripSpyMonitorModule.cc.

References i, outfile_, outfileMap_, and outfileNames_.

 { 
   const unsigned int nFiles = outfileNames_.size();
   for (unsigned int i(0); i<nFiles; i++){
     outfile_[i].close();
   }
 
   outfileMap_.clear();
   outfileNames_.clear();
 
 }

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 229 of file SiStripSpyMonitorModule.cc.

 {
 
   static bool firstEvent = true;
   //update cabling and pedestals
   const SiStripFedCabling* lCabling = utility_.getCabling( iSetup );
   if (firstEvent && writeCabling_){
     std::ofstream lOutCabling;
     lOutCabling.open("trackerDetId_FEDIdChNum_list.txt",std::ios::out);
     for (uint16_t lFedId = sistrip::FED_ID_MIN; lFedId <= sistrip::FED_ID_MAX; ++lFedId) {//loop on feds
       for (uint16_t lFedChannel = 0; lFedChannel < sistrip::FEDCH_PER_FED; lFedChannel++){//loop on channels
     const FedChannelConnection & lConnection = lCabling->fedConnection(lFedId,lFedChannel);
     if (!lConnection.isConnected()) continue;
     uint32_t lDetId = lConnection.detId();
     lOutCabling << "FED ID = " << lFedId 
             << ", Channel = " << lFedChannel 
             << ",fedkey = " << sistrip::FEDCH_PER_FED*lFedId + lFedChannel
             << ", detId = " << lDetId 
             << std::endl;
       }
     }
     lOutCabling.close();
   }
 
   //For spy data
   //get map of TotalEventCount and L1ID, indexed by fedId, and APVaddress indexed by fedIndex.
   edm::Handle<std::vector<uint32_t> > lSpyL1IDHandle,lSpyTotCountHandle,lSpyAPVeHandle;
   try {
     /*
     iEvent.getByLabel(spyL1Tag_,lSpyL1IDHandle);
     iEvent.getByLabel(spyTotCountTag_,lSpyTotCountHandle);
     iEvent.getByLabel(spyAPVeTag_,lSpyAPVeHandle);
     */
     iEvent.getByToken(spyL1Token_,lSpyL1IDHandle);
     iEvent.getByToken(spyTotCountToken_,lSpyTotCountHandle);
     iEvent.getByToken(spyAPVeToken_,lSpyAPVeHandle);
   }
   catch (const cms::Exception& e) {
     LogError("SiStripSpyMonitorModule") << e.what() ;
     return;
   }
   //const std::map<uint32_t,uint32_t> & lSpyMaxCountMap = *lSpyL1IDHandle;
   //const std::map<uint32_t,uint32_t> & lSpyMinCountMap = *lSpyTotCountHandle;
   const std::vector<uint32_t> & lSpyAPVeVec = *lSpyAPVeHandle;
 
   //retrieve the scope digis
   edm::Handle<edm::DetSetVector<SiStripRawDigi> > digisHandle;
   try {
     //    iEvent.getByLabel(spyScopeRawDigisTag_, digisHandle);
     iEvent.getByToken(spyScopeRawDigisToken_, digisHandle);
   }
   catch (const cms::Exception& e) {
     LogError("SiStripSpyMonitorModule") << e.what() ;
     return;
   }
   const edm::DetSetVector<SiStripRawDigi> * lInputDigis = digisHandle.product();
 
   //retrieve the reordered payload digis
   edm::Handle<edm::DetSetVector<SiStripRawDigi> > payloadHandle;
   try {
     //    iEvent.getByLabel(spyPedSubtrDigisTag_, payloadHandle);
     iEvent.getByToken(spyPedSubtrDigisToken_, payloadHandle);
   }
   catch (const cms::Exception& e) {
     LogError("SiStripSpyMonitorModule") << e.what() ;
     return;
   }
   const edm::DetSetVector<SiStripRawDigi> * lPayloadDigis = payloadHandle.product();
 
 
   //for first event only
   //loop first on channels to calculate majority value of when the first header bit is found.
   //output info message to give the value found
   //should currently be 6 but may vary in the futur
   //then we can check firstTrailerBit is +256+24 after
 
   if (firstEvent){
     sistrip::SpyUtilities::getMajorityHeader(lInputDigis,firstHeaderBit_);
     firstTrailerBit_ = firstHeaderBit_+24+sistrip::STRIPS_PER_FEDCH;
   }
 
   //initialise some counters, filled in histos eventually
   SPYHistograms::ErrorCounters lCounters;
   lCounters.nNoData = 0;
   lCounters.nLowRange = 0;
   lCounters.nHighRange = 0;
   lCounters.nMinZero = 0;
   lCounters.nMaxSat = 0;
   lCounters.nLowPb = 0;
   lCounters.nHighPb = 0;
   lCounters.nOOS = 0;
   lCounters.nOtherPbs = 0;
   lCounters.nAPVError = 0;
   lCounters.nAPVAddressError = 0;
   lCounters.nNegPeds = 0;
 
 
 
   //fill event number for output text files
   const unsigned int nFiles = outfileNames_.size();
   for (unsigned int i(0); i<nFiles; i++){
     outfile_[i] << "**** evt "<< iEvent.id().event() << " ****" << std::endl ;
   }
 
 
   //loop over all FEDs and channels
   
   for (uint16_t lFedId = sistrip::FED_ID_MIN; lFedId <= sistrip::FED_ID_MAX; ++lFedId) {//loop on feds
 
     SPYHistograms::Errors lFEDErrors;
     lFEDErrors.hasNoData = false;
     lFEDErrors.hasLowRange = false; 
     lFEDErrors.hasHighRange = false; 
     lFEDErrors.hasMinZero = false;
     lFEDErrors.hasMaxSat = false;
     lFEDErrors.hasLowPb = false;
     lFEDErrors.hasHighPb = false;
     lFEDErrors.hasOOS = false;
     lFEDErrors.hasOtherPbs = false;
     lFEDErrors.hasErrorBit0 = false;
     lFEDErrors.hasErrorBit1 = false;
     lFEDErrors.hasAPVAddressError0 = false;
     lFEDErrors.hasAPVAddressError1 = false;
     lFEDErrors.hasNegPeds = false;
  
     uint32_t lAPVAddrRef = lSpyAPVeVec.at(lFedId);
 
     for (uint16_t lFedChannel = 0; lFedChannel < sistrip::FEDCH_PER_FED; lFedChannel++){//loop on channels
       
       uint32_t lFedIndex = sistrip::FEDCH_PER_FED*lFedId + lFedChannel;
       
       const FedChannelConnection & lConnection = lCabling->fedConnection(lFedId,lFedChannel);
 
       if (!lConnection.isConnected()) continue;
 
       uint32_t lDetId = lConnection.detId();
       //uint16_t lNPairs = lConnection.nApvPairs();
       uint16_t lPair = lConnection.apvPairNumber();
 
       edm::DetSetVector<SiStripRawDigi>::const_iterator lDigis = lInputDigis->find(lFedIndex);
 
       //pedsubtr digis
       edm::DetSetVector<SiStripRawDigi>::const_iterator lPayload = lPayloadDigis->find(lDetId);
        
       //no digis found, continue.
       if (lDigis == lInputDigis->end()) {
     LogDebug("SiStripSpyMonitorModule") << " -- digis not found in ScopeRawDigis map for FEDID " 
                           << lFedId << " and FED channel " << lFedChannel << std::endl;
     continue;
       }
 
       sistrip::SpyUtilities::Frame lFrame = sistrip::SpyUtilities::extractFrameInfo(*lDigis);
       
       SPYHistograms::Errors lErrors;
       lErrors.hasNoData = false;
       lErrors.hasLowRange = false; 
       lErrors.hasHighRange = false; 
       lErrors.hasMinZero = false;
       lErrors.hasMaxSat = false;
       lErrors.hasLowPb = false;
       lErrors.hasHighPb = false;
       lErrors.hasOOS = false;
       lErrors.hasOtherPbs = false;
       lErrors.hasErrorBit0 = false;
       lErrors.hasErrorBit1 = false;
       lErrors.hasAPVAddressError0 = false;
       lErrors.hasAPVAddressError1 = false;
       lErrors.hasNegPeds = false;
       
       uint16_t lRange = sistrip::SpyUtilities::range(lFrame);
       uint16_t lThreshold = sistrip::SpyUtilities::threshold(lFrame);
 
       if (lRange == 0) {
     lCounters.nNoData++;
     lErrors.hasNoData = true;
     lFEDErrors.hasNoData = true;
     if (outfileMap_.find("NoData") != outfileMap_.end()) outfile_[outfileMap_["NoData"]] << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lFrame.digitalLow == 0 && lRange > 0) {
     lCounters.nMinZero++;
     lErrors.hasMinZero = true;
     lFEDErrors.hasMinZero = true;
     if (outfileMap_.find("MinZero") != outfileMap_.end()) outfile_[outfileMap_["MinZero"]] << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lFrame.digitalHigh >= 1023){
     lCounters.nMaxSat++;
     lErrors.hasMaxSat = true;
     lFEDErrors.hasMaxSat = true;
     if (outfileMap_.find("MaxSat") != outfileMap_.end()) outfile_[outfileMap_["MaxSat"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lRange > 0 && lRange < frameQuality_.minDigiRange) {
     lCounters.nLowRange++;
     lErrors.hasLowRange = true;
     lFEDErrors.hasLowRange = true;
     if (outfileMap_.find("LowRange") != outfileMap_.end()) outfile_[outfileMap_["LowRange"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lRange > frameQuality_.maxDigiRange) {
     lCounters.nHighRange++;
     lErrors.hasHighRange = true;
     lFEDErrors.hasHighRange = true;
     if (outfileMap_.find("HighRange") != outfileMap_.end()) outfile_[outfileMap_["HighRange"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lFrame.digitalLow < frameQuality_.minZeroLight || lFrame.digitalLow > frameQuality_.maxZeroLight) {
     lCounters.nLowPb++;
     lErrors.hasLowPb = true;
     lFEDErrors.hasLowPb = true;
     if (outfileMap_.find("LowPb") != outfileMap_.end()) outfile_[outfileMap_["LowPb"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lFrame.digitalHigh < frameQuality_.minTickHeight || lFrame.digitalHigh > frameQuality_.maxTickHeight){
     lCounters.nHighPb++;
     lErrors.hasHighPb = true;
     lFEDErrors.hasHighPb = true;
     if (outfileMap_.find("HighPb") != outfileMap_.end()) outfile_[outfileMap_["HighPb"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if ( lFrame.firstHeaderBit != firstHeaderBit_ && //header in wrong position
         ( (lFrame.firstHeaderBit != sistrip::SPY_SAMPLES_PER_CHANNEL && //header and 
            lFrame.firstTrailerBit != sistrip::SPY_SAMPLES_PER_CHANNEL && //trailer found
            lFrame.firstTrailerBit-lFrame.firstHeaderBit == 280) ||      //+ right distance between them
           (lFrame.firstHeaderBit != sistrip::SPY_SAMPLES_PER_CHANNEL && // or header found
            lFrame.firstTrailerBit == sistrip::SPY_SAMPLES_PER_CHANNEL && // and trailer not found
            lFrame.firstHeaderBit > 16 ) || // corresponding to back-to-back frame late enough
           (lFrame.firstHeaderBit == sistrip::SPY_SAMPLES_PER_CHANNEL && // or header not found
            identifyTickmarks(*lDigis,lThreshold) ) // but such that tickmark compatible with OOS frame
           )
         ){
     lCounters.nOOS++;
     lErrors.hasOOS = true;
     lFEDErrors.hasOOS = true;
     if (outfileMap_.find("OOS") != outfileMap_.end()) outfile_[outfileMap_["OOS"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if ( !(lFrame.firstHeaderBit == firstHeaderBit_ && lFrame.firstTrailerBit == firstTrailerBit_) ){
     lCounters.nOtherPbs++;
     lErrors.hasOtherPbs = true;
     lFEDErrors.hasOtherPbs = true;
     if (outfileMap_.find("OtherPbs") != outfileMap_.end()) outfile_[outfileMap_["OtherPbs"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
       else if (lFrame.apvErrorBit.first || lFrame.apvErrorBit.second) {
     if (lFrame.apvErrorBit.first) {
       lCounters.nAPVError++;
       lErrors.hasErrorBit0 = true;
       lFEDErrors.hasErrorBit0 = true;
     }
     if (lFrame.apvErrorBit.second) {
       lCounters.nAPVError++;
       lErrors.hasErrorBit1 = true;
       lFEDErrors.hasErrorBit1 = true;
     }
     if (outfileMap_.find("APVError") != outfileMap_.end()) {
       outfile_[outfileMap_["APVError"]]  << lFedId << " " << lFedChannel << " " << lDetId ;
       if (lFrame.apvErrorBit.first) outfile_[outfileMap_["APVError"]] << " APV0" << std::endl;
       if (lFrame.apvErrorBit.second) outfile_[outfileMap_["APVError"]] << " APV1" << std::endl;
     }
       }
       else if ( lFrame.apvAddress.first != lAPVAddrRef || 
         lFrame.apvAddress.second != lAPVAddrRef ) {
     if (lFrame.apvAddress.first != lAPVAddrRef){
       lCounters.nAPVAddressError++;
       lErrors.hasAPVAddressError0 = true;
       lFEDErrors.hasAPVAddressError0 = true;
     }
     if (lFrame.apvAddress.second != lAPVAddrRef){
       lCounters.nAPVAddressError++;
       lErrors.hasAPVAddressError1 = true;
       lFEDErrors.hasAPVAddressError1 = true;
     }
     if (outfileMap_.find("APVAddressError") != outfileMap_.end()){
       outfile_[outfileMap_["APVAddressError"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       if (lFrame.apvAddress.first != lAPVAddrRef) outfile_[outfileMap_["APVAddressError"]]  << " APV0" << std::endl;
       if (lFrame.apvAddress.second != lAPVAddrRef) outfile_[outfileMap_["APVAddressError"]]  << " APV1" << std::endl;
     }
       }
       else if (lPayload != lPayloadDigis->end() && hasNegativePedSubtr(*lPayload,lPair)){
     lCounters.nNegPeds++;
     lErrors.hasNegPeds = true;
     lFEDErrors.hasNegPeds = true;
     if (outfileMap_.find("NegPeds") != outfileMap_.end()) outfile_[outfileMap_["NegPeds"]]  << lFedId << " " << lFedChannel << " " << lDetId << std::endl;
       }
 
       histManager_.fillDetailedHistograms(lErrors,lFrame,lFedId,lFedChannel);
 
     }//loop on channels
 
 
     histManager_.fillFEDHistograms(lFEDErrors,lFedId);
 
   }//loop on feds
 
   double lTime;
   //if (fillWithEvtNum_) 
   //lTime = iEvent.id().event();
   //else if (fillWithLocalEvtNum_) lTime = evt_;
   //no orbit number for spy data !!
   //else lTime = iEvent.orbitNumber()/11223.;
   lTime = iEvent.id().event();
   if (fillWithLocalEvtNum_) lTime = evt_;
   
   histManager_.fillCountersHistograms(lCounters,lTime);
  
 
   //used to fill histo vs time with local event number....
   evt_++;
   firstEvent = false;
 
 }//analyze method

void SiStripSpyMonitorModule::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 202 of file SiStripSpyMonitorModule.cc.

References SPYHistograms::bookAllFEDHistograms(), SPYHistograms::bookTopLevelHistograms(), dqm_, evt_, fillAllDetailedHistograms_, firstHeaderBit_, firstTrailerBit_, folderName_, histManager_, DQMStore::pwd(), and DQMStore::setCurrentFolder().

 {
 
   //get DQM store
   dqm_ = &(*edm::Service<DQMStore>());
   dqm_->setCurrentFolder(folderName_);
 
   LogInfo("SiStripSpyMonitorModule") << " Histograms will be written in " 
                      << folderName_ 
                      << ". Current folder is : " 
                      << dqm_->pwd() 
                      << std::endl;
   
   //this propagates dqm_ to the histoclass, must be called !
   histManager_.bookTopLevelHistograms(dqm_);
   
   if (fillAllDetailedHistograms_) histManager_.bookAllFEDHistograms();
 
   evt_ = 0;
   firstHeaderBit_ = 0;
   firstTrailerBit_ = 0;
 }

void SiStripSpyMonitorModule::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 539 of file SiStripSpyMonitorModule.cc.

References dqm_, dqmStoreFileName_, i, outfile_, outfileMap_, outfileNames_, DQMStore::save(), and writeDQMStore_.

 {
 
   LogInfo("SiStripSpyMonitorModule") << "WriteDQMStore ? " << writeDQMStore_ 
                      << " file name = " << dqmStoreFileName_
                      << std::endl;
   if (writeDQMStore_) dqm_->save(dqmStoreFileName_);
 
   const unsigned int nFiles = outfileNames_.size();
   for (unsigned int i(0); i<nFiles; i++){
     outfile_[i].close();
   }
 
   outfileMap_.clear();
   outfileNames_.clear();
 
 }

edm::DetSetVector< SiStripRawDigi >::detset::const_iterator SiStripSpyMonitorModule::findTwoConsecutive	(	const edm::DetSetVector< SiStripRawDigi >::detset &	channelDigis,
		const uint16_t	threshold,
		uint16_t &	aCounter
	)

private

Definition at line 611 of file SiStripSpyMonitorModule.cc.

References edm::DetSet< T >::begin(), edm::DetSet< T >::end(), sistrip::SPY_SAMPLES_PER_CHANNEL, and dtDQMClient_cfg::threshold.

Referenced by identifyTickmarks().

 {
   const edm::DetSetVector<SiStripRawDigi>::detset::const_iterator endChannelDigis = channelDigis.end();
   edm::DetSetVector<SiStripRawDigi>::detset::const_iterator lStart = channelDigis.begin()+aCounter;
 
   bool foundTrailer = false;
   // Loop over digis looking for last two above threshold
   uint8_t aboveThreshold = 0;
 
   for (; lStart != endChannelDigis; ++lStart) {
     if ( lStart->adc() > threshold) {
       aboveThreshold++;
     }
     else {
       aboveThreshold = 0;
     }
     if (aboveThreshold == 2) {foundTrailer = true; break; }
     aCounter++;
   }//end of loop over digis
   
   if (foundTrailer) return lStart;
   else {
     aCounter = sistrip::SPY_SAMPLES_PER_CHANNEL;
     return endChannelDigis;
   }
 }

bool SiStripSpyMonitorModule::hasNegativePedSubtr	(	const edm::DetSetVector< SiStripRawDigi >::detset &	channelDigis,
		uint16_t	aPair
	)

private

Definition at line 559 of file SiStripSpyMonitorModule.cc.

References edm::DetSet< T >::begin(), prof2calltree::count, edm::DetSet< T >::end(), and sistrip::STRIPS_PER_FEDCH.

Referenced by analyze().

 {
   edm::DetSetVector<SiStripRawDigi>::detset::const_iterator iDigi = channelDigis.begin();
   const edm::DetSetVector<SiStripRawDigi>::detset::const_iterator endChannelDigis = channelDigis.end();
 
   uint32_t count = 0;
   for (; iDigi != endChannelDigis; ++iDigi) {
     const uint16_t val = iDigi->adc();
     uint16_t lPair = static_cast<uint16_t>(count/sistrip::STRIPS_PER_FEDCH);
     if (val == 0 && lPair == aPair) return true;
     count++;
   }
 
   return false;
 
 }

bool SiStripSpyMonitorModule::identifyTickmarks	(	const edm::DetSetVector< SiStripRawDigi >::detset &	channelDigis,
		const uint16_t	threshold
	)

private

Definition at line 577 of file SiStripSpyMonitorModule.cc.

References prof2calltree::count, edm::DetSet< T >::detId(), edm::DetSet< T >::end(), evt_, findTwoConsecutive(), LogDebug, and sistrip::SPY_SAMPLES_PER_CHANNEL.

Referenced by analyze().

 {
 
   //start from the end
   uint16_t count = sistrip::SPY_SAMPLES_PER_CHANNEL-3;
   uint16_t lastPos = sistrip::SPY_SAMPLES_PER_CHANNEL;
   uint16_t nTrailers = 0;
   edm::DetSetVector<SiStripRawDigi>::detset::const_iterator iDigi = channelDigis.end();
 
   for ( ; count == 0; count--){
     iDigi = findTwoConsecutive(channelDigis,threshold,count);
     //if found, in different position = 70 before than previous value, go and look 70 before
     if (iDigi != channelDigis.end() && 
     (lastPos == sistrip::SPY_SAMPLES_PER_CHANNEL || count == lastPos+1-70)) {
       nTrailers++;
       lastPos = count-1;
       count -= 70;
     }
     //else keep looking
     else count--;
   }
 
   if (nTrailers > 1) LogDebug("SiStripSpyMonitorModule") << " -- Found "
                              << nTrailers << " trailers every 70 clock cycles for channel " 
                              << channelDigis.detId() 
                              << ", evt " << evt_
                              << std::endl;
   //if only one found, should be < 280 otherwise header should have been found and this method would not be called
   return (nTrailers > 1) || (nTrailers == 1 && lastPos < 280);
 
 }