#include <HcalDetDiagTimingMonitor.h>

Inheritance diagram for HcalDetDiagTimingMonitor:

Public Member Functions
void	analyze (const edm::Event &iEvent, const edm::EventSetup &iSetup)

void	bookHistograms (DQMStore::IBooker &ib, const edm::Run &run, const edm::EventSetup &c)

void	done ()

double	get_ped_hbhe (int eta, int phi, int depth, int cup)

double	get_ped_hf (int eta, int phi, int depth, int cup)

double	get_ped_ho (int eta, int phi, int depth, int cup)

double	GetTime (double *data, int n)

	HcalDetDiagTimingMonitor (const edm::ParameterSet &ps)

bool	isSignal (double *data, int n)

void	reset ()

void	set_hbhe (int eta, int phi, int depth, int cap, float val)

void	set_hf (int eta, int phi, int depth, int cap, float val)

void	set_ho (int eta, int phi, int depth, int cap, float val)

void	setup (DQMStore::IBooker &)

	~HcalDetDiagTimingMonitor ()

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

	HcalBaseDQMonitor ()

virtual	~HcalBaseDQMonitor ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Public Attributes
double	HBHE [100][73][5][4]

double	HF [100][73][5][4]

double	HO [100][73][5][4]

double	nHBHE [100][73][5][4]

double	nHF [100][73][5][4]

double	nHO [100][73][5][4]

double	occHBHE [100][73][5]

double	occHF [100][73][5]

double	occHO [100][73][5]

double	occSum

Private Member Functions
void	CheckTiming ()

Private Attributes
bool	CosmicsCorr_

int	GCTTriggerBit1_

int	GCTTriggerBit2_

int	GCTTriggerBit3_

int	GCTTriggerBit4_

int	GCTTriggerBit5_

MonitorElement *	HBHEShape

MonitorElement *	HBTimeDT

MonitorElement *	HBTimeGCT

MonitorElement *	HBTimeHO

MonitorElement *	HBTimeRPC

MonitorElement *	HETimeCSCm

MonitorElement *	HETimeCSCp

MonitorElement *	HETimeRPCm

MonitorElement *	HETimeRPCp

MonitorElement *	HFTimeCSCm

MonitorElement *	HFTimeCSCp

MonitorElement *	HOShape

MonitorElement *	HOTimeDT

MonitorElement *	HOTimeGCT

MonitorElement *	HOTimeHO

MonitorElement *	HOTimeRPC

edm::InputTag	inputLabelDigi_

edm::InputTag	L1ADataLabel_

MonitorElement *	Summary

edm::EDGetTokenT < HBHEDigiCollection >	tok_hbhe_

edm::EDGetTokenT < HFDigiCollection >	tok_hf_

edm::EDGetTokenT < HODigiCollection >	tok_ho_

edm::EDGetTokenT < L1GlobalTriggerReadoutRecord >	tok_l1_

edm::EDGetTokenT < L1MuGMTReadoutCollection >	tok_l1mu_

edm::EDGetTokenT < FEDRawDataCollection >	tok_raw_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from HcalBaseDQMonitor
virtual void	beginLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)

virtual void	cleanup (void)

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

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 (DQMStore::IBooker &ib, 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_

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

Author: D. Vishnevskiy

Definition at line 28 of file HcalDetDiagTimingMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 38 of file HcalDetDiagTimingMonitor.cc.

                                                                             :HcalBaseDQMonitor(ps)
 {
   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","DetDiagTimingMonitor_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);
 
   GCTTriggerBit1_= ps.getUntrackedParameter<int>("GCTTriggerBit1", 15);         
   GCTTriggerBit2_= ps.getUntrackedParameter<int>("GCTTriggerBit2", 16);         
   GCTTriggerBit3_= ps.getUntrackedParameter<int>("GCTTriggerBit3", 17);         
   GCTTriggerBit4_= ps.getUntrackedParameter<int>("GCTTriggerBit4", 18);         
   GCTTriggerBit5_= ps.getUntrackedParameter<int>("GCTTriggerBit5", 16);         
   CosmicsCorr_   = ps.getUntrackedParameter<bool>("CosmicsCorr", true); 
   
   L1ADataLabel_  = ps.getUntrackedParameter<edm::InputTag>("gtLabel");
   inputLabelDigi_= ps.getUntrackedParameter<edm::InputTag>("digiLabel");
 
   // register for data access
   tok_raw_ = consumes<FEDRawDataCollection>(ps.getUntrackedParameter<edm::InputTag>("FEDRawDataCollection",edm::InputTag("source","")));
   tok_l1_ = consumes<L1GlobalTriggerReadoutRecord>(L1ADataLabel_);
   tok_l1mu_ = consumes<L1MuGMTReadoutCollection>(L1ADataLabel_);
   tok_hbhe_ = consumes<HBHEDigiCollection>(inputLabelDigi_);
   tok_ho_ = consumes<HODigiCollection>(inputLabelDigi_);
   tok_hf_ = consumes<HFDigiCollection>(inputLabelDigi_);
 }

