#include <HcalNZSMonitor.h>

Inheritance diagram for HcalNZSMonitor:

Public Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &s)

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

void	cleanup ()

	HcalNZSMonitor (const edm::ParameterSet &ps)

void	processEvent (const FEDRawDataCollection &rawraw, const edm::TriggerResults &, int bxNum, const edm::TriggerNames &triggerNames)

void	reset ()

void	setup ()

void	unpack (const FEDRawData &raw, const HcalElectronicsMap &emap)

	~HcalNZSMonitor ()

Public Member Functions inherited from HcalBaseDQMonitor
	HcalBaseDQMonitor (const edm::ParameterSet &ps)

	HcalBaseDQMonitor ()

virtual	~HcalBaseDQMonitor ()

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
bool	isUnsuppressed (HcalHTRData &payload)

Private Attributes
edm::InputTag	hltresultsLabel_

MonitorElement *	meBXtriggered_

MonitorElement *	meFEDsizesNZS_

MonitorElement *	meFEDsizeVsLumi_

MonitorElement *	meFullCMSdataSize_

MonitorElement *	meIsUS_

MonitorElement *	meL1evtNumber_

MonitorElement *	meTrigFrac_

std::vector< int >	nAcc

int	nAcc_Total

int	nAndAcc

int	period_

edm::InputTag	rawdataLabel_

std::vector< int >	selFEDs_

edm::EDGetTokenT < FEDRawDataCollection >	tok_raw_

edm::EDGetTokenT < edm::TriggerResults >	tok_res_

std::vector< std::string >	triggers_

uint64_t	UScount [32][15]

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 HcalBaseDQMonitor
virtual void	beginJob ()

virtual void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

virtual void	endJob (void)

virtual void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

virtual void	endRun (const edm::Run &run, const edm::EventSetup &c)

void	getLogicalMap (const edm::EventSetup &c)

bool	IsAllowedCalibType ()

bool	LumiInOrder (int lumisec)

void	SetupEtaPhiHists (EtaPhiHists &hh, std::string Name, std::string Units)

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 HcalBaseDQMonitor
std::vector< int >	AllowedCalibTypes_

int	badChannelStatusMask_

int	currentLS

int	currenttype_

DQMStore *	dbe_

int	debug_

bool	enableCleanup_

bool	eventAllowed_

bool	HBpresent_

bool	HEpresent_

bool	HFpresent_

bool	HOpresent_

int	ievt_

std::map< unsigned int, int >	KnownBadCells_

int	levt_

HcalLogicalMap *	logicalMap_

bool	makeDiagnostics_

MonitorElement *	meIevt_

MonitorElement *	meLevt_

bool	mergeRuns_

MonitorElement *	meTevt_

MonitorElement *	meTevtHist_

bool	needLogicalMap_

int	NLumiBlocks_

bool	Online_

std::string	prefixME_

MonitorElement *	ProblemsCurrentLB

MonitorElement *	ProblemsVsLB

MonitorElement *	ProblemsVsLB_HB

MonitorElement *	ProblemsVsLB_HBHEHF

MonitorElement *	ProblemsVsLB_HE

MonitorElement *	ProblemsVsLB_HF

MonitorElement *	ProblemsVsLB_HO

bool	skipOutOfOrderLS_

std::string	subdir_

int	tevt_

Detailed Description

Definition at line 10 of file HcalNZSMonitor.h.

Constructor & Destructor Documentation

HcalNZSMonitor::HcalNZSMonitor ( const edm::ParameterSet & ps )

Definition at line 11 of file HcalNZSMonitor.cc.