HcalDetDiagTimingMonitor::~HcalDetDiagTimingMonitor ( )

Definition at line 75 of file HcalDetDiagTimingMonitor.cc.

75 {}

Member Function Documentation

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 138 of file HcalDetDiagTimingMonitor.cc.

References hcaldqm::constants::adc2fC, HcalBaseDQMonitor::analyze(), CheckTiming(), CosmicsCorr_, FEDRawData::data(), data, HLT_25ns10e33_v2_cff::depth, relval_parameters_module::energy, eta, MonitorElement::Fill(), GCTTriggerBit1_, GCTTriggerBit2_, GCTTriggerBit3_, GCTTriggerBit4_, GCTTriggerBit5_, get_ped_hbhe(), get_ped_hf(), get_ped_ho(), edm::Event::getByToken(), L1MuGMTReadoutCollection::getRecords(), GetTime(), HBHEShape, HBTimeDT, HBTimeGCT, HBTimeHO, HBTimeRPC, hc_Null, HcalBarrel, HcalEndcap, HETimeCSCm, HETimeCSCp, HETimeRPCm, HETimeRPCp, HFTimeCSCm, HFTimeCSCp, HOShape, HOTimeDT, HOTimeGCT, HOTimeHO, HOTimeRPC, i, HcalBaseDQMonitor::ievt_, HcalBaseDQMonitor::IsAllowedCalibType(), isSignal(), edm::HandleBase::isValid(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), MAXCSC, MAXDTBX, FEDNumbering::MAXHCALFEDID, FEDNumbering::MAXHCALuTCAFEDID, MAXRPC, FEDNumbering::MINHCALFEDID, FEDNumbering::MINHCALuTCAFEDID, N, occHBHE, occHF, occHO, occSum, phi, edm::Handle< T >::product(), lumiPlot::rawdata, set_hbhe(), set_hf(), set_ho(), funct::sin(), FEDRawData::size(), tok_hbhe_, tok_hf_, tok_ho_, tok_l1_, tok_l1mu_, tok_raw_, TRIG_CSC, TRIG_DT, TRIG_GCT, TRIG_RPC, and TRIG_RPCF.

 {
   HcalBaseDQMonitor::analyze(iEvent, iSetup);
   if (!IsAllowedCalibType()) return;
   if (LumiInOrder(iEvent.luminosityBlock())==false) return;
 //  HcalBaseDQMonitor::analyze(iEvent, iSetup);
   
   int eta,phi,depth,nTS,BXinEVENT=1,TRIGGER=0;
   
   // We do not want to look at Abort Gap events
   edm::Handle<FEDRawDataCollection> rawdata;
   iEvent.getByToken(tok_raw_,rawdata);
   //checking FEDs for calibration information
   if(!rawdata.isValid()) return;
   for(int i=FEDNumbering::MINHCALFEDID;
           i<=FEDNumbering::MAXHCALuTCAFEDID; i++) 
   {
     if (i>FEDNumbering::MAXHCALFEDID && i<FEDNumbering::MINHCALuTCAFEDID)
         continue;
 
     const FEDRawData& fedData = rawdata->FEDData(i) ;
     if ( fedData.size() < 24 ) continue ;
     if(((const HcalDCCHeader*)(fedData.data()))->getCalibType()!=hc_Null) return;
   }
   bool GCTTrigger1=false,GCTTrigger2=false,GCTTrigger3=false,GCTTrigger4=false,GCTTrigger5=false,HOselfTrigger=false;
   // Check GCT trigger bits
   edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
   iEvent.getByToken(tok_l1_, gtRecord);
   if(gtRecord.isValid()){
 
     const TechnicalTriggerWord tWord = gtRecord->technicalTriggerWord();
     const DecisionWord         dWord = gtRecord->decisionWord();
     //bool HFselfTrigger   = tWord.at(9);
     if (!tWord.empty()) HOselfTrigger    = tWord.at(11);
 
     if (!dWord.empty()) 
       {
     GCTTrigger1      = dWord.at(GCTTriggerBit1_);     
     GCTTrigger2      = dWord.at(GCTTriggerBit2_);     
     GCTTrigger3      = dWord.at(GCTTriggerBit3_);     
     GCTTrigger4      = dWord.at(GCTTriggerBit4_);     
     GCTTrigger5      = dWord.at(GCTTriggerBit5_);     
       }
 
     // define trigger trigger source (example from GMT group)
     edm::Handle<L1MuGMTReadoutCollection> gmtrc_handle; 
     iEvent.getByToken(tok_l1mu_,gmtrc_handle);
     if(!gmtrc_handle.isValid()) return;
     L1MuGMTReadoutCollection const* gmtrc = gmtrc_handle.product();
    
     int idt   =0;
     int icsc  =0;
     int irpcb =0;
     int irpcf =0;
     int ndt[5]   = {0,0,0,0,0};
     int ncsc[5]  = {0,0,0,0,0};
     int nrpcb[5] = {0,0,0,0,0};
     int nrpcf[5] = {0,0,0,0,0};
     int N;  
     std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
     std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;
     N=0;
     int NN=0;
     for(igmtrr=gmt_records.begin(); igmtrr!=gmt_records.end(); igmtrr++) {
       if(igmtrr->getBxInEvent()==0) BXinEVENT=NN;
       NN++;
       std::vector<L1MuRegionalCand>::const_iterator iter1;
       std::vector<L1MuRegionalCand> rmc;
       // DTBX Trigger
       rmc = igmtrr->getDTBXCands(); 
       for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
     if ( idt < MAXDTBX && !(*iter1).empty() ) {
       idt++; 
       if(N<5) ndt[N]++; 
                  
     }    
       }
       // CSC Trigger
       rmc = igmtrr->getCSCCands(); 
       for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
     if ( icsc < MAXCSC && !(*iter1).empty() ) {
       icsc++; 
       if(N<5) ncsc[N]++; 
     } 
       }
       // RPCb Trigger
       rmc = igmtrr->getBrlRPCCands();
       for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
     if ( irpcb < MAXRPC && !(*iter1).empty() ) {
       irpcb++;
       if(N<5) nrpcb[N]++;
                 
     }  
       }
       // RPCfwd Trigger
       rmc = igmtrr->getFwdRPCCands();
       for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
     if ( irpcf < MAXRPC && !(*iter1).empty() ) {
       irpcf++;
       if(N<5) nrpcf[N]++;
                 
     }  
       }
         
       N++;
     }
     if(ncsc[BXinEVENT]>0 ) {   TRIGGER=+TRIG_CSC;  }
     if(ndt[BXinEVENT]>0  ) {   TRIGGER=+TRIG_DT;   }
     if(nrpcb[BXinEVENT]>0) {   TRIGGER=+TRIG_RPC;  }
     if(nrpcf[BXinEVENT]>0) {   TRIGGER=+TRIG_RPCF; }
     if(GCTTrigger1 || GCTTrigger2 || GCTTrigger3 || GCTTrigger4 || GCTTrigger5) {  TRIGGER=+TRIG_GCT; }
   }
   if(ievt_<100){
     edm::Handle<HBHEDigiCollection> hbhe; 
     iEvent.getByToken(tok_hbhe_,hbhe);
     if(hbhe.isValid()){   
       for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
     eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
     for(int i=0;i<nTS;i++) if(digi->sample(i).adc()<20) set_hbhe(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
       }   
     }    
     edm::Handle<HODigiCollection> ho; 
     iEvent.getByToken(tok_ho_,ho);
     if(ho.isValid()){
       for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
     eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
     for(int i=0;i<nTS;i++) if(digi->sample(i).adc()<20) set_ho(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
       }   
     }
     edm::Handle<HFDigiCollection> hf;
     iEvent.getByToken(tok_hf_,hf);
     if(hf.isValid()){
       for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
     eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
     for(int i=0;i<nTS;i++) if(digi->sample(i).adc()<20) set_hf(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
       }   
     }
     return;   
   }
   double data[20];
   edm::Handle<HBHEDigiCollection> hbhe; 
   iEvent.getByToken(tok_hbhe_,hbhe);
   if(hbhe.isValid()){ 
     for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
       eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
       for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()]-get_ped_hbhe(eta,phi,depth,digi->sample(i).capid());
 
       double energy=0;
       for(int i=0;i<nTS;i++) energy+=data[i]; 
       if(energy>20) for(int i=0;i<nTS;i++) HBHEShape->Fill(i,data[i]);
 
       if(!isSignal(data,nTS)) continue;
       
       occHBHE[eta+50][phi][depth]+=1.0; occSum+=1.0;
       if((occHBHE[eta+50][phi][depth]/(double)(ievt_))>0.001) continue;
          
       double Time=GetTime(data,nTS);
       if(digi->id().subdet()==HcalBarrel){
     if(TRIGGER==TRIG_GCT) HBTimeGCT->Fill(Time);
     if(CosmicsCorr_) Time+=(7.5*sin((phi*5.0)/180.0*3.14159))/25.0;
     if(TRIGGER==TRIG_DT)  HBTimeDT->Fill(Time);
     if(HOselfTrigger)     HBTimeHO->Fill(Time);
     if(TRIGGER==TRIG_RPC) HBTimeRPC->Fill(Time);
       }
       if(digi->id().subdet()==HcalEndcap){
     if(CosmicsCorr_) Time+=(3.5*sin((phi*5.0)/180.0*3.14159))/25.0; 
     if(TRIGGER==TRIG_CSC && eta>0)  HETimeCSCp->Fill(Time);
     if(TRIGGER==TRIG_CSC && eta<0)  HETimeCSCm->Fill(Time);  
     if(TRIGGER==TRIG_RPCF && eta>0) HETimeRPCp->Fill(Time);
     if(TRIGGER==TRIG_RPCF && eta<0) HETimeRPCm->Fill(Time);          
       }
      
     }    
   }
   edm::Handle<HODigiCollection> ho; 
   iEvent.getByToken(tok_ho_,ho);
   if(ho.isValid()){ 
     for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
       eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
       for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()]-get_ped_ho(eta,phi,depth,digi->sample(i).capid());
 
       double energy=0;
       for(int i=0;i<nTS;i++) energy+=data[i]; 
       if(energy>100) for(int i=0;i<nTS;i++) HOShape->Fill(i,data[i]);
       if(energy<100) continue;
 
       occHO[eta+50][phi][depth]+=1.0;
       occSum+=1.0;
       if((occHO[eta+50][phi][depth]/(double)(ievt_))>0.001) continue;
          
       double Time=GetTime(data,nTS);
       if(CosmicsCorr_) Time+=(12.0*sin((phi*5.0)/180.0*3.14159))/25.0;  
       if(TRIGGER==TRIG_DT)  HOTimeDT ->Fill(Time);
       if(HOselfTrigger)     HOTimeHO ->Fill(Time);
       if(TRIGGER==TRIG_RPC) HOTimeRPC->Fill(Time);
       if(TRIGGER==TRIG_GCT) HOTimeGCT->Fill(Time);
     }   
   }
   edm::Handle<HFDigiCollection> hf; 
   iEvent.getByToken(tok_hf_,hf);
   if(hf.isValid()){ 
     for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
       eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
       double energy=0;
       for(int i=0;i<nTS;i++){
     data[i]=adc2fC[digi->sample(i).adc()]-get_ped_hf(eta,phi,depth,digi->sample(i).capid());
     energy+=data[i]; 
       }
         
       if(energy<25.0) continue;
       occHF[eta+50][phi][depth]+=1.0;
       occSum+=1.0;
         
       double Time=GetTime(data,nTS);
       if((occHF[eta+50][phi][depth]/(double)(ievt_))>0.01) continue;
         
       if(TRIGGER==TRIG_CSC && eta>0) HFTimeCSCp->Fill(Time); 
       if(TRIGGER==TRIG_CSC && eta<0) HFTimeCSCm->Fill(Time); 
     }   
   }
   if((ievt_%500)==0){
     CheckTiming();
     //printf("Run: %i, Evants processed: %i\n",iEvent.run(),ievt_);       
   }
 }

void HcalDetDiagTimingMonitor::bookHistograms	(	DQMStore::IBooker &	ib,
		const edm::Run &	run,
		const edm::EventSetup &	c
	)

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 79 of file HcalDetDiagTimingMonitor.cc.

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

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

void HcalDetDiagTimingMonitor::CheckTiming ( )

private

Definition at line 370 of file HcalDetDiagTimingMonitor.cc.

References MonitorElement::getEntries(), HBTimeDT, HBTimeGCT, HBTimeHO, HBTimeRPC, HETimeCSCm, HETimeCSCp, HETimeRPCm, HETimeRPCp, HFTimeCSCm, HFTimeCSCp, HOTimeDT, HOTimeGCT, HOTimeHO, HOTimeRPC, MonitorElement::setBinContent(), and Summary.

Referenced by analyze().

                                           {
   if(HBTimeDT->getEntries()>10){
     Summary->setBinContent(1,1,1); 
   }
   if(HBTimeRPC->getEntries()>10){
     Summary->setBinContent(2,1,1); 
   } 
   if(HBTimeGCT->getEntries()>10){
     Summary->setBinContent(3,1,1); 
   } 
   if(HBTimeHO->getEntries()>10){
     Summary->setBinContent(6,1,1); 
   } 
   if(HOTimeDT->getEntries()>10){
     Summary->setBinContent(1,2,1); 
   }
   if(HOTimeRPC->getEntries()>10){
     Summary->setBinContent(2,2,1); 
   } 
   if(HOTimeGCT->getEntries()>10){
     Summary->setBinContent(3,2,1); 
   } 
   if(HOTimeHO->getEntries()>10){
     Summary->setBinContent(6,2,1); 
   } 
   if(HETimeCSCp->getEntries()>10){
     Summary->setBinContent(4,4,1); 
   } 
   if(HETimeCSCm->getEntries()>10){
     Summary->setBinContent(4,3,1); 
   } 
   if(HETimeRPCp->getEntries()>10){
     Summary->setBinContent(5,4,1); 
   } 
   if(HETimeRPCm->getEntries()>10){
     Summary->setBinContent(5,3,1); 
   } 
   if(HFTimeCSCp->getEntries()>10){
     Summary->setBinContent(4,6,1); 
   } 
   if(HFTimeCSCm->getEntries()>10){
     Summary->setBinContent(4,4,1); 
   } 
 }