References HcalBaseDQMonitor::AllowedCalibTypes_, HcalBaseDQMonitor::debug_, HcalBaseDQMonitor::enableCleanup_, edm::ParameterSet::getUntrackedParameter(), hltresultsLabel_, HcalBaseDQMonitor::makeDiagnostics_, HcalBaseDQMonitor::mergeRuns_, HcalBaseDQMonitor::NLumiBlocks_, HcalBaseDQMonitor::Online_, period_, HcalBaseDQMonitor::prefixME_, rawdataLabel_, HcalBaseDQMonitor::skipOutOfOrderLS_, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, tok_raw_, tok_res_, and triggers_.

 {
   Online_                = ps.getUntrackedParameter<bool>("online",false);
   mergeRuns_             = ps.getUntrackedParameter<bool>("mergeRuns",false);
   enableCleanup_         = ps.getUntrackedParameter<bool>("enableCleanup",false);
   debug_                 = ps.getUntrackedParameter<int>("debug",0);
   prefixME_              = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
   if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
     prefixME_.append("/");
   subdir_                = ps.getUntrackedParameter<std::string>("TaskFolder","NZSMonitor_Hcal"); 
   if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
     subdir_.append("/");
   subdir_=prefixME_+subdir_;
   AllowedCalibTypes_     = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
   skipOutOfOrderLS_      = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",false);
   NLumiBlocks_           = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
   makeDiagnostics_       = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
 
   triggers_=ps.getUntrackedParameter<std::vector<std::string> >("nzsHLTnames"); //["HLT_HcalPhiSym","HLT_HcalNZS_8E29]
   period_=ps.getUntrackedParameter<int>("NZSeventPeriod",4096); //4096
   rawdataLabel_          = ps.getUntrackedParameter<edm::InputTag>("RawDataLabel");
   hltresultsLabel_       = ps.getUntrackedParameter<edm::InputTag>("HLTResultsLabel");
 
   // register for data access
   tok_raw_ = consumes<FEDRawDataCollection>(rawdataLabel_);
   tok_res_ = consumes<edm::TriggerResults>(hltresultsLabel_);
 
 } 

HcalNZSMonitor::~HcalNZSMonitor ( )

Definition at line 40 of file HcalNZSMonitor.cc.

40 {}

Member Function Documentation

void HcalNZSMonitor::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	s
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 136 of file HcalNZSMonitor.cc.

References HcalBaseDQMonitor::analyze(), edm::EventBase::bunchCrossing(), HcalBaseDQMonitor::debug_, edm::Event::getByToken(), hltresultsLabel_, HcalBaseDQMonitor::IsAllowedCalibType(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), processEvent(), rawdataLabel_, tok_raw_, tok_res_, and edm::Event::triggerNames().

 {
   if (!IsAllowedCalibType()) return;
   if (LumiInOrder(e.luminosityBlock())==false) return;
   
   edm::Handle<FEDRawDataCollection> rawraw;
 
   if (!(e.getByToken(tok_raw_,rawraw)))
     {
       edm::LogWarning("HcalNZSMonitor")<<" raw data with label "<<rawdataLabel_<<" not available";
       return;
     }
 
   edm::Handle<edm::TriggerResults> hltRes;
   if (!(e.getByToken(tok_res_,hltRes)))
     {
       if (debug_>0) edm::LogWarning("HcalNZSMonitor")<<" Could not get HLT results with tag "<<hltresultsLabel_<<std::endl;
       return;
     }
 
   const edm::TriggerNames & triggerNames = e.triggerNames(*hltRes);
   // Collections were found; increment counters
   HcalBaseDQMonitor::analyze(e,s);
 
   processEvent(*rawraw, *hltRes, e.bunchCrossing(), triggerNames);
 
 } // void HcalNZSMonitor::analyze(...)