void HcalDetDiagTimingMonitor::done ( )

Definition at line 415 of file HcalDetDiagTimingMonitor.cc.

415 { }

double HcalDetDiagTimingMonitor::get_ped_hbhe	(	int	eta,
		int	phi,
		int	depth,
		int	cup
	)

inline

Definition at line 67 of file HcalDetDiagTimingMonitor.h.

References HLT_25ns10e33_v2_cff::depth, HBHE, nHBHE, and phi.

Referenced by analyze().

                                                         {
       if(nHBHE[eta+50][phi][depth][cup]<10) return 2.5; 
       if(nHBHE[eta+50][phi][depth][cup]!=0){
          double ped=HBHE[eta+50][phi][depth][cup]/nHBHE[eta+50][phi][depth][cup];
          if(ped>1.5 && ped<4.5) return ped;
       } 
       return 9999; 
   }   

double HcalDetDiagTimingMonitor::get_ped_hf	(	int	eta,
		int	phi,
		int	depth,
		int	cup
	)

inline

Definition at line 83 of file HcalDetDiagTimingMonitor.h.

References HLT_25ns10e33_v2_cff::depth, HF, nHF, and phi.

Referenced by analyze().

                                                       {
       if(nHF[eta+50][phi][depth][cup]<10) return 2.5; 
       if(nHF[eta+50][phi][depth][cup]!=0){
          double ped=HF[eta+50][phi][depth][cup]/nHF[eta+50][phi][depth][cup];
          if(ped>1.5 && ped<4.5) return ped;
       }
       return 9999; 
   }   

double HcalDetDiagTimingMonitor::get_ped_ho	(	int	eta,
		int	phi,
		int	depth,
		int	cup
	)

inline

Definition at line 75 of file HcalDetDiagTimingMonitor.h.

References HLT_25ns10e33_v2_cff::depth, HO, nHO, and phi.

Referenced by analyze().

                                                       {
       if(nHO[eta+50][phi][depth][cup]<10) return 8.5; 
       if(nHO[eta+50][phi][depth][cup]!=0){
          double ped=HO[eta+50][phi][depth][cup]/nHO[eta+50][phi][depth][cup];
          if(ped>5 && ped<15) return ped;
       }
       return 9999; 
   }   

double HcalDetDiagTimingMonitor::GetTime	(	double *	data,
		int	n
	)

inline

Definition at line 33 of file HcalDetDiagTimingMonitor.h.

References j, and gen::n.

Referenced by analyze().

                                     {
         int MaxI=-100; double Time=0,SumT=0,MaxT=-10;
         for(int j=0;j<n;++j) if(MaxT<data[j]){ MaxT=data[j]; MaxI=j; }
     if(MaxI>=0){
            Time=MaxI*data[MaxI];
            SumT=data[MaxI];
        if(MaxI>0){ Time+=(MaxI-1)*data[MaxI-1]; SumT+=data[MaxI-1]; }
        if(MaxI<(n-1)){ Time+=(MaxI+1)*data[MaxI+1]; SumT+=data[MaxI+1]; }
        Time=Time/SumT;
     }
         return Time;
   }