void HcalNZSMonitor::beginRun	(	const edm::Run &	run,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 63 of file HcalNZSMonitor.cc.

References HcalBaseDQMonitor::beginRun(), gather_cfg::cout, HcalBaseDQMonitor::debug_, HcalBaseDQMonitor::mergeRuns_, reset(), setup(), and HcalBaseDQMonitor::tevt_.

 {
   if (debug_>1) std::cout <<"HcalNZSMonitor::beginRun"<<std::endl;
   HcalBaseDQMonitor::beginRun(run,c);
 
   if (tevt_==0) this->setup(); // set up histograms if they have not been created before
   if (mergeRuns_==false)
     this->reset();
 
   return;
 
 } // void HcalNZSMonitor::beginRun(...)

void HcalNZSMonitor::cleanup ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 53 of file HcalNZSMonitor.cc.

References HcalBaseDQMonitor::dbe_, DQMStore::removeContents(), DQMStore::setCurrentFolder(), and HcalBaseDQMonitor::subdir_.

 {
   if(dbe_)
     {
       dbe_->setCurrentFolder(subdir_);
       dbe_->removeContents();
     }
 }//void HcalNZSMonitor::cleanup()

bool HcalNZSMonitor::isUnsuppressed ( HcalHTRData & payload )

private

void HcalNZSMonitor::processEvent	(	const FEDRawDataCollection &	rawraw,
		const edm::TriggerResults &	trigRes,
		int	bxNum,
		const edm::TriggerNames &	triggerNames
	)

Level-1 event number generated by the TTC system

Definition at line 165 of file HcalNZSMonitor.cc.

References edm::HLTGlobalStatus::accept(), gather_cfg::cout, HcalBaseDQMonitor::currentLS, FEDRawData::data(), HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::debug_, FEDRawDataCollection::FEDData(), MonitorElement::Fill(), HcalDCCHeader::getSpigotData(), HcalDCCHeader::getSpigotPresent(), i, HcalHTRData::isUnsuppressed(), j, gen::k, ConfigFiles::l, FEDNumbering::MAXFEDID, meBXtriggered_, meFEDsizesNZS_, meFEDsizeVsLumi_, meFullCMSdataSize_, meIsUS_, meL1evtNumber_, meTrigFrac_, FEDNumbering::MINTriggerGTPFEDID, nAcc, nAcc_Total, nAndAcc, AlCaHLTBitMon_ParallelJobs::p, period_, selFEDs_, MonitorElement::setBinContent(), FEDRawData::size(), edm::TriggerNames::size(), HcalDCCHeader::SPIGOT_COUNT, edm::TriggerNames::triggerName(), and triggers_.

Referenced by analyze().

 {
   if(!dbe_) 
     { 
       if (debug_>0) std::cout <<"HcalNZSMonitor::processEvent DQMStore not instantiated!!!"<<std::endl;  
       return;
     }
 
   const unsigned int nTrig(triggerNames.size());
  
   std::vector<bool> trigAcc;
   for (unsigned int i=0; i<triggers_.size(); i++) trigAcc.push_back(false);
   
    for (unsigned int k=0; k<nTrig; k++)
      {
        for (unsigned int i=0; i<triggers_.size(); i++)
      {
        if (triggerNames.triggerName(k) == triggers_[i] && trigRes.accept(k)) trigAcc[i]=true;
      }
      }
    bool andAcc=true;
    bool orAcc=false;
    for (unsigned int p=0; p<triggers_.size(); p++)
      {
        if (!trigAcc[p]) andAcc=false;
        if (trigAcc[p]) 
      {
        orAcc=true;
        nAcc[p]++;
      }
      }
    
    if (orAcc)
      {
        nAcc_Total++;
        for (unsigned int l=0; l<triggers_.size(); l++)
      {
        meTrigFrac_->setBinContent(l+1,nAcc[l]/(float)nAcc_Total);
      }
      }
 
    if (andAcc) 
      {
        nAndAcc++;
        meTrigFrac_->setBinContent(triggers_.size()+1,nAndAcc/(float)nAcc_Total);
      }
    
   bool processevent=false;
   if (orAcc) processevent=true;
 
   if (!processevent) return;
 
   meBXtriggered_->Fill(bxNum+0.001,1);
 
   //calculate full HCAL data size:
   size_t hcalSize=0;
   bool hcalIsZS = false;
   for (unsigned int k=0; k<selFEDs_.size(); k++)
     {
       const FEDRawData & fedData = rawraw.FEDData(selFEDs_[k]);
       hcalSize+=fedData.size();
       meFEDsizesNZS_->Fill(selFEDs_[k]+0.001,fedData.size()/1024);
 
       const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(fedData.data());
       if (dccHeader==0) continue;  // protection against bad data -- saw this happen in file /store/streamer/Data/A/000/131/540/Data.00131540.0200.A.storageManager.00.0000.dat; not yet sure why -- Jeff, 22 March 2010; this was due to empty (masked?) HO FEDs 724 and 727 -- Grigory, 25/03/2010 
 
       HcalHTRData htr;
       int nspigot =0; 
       for (int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) 
     {    
       nspigot++;
       if (!dccHeader->getSpigotPresent(spigot)) continue;
       // Load the given decoder with the pointer and length from this spigot.
       dccHeader->getSpigotData(spigot,htr, fedData.size()); 
       if(k != 20 && nspigot !=14 ) 
         {      
           if ( !htr.isUnsuppressed() ) hcalIsZS = true; 
         }
     }
     }
   
   if (hcalIsZS) meIsUS_->Fill(0,1);
   else meIsUS_->Fill(1,1);
 
   meFEDsizeVsLumi_->Fill(currentLS+0.001, hcalSize/1024);
   
   size_t fullSize=0;
   for (int j=0; j<=FEDNumbering::MAXFEDID; ++j )
     {
       const FEDRawData & fedData = rawraw.FEDData(j);
       fullSize+=fedData.size();
     }
 
   meFullCMSdataSize_->Fill(currentLS+0.001,fullSize/1024);
   
   // get Trigger FED-Id
   const FEDRawData& fedData = rawraw.FEDData(FEDNumbering::MINTriggerGTPFEDID) ;
   FEDHeader header(fedData.data()) ;
   
   if (header.lvl1ID()%period_==0) meL1evtNumber_->Fill(1,1);
   else meL1evtNumber_->Fill(0,1);
   return;
 
 } //void HcalNZSMonitor::processEvent(...)

void HcalNZSMonitor::reset ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 42 of file HcalNZSMonitor.cc.

References meBXtriggered_, meFEDsizesNZS_, meFEDsizeVsLumi_, meFullCMSdataSize_, meIsUS_, meL1evtNumber_, meTrigFrac_, and MonitorElement::Reset().

Referenced by beginRun().

 {
   meFEDsizeVsLumi_->Reset();
   meFEDsizesNZS_->Reset();
   meL1evtNumber_->Reset();
   meIsUS_->Reset();
   meBXtriggered_->Reset();
   meTrigFrac_->Reset();
   meFullCMSdataSize_->Reset();
 } // void HcalNZSMonitor::reset()

void HcalNZSMonitor::setup ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 77 of file HcalNZSMonitor.cc.

References DQMStore::book1D(), DQMStore::bookProfile(), gather_cfg::cout, HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::debug_, MonitorElement::getTProfile(), i, gen::k, FEDNumbering::MAXHCALFEDID, meBXtriggered_, meFEDsizesNZS_, meFEDsizeVsLumi_, meFullCMSdataSize_, meIsUS_, meL1evtNumber_, meTrigFrac_, FEDNumbering::MINHCALFEDID, nAcc, nAcc_Total, nAndAcc, HcalBaseDQMonitor::NLumiBlocks_, selFEDs_, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), HcalBaseDQMonitor::setup(), HcalBaseDQMonitor::subdir_, and triggers_.