bool HcalDetDiagTimingMonitor::isSignal	(	double *	data,
		int	n
	)

inline

Definition at line 45 of file HcalDetDiagTimingMonitor.h.

References i, and bookConverter::max.

Referenced by analyze().

                                    {
         int Imax=-1; double max=-100;
         for(int i=0;i<n;i++) if(data[i]>max){max=data[i]; Imax=i;}
         if(Imax==0 && Imax==(n-1)) return false;
         float sum=data[Imax-1]+data[Imax+1];
         if(data[Imax]>5.5 && sum>(data[Imax]*0.25)) return true;
         return false;
   }

void HcalDetDiagTimingMonitor::reset ( void )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 77 of file HcalDetDiagTimingMonitor.cc.

Referenced by bookHistograms().

77 {}

void HcalDetDiagTimingMonitor::set_hbhe	(	int	eta,
		int	phi,
		int	depth,
		int	cap,
		float	val
	)

inline

Definition at line 55 of file HcalDetDiagTimingMonitor.h.

References cap, HLT_25ns10e33_v2_cff::depth, HBHE, nHBHE, and phi.

Referenced by analyze().

                                                             {
        HBHE[eta+50][phi][depth][cap]+=val;
        nHBHE[eta+50][phi][depth][cap]+=1.0;
   }   

void HcalDetDiagTimingMonitor::set_hf	(	int	eta,
		int	phi,
		int	depth,
		int	cap,
		float	val
	)

inline

Definition at line 63 of file HcalDetDiagTimingMonitor.h.

References cap, HLT_25ns10e33_v2_cff::depth, HF, nHF, and phi.

Referenced by analyze().

                                                           {
        HF[eta+50][phi][depth][cap]+=val;
        nHF[eta+50][phi][depth][cap]+=1.0;
   }

void HcalDetDiagTimingMonitor::set_ho	(	int	eta,
		int	phi,
		int	depth,
		int	cap,
		float	val
	)

inline

Definition at line 59 of file HcalDetDiagTimingMonitor.h.

References cap, HLT_25ns10e33_v2_cff::depth, HO, nHO, and phi.

Referenced by analyze().

                                                           {
        HO[eta+50][phi][depth][cap]+=val;
        nHO[eta+50][phi][depth][cap]+=1.0;
   }   

void HcalDetDiagTimingMonitor::setup ( DQMStore::IBooker & ib )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 94 of file HcalDetDiagTimingMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), HBHEShape, HBTimeDT, HBTimeGCT, HBTimeHO, HBTimeRPC, HETimeCSCm, HETimeCSCp, HETimeRPCm, HETimeRPCp, HFTimeCSCm, HFTimeCSCp, HOShape, HOTimeDT, HOTimeGCT, HOTimeHO, HOTimeRPC, i, j, MonitorElement::setBinContent(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), HcalBaseDQMonitor::setup(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, and Summary.

Referenced by bookHistograms().

 {
   
   HcalBaseDQMonitor::setup(ib);
 
   std::string str;
      ib.setCurrentFolder(subdir_);   
      str="Hcal Timing summary"; Summary = ib.book2D(str,str,6,0,6,6,0,6); 
      Summary->setBinLabel(1,"DT",1);
      Summary->setBinLabel(2,"RPC",1);
      Summary->setBinLabel(3,"GCT",1);
      Summary->setBinLabel(4,"CSC",1);
      Summary->setBinLabel(5,"RPCf",1);
      Summary->setBinLabel(6,"bit11",1);
      Summary->setBinLabel(1,"HB",2);
      Summary->setBinLabel(2,"HO",2);
      Summary->setBinLabel(3,"HEM",2);
      Summary->setBinLabel(4,"HEP",2);
      Summary->setBinLabel(5,"HFM",2);
      Summary->setBinLabel(6,"HFP",2);
      for(int i=1;i<=6;i++) for(int j=1;j<=6;j++) Summary->setBinContent(i,j,-1);
      
      ib.setCurrentFolder(subdir_+"Timing Plots");
      str="HB Timing (DT Trigger)";                      HBTimeDT  = ib.book1D(str,str,100,0,10); 
      str="HO Timing (DT Trigger)";                      HOTimeDT  = ib.book1D(str,str,100,0,10); 
      str="HB Timing (RPC Trigger)";                     HBTimeRPC = ib.book1D(str,str,100,0,10); 
      str="HO Timing (RPC Trigger)";                     HOTimeRPC = ib.book1D(str,str,100,0,10); 
      str="HB Timing (HO SelfTrigger tech bit 11)";      HBTimeHO  = ib.book1D(str,str,100,0,10); 
      str="HO Timing (HO SelfTrigger tech bit 11)";      HOTimeHO  = ib.book1D(str,str,100,0,10); 
      
      str="HB Timing (GCT Trigger alg bit 15 16 17 18)"; HBTimeGCT  =ib.book1D(str,str,100,0,10); 
      str="HO Timing (GCT Trigger alg bit 15 16 17 18)"; HOTimeGCT  =ib.book1D(str,str,100,0,10); 
      
      str="HEP Timing (CSC Trigger)";                    HETimeCSCp =ib.book1D(str,str,100,0,10); 
      str="HEM Timing (CSC Trigger)";                    HETimeCSCm =ib.book1D(str,str,100,0,10);
      str="HEP Timing (RPCf Trigger)";                   HETimeRPCp =ib.book1D(str,str,100,0,10); 
      str="HEM Timing (RPCf Trigger)";                   HETimeRPCm =ib.book1D(str,str,100,0,10);
      str="HFP Timing (CSC Trigger)";                    HFTimeCSCp =ib.book1D(str,str,100,0,10); 
      str="HFM Timing (CSC Trigger)";                    HFTimeCSCm =ib.book1D(str,str,100,0,10);     
      str="HBHE Shape";                                  HBHEShape  =ib.book1D(str,str,10,-0.5,9.5);
      str="HO Shape";                                    HOShape    =ib.book1D(str,str,10,-0.5,9.5);
 
 } 