Referenced by beginRun().

 {
   HcalBaseDQMonitor::setup();
   
   if(debug_>1) std::cout << "<HcalNZSMonitor::setup> About to pushback fedUnpackList_" << std::endl;
 
   selFEDs_.clear();
   for (int i=FEDNumbering::MINHCALFEDID; i<=FEDNumbering::MAXHCALFEDID; i++)
     {
       selFEDs_.push_back(i);
     }
 
   nAcc.clear();
   for (unsigned int i=0; i<triggers_.size(); i++) nAcc.push_back(0);
 
   nAndAcc=0;
   nAcc_Total=0;
   
   if (debug_>1) std::cout <<"<HcalNZSMonitor::setup>  Creating histograms"<<std::endl;
   if (dbe_)
     {
       dbe_->setCurrentFolder(subdir_);
       
       meFEDsizesNZS_=dbe_->bookProfile("FED sizes","FED sizes",32,699.5,731.5,100,-1000.0,12000.0,"");
       meFEDsizesNZS_->setAxisTitle("FED number",1);
       meFEDsizesNZS_->setAxisTitle("average size (KB)",2);
       meFEDsizesNZS_->getTProfile()->SetMarkerStyle(22);
       
       meFEDsizeVsLumi_=dbe_->bookProfile("FED_size_Vs_lumi_block_number",
                      "FED size Vs lumi block number;lumiblock number;average HCAL FED size (kB)",
                      NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,10000,"");
       meFEDsizeVsLumi_->getTProfile()->SetMarkerStyle(22);
 
       meL1evtNumber_=dbe_->book1D("Is_L1_event_number_multiple_of_NZS_period",
                   "Is L1 event number multiple of NZS period",2,0,2);
       meL1evtNumber_->setBinLabel(1, "NO", 1);
       meL1evtNumber_->setBinLabel(2, "YES", 1);
 
       meIsUS_=dbe_->book1D("IsUnsuppressed_bit","IsUnsuppressed bit",2,0,2);
       meIsUS_->setBinLabel(1,"NO",1);
       meIsUS_->setBinLabel(2,"YES",1);
 
       meBXtriggered_=dbe_->book1D("Triggered_BX_number","Triggered BX number",3850,0,3850);
       meBXtriggered_->setAxisTitle("BX number",1);
 
       meTrigFrac_=dbe_->book1D("HLT_accept_fractions","HLT accept fractions",triggers_.size()+1,0,triggers_.size()+1);
       for (unsigned int k=0; k<triggers_.size(); k++) meTrigFrac_->setBinLabel(k+1,triggers_[k].c_str(),1);
       meTrigFrac_->setBinLabel(triggers_.size()+1,"AND",1);
 
       meFullCMSdataSize_=dbe_->bookProfile("full_CMS_datasize",
                        "full CMS data size;lumiblock number;average FEDRawDataCollection size (kB)",
                        NLumiBlocks_,0.5,NLumiBlocks_+0.5,100,0,10000,"");
       meFullCMSdataSize_->getTProfile()->SetMarkerStyle(22);
 
     } // if (dbe_)
   return;
 } // void HcalNZSMonitor::setup()