HcalRecHitMonitor Class Reference

#include <HcalRecHitMonitor.h>

Inheritance diagram for HcalRecHitMonitor:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
void	beginRun (const edm::Run &run, const edm::EventSetup &c)
void	cleanup ()
void	endJob ()
void	endLuminosityBlock (const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
void	endRun (const edm::Run &run, const edm::EventSetup &c)
	HcalRecHitMonitor (const edm::ParameterSet &ps)
void	processEvent (const HBHERecHitCollection &hbHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, int BCN, const edm::Event &iEvent)
void	processEvent_rechit (const HBHERecHitCollection &hbheHits, const HORecHitCollection &hoHits, const HFRecHitCollection &hfHits, bool passedHcalHLT, bool passedMinBiasHLT, int BCN)
void	reset ()
void	setup ()
void	zeroCounters ()
	~HcalRecHitMonitor ()
Public Member Functions inherited from HcalBaseDQMonitor
	HcalBaseDQMonitor (const edm::ParameterSet &ps)
	HcalBaseDQMonitor ()
virtual	~HcalBaseDQMonitor ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
void	fill_Nevents ()

Private Attributes
double	energy2_ [85][72][4]
double	energy2_thresh_ [85][72][4]
double	energy_ [85][72][4]
double	energy_thresh_ [85][72][4]
double	energyThreshold_
double	ETThreshold_
MonitorElement *	h_FlagMap_DIGITIME
MonitorElement *	h_FlagMap_HPDMULT
MonitorElement *	h_FlagMap_LONGSHORT
MonitorElement *	h_FlagMap_PULSESHAPE
MonitorElement *	h_FlagMap_TIMEADD
MonitorElement *	h_FlagMap_TIMEERROR
MonitorElement *	h_FlagMap_TIMESUBTRACT
MonitorElement *	h_HBflagcounter
MonitorElement *	h_HBHE_FlagCorr
MonitorElement *	h_HBHEHPDMult_vs_LS
MonitorElement *	h_HBHEPulseShape_vs_LS
MonitorElement *	h_HBM_weightedTime
MonitorElement *	h_HBOccupancy
MonitorElement *	h_HBP_weightedTime
MonitorElement *	h_HBsizeVsLS
MonitorElement *	h_HBThreshOccupancy
MonitorElement *	h_HBThreshTime
MonitorElement *	h_HBTime
MonitorElement *	h_HBTimeVsEnergy
MonitorElement *	h_HE_HcalHLT_energydifference
MonitorElement *	h_HE_HcalHLT_weightedtimedifference
MonitorElement *	h_HEenergydifference
MonitorElement *	h_HEflagcounter
MonitorElement *	h_HEM_weightedTime
MonitorElement *	h_HEOccupancy
MonitorElement *	h_HEP_weightedTime
MonitorElement *	h_HEsizeVsLS
MonitorElement *	h_HEThreshOccupancy
MonitorElement *	h_HEThreshTime
MonitorElement *	h_HETime
MonitorElement *	h_HEtimedifference
MonitorElement *	h_HETimeVsEnergy
MonitorElement *	h_HF_FlagCorr
MonitorElement *	h_HF_HcalHLT_energydifference
MonitorElement *	h_HF_HcalHLT_weightedtimedifference
MonitorElement *	h_HFDigiTime_vs_LS
MonitorElement *	h_HFenergydifference
MonitorElement *	h_HFflagcounter
MonitorElement *	h_HFLongShort_vs_LS
MonitorElement *	h_HFM_weightedTime
MonitorElement *	h_HFOccupancy
MonitorElement *	h_HFP_weightedTime
MonitorElement *	h_HFsizeVsLS
MonitorElement *	h_HFThreshOccupancy
MonitorElement *	h_HFThreshTime
MonitorElement *	h_HFTime
MonitorElement *	h_HFtimedifference
MonitorElement *	h_HFTimeVsEnergy
MonitorElement *	h_HOflagcounter
MonitorElement *	h_HOOccupancy
MonitorElement *	h_HOsizeVsLS
MonitorElement *	h_HOThreshOccupancy
MonitorElement *	h_HOThreshTime
MonitorElement *	h_HOTime
MonitorElement *	h_HOTimeVsEnergy
MonitorElement *	h_LumiPlot_BX_allevents
MonitorElement *	h_LumiPlot_BX_HcalHLTEvents
MonitorElement *	h_LumiPlot_BX_HcalHLTEvents_notimecut
MonitorElement *	h_LumiPlot_BX_MinBiasEvents
MonitorElement *	h_LumiPlot_BX_MinBiasEvents_notimecut
MonitorElement *	h_LumiPlot_LS_allevents
MonitorElement *	h_LumiPlot_LS_HcalHLTEvents
MonitorElement *	h_LumiPlot_LS_HcalHLTEvents_notimecut
MonitorElement *	h_LumiPlot_LS_MinBiasEvents
MonitorElement *	h_LumiPlot_LS_MinBiasEvents_notimecut
MonitorElement *	h_LumiPlot_MinTime_vs_MinHT
MonitorElement *	h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus
MonitorElement *	h_LumiPlot_SumHT_HFPlus_vs_HFMinus
MonitorElement *	h_LumiPlot_timeHFPlus_vs_timeHFMinus
MonitorElement *	h_LumiPlot_timeHT_HFM
MonitorElement *	h_LumiPlot_timeHT_HFP
MonitorElement *	h_rechitieta
MonitorElement *	h_rechitieta_05
MonitorElement *	h_rechitieta_10
MonitorElement *	h_rechitieta_100
MonitorElement *	h_rechitieta_25
MonitorElement *	h_rechitieta_thresh
MonitorElement *	h_rechitiphi
MonitorElement *	h_rechitiphi_05
MonitorElement *	h_rechitiphi_10
MonitorElement *	h_rechitiphi_100
MonitorElement *	h_rechitiphi_25
MonitorElement *	h_rechitiphi_thresh
MonitorElement *	h_TriggeredEvents
double	HB_occupancy_ [260]
double	HB_occupancy_thresh_ [260]
double	HBenergyThreshold_
double	HBETThreshold_
int	HBflagcounter_ [32]
edm::InputTag	hbheRechitLabel_
bool	HBpresent_
double	HBtime_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HBtime_thresh_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
std::vector< std::string >	HcalHLTBits_
double	HE_occupancy_ [260]
double	HE_occupancy_thresh_ [260]
double	HEenergyThreshold_
double	HEETThreshold_
int	HEflagcounter_ [32]
bool	HEpresent_
double	HEtime_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HEtime_thresh_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HF_occupancy_ [174]
double	HF_occupancy_thresh_ [174]
double	HFenergyLong_ [200]
double	HFenergyLong_thresh_ [200]
double	HFenergyShort_ [200]
double	HFenergyShort_thresh_ [200]
double	HFenergyThreshold_
double	HFETThreshold_
int	HFflagcounter_ [32]
double	HFlong_occupancy_ [865]
double	HFlong_occupancy_thresh_ [865]
MonitorElement *	HFP_HFM_Energy
bool	HFpresent_
edm::InputTag	hfRechitLabel_
double	HFshort_occupancy_ [865]
double	HFshort_occupancy_thresh_ [865]
double	HFtime_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HFtime_thresh_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HFtimeLong_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HFtimeLong_thresh_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HFtimeShort_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HFtimeShort_thresh_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
edm::InputTag	hltresultsLabel_
double	HO_occupancy_ [218]
double	HO_occupancy_thresh_ [218]
double	HOenergyThreshold_
double	HOETThreshold_
int	HOflagcounter_ [32]
bool	HOpresent_
edm::InputTag	hoRechitLabel_
double	HOtime_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
double	HOtime_thresh_ [RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]
edm::InputTag	l1gtLabel_
std::vector< std::string >	MinBiasHLTBits_
unsigned int	occupancy_ [85][72][4]
unsigned int	occupancy_thresh_ [85][72][4]
EtaPhiHists	OccupancyByDepth
EtaPhiHists	OccupancyThreshByDepth
bool	setupDone_
EtaPhiHists	SqrtSumEnergy2ByDepth
EtaPhiHists	SqrtSumEnergy2ThreshByDepth
EtaPhiHists	SumEnergyByDepth
EtaPhiHists	SumEnergyThreshByDepth
EtaPhiHists	SumTimeByDepth
EtaPhiHists	SumTimeThreshByDepth
double	time_ [85][72][4]
double	time_thresh_ [85][72][4]
double	timediffThresh_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
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)
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::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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

Date:

2012/11/12 20:52:11

Revision:

1.53

Author

J. Temple - Univ. of Maryland

Definition at line 17 of file HcalRecHitMonitor.h.

Constructor & Destructor Documentation

HcalRecHitMonitor::HcalRecHitMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 16 of file HcalRecHitMonitor.cc.

References HcalBaseDQMonitor::AllowedCalibTypes_, HcalBaseDQMonitor::debug_, HcalBaseDQMonitor::enableCleanup_, energyThreshold_, ETThreshold_, edm::ParameterSet::getUntrackedParameter(), HBenergyThreshold_, HBETThreshold_, hbheRechitLabel_, HcalHLTBits_, HEenergyThreshold_, HEETThreshold_, HFenergyThreshold_, HFETThreshold_, hfRechitLabel_, hltresultsLabel_, HOenergyThreshold_, HOETThreshold_, hoRechitLabel_, l1gtLabel_, HcalBaseDQMonitor::makeDiagnostics_, HcalBaseDQMonitor::mergeRuns_, MinBiasHLTBits_, HcalBaseDQMonitor::needLogicalMap_, HcalBaseDQMonitor::NLumiBlocks_, HcalBaseDQMonitor::Online_, HcalBaseDQMonitor::prefixME_, setupDone_, HcalBaseDQMonitor::skipOutOfOrderLS_, AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, and timediffThresh_.

{
  // Common Base Class parameters
  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","RecHitMonitor_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);

  hbheRechitLabel_       = ps.getUntrackedParameter<edm::InputTag>("hbheRechitLabel");
  hoRechitLabel_         = ps.getUntrackedParameter<edm::InputTag>("hoRechitLabel");
  hfRechitLabel_         = ps.getUntrackedParameter<edm::InputTag>("hfRechitLabel");
  l1gtLabel_             = ps.getUntrackedParameter<edm::InputTag>("L1GTLabel"); // should be l1GtUnpack

  hltresultsLabel_       = ps.getUntrackedParameter<edm::InputTag>("HLTResultsLabel");
  HcalHLTBits_           = ps.getUntrackedParameter<std::vector<std::string> >("HcalHLTBits");
  MinBiasHLTBits_        = ps.getUntrackedParameter<std::vector<std::string> >("MinBiasHLTBits");

  // energy/ET threshold plots also require that at least one MinBias trigger bit fires
  energyThreshold_       = ps.getUntrackedParameter<double>("energyThreshold",2);
  HBenergyThreshold_     = ps.getUntrackedParameter<double>("HB_energyThreshold",energyThreshold_);
  HEenergyThreshold_     = ps.getUntrackedParameter<double>("HE_energyThreshold",energyThreshold_);
  HOenergyThreshold_     = ps.getUntrackedParameter<double>("HO_energyThreshold",energyThreshold_);
  HFenergyThreshold_     = ps.getUntrackedParameter<double>("HF_energyThreshold",energyThreshold_);
  
  ETThreshold_       = ps.getUntrackedParameter<double>("ETThreshold",0);
  HBETThreshold_     = ps.getUntrackedParameter<double>("HB_ETThreshold",ETThreshold_);
  HEETThreshold_     = ps.getUntrackedParameter<double>("HE_ETThreshold",ETThreshold_);
  HOETThreshold_     = ps.getUntrackedParameter<double>("HO_ETThreshold",ETThreshold_);
  HFETThreshold_     = ps.getUntrackedParameter<double>("HF_ETThreshold",ETThreshold_);

  timediffThresh_    = ps.getUntrackedParameter<double>("collisiontimediffThresh",10.);
  setupDone_         = false;
  needLogicalMap_    = true;
} //constructor

HcalRecHitMonitor::~HcalRecHitMonitor

(

)

Definition at line 62 of file HcalRecHitMonitor.cc.

{
} //destructor

Member Function Documentation

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 583 of file HcalRecHitMonitor.cc.

References HcalBaseDQMonitor::analyze(), edm::EventBase::bunchCrossing(), gather_cfg::cout, HcalBaseDQMonitor::currentLS, HcalBaseDQMonitor::debug_, alignCSCRings::e, MonitorElement::Fill(), edm::Event::getByLabel(), HcalBaseDQMonitor::getLogicalMap(), h_LumiPlot_BX_allevents, h_LumiPlot_LS_allevents, hbheRechitLabel_, hfRechitLabel_, hoRechitLabel_, HcalBaseDQMonitor::IsAllowedCalibType(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), and processEvent().

{
  getLogicalMap(s);
  if (debug_>0)  std::cout <<"HcalRecHitMonitor::analyze; debug = "<<debug_<<std::endl;

  if (!IsAllowedCalibType()) return;
  if (LumiInOrder(e.luminosityBlock())==false) return;

  // Get objects
  edm::Handle<HBHERecHitCollection> hbhe_rechit;
  edm::Handle<HORecHitCollection> ho_rechit;
  edm::Handle<HFRecHitCollection> hf_rechit;

  if (!(e.getByLabel(hbheRechitLabel_,hbhe_rechit)))
    {
      edm::LogWarning("HcalHotCellMonitor")<< hbheRechitLabel_<<" hbhe_rechit not available";
      return;
    }

  if (!(e.getByLabel(hfRechitLabel_,hf_rechit)))
    {
      edm::LogWarning("HcalHotCellMonitor")<< hfRechitLabel_<<" hf_rechit not available";
      return;
    }
  
  if (!(e.getByLabel(hoRechitLabel_,ho_rechit)))
    {
      edm::LogWarning("HcalHotCellMonitor")<< hoRechitLabel_<<" ho_rechit not available";
      return;
    }


  h_LumiPlot_LS_allevents->Fill(currentLS);
  h_LumiPlot_BX_allevents->Fill(e.bunchCrossing());
  processEvent(*hbhe_rechit, *ho_rechit, *hf_rechit, e.bunchCrossing(), e);

  HcalBaseDQMonitor::analyze(e,s);
} // void HcalRecHitMonitor::analyze()

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 493 of file HcalRecHitMonitor.cc.

References HcalBaseDQMonitor::beginRun(), DQMStore::bookString(), gather_cfg::cout, HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::debug_, HcalHLTBits_, i, HcalBaseDQMonitor::mergeRuns_, MinBiasHLTBits_, reset(), DQMStore::setCurrentFolder(), setup(), AlCaHLTBitMon_QueryRunRegistry::string, HcalBaseDQMonitor::subdir_, and HcalBaseDQMonitor::tevt_.

{

  if (debug_>0) std::cout <<"HcalRecHitMonitor::beginRun():  task =  '"<<subdir_<<"'"<<std::endl;
  HcalBaseDQMonitor::beginRun(run, c);
  if (tevt_==0) // create histograms, if they haven't been created already
    this->setup();
  // Clear histograms at the start of each run if not merging runs
  if (mergeRuns_==false)
    this->reset();

  if (tevt_!=0) return;
  // create histograms displaying trigger parameters?  Specify names?
  dbe_->setCurrentFolder(subdir_+"rechit_parameters");
  std::string tnames="";
  if (HcalHLTBits_.size()>0)
    tnames=HcalHLTBits_[0];
  for (unsigned int i=1;i<HcalHLTBits_.size();++i)
    tnames=tnames + " OR " + HcalHLTBits_[i];
  dbe_->bookString("HcalHLTriggerRequirements",tnames);
  tnames="";
  if (MinBiasHLTBits_.size()>0)
    tnames=MinBiasHLTBits_[0];
  for (unsigned int i=1;i<MinBiasHLTBits_.size();++i)
    tnames=tnames + " OR " + MinBiasHLTBits_[i];
  dbe_->bookString("MinBiasHLTriggerRequirements",tnames);
  return;
  
} //void HcalRecHitMonitor::beginRun(...)

void HcalRecHitMonitor::cleanup ( void  )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 556 of file HcalRecHitMonitor.cc.

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

Referenced by endJob().

{
  //Add code to clean out subdirectories
  if (!enableCleanup_) return;
  if (dbe_)
    {
      dbe_->setCurrentFolder(subdir_); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"rechit_parameters"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"rechit_parameters/thresholds"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits/sumplots"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/sumplots"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers/passedTechTriggers/"); dbe_->removeContents();

      dbe_->setCurrentFolder(subdir_+"AnomalousCellFlags"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"diagnostics/hb"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"diagnostics/he"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"diagnostics/ho"); dbe_->removeContents();
      dbe_->setCurrentFolder(subdir_+"diagnostics/hf"); dbe_->removeContents();
    }
  return;
} // void HcalRecHitMonitor::cleanup()

void HcalRecHitMonitor::endJob ( void  )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 546 of file HcalRecHitMonitor.cc.

References cleanup(), HcalBaseDQMonitor::cleanup(), and HcalBaseDQMonitor::enableCleanup_.

{
  if (!enableCleanup_) return;
  HcalBaseDQMonitor::cleanup();
  this->cleanup();
}

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 1233 of file HcalRecHitMonitor.cc.

References fill_Nevents(), HcalBaseDQMonitor::LumiInOrder(), and edm::LuminosityBlockBase::luminosityBlock().

{
  // don't fill lumi block information if it's already been filled
  if (LumiInOrder(lumiSeg.luminosityBlock())==false) return;
  fill_Nevents();
  return;
} //endLuminosityBlock

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

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 524 of file HcalRecHitMonitor.cc.

References gather_cfg::cout, HcalBaseDQMonitor::debug_, and HcalBaseDQMonitor::subdir_.

{
  if (debug_>0) std::cout <<"HcalRecHitMonitor::endRun():  task =  '"<<subdir_<<"'"<<std::endl;

  //Any special fill calls needed?  Shouldn't be necessary; last endLuminosityBlock should do necessary fills
} // void HcalRecHitMonitor::endrun(...)

void HcalRecHitMonitor::fill_Nevents ( void  )

private

Definition at line 1244 of file HcalRecHitMonitor.cc.

References HcalObjRepresent::CalcIeta(), gather_cfg::cout, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, energy2_, energy2_thresh_, energy_, energy_thresh_, eta(), MonitorElement::Fill(), HcalObjRepresent::FillUnphysicalHEHFBins(), h_HBflagcounter, h_HBOccupancy, h_HBThreshOccupancy, h_HBThreshTime, h_HBTime, h_HEflagcounter, h_HEOccupancy, h_HEThreshOccupancy, h_HEThreshTime, h_HETime, h_HFflagcounter, h_HFOccupancy, h_HFThreshOccupancy, h_HFThreshTime, h_HFTime, h_HOflagcounter, h_HOOccupancy, h_HOThreshOccupancy, h_HOThreshTime, h_HOTime, h_rechitieta, h_rechitieta_thresh, h_rechitiphi, h_rechitiphi_thresh, HB_occupancy_, HB_occupancy_thresh_, HBflagcounter_, HBtime_, HBtime_thresh_, HE_occupancy_, HE_occupancy_thresh_, HEflagcounter_, HEtime_, HEtime_thresh_, HF_occupancy_, HF_occupancy_thresh_, HFflagcounter_, HFtime_, HFtime_thresh_, HO_occupancy_, HO_occupancy_thresh_, HOflagcounter_, HOtime_, HOtime_thresh_, i, HcalBaseDQMonitor::ievt_, gen::k, occupancy_, occupancy_thresh_, OccupancyByDepth, OccupancyThreshByDepth, phi, RECHITMON_TIME_MAX, RECHITMON_TIME_MIN, MonitorElement::setBinContent(), mathSSE::sqrt(), SqrtSumEnergy2ByDepth, SqrtSumEnergy2ThreshByDepth, SumEnergyByDepth, SumEnergyThreshByDepth, SumTimeByDepth, SumTimeThreshByDepth, time_, and time_thresh_.

Referenced by endLuminosityBlock().

{
  // looking at the contents of HbFlagcounters
  if (debug_>0)
    {
      for (int k = 0; k < 32; k++){
    std::cout << "<HcalRecHitMonitor::fill_Nevents>  HF Flag counter:  Bin #" << k+1 << " = "<< HFflagcounter_[k] << std::endl;
      }
    }

  for (int i=0;i<32;i++)
    {
      h_HBflagcounter->Fill(i,HBflagcounter_[i]);
      h_HEflagcounter->Fill(i,HEflagcounter_[i]);
      h_HOflagcounter->Fill(i,HOflagcounter_[i]);
      h_HFflagcounter->Fill(i,HFflagcounter_[i]);
      HBflagcounter_[i]=0;
      HEflagcounter_[i]=0;
      HOflagcounter_[i]=0;
      HFflagcounter_[i]=0;
    }

  // Fill Occupancy & Sum Energy, Time plots
  int myieta=-1;
  if (ievt_>0)
    {
      for (int mydepth=0;mydepth<4;++mydepth)
    {
      for (int eta=0;eta<OccupancyByDepth.depth[mydepth]->getNbinsX();++eta)
        {
          myieta=CalcIeta(eta,mydepth+1);

          for (int phi=0;phi<72;++phi)
        {
          if (occupancy_[eta][phi][mydepth]>0)
            {
              h_rechitieta->Fill(myieta,occupancy_[eta][phi][mydepth]);
              h_rechitiphi->Fill(phi+1,occupancy_[eta][phi][mydepth]);
            }
          if (occupancy_thresh_[eta][phi][mydepth]>0)
            {
              h_rechitieta_thresh->Fill(myieta,occupancy_thresh_[eta][phi][mydepth]);
              h_rechitiphi_thresh->Fill(phi+1,occupancy_thresh_[eta][phi][mydepth]);
            }
          OccupancyByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,occupancy_[eta][phi][mydepth]);
          SumEnergyByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,energy_[eta][phi][mydepth]);
                  SqrtSumEnergy2ByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,sqrt(energy2_[eta][phi][mydepth]));
          SumTimeByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,time_[eta][phi][mydepth]);

          OccupancyThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,occupancy_thresh_[eta][phi][mydepth]);
          SumEnergyThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,energy_thresh_[eta][phi][mydepth]);
          SqrtSumEnergy2ThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,sqrt(energy2_thresh_[eta][phi][mydepth]));
          SumTimeThreshByDepth.depth[mydepth]->setBinContent(eta+1,phi+1,time_thresh_[eta][phi][mydepth]);
        } // for (int phi=0;phi<72;++phi)
        } // for (int eta=0;eta<OccupancyByDepth...;++eta)
    } // for (int mydepth=0;...)

      FillUnphysicalHEHFBins(OccupancyByDepth);
      FillUnphysicalHEHFBins(OccupancyThreshByDepth);
      FillUnphysicalHEHFBins(SumEnergyByDepth);
      FillUnphysicalHEHFBins(SqrtSumEnergy2ByDepth);
      FillUnphysicalHEHFBins(SumEnergyThreshByDepth);
      FillUnphysicalHEHFBins(SqrtSumEnergy2ThreshByDepth);
      FillUnphysicalHEHFBins(SumTimeByDepth);
      FillUnphysicalHEHFBins(SumTimeThreshByDepth);

    } // if (ievt_>0)

  // Fill subdet plots

  for (int i=0;i<(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN);++i)
    {
      if (HBtime_[i]!=0)
    {
      h_HBTime->setBinContent(i+1,HBtime_[i]);
    }
      if (HBtime_thresh_[i]!=0)
    {
      h_HBThreshTime->setBinContent(i+1,HBtime_thresh_[i]);
    }
      if (HEtime_[i]!=0)
    {

      h_HETime->setBinContent(i+1,HEtime_[i]);
    }
      if (HEtime_thresh_[i]!=0)
    {
      h_HEThreshTime->setBinContent(i+1,HEtime_thresh_[i]);
    }
      if (HOtime_[i]!=0)
    {
      h_HOTime->setBinContent(i+1,HOtime_[i]);
    }
      if (HOtime_thresh_[i]!=0)
    {
      h_HOThreshTime->setBinContent(i+1,HOtime_thresh_[i]);
    }
      if (HFtime_[i]!=0)
    {
      h_HFTime->setBinContent(i+1,HFtime_[i]);
    }
      if (HFtime_thresh_[i]!=0)
    {
      h_HFThreshTime->setBinContent(i+1,HFtime_thresh_[i]);
    }
    } // for (int  i=0;i<(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN);++i)

  for (int i=0;i<260;++i)
    {
      if (HB_occupancy_[i]>0)
    {
      h_HBOccupancy->setBinContent(i+1,HB_occupancy_[i]);
    }
      if (HB_occupancy_thresh_[i]>0)
    {
      h_HBThreshOccupancy->setBinContent(i+1,HB_occupancy_thresh_[i]);
    }
      if (HE_occupancy_[i]>0)
    {
      h_HEOccupancy->setBinContent(i+1,HE_occupancy_[i]);
    }
      if (HE_occupancy_thresh_[i]>0)
    {
      h_HEThreshOccupancy->setBinContent(i+1,HE_occupancy_thresh_[i]);
    }
    }//for (int i=0;i<260;++i)

  for (int i=0;i<217;++i)
    {
      if (HO_occupancy_[i]>0)
    {
      h_HOOccupancy->setBinContent(i+1,HO_occupancy_[i]);
    }
      if (HO_occupancy_thresh_[i]>0)
    {
      h_HOThreshOccupancy->setBinContent(i+1,HO_occupancy_thresh_[i]);
    }
    }//  for (int i=0;i<217;++i)

  for (int i=0;i<173;++i)
    {
      if (HF_occupancy_[i]>0)
    {
      h_HFOccupancy->setBinContent(i+1,HF_occupancy_[i]);
    }
      if (HF_occupancy_thresh_[i]>0)
    {
      h_HFThreshOccupancy->setBinContent(i+1,HF_occupancy_thresh_[i]);
    }
    }//  for (int i=0;i<173;++i)

  //zeroCounters();

  if (debug_>0)
    std::cout <<"<HcalRecHitMonitor::fill_Nevents> FILLED REC HIT CELL PLOTS"<<std::endl;

} // void HcalRecHitMonitor::fill_Nevents(void)

void HcalRecHitMonitor::processEvent	(	const HBHERecHitCollection &	hbHits,
		const HORecHitCollection &	hoHits,
		const HFRecHitCollection &	hfHits,
		int	BCN,
		const edm::Event &	iEvent
	)

Definition at line 623 of file HcalRecHitMonitor.cc.

References gather_cfg::cout, HcalBaseDQMonitor::debug_, MonitorElement::Fill(), edm::Event::getByLabel(), h_TriggeredEvents, HcalHLTBits_, hltresultsLabel_, i, gen::k, MinBiasHLTBits_, processEvent_rechit(), edm::TriggerNames::size(), edm::TriggerNames::triggerName(), and edm::Event::triggerNames().

Referenced by analyze().

{

  
  if (debug_>1) std::cout <<"<HcalRecHitMonitor::processEvent> Processing event..."<<std::endl;


  bool passedHcalHLT=false;
  bool passedMinBiasHLT=false;

  edm::Handle<edm::TriggerResults> hltRes;
  if (!(iEvent.getByLabel(hltresultsLabel_,hltRes)))
    {
      if (debug_>0) edm::LogWarning("HcalRecHitMonitor")<<" Could not get HLT results with tag "<<hltresultsLabel_<<std::endl;
    }
  else
    {
      const edm::TriggerNames & triggerNames = iEvent.triggerNames(*hltRes);
      const unsigned int nTrig(triggerNames.size());
      for (unsigned int i=0;i<nTrig;++i)
    {
      // trigger decision is based on 'OR' of any specified trigger names
      for (unsigned int k=0;k<HcalHLTBits_.size();++k)
        {
          // if (triggerNames.triggerName(i)==HcalHLTBits_[k] && hltRes->accept(i))
          if (triggerNames.triggerName(i).find(HcalHLTBits_[k])!=std::string::npos && hltRes->accept(i))
        { 
          passedHcalHLT=true;
          break;
        }
        }
      // repeat for minbias triggers
      for (unsigned int k=0;k<MinBiasHLTBits_.size();++k)
        {
          // if (triggerNames.triggerName(i)==MinBiasHLTBits_[k] && hltRes->accept(i))      
          if (triggerNames.triggerName(i).find(MinBiasHLTBits_[k])!=std::string::npos && hltRes->accept(i))
        { 
          passedMinBiasHLT=true;
          break;
        }
        }
    }
    } //else

  if (debug_>2 && passedHcalHLT)  std::cout <<"\t<HcalRecHitMonitor::processEvent>  Passed Hcal HLT trigger  "<<std::endl;
  if (debug_>2 && passedMinBiasHLT)  std::cout <<"\t<HcalRecHitMonitor::processEvent>  Passed MinBias HLT trigger  "<<std::endl;
  
  h_TriggeredEvents->Fill(0); // all events
  if (passedMinBiasHLT) h_TriggeredEvents->Fill(1); // Minbias;
  if (passedHcalHLT)    h_TriggeredEvents->Fill(2); // hcal HLT
  processEvent_rechit(hbHits, hoHits, hfHits,passedHcalHLT,passedMinBiasHLT,BCN);

  return;
} // void HcalRecHitMonitor::processEvent(...)

void HcalRecHitMonitor::processEvent_rechit	(	const HBHERecHitCollection &	hbheHits,
		const HORecHitCollection &	hoHits,
		const HFRecHitCollection &	hfHits,
		bool	passedHcalHLT,
		bool	passedMinBiasHLT,
		int	BCN
	)

Definition at line 687 of file HcalRecHitMonitor.cc.

References abs, edm::SortedCollection< T, SORT >::begin(), CalcEtaBin(), gather_cfg::cout, HcalBaseDQMonitor::currentLS, HcalBaseDQMonitor::debug_, EtaPhiHists::depth, edm::SortedCollection< T, SORT >::end(), energy2_, energy2_thresh_, energy_, energy_thresh_, f, MonitorElement::Fill(), HcalLogicalMap::getHcalFrontEndId(), h_FlagMap_DIGITIME, h_FlagMap_HPDMULT, h_FlagMap_LONGSHORT, h_FlagMap_PULSESHAPE, h_FlagMap_TIMEADD, h_FlagMap_TIMEERROR, h_FlagMap_TIMESUBTRACT, h_HBflagcounter, h_HBHE_FlagCorr, h_HBHEHPDMult_vs_LS, h_HBHEPulseShape_vs_LS, h_HBM_weightedTime, h_HBP_weightedTime, h_HBsizeVsLS, h_HBThreshTime, h_HBTime, h_HBTimeVsEnergy, h_HE_HcalHLT_energydifference, h_HE_HcalHLT_weightedtimedifference, h_HEenergydifference, h_HEflagcounter, h_HEM_weightedTime, h_HEP_weightedTime, h_HEsizeVsLS, h_HEThreshTime, h_HETime, h_HEtimedifference, h_HETimeVsEnergy, h_HF_FlagCorr, h_HF_HcalHLT_energydifference, h_HF_HcalHLT_weightedtimedifference, h_HFDigiTime_vs_LS, h_HFenergydifference, h_HFflagcounter, h_HFLongShort_vs_LS, h_HFM_weightedTime, h_HFP_weightedTime, h_HFsizeVsLS, h_HFThreshTime, h_HFTime, h_HFtimedifference, h_HFTimeVsEnergy, h_HOflagcounter, h_HOsizeVsLS, h_HOThreshTime, h_HOTime, h_HOTimeVsEnergy, h_LumiPlot_BX_HcalHLTEvents, h_LumiPlot_BX_HcalHLTEvents_notimecut, h_LumiPlot_BX_MinBiasEvents, h_LumiPlot_BX_MinBiasEvents_notimecut, h_LumiPlot_LS_HcalHLTEvents, h_LumiPlot_LS_HcalHLTEvents_notimecut, h_LumiPlot_LS_MinBiasEvents, h_LumiPlot_LS_MinBiasEvents_notimecut, h_LumiPlot_MinTime_vs_MinHT, h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus, h_LumiPlot_SumHT_HFPlus_vs_HFMinus, h_LumiPlot_timeHFPlus_vs_timeHFMinus, h_LumiPlot_timeHT_HFM, h_LumiPlot_timeHT_HFP, h_rechitieta_05, h_rechitieta_10, h_rechitieta_100, h_rechitieta_25, h_rechitiphi_05, h_rechitiphi_10, h_rechitiphi_100, h_rechitiphi_25, HB_occupancy_, HB_occupancy_thresh_, HBenergyThreshold_, HBETThreshold_, HBflagcounter_, HcalCaloFlagLabels::HBHEHpdHitMultiplicity, HcalCaloFlagLabels::HBHEPulseShape, HBtime_, HBtime_thresh_, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, HE_occupancy_, HE_occupancy_thresh_, HEenergyThreshold_, HEETThreshold_, HEflagcounter_, HEtime_, HEtime_thresh_, HF_occupancy_, HF_occupancy_thresh_, HcalCaloFlagLabels::HFDigiTime, HFenergyThreshold_, HFETThreshold_, HFflagcounter_, HcalCaloFlagLabels::HFLongShort, HFP_HFM_Energy, HFtime_, HFtime_thresh_, HO_occupancy_, HO_occupancy_thresh_, HOenergyThreshold_, HOETThreshold_, HOflagcounter_, HOtime_, HOtime_thresh_, i, HcalDetId::ieta(), HcalBaseDQMonitor::logicalMap_, min, occupancy_, occupancy_thresh_, OccupancyByDepth, OccupancyThreshByDepth, funct::pow(), RECHITMON_TIME_MAX, RECHITMON_TIME_MIN, submit::rm, edm::SortedCollection< T, SORT >::size(), SqrtSumEnergy2ByDepth, SumEnergyByDepth, SumTimeByDepth, theHBHEEtaBounds, theHFEtaBounds, time_, time_thresh_, timediffThresh_, HcalCaloFlagLabels::TimingAddedBit, HcalCaloFlagLabels::TimingErrorBit, HcalCaloFlagLabels::TimingSubtractedBit, and MonitorElement::update().

Referenced by processEvent().

{
  // Gather rechit info
  
  //const float area[]={0.111,0.175,0.175,0.175,0.175,0.175,0.174,0.178,0.172,0.175,0.178,0.346,0.604};

  if (debug_>1) std::cout <<"<HcalRecHitMonitor::processEvent_rechitenergy> Processing rechits..."<<std::endl;
  
  // loop over HBHE
  
  int     hbocc=0;
  int     heocc=0;
  int     hboccthresh=0;
  int     heoccthresh=0;

  double HtPlus =0, HtMinus=0;
  double HFePlus=0, HFeMinus=0;
  double HBePlus=0, HBeMinus=0;
  double HEePlus=0, HEeMinus=0;
  double HFtPlus=0, HFtMinus=0;
  double HBtPlus=0, HBtMinus=0;
  double HEtPlus=0, HEtMinus=0;

  int hbpocc=0, hbmocc=0, hepocc=0, hemocc=0, hfpocc=0, hfmocc=0;

  for (unsigned int i=0;i<4;++i)
    {
      OccupancyByDepth.depth[i]->update();
      OccupancyThreshByDepth.depth[i]->update();
      SumEnergyByDepth.depth[i]->update();
      SqrtSumEnergy2ByDepth.depth[i]->update();
      SumTimeByDepth.depth[i]->update();
    }
    
  h_HBflagcounter->update();
  h_HEflagcounter->update();
  h_HFflagcounter->update();
  h_HOflagcounter->update();
  
  
  for (HBHERecHitCollection::const_iterator HBHEiter=hbheHits.begin(); HBHEiter!=hbheHits.end(); ++HBHEiter) 
    { // loop over all hits
      float en = HBHEiter->energy();
      float ti = HBHEiter->time();
      HcalDetId id(HBHEiter->detid().rawId());
      int ieta = id.ieta();
      int iphi = id.iphi();
      int depth = id.depth();

      if (en>0.5)
        {
          h_rechitieta_05->Fill(ieta);
          h_rechitiphi_05->Fill(iphi);
          if (en>1.)
            {
              h_rechitieta_10->Fill(ieta);
              h_rechitiphi_10->Fill(iphi);
              if (en>2.5)
                {
                  h_rechitieta_25->Fill(ieta);
                  h_rechitiphi_25->Fill(iphi);
                  if (en>10.)
                    {
                      h_rechitieta_100->Fill(ieta);
                      h_rechitiphi_100->Fill(iphi);
                    }
                }
            }
        }



      HcalSubdetector subdet = id.subdet();
      double fEta=fabs(0.5*(theHBHEEtaBounds[abs(ieta)-1]+theHBHEEtaBounds[abs(ieta)]));

      int calcEta = CalcEtaBin(subdet,ieta,depth);
      int rbxindex=logicalMap_->getHcalFrontEndId(HBHEiter->detid()).rbxIndex();
      int rm= logicalMap_->getHcalFrontEndId(HBHEiter->detid()).rm();

      // Fill HBHE flag plots
      h_HBHE_FlagCorr->Fill(HBHEiter->flagField(HcalCaloFlagLabels::HBHEPulseShape),
                HBHEiter->flagField(HcalCaloFlagLabels::HBHEHpdHitMultiplicity)); 

      if (HBHEiter->flagField(HcalCaloFlagLabels::HBHEHpdHitMultiplicity))
    {
      h_FlagMap_HPDMULT->Fill(rbxindex,rm);
      h_HBHEHPDMult_vs_LS->Fill(currentLS,1);
    }
      if (HBHEiter->flagField(HcalCaloFlagLabels::HBHEPulseShape))
    {
      h_FlagMap_PULSESHAPE->Fill(rbxindex,rm);
      h_HBHEPulseShape_vs_LS->Fill(currentLS,1);
    }
      if (HBHEiter->flagField(HcalCaloFlagLabels::TimingSubtractedBit))
    h_FlagMap_TIMESUBTRACT->Fill(rbxindex,rm);
      else if (HBHEiter->flagField(HcalCaloFlagLabels::TimingAddedBit))
    h_FlagMap_TIMEADD->Fill(rbxindex,rm);
      else if (HBHEiter->flagField(HcalCaloFlagLabels::TimingErrorBit))
    h_FlagMap_TIMEERROR->Fill(rbxindex,rm);

      if (subdet==HcalBarrel)
    {
      if (en>HBenergyThreshold_)
        h_HBTimeVsEnergy->Fill(en,ti);
      //Looping over HB searching for flags --- cris
      for (int f=0;f<32;f++)
        {
          // Let's display HSCP just to see if these bits are set
          /*
           if (f == HcalCaloFlagLabels::HSCP_R1R2)         continue;
               if (f == HcalCaloFlagLabels::HSCP_FracLeader)   continue;
               if (f == HcalCaloFlagLabels::HSCP_OuterEnergy)  continue;
               if (f == HcalCaloFlagLabels::HSCP_ExpFit)       continue;
          */
          if (HBHEiter->flagField(f))
        ++HBflagcounter_[f];
        }
      ++occupancy_[calcEta][iphi-1][depth-1];
      energy_[calcEta][iphi-1][depth-1]+=en;
          energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
      time_[calcEta][iphi-1][depth-1]+=ti;
      if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
        h_HBTime->Fill(ti);
      else
        ++HBtime_[int(ti-RECHITMON_TIME_MIN)];
      ++hbocc; 

      // Threshold plots;  require E> threshold and minbias trigger
      if (
          en>=HBenergyThreshold_ && 
          en/cosh(fEta)>=HBETThreshold_ 
          ) 
        {
          if (passedMinBiasHLT==true)
        {
          ++occupancy_thresh_[calcEta][iphi-1][depth-1];
          energy_thresh_[calcEta][iphi-1][depth-1]+=en;
          energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
          time_thresh_[calcEta][iphi-1][depth-1]+=ti;
          
          ++hboccthresh;
          if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
            h_HBThreshTime->Fill(ti);
          else
            ++HBtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
        }

          if (ieta>0)
        {
          HBePlus+=en;
          HBtPlus+=ti*en;
          hbpocc++;
        }
          else
        {
          HBeMinus+=en;
          HBtMinus+=ti*en;
          hbmocc++;
        }
        } // if (HB  en>thresh, ET>thresh)
    } // if (id.subdet()==HcalBarrel)

      else if (subdet==HcalEndcap)
    {
      if (en>HEenergyThreshold_)
        h_HETimeVsEnergy->Fill(en,ti);
      //Looping over HE searching for flags --- cris
      for (int f=0;f<32;f++)
            {
              if (HBHEiter->flagField(f))
                ++HEflagcounter_[f];
            }

      ++occupancy_[calcEta][iphi-1][depth-1];
      energy_[calcEta][iphi-1][depth-1]+=en;
          energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
      time_[calcEta][iphi-1][depth-1]+=ti;

      ++heocc;
      if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
        h_HETime->Fill(ti);
      else
        ++HEtime_[int(ti-RECHITMON_TIME_MIN)];

      // Threshold plots require e>E_thresh, ET>ET_thresh
      if (en>=HEenergyThreshold_
          && en/cosh(fEta)>=HEETThreshold_
          )
        {
          // occupancy plots also require passedMinBiasHLT
          if (passedMinBiasHLT==true)
        {
          ++occupancy_thresh_[calcEta][iphi-1][depth-1];
          energy_thresh_[calcEta][iphi-1][depth-1]+=en;
          energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
          time_thresh_[calcEta][iphi-1][depth-1]+=ti;
          ++heoccthresh;
          if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
              h_HEThreshTime->Fill(ti);
          else
            ++HEtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
        }
          // ePlus, tPlus calculated regardless of trigger
          if (ieta>0)
        {
          HEePlus+=en;
          HEtPlus+=ti*en;
          hepocc++;
        }
          else
        {
          HEeMinus+=en;
          HEtMinus+=ti*en;
          hemocc++;
        }
        } // if (en>=HEenergyThreshold_ && ET>threshold)

    } // else if (id.subdet()==HcalEndcap)
     
    } //for (HBHERecHitCollection::const_iterator HBHEiter=...)
  
  // Calculate normalized time
  HEePlus>0  ?  HEtPlus/=HEePlus   :  HEtPlus=10000;
  HEeMinus>0 ?  HEtMinus/=HEeMinus :  HEtMinus=-10000;
  HBePlus>0  ?  HBtPlus/=HBePlus   :  HBtPlus=10000;
  HBeMinus>0 ?  HBtMinus/=HBeMinus :  HBtMinus=-10000;
  
  ++HB_occupancy_[hbocc/10];
  ++HE_occupancy_[heocc/10];
  ++HB_occupancy_thresh_[hboccthresh/10];
  ++HE_occupancy_thresh_[heoccthresh/10];
  h_HBsizeVsLS->Fill(currentLS,hbocc);
  h_HEsizeVsLS->Fill(currentLS,heocc);

  // loop over HO

  h_HOsizeVsLS->Fill(currentLS,hoHits.size());
  int hoocc=0;
  int hooccthresh=0;
  for (HORecHitCollection::const_iterator HOiter=hoHits.begin(); HOiter!=hoHits.end(); ++HOiter) 
    { // loop over all hits
      float en = HOiter->energy();
      float ti = HOiter->time();
      if (en>HOenergyThreshold_)
    h_HOTimeVsEnergy->Fill(en,ti);

      HcalDetId id(HOiter->detid().rawId());
      int ieta = id.ieta();
      int iphi = id.iphi();
      int depth = id.depth();

      if (en>0.5)
        {
          h_rechitieta_05->Fill(ieta);
          h_rechitiphi_05->Fill(iphi);
          if (en>1.)
            {
              h_rechitieta_10->Fill(ieta);
              h_rechitiphi_10->Fill(iphi);
              if (en>2.5)
                {
                  h_rechitieta_25->Fill(ieta);
                  h_rechitiphi_25->Fill(iphi);
                  if (en>10.)
                    {
                      h_rechitieta_100->Fill(ieta);
                      h_rechitiphi_100->Fill(iphi);
                    }
                }
            }
        }



      int calcEta = CalcEtaBin(HcalOuter,ieta,depth);
      double fEta=fabs(0.5*(theHBHEEtaBounds[abs(ieta)-1]+theHBHEEtaBounds[abs(ieta)]));
      
      int rbxindex=logicalMap_->getHcalFrontEndId(HOiter->detid()).rbxIndex();
      int rm= logicalMap_->getHcalFrontEndId(HOiter->detid()).rm();
      
      if (HOiter->flagField(HcalCaloFlagLabels::TimingSubtractedBit))
    h_FlagMap_TIMESUBTRACT->Fill(rbxindex,rm);
      else if (HOiter->flagField(HcalCaloFlagLabels::TimingAddedBit))
    h_FlagMap_TIMEADD->Fill(rbxindex,rm);
      else if (HOiter->flagField(HcalCaloFlagLabels::TimingErrorBit))
    h_FlagMap_TIMEERROR->Fill(rbxindex,rm);
      
      
      //Looping over HO searching for flags --- cris
      for (int f=0;f<32;f++)
    {
      if (HOiter->flagField(f))
        HOflagcounter_[f]++;
    }
      
      ++occupancy_[calcEta][iphi-1][depth-1];
      energy_[calcEta][iphi-1][depth-1]+=en;
      energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
      time_[calcEta][iphi-1][depth-1]+=ti;
      ++hoocc;
      if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
    h_HOTime->Fill(ti);
      else
    ++HOtime_[int(ti-RECHITMON_TIME_MIN)];

      // We don't calculate HOplus/HOminus values (independent of trigger), so require min bias trigger 
      // along with E, ET thresholds directly in this HO loop:

      if (en>=HOenergyThreshold_  
      && en/cosh(fEta)>=HOETThreshold_
      && passedMinBiasHLT==true   
      )
    {
      ++occupancy_thresh_[calcEta][iphi-1][depth-1];
      energy_thresh_[calcEta][iphi-1][depth-1]+=en;
      energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
      time_thresh_[calcEta][iphi-1][depth-1]+=ti;

      ++hooccthresh;
      if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
        h_HOThreshTime->Fill(ti);
      else
        ++HOtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
    } 
    } // loop over all HO hits

  ++HO_occupancy_[hoocc/10];
  ++HO_occupancy_thresh_[hooccthresh/10];
   
  // loop over HF
  h_HFsizeVsLS->Fill(currentLS,hfHits.size());

  HtPlus=0; HtMinus=0;
  
  int hfocc=0;
  int hfoccthresh=0;
  for (HFRecHitCollection::const_iterator HFiter=hfHits.begin(); HFiter!=hfHits.end(); ++HFiter) 
    { // loop over all hits
      float en = HFiter->energy();
      float ti = HFiter->time();
      if (en> HFenergyThreshold_)
    h_HFTimeVsEnergy->Fill(en,ti);

      HcalDetId id(HFiter->detid().rawId());
      int ieta = id.ieta();
      int iphi = id.iphi();
      int depth = id.depth();

      if (en>0.5)
    {
      h_rechitieta_05->Fill(ieta);
      h_rechitiphi_05->Fill(iphi);
      if (en>1.)
        {
          h_rechitieta_10->Fill(ieta);
          h_rechitiphi_10->Fill(iphi);
          if (en>2.5)
        {
          h_rechitieta_25->Fill(ieta);
          h_rechitiphi_25->Fill(iphi);
          if (en>10.)
            {
              h_rechitieta_100->Fill(ieta);
              h_rechitiphi_100->Fill(iphi);
            }
        }
        }
    }

      double fEta=fabs(0.5*(theHFEtaBounds[abs(ieta)-29]+theHFEtaBounds[abs(ieta)-28]));
      int calcEta = CalcEtaBin(HcalForward,ieta,depth);

      int rbxindex=logicalMap_->getHcalFrontEndId(HFiter->detid()).rbxIndex();
      int rm= logicalMap_->getHcalFrontEndId(HFiter->detid()).rm(); 
     
      h_HF_FlagCorr->Fill(HFiter->flagField(HcalCaloFlagLabels::HFDigiTime),HFiter->flagField(HcalCaloFlagLabels::HFLongShort)); 
      if (HFiter->flagField(HcalCaloFlagLabels::TimingSubtractedBit))
    h_FlagMap_TIMESUBTRACT->Fill(rbxindex,rm);
      else if (HFiter->flagField(HcalCaloFlagLabels::TimingAddedBit))
    h_FlagMap_TIMEADD->Fill(rbxindex,rm);
      else if (HFiter->flagField(HcalCaloFlagLabels::TimingErrorBit))
    h_FlagMap_TIMEERROR->Fill(rbxindex,rm);

      if (HFiter->flagField(HcalCaloFlagLabels::HFDigiTime))
    {
      h_FlagMap_DIGITIME->Fill(rbxindex,rm);
      h_HFDigiTime_vs_LS->Fill(currentLS,1);
    }
      if (HFiter->flagField(HcalCaloFlagLabels::HFLongShort))
    {
      h_FlagMap_LONGSHORT->Fill(rbxindex,rm);
      h_HFLongShort_vs_LS->Fill(currentLS,1);
    }
      //Looping over HF searching for flags --- cris
      for (int f=0;f<32;f++)
    {
      if (HFiter->flagField(f))
        HFflagcounter_[f]++;
    }

      // Occupancy plots, without threshold
      ++occupancy_[calcEta][iphi-1][depth-1];
      energy_[calcEta][iphi-1][depth-1]+=en;
      energy2_[calcEta][iphi-1][depth-1]+=pow(en,2);
      time_[calcEta][iphi-1][depth-1]+=ti;
      ++hfocc;
      if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
    h_HFTime->Fill(ti);
      else
    ++HFtime_[int(ti-RECHITMON_TIME_MIN)];

      ieta>0 ?  HtPlus+=en/cosh(fEta)  : HtMinus+=en/cosh(fEta);  // add energy from all cells, or only those > threshold?

      if (en>=HFenergyThreshold_ && 
      en/cosh(fEta)>=HFETThreshold_ 
      )
    {
      // Occupancy plots require min bias trigger, along with thresholds exceeded
      if (passedMinBiasHLT)
        {
          ++occupancy_thresh_[calcEta][iphi-1][depth-1];
          energy_thresh_[calcEta][iphi-1][depth-1]+=en;
          energy2_thresh_[calcEta][iphi-1][depth-1]+=pow(en,2);
          time_thresh_[calcEta][iphi-1][depth-1]+=ti;
          
          ++hfoccthresh;
          if (ti<RECHITMON_TIME_MIN || ti>RECHITMON_TIME_MAX)
        h_HFThreshTime->Fill(ti);
          else
        ++HFtime_thresh_[int(ti-RECHITMON_TIME_MIN)];
        }

      if (ieta>0)
        {
          HFtPlus+=en*ti;
          HFePlus+=en;
          hfpocc++;
        }
      else if (ieta<0)
        {
          HFtMinus+=en*ti;
          HFeMinus+=en;
          hfmocc++;
        }
    } // if (en>thresh, ET>thresh)
    } // loop over all HF hits
     
  ++HF_occupancy_[hfocc/10];
  ++HF_occupancy_thresh_[hfoccthresh/10];

  
  // Form event-wide variables (time averages, etc.), and plot them

  // Calculate weighted times.  (Set tPlus, tMinus to overflow in case where total energy < 0)
  HFePlus>0  ? HFtPlus/=HFePlus    : HFtPlus  =  10000;
  HFeMinus>0 ? HFtMinus/=HFeMinus  : HFtMinus = -10000;
     
  double mintime=99;  // used to be min(tPlus,tMinus);
  double minHT=std::min(HtMinus,HtPlus);
  minHT==HtMinus ?  mintime=HFtMinus : mintime = HFtPlus;
  //mintime = min(HFtPlus,HFtMinus); // I think we might want to use this value for mintime?


  h_LumiPlot_MinTime_vs_MinHT->Fill(minHT, mintime);
  h_LumiPlot_timeHT_HFM->Fill(HtMinus,HFtMinus);
  h_LumiPlot_timeHT_HFP->Fill(HtPlus,HFtPlus);

  if (passedMinBiasHLT==true)
    {
      h_LumiPlot_SumHT_HFPlus_vs_HFMinus->Fill(HtMinus,HtPlus);
      // HtMinus, HtPlus require no energy cuts for their contributing cells
      // HFeMinus, HFePlus require that cells be > threshold cut
      
      if (HtMinus>1 && HtPlus > 1) // is this the condition we want, or do we want hfmocc>0 && hfpocc >0?
        {
      h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus->Fill(HFeMinus,HFePlus);
      h_LumiPlot_timeHFPlus_vs_timeHFMinus->Fill(HFtMinus,HFtPlus);

      h_HFP_weightedTime->Fill(HFtPlus);
      h_HFM_weightedTime->Fill(HFtMinus);
      h_HBP_weightedTime->Fill(HBtPlus);
      h_HBM_weightedTime->Fill(HBtMinus);
      
      h_HEP_weightedTime->Fill(HEtPlus);
      h_HEM_weightedTime->Fill(HEtMinus);

      if (hepocc>0 && hemocc>0)
        {
          h_HEtimedifference->Fill(HEtPlus-HEtMinus);
          if (HEePlus-HEeMinus!=0) h_HEenergydifference->Fill((HEePlus-HEeMinus)/(HEePlus+HEeMinus));
        }
      if (hfpocc>0 && hfmocc>0)  // which condition do we want?
        {
          h_HFtimedifference->Fill((HFtPlus)-(HFtMinus));
          if (HFePlus+HFeMinus!=0) h_HFenergydifference->Fill((HFePlus-HFeMinus)/(HFePlus+HFeMinus));
        }

      h_LumiPlot_LS_MinBiasEvents_notimecut->Fill(currentLS);
      h_LumiPlot_BX_MinBiasEvents_notimecut->Fill(BCN);
      if (fabs(HFtPlus-HFtMinus)<timediffThresh_)
        {
          h_LumiPlot_LS_MinBiasEvents->Fill(currentLS);
          h_LumiPlot_BX_MinBiasEvents->Fill(BCN);
        }
      
      HFP_HFM_Energy->Fill(HFeMinus/1000., HFePlus/1000.);
    }
 
      if (debug_>1) std::cout <<"\t<HcalRecHitMonitor:: HF averages>  TPLUS = "<<HFtPlus<<"  EPLUS = "<<HFePlus<<"  TMINUS = "<<HFtMinus<<"  EMINUS = "<<HFeMinus<<"  Weighted Time Diff = "<<((HFtPlus)-(HFtMinus))<<std::endl;
      

    } // if (passedMinBiasHLT)

  if (passedHcalHLT && HtMinus>1 && HtPlus> 1 )
    {
      if (hfpocc>0 && hfmocc>0)
    {
      h_HF_HcalHLT_weightedtimedifference->Fill(HFtPlus-HFtMinus);
      if (HFePlus+HFeMinus!=0) h_HF_HcalHLT_energydifference->Fill((HFePlus-HFeMinus)/(HFePlus+HFeMinus));
    }
      if  (hepocc>0 && hemocc>0)
    {
      h_HE_HcalHLT_weightedtimedifference->Fill(HEtPlus-HEtMinus);
      if (HEePlus-HEeMinus!=0) h_HE_HcalHLT_energydifference->Fill((HEePlus-HEeMinus)/(HEePlus+HEeMinus));
    }

      h_LumiPlot_LS_HcalHLTEvents_notimecut->Fill(currentLS);
      h_LumiPlot_BX_HcalHLTEvents_notimecut->Fill(BCN);
      if (fabs(HFtPlus-HFtMinus)<timediffThresh_)
    {
      h_LumiPlot_LS_HcalHLTEvents->Fill(currentLS);
      h_LumiPlot_BX_HcalHLTEvents->Fill(BCN);
    }
     } // passsed Hcal HLT
   
 return;
} // void HcalRecHitMonitor::processEvent_rechitenergy

void HcalRecHitMonitor::reset ( void  )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 534 of file HcalRecHitMonitor.cc.

References HcalBaseDQMonitor::dbe_, MonitorElement::DQM_KIND_TH1F, MonitorElement::DQM_KIND_TH2F, MonitorElement::DQM_KIND_TPROFILE, DQMStore::getAllContents(), i, HcalObjRepresent::Reset(), and HcalBaseDQMonitor::subdir_.

Referenced by beginRun(), and setup().

{
  std::vector<MonitorElement*> hists = dbe_->getAllContents(subdir_);
  for (unsigned int i=0;i<hists.size();++i)
    {
      if (hists[i]->kind()==MonitorElement::DQM_KIND_TH1F ||
      hists[i]->kind()==MonitorElement::DQM_KIND_TH2F ||
      hists[i]->kind()==MonitorElement::DQM_KIND_TPROFILE)
    hists[i]->Reset();
    }
}  

void HcalRecHitMonitor::setup ( void  )

virtual

Reimplemented from HcalBaseDQMonitor.

Definition at line 70 of file HcalRecHitMonitor.cc.

References HcalCaloFlagLabels::ADCSaturationBit, DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookFloat(), DQMStore::bookProfile(), gather_cfg::cout, HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::debug_, MonitorElement::Fill(), HcalCaloFlagLabels::Fraction2TS, MonitorElement::getTH1F(), h_FlagMap_DIGITIME, h_FlagMap_HPDMULT, h_FlagMap_LONGSHORT, h_FlagMap_PULSESHAPE, h_FlagMap_TIMEADD, h_FlagMap_TIMEERROR, h_FlagMap_TIMESUBTRACT, h_HBflagcounter, h_HBHE_FlagCorr, h_HBHEHPDMult_vs_LS, h_HBHEPulseShape_vs_LS, h_HBM_weightedTime, h_HBOccupancy, h_HBP_weightedTime, h_HBsizeVsLS, h_HBThreshOccupancy, h_HBThreshTime, h_HBTime, h_HBTimeVsEnergy, h_HE_HcalHLT_energydifference, h_HE_HcalHLT_weightedtimedifference, h_HEenergydifference, h_HEflagcounter, h_HEM_weightedTime, h_HEOccupancy, h_HEP_weightedTime, h_HEsizeVsLS, h_HEThreshOccupancy, h_HEThreshTime, h_HETime, h_HEtimedifference, h_HETimeVsEnergy, h_HF_FlagCorr, h_HF_HcalHLT_energydifference, h_HF_HcalHLT_weightedtimedifference, h_HFDigiTime_vs_LS, h_HFenergydifference, h_HFflagcounter, h_HFLongShort_vs_LS, h_HFM_weightedTime, h_HFOccupancy, h_HFP_weightedTime, h_HFsizeVsLS, h_HFThreshOccupancy, h_HFThreshTime, h_HFTime, h_HFtimedifference, h_HFTimeVsEnergy, h_HOflagcounter, h_HOOccupancy, h_HOsizeVsLS, h_HOThreshOccupancy, h_HOThreshTime, h_HOTime, h_HOTimeVsEnergy, h_LumiPlot_BX_allevents, h_LumiPlot_BX_HcalHLTEvents, h_LumiPlot_BX_HcalHLTEvents_notimecut, h_LumiPlot_BX_MinBiasEvents, h_LumiPlot_BX_MinBiasEvents_notimecut, h_LumiPlot_LS_allevents, h_LumiPlot_LS_HcalHLTEvents, h_LumiPlot_LS_HcalHLTEvents_notimecut, h_LumiPlot_LS_MinBiasEvents, h_LumiPlot_LS_MinBiasEvents_notimecut, h_LumiPlot_MinTime_vs_MinHT, h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus, h_LumiPlot_SumHT_HFPlus_vs_HFMinus, h_LumiPlot_timeHFPlus_vs_timeHFMinus, h_LumiPlot_timeHT_HFM, h_LumiPlot_timeHT_HFP, h_rechitieta, h_rechitieta_05, h_rechitieta_10, h_rechitieta_100, h_rechitieta_25, h_rechitieta_thresh, h_rechitiphi, h_rechitiphi_05, h_rechitiphi_10, h_rechitiphi_100, h_rechitiphi_25, h_rechitiphi_thresh, h_TriggeredEvents, HBenergyThreshold_, HBETThreshold_, HcalCaloFlagLabels::HBHEHpdHitMultiplicity, HcalCaloFlagLabels::HBHEPulseShape, HcalCaloFlagLabels::HBHETimingShapedCutsBits, HcalCaloFlagLabels::HBHETimingTrustBits, HEenergyThreshold_, HEETThreshold_, HcalCaloFlagLabels::HFDigiTime, HFenergyThreshold_, HFETThreshold_, HcalCaloFlagLabels::HFLongShort, HFP_HFM_Energy, HcalCaloFlagLabels::HFTimingTrustBits, HOenergyThreshold_, HOETThreshold_, HcalCaloFlagLabels::HSCP_ExpFit, HcalCaloFlagLabels::HSCP_FracLeader, HcalCaloFlagLabels::HSCP_OuterEnergy, HcalCaloFlagLabels::HSCP_R1R2, HcalBaseDQMonitor::NLumiBlocks_, OccupancyByDepth, OccupancyThreshByDepth, RECHITMON_TIME_MAX, RECHITMON_TIME_MIN, reset(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), HcalBaseDQMonitor::setup(), setupDone_, HcalBaseDQMonitor::SetupEtaPhiHists(), SqrtSumEnergy2ByDepth, SqrtSumEnergy2ThreshByDepth, HcalBaseDQMonitor::subdir_, SumEnergyByDepth, SumEnergyThreshByDepth, SumTimeByDepth, SumTimeThreshByDepth, timediffThresh_, HcalCaloFlagLabels::TimingAddedBit, HcalCaloFlagLabels::TimingErrorBit, HcalCaloFlagLabels::TimingSubtractedBit, and zeroCounters().

Referenced by beginRun().

{

  if (setupDone_)
  {
    // Always do a zeroing/resetting so that empty histograms/counter
    // will always appear.
  
    // clear all counters, reset histograms
    this->zeroCounters();
    this->reset();
    return;
  }
  else
    setupDone_=true;

  HcalBaseDQMonitor::setup();


  if (debug_>0)
    std::cout <<"<HcalRecHitMonitor::setup>  Setting up histograms"<<std::endl;

  // RecHit Monitor - specific cfg variables

  if (debug_>1)
    std::cout <<"<HcalRecHitMonitor::setup>  Creating Histograms"<<std::endl;

  dbe_->setCurrentFolder(subdir_);
  h_TriggeredEvents=dbe_->book1D("EventTriggers","EventTriggers",3,-0.5,2.5);
  h_TriggeredEvents->setBinLabel(1,"AllEvents");
  h_TriggeredEvents->setBinLabel(2,"HLT_Minbias");
  h_TriggeredEvents->setBinLabel(3,"HLT_Hcal");

  dbe_->setCurrentFolder(subdir_+"rechit_parameters");
  MonitorElement* THR;
  dbe_->setCurrentFolder(subdir_+"rechit_parameters/thresholds");
  THR=dbe_->bookFloat("HB_Rechit_Energy_Threshold");
  THR->Fill(HBenergyThreshold_);
  THR=dbe_->bookFloat("HE_Rechit_Energy_Threshold");
  THR->Fill(HEenergyThreshold_);
  THR=dbe_->bookFloat("HO_Rechit_Energy_Threshold");
  THR->Fill(HOenergyThreshold_);
  THR=dbe_->bookFloat("HF_Rechit_Energy_Threshold");
  THR->Fill(HFenergyThreshold_);
  THR=dbe_->bookFloat("HB_Rechit_ET_Threshold");
  THR->Fill(HBETThreshold_);
  THR=dbe_->bookFloat("HE_Rechit_ET_Threshold");
  THR->Fill(HEETThreshold_);
  THR=dbe_->bookFloat("HO_Rechit_ET_Threshold");
  THR->Fill(HOETThreshold_);
  THR=dbe_->bookFloat("HF_Rechit_ET_Threshold");
  THR->Fill(HFETThreshold_);
  THR=dbe_->bookFloat("Maximum_HFM_HFP_time_difference_for_luminosityplots");
  THR->Fill(timediffThresh_);

  
  // Set up histograms that are filled by all rechits
  dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits");
  SetupEtaPhiHists(OccupancyByDepth,"RecHit Occupancy","");
  h_rechitieta = dbe_->book1D("HcalRecHitIeta",
                  "Hcal RecHit ieta",
                  83,-41.5,41.5);
  h_rechitiphi = dbe_->book1D("HcalRecHitIphi",
                  "Hcal RecHit iphi",
                  72,0.5,72.5);

  h_rechitieta_05 = dbe_->book1D("HcalRecHitIeta05",
                 "Hcal RecHit ieta E>0.5 GeV",
                 83,-41.5,41.5);
  h_rechitiphi_05 = dbe_->book1D("HcalRecHitIphi05",
                 "Hcal RecHit iphi E>0.5 GeV",
                 72,0.5,72.5);
  h_rechitieta_10 = dbe_->book1D("HcalRecHitIeta10",
                 "Hcal RecHit ieta E>1.0 GeV",
                 83,-41.5,41.5);
  h_rechitiphi_10 = dbe_->book1D("HcalRecHitIphi10",
                 "Hcal RecHit iphi E>1.0 GeV",
                 72,0.5,72.5);
  h_rechitieta_25 = dbe_->book1D("HcalRecHitIeta25",
                 "Hcal RecHit ieta E>2.5 GeV",
                 83,-41.5,41.5);
  h_rechitiphi_25 = dbe_->book1D("HcalRecHitIphi25",
                 "Hcal RecHit iphi E>2.5 GeV",
                 72,0.5,72.5);
  h_rechitieta_100 = dbe_->book1D("HcalRecHitIeta100",
                  "Hcal RecHit ieta E>10.0 GeV",
                  83,-41.5,41.5);
  h_rechitiphi_100 = dbe_->book1D("HcalRecHitIphi100",
                  "Hcal RecHit iphi E>10.0 GeV",
                  72,0.5,72.5);



  h_LumiPlot_LS_allevents = dbe_->book1D("AllEventsPerLS",
                     "LS # of all events",
                     NLumiBlocks_,0.5,NLumiBlocks_+0.5);
  h_LumiPlot_BX_allevents = dbe_->book1D("BX_allevents",
                     "BX # of all events",
                     3600,0,3600);
  h_LumiPlot_MinTime_vs_MinHT = dbe_->book2D("MinTime_vs_MinSumET",
                         "Energy-Weighted Time vs Min (HF+,HF-) Scalar Sum ET;min Sum ET(GeV);time(ns)",
                         100,0,10,80,-40,40);

  h_LumiPlot_timeHT_HFM = dbe_->book2D("HFM_Time_vs_SumET",
                       "Energy-Weighted Time vs HFMinus Scalar Sum ET;Sum ET(GeV);time(ns)",
                       100,0,10,80,-40,40);
  
  h_LumiPlot_timeHT_HFP = dbe_->book2D("HFP_Time_vs_SumET",
                       "Energy-Weighted Time vs HFPlus Scalar Sum ET;Sum ET(GeV);time(ns)",
                       100,0,10,80,-40,40);


  dbe_->setCurrentFolder(subdir_+"Distributions_AllRecHits/sumplots");
  SetupEtaPhiHists(SumEnergyByDepth,"RecHit Summed Energy","GeV");
  SetupEtaPhiHists(SqrtSumEnergy2ByDepth,"RecHit Sqrt Summed Energy2","GeV");
  SetupEtaPhiHists(SumTimeByDepth,"RecHit Summed Time","nS");

  // Histograms for events that passed MinBias triggers
  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias");

  h_HBP_weightedTime = dbe_->book1D("WeightedTime_HBP","Weighted Time for HBP",
                    300,-150,150);
  h_HBM_weightedTime = dbe_->book1D("WeightedTime_HBM","Weighted Time for HBM",
                    300,-150,150);
  h_HEP_weightedTime = dbe_->book1D("WeightedTime_HEP","Weighted Time for HEP",
                    300,-150,150);
  h_HEM_weightedTime = dbe_->book1D("WeightedTime_HEM","Weighted Time for HEM",
                    300,-150,150);
  h_HFP_weightedTime = dbe_->book1D("WeightedTime_HFP","Weighted Time for HFP",
                    300,-150,150);
  h_HFM_weightedTime = dbe_->book1D("WeightedTime_HFM","Weighted Time for HFM",
                    300,-150,150);

  h_HFtimedifference = dbe_->book1D("HFweightedtimeDifference",
                    "Energy-Weighted time difference between HF+ and HF- passing MinBias (no HT cut)",
                    251,-250.5,250.5);
  h_HEtimedifference = dbe_->book1D("HEweightedtimeDifference",
                    "Energy-Weighted time difference between HE+ and HE- passing MinBias (no HT cut)",
                    251,-250.5,250.5);
  
  HFP_HFM_Energy = dbe_->book2D("HFP_HFM_Energy",
                "HFP VS HFM Energy; Total Energy in HFMinus (TeV); Total Energy in HFPlus (TeV)",
                100,0,100, 100,0,100);
  
  // Would these work better as 2D plots?
  h_HFenergydifference = dbe_->book1D("HFenergyDifference",
                      "Sum(E_HFPlus - E_HFMinus)/Sum(E_HFPlus + E_HFMinus)",
                      200,-1,1);
  h_HEenergydifference = dbe_->book1D("HEenergyDifference",
                      "Sum(E_HEPlus - E_HEMinus)/Sum(E_HEPlus + E_HEMinus)",
                      200,-1,1);

  h_LumiPlot_LS_MinBiasEvents=dbe_->book1D("MinBiasEventsPerLS",
                       "Number of MinBias Events vs LS (HT cut and HFM-HFP time cut)",
                       NLumiBlocks_/10,0.5,NLumiBlocks_+0.5); 
  h_LumiPlot_LS_MinBiasEvents_notimecut=dbe_->book1D("MinBiasEventsPerLS_notimecut",
                             "Number of Events with MinBias vs LS (HFM,HFP HT>1,no time cut)",
                             NLumiBlocks_/10,0.5,NLumiBlocks_+0.5); 

  h_LumiPlot_SumHT_HFPlus_vs_HFMinus = dbe_->book2D("SumHT_plus_minus",
                            "HF+ Sum HT vs HF- Sum HT",60,0,30,60,0,30);
  h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus = dbe_->book2D("SumEnergy_plus_minus",
                            "HF+ Sum Energy  vs HF- Sum Energy",
                            60,0,150,60,0,150);
  h_LumiPlot_timeHFPlus_vs_timeHFMinus = dbe_->book2D("timeHFplus_vs_timeHFminus",
                              "Energy-weighted time average of HF+ vs HF-",
                              60,-60,60,60,-60,60);
  h_LumiPlot_BX_MinBiasEvents = dbe_->book1D("BX_MinBias_Events_TimeCut",
                      "BX # of MinBias events (HFM & HFP HT>1 & HFM-HFP time cut)",
                      3600,0,3600);
  h_LumiPlot_BX_MinBiasEvents_notimecut = dbe_->book1D("BX_MinBias_Events_notimecut",
                               "BX # of MinBias events (HFM,HFP HT>1, no time cut)",
                               3600,0,3600);
  // threshold plots must pass MinBias Trigger
  SetupEtaPhiHists(OccupancyThreshByDepth,"Above Threshold RecHit Occupancy","");
  h_rechitieta_thresh = dbe_->book1D("HcalRecHitIeta_thresh",
                  "Hcal RecHit ieta above energy and ET threshold",
                  83,-41.5,41.5);
  h_rechitiphi_thresh = dbe_->book1D("HcalRecHitIphi_thresh",
                  "Hcal RecHit iphi above energy and ET threshold",
                  72,0.5,72.5);

  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/sumplots");
  SetupEtaPhiHists(SumEnergyThreshByDepth,"Above Threshold RecHit Summed Energy","GeV");
  SetupEtaPhiHists(SumTimeThreshByDepth,"Above Threshold RecHit Summed Time","nS");
  SetupEtaPhiHists(SqrtSumEnergy2ThreshByDepth,"Above Threshold RecHit Sqrt Summed Energy2","GeV");

  dbe_->setCurrentFolder(subdir_+"Distributions_PassedMinBias/rechit_1D_plots");
  h_HBThreshTime=dbe_->book1D("HB_time_thresh", 
                  "HB RecHit Time Above Threshold",
                  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HBThreshOccupancy=dbe_->book1D("HB_occupancy_thresh",
                   "HB RecHit Occupancy Above Threshold",260,-0.5,2599.5);
  h_HEThreshTime=dbe_->book1D("HE_time_thresh", 
                  "HE RecHit Time Above Threshold",
                  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HEThreshOccupancy=dbe_->book1D("HE_occupancy_thresh",
                   "HE RecHit Occupancy Above Threshold",260,-0.5,2599.5);
  h_HOThreshTime=dbe_->book1D("HO_time_thresh", 
                  "HO RecHit Time Above Threshold",
                  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HOThreshOccupancy=dbe_->book1D("HO_occupancy_thresh",
                   "HO RecHit Occupancy Above Threshold",217,-0.5,2169.5);
  h_HFThreshTime=dbe_->book1D("HF_time_thresh", 
                  "HF RecHit Time Above Threshold",
                  int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HFThreshOccupancy=dbe_->book1D("HF_occupancy_thresh",
                   "HF RecHit Occupancy Above Threshold",
                   173,-0.5,1729.5);

  // Histograms for events that did passed Hcal-specfied HLT triggers
  dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers");
  
  h_LumiPlot_BX_HcalHLTEvents = dbe_->book1D("BX_HcalHLT_Events_TimeCut",
                      "BX # of HcalHLT events (HFM & HFP HT>1 & HFM-HFP time cut)",
                      3600,0,3600);
  h_LumiPlot_BX_HcalHLTEvents_notimecut = dbe_->book1D("BX_HcalHLT_Events_notimecut",
                               "BX # of HcalHLT events (HFM,HFP HT>1, no time cut)",
                               3600,0,3600);
  h_LumiPlot_LS_HcalHLTEvents=dbe_->book1D("HcalHLTEventsPerLS",
                       "Number of HcalHLT Events vs LS (HT cut and HFM-HFP time cut)",
                       NLumiBlocks_/10,0.5,NLumiBlocks_+0.5); 
  h_LumiPlot_LS_HcalHLTEvents_notimecut=dbe_->book1D("HcalHLTEventsPerLS_notimecut",
                             "Number of Events with HcalHLT vs LS (HFM,HFP HT>1,no time cut)",
                             NLumiBlocks_/10,0.5,NLumiBlocks_+0.5); 
  

  dbe_->setCurrentFolder(subdir_+"Distributions_PassedHcalHLTriggers/");
  h_HF_HcalHLT_weightedtimedifference = dbe_->book1D("HF_HcalHLT_weightedtimeDifference",
                         "Energy-Weighted time difference between HF+ and HF- Hcal HLT",
                         251,-250.5,250.5);
  h_HE_HcalHLT_weightedtimedifference = dbe_->book1D("HE_HcalHLT_weightedtimeDifference",
                         "Energy-Weighted time difference between HE+ and HE- Hcal HLT",
                         251,-250.5,250.5);
  h_HF_HcalHLT_energydifference = dbe_->book1D("HF_HcalHLT_energyDifference",
                         "Sum(E_HFPlus - E_HFMinus)/Sum(E_HFPlus + E_HFMinus)",
                         200,-1,1);
  h_HE_HcalHLT_energydifference = dbe_->book1D("HE_HcalHLT_energyDifference",
                         "Sum(E_HEPlus - E_HEMinus)/Sum(E_HEPlus + E_HEMinus)",
                         200,-1,1);

  // Do we want separate directories for Minbias, other flags at some point?
  dbe_->setCurrentFolder(subdir_+"AnomalousCellFlags");// HB Flag Histograms


  h_HFLongShort_vs_LS=dbe_->book1D("HFLongShort_vs_LS",
                   "HFLongShort Flags vs Lumi Section",
                   NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
  h_HFDigiTime_vs_LS=dbe_->book1D("HFDigiTime_vs_LS",
                  "HFDigiTime Flags vs Lumi Section",
                  NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
  h_HBHEHPDMult_vs_LS=dbe_->book1D("HBHEHPDMult_vs_LS",
                   "HBHEHPDMult Flags vs Lumi Section",
                   NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);
  h_HBHEPulseShape_vs_LS=dbe_->book1D("HBHEPulseShape_vs_LS",
                      "HBHEPulseShape Flags vs Lumi Section",
                      NLumiBlocks_/10,0.5,0.5+NLumiBlocks_);

  h_HF_FlagCorr=dbe_->book2D("HF_FlagCorrelation",
                 "HF LongShort vs. DigiTime flags; DigiTime; LongShort", 
                 2,-0.5,1.5,2,-0.5,1.5);
  h_HF_FlagCorr->setBinLabel(1,"OFF",1);
  h_HF_FlagCorr->setBinLabel(2,"ON",1);
  h_HF_FlagCorr->setBinLabel(1,"OFF",2);
  h_HF_FlagCorr->setBinLabel(2,"ON",2);

  h_HBHE_FlagCorr=dbe_->book2D("HBHE_FlagCorrelation",
                   "HBHE HpdHitMultiplicity vs. PulseShape flags; PulseShape; HpdHitMultiplicity", 
                   2,-0.5,1.5,2,-0.5,1.5);
  h_HBHE_FlagCorr->setBinLabel(1,"OFF",1);
  h_HBHE_FlagCorr->setBinLabel(2,"ON",1);
  h_HBHE_FlagCorr->setBinLabel(1,"OFF",2);
  h_HBHE_FlagCorr->setBinLabel(2,"ON",2);

  h_FlagMap_HPDMULT=dbe_->book2D("FlagMap_HPDMULT",
                 "RBX Map of HBHEHpdHitMultiplicity Flags;RBX;RM",
                  72,-0.5,71.5,4,0.5,4.5);
  h_FlagMap_PULSESHAPE=dbe_->book2D("FlagMap_PULSESHAPE",
                    "RBX Map of HBHEPulseShape Flags;RBX;RM",
                  72,-0.5,71.5,4,0.5,4.5);
  h_FlagMap_DIGITIME=dbe_->book2D("FlagMap_DIGITIME",
                  "RBX Map of HFDigiTime Flags;RBX;RM",
                  24,131.5,155.5,4,0.5,4.5);
  h_FlagMap_LONGSHORT=dbe_->book2D("FlagMap_LONGSHORT",
                   "RBX Map of HFLongShort Flags;RBX;RM",
                   24,131.5,155.5,4,0.5,4.5);

  h_FlagMap_TIMEADD=dbe_->book2D("FlagMap_TIMEADD",
                 "RBX Map of Timing Added Flags;RBX;RM",
                   156,-0.5,155.5,4,0.5,4.5);
  h_FlagMap_TIMESUBTRACT=dbe_->book2D("FlagMap_TIMESUBTRACT",
                      "RBX Map of Timing Subtracted Flags;RBX;RM",
                   156,-0.5,155.5,4,0.5,4.5);
  h_FlagMap_TIMEERROR=dbe_->book2D("FlagMap_TIMEERROR",
                   "RBX Map of Timing Error Flags;RBX;RM",
                   156,-0.5,155.5,4,0.5,4.5);

  h_HBflagcounter=dbe_->book1D("HBflags","HB flags",32,-0.5,31.5);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEHpdHitMultiplicity, "HpdHitMult",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEPulseShape, "PulseShape",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_R1R2, "HSCP R1R2",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_FracLeader, "HSCP FracLeader",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_OuterEnergy, "HSCP OuterEnergy",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_ExpFit, "HSCP ExpFit",1);
  // 2-bit timing counter
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust1",1);
  h_HBflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust2",1);
  //3-bit timing shape cut
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape1",1);
  h_HBflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape2",1);
  h_HBflagcounter->setBinLabel(3+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape3",1);

  // common flags
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
  h_HBflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);

  // HE Flag Histograms
  h_HEflagcounter=dbe_->book1D("HEflags","HE flags",32,-0.5,31.5);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEHpdHitMultiplicity, "HpdHitMult",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHEPulseShape, "PulseShape",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_R1R2, "HSCP R1R2",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_FracLeader, "HSCP FracLeader",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_OuterEnergy, "HSCP OuterEnergy",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HSCP_ExpFit, "HSCP ExpFit",1);
  // 2-bit timing counter
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust1",1);
  h_HEflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingTrustBits,"TimingTrust2",1);
  //3-bit timing shape cut
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape1",1);
  h_HEflagcounter->setBinLabel(2+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape2",1);
  h_HEflagcounter->setBinLabel(3+HcalCaloFlagLabels::HBHETimingShapedCutsBits,"TimingShape3",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
  h_HEflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);

  // HO Flag Histograms
  h_HOflagcounter=dbe_->book1D("HOflags","HO flags",32,-0.5,31.5);
  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
  h_HOflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);

  // HF Flag Histograms
  h_HFflagcounter=dbe_->book1D("HFflags","HF flags",32,-0.5,31.5);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::HFLongShort, "LongShort",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::HFDigiTime, "DigiTime",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::HFTimingTrustBits,"TimingTrust1",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingSubtractedBit, "Subtracted",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingAddedBit, "Added",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::TimingErrorBit, "TimingError",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::ADCSaturationBit, "Saturation",1);
  h_HFflagcounter->setBinLabel(1+HcalCaloFlagLabels::Fraction2TS,"Frac2TS",1);

  h_HBflagcounter->getTH1F()->LabelsOption("v");
  h_HEflagcounter->getTH1F()->LabelsOption("v");
  h_HOflagcounter->getTH1F()->LabelsOption("v");
  h_HFflagcounter->getTH1F()->LabelsOption("v");


  // Diagnostic plots are currently filled for all rechits (no trigger/threshold requirement)
  // hb
  dbe_->setCurrentFolder(subdir_+"diagnostics/hb");
  
  h_HBTimeVsEnergy=dbe_->book2D("HBTimeVsEnergy","HB Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);

  h_HBsizeVsLS=dbe_->bookProfile("HBRecHitsVsLB","HB RecHits vs Luminosity Block",
                 NLumiBlocks_,0.5,NLumiBlocks_+0.5,
                 100,0,10000);

  h_HBTime=dbe_->book1D("HB_time","HB RecHit Time",
            int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HBOccupancy=dbe_->book1D("HB_occupancy",
                 "HB RecHit Occupancy",260,-0.5,2599.5);
  
  //he
  dbe_->setCurrentFolder(subdir_+"diagnostics/he");

  h_HETimeVsEnergy=dbe_->book2D("HETimeVsEnergy","HE Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);

  h_HEsizeVsLS=dbe_->bookProfile("HERecHitsVsLB","HE RecHits vs Luminosity Block",
                 NLumiBlocks_,0.5,NLumiBlocks_+0.5,
                 100,0,10000);
    
  h_HETime=dbe_->book1D("HE_time","HE RecHit Time",
            int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HEOccupancy=dbe_->book1D("HE_occupancy","HE RecHit Occupancy",260,-0.5,2599.5);
   
  // ho
  dbe_->setCurrentFolder(subdir_+"diagnostics/ho"); 

  h_HOTimeVsEnergy=dbe_->book2D("HOTimeVsEnergy","HO Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);

  h_HOsizeVsLS=dbe_->bookProfile("HORecHitsVsLB","HO RecHits vs Luminosity Block",
                 NLumiBlocks_,0.5,NLumiBlocks_+0.5,
                 100,0,10000);
  h_HOTime=dbe_->book1D("HO_time",
            "HO RecHit Time",
            int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HOOccupancy=dbe_->book1D("HO_occupancy",
                 "HO RecHit Occupancy",217,-0.5,2169.5);
  
  // hf
  dbe_->setCurrentFolder(subdir_+"diagnostics/hf");

  h_HFTimeVsEnergy=dbe_->book2D("HFTimeVsEnergy","HF Time Vs Energy (All RecHits);Energy (GeV); time(nS)",100,0,500,40,-100,100);

  h_HFsizeVsLS=dbe_->bookProfile("HFRecHitsVsLB",
                 "HF RecHits vs Luminosity Block",
                 NLumiBlocks_,0.5,NLumiBlocks_+0.5,
                 100, 0,10000); 
  h_HFTime=dbe_->book1D("HF_time","HF RecHit Time",
            int(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN),RECHITMON_TIME_MIN,RECHITMON_TIME_MAX);
  h_HFOccupancy=dbe_->book1D("HF_occupancy","HF RecHit Occupancy",173,-0.5,1729.5);

  return;
} //void HcalRecHitMonitor::setup(...)

void HcalRecHitMonitor::zeroCounters

(

void

)

Definition at line 1403 of file HcalRecHitMonitor.cc.

References energy2_, energy2_thresh_, energy_, energy_thresh_, HB_occupancy_, HB_occupancy_thresh_, HBflagcounter_, HBtime_, HBtime_thresh_, HE_occupancy_, HE_occupancy_thresh_, HEflagcounter_, HEtime_, HEtime_thresh_, HF_occupancy_, HF_occupancy_thresh_, HFenergyLong_, HFenergyLong_thresh_, HFenergyShort_, HFenergyShort_thresh_, HFflagcounter_, HFlong_occupancy_, HFlong_occupancy_thresh_, HFshort_occupancy_, HFshort_occupancy_thresh_, HFtime_, HFtime_thresh_, HFtimeLong_, HFtimeLong_thresh_, HFtimeShort_, HFtimeShort_thresh_, HO_occupancy_, HO_occupancy_thresh_, HOflagcounter_, HOtime_, HOtime_thresh_, i, j, gen::k, occupancy_, occupancy_thresh_, RECHITMON_TIME_MAX, RECHITMON_TIME_MIN, time_, and time_thresh_.

Referenced by setup().

{
  // Set all histogram counters back to zero

  for (int i=0;i<32;++i)
    {
      HBflagcounter_[i]=0;
      HEflagcounter_[i]=0;
      HOflagcounter_[i]=0;
      HFflagcounter_[i]=0;

    }
  // TH2F counters
  for (int i=0;i<85;++i)
    {
      for (int j=0;j<72;++j)
    {
      for (int k=0;k<4;++k)
        {
          occupancy_[i][j][k]=0;
          occupancy_thresh_[i][j][k]=0;
          energy_[i][j][k]=0;
              energy2_[i][j][k]=0;
          energy_thresh_[i][j][k]=0;
          energy2_thresh_[i][j][k]=0;
          time_[i][j][k]=0;
          time_thresh_[i][j][k]=0;
        }
    } // for (int j=0;j<PHIBINS;++j)
    } // for (int i=0;i<87;++i)

  // TH1F counters
  
  for (int i=0;i<200;++i)
    {
      HFenergyLong_[i]=0;
      HFenergyLong_thresh_[i]=0;
      HFenergyShort_[i]=0;
      HFenergyShort_thresh_[i]=0;
    }

  // time
  for (int i=0;i<(RECHITMON_TIME_MAX-RECHITMON_TIME_MIN);++i)
    {
      HBtime_[i]=0;
      HBtime_thresh_[i]=0;
      HEtime_[i]=0;
      HEtime_thresh_[i]=0;
      HOtime_[i]=0;
      HOtime_thresh_[i]=0;
      HFtime_[i]=0;
      HFtime_thresh_[i]=0;
      HFtimeLong_[i]=0;
      HFtimeLong_thresh_[i]=0;
      HFtimeShort_[i]=0;
      HFtimeShort_thresh_[i]=0;
    }

  // occupancy
  for (int i=0;i<865;++i)
    {
      if (i<260)
    {
      HB_occupancy_[i]=0;
      HE_occupancy_[i]=0;
      HB_occupancy_thresh_[i]=0;
      HE_occupancy_thresh_[i]=0;
    }
      if (i<218)
    {
      HO_occupancy_[i]=0;
      HO_occupancy_thresh_[i]=0;
    }
      if (i<174)
    {
      HF_occupancy_[i]=0;
      HF_occupancy_thresh_[i]=0;
    }

      HFlong_occupancy_[i] =0;
      HFshort_occupancy_[i]=0;
      HFlong_occupancy_thresh_[i] =0;
      HFshort_occupancy_thresh_[i]=0;
    } // for (int i=0;i<865;++i)

  return;
} //void HcalRecHitMonitor::zeroCounters(void)

Member Data Documentation

double HcalRecHitMonitor::energy2_[85][72][4]

private

Definition at line 79 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::energy2_thresh_[85][72][4]

private

Definition at line 81 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::energy_[85][72][4]

private

Definition at line 78 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::energy_thresh_[85][72][4]

private

Definition at line 80 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::energyThreshold_

private

Definition at line 53 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor().

double HcalRecHitMonitor::ETThreshold_

private

Definition at line 59 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor().

MonitorElement* HcalRecHitMonitor::h_FlagMap_DIGITIME

private

Definition at line 168 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_FlagMap_HPDMULT

private

Definition at line 166 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_FlagMap_LONGSHORT

private

Definition at line 169 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_FlagMap_PULSESHAPE

private

Definition at line 167 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_FlagMap_TIMEADD

private

Definition at line 170 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_FlagMap_TIMEERROR

private

Definition at line 172 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_FlagMap_TIMESUBTRACT

private

Definition at line 171 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBflagcounter

private

Definition at line 161 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBHE_FlagCorr

private

Definition at line 180 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBHEHPDMult_vs_LS

private

Definition at line 176 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBHEPulseShape_vs_LS

private

Definition at line 177 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBM_weightedTime

private

Definition at line 221 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBOccupancy

private

Definition at line 143 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBP_weightedTime

private

Definition at line 220 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBsizeVsLS

private

Definition at line 136 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBThreshOccupancy

private

Definition at line 144 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBThreshTime

private

Definition at line 142 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBTime

private

Definition at line 141 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HBTimeVsEnergy

private

Definition at line 223 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HE_HcalHLT_energydifference

private

Definition at line 192 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HE_HcalHLT_weightedtimedifference

private

Definition at line 191 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEenergydifference

private

Definition at line 187 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEflagcounter

private

Definition at line 162 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEM_weightedTime

private

Definition at line 219 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEOccupancy

private

Definition at line 148 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEP_weightedTime

private

Definition at line 218 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEsizeVsLS

private

Definition at line 137 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEThreshOccupancy

private

Definition at line 149 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEThreshTime

private

Definition at line 147 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HETime

private

Definition at line 146 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HEtimedifference

private

Definition at line 186 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HETimeVsEnergy

private

Definition at line 224 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HF_FlagCorr

private

Definition at line 179 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HF_HcalHLT_energydifference

private

Definition at line 190 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HF_HcalHLT_weightedtimedifference

private

Definition at line 189 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFDigiTime_vs_LS

private

Definition at line 175 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFenergydifference

private

Definition at line 185 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFflagcounter

private

Definition at line 164 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFLongShort_vs_LS

private

Definition at line 174 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFM_weightedTime

private

Definition at line 217 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFOccupancy

private

Definition at line 158 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFP_weightedTime

private

Definition at line 216 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFsizeVsLS

private

Definition at line 139 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFThreshOccupancy

private

Definition at line 159 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFThreshTime

private

Definition at line 157 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFTime

private

Definition at line 156 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFtimedifference

private

Definition at line 184 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HFTimeVsEnergy

private

Definition at line 226 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOflagcounter

private

Definition at line 163 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOOccupancy

private

Definition at line 153 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOsizeVsLS

private

Definition at line 138 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOThreshOccupancy

private

Definition at line 154 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOThreshTime

private

Definition at line 152 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOTime

private

Definition at line 151 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_HOTimeVsEnergy

private

Definition at line 225 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_BX_allevents

private

Definition at line 205 of file HcalRecHitMonitor.h.

Referenced by analyze(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_BX_HcalHLTEvents

private

Definition at line 208 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_BX_HcalHLTEvents_notimecut

private

Definition at line 209 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_BX_MinBiasEvents

private

Definition at line 206 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_BX_MinBiasEvents_notimecut

private

Definition at line 207 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_LS_allevents

private

Definition at line 194 of file HcalRecHitMonitor.h.

Referenced by analyze(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_LS_HcalHLTEvents

private

Definition at line 197 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_LS_HcalHLTEvents_notimecut

private

Definition at line 198 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_LS_MinBiasEvents

private

Definition at line 195 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_LS_MinBiasEvents_notimecut

private

Definition at line 196 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_MinTime_vs_MinHT

private

Definition at line 211 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_SumEnergy_HFPlus_vs_HFMinus

private

Definition at line 203 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_SumHT_HFPlus_vs_HFMinus

private

Definition at line 200 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_timeHFPlus_vs_timeHFMinus

private

Definition at line 201 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_timeHT_HFM

private

Definition at line 212 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_LumiPlot_timeHT_HFP

private

Definition at line 213 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitieta

private

Definition at line 121 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitieta_05

private

Definition at line 124 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitieta_10

private

Definition at line 125 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitieta_100

private

Definition at line 127 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitieta_25

private

Definition at line 126 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitieta_thresh

private

Definition at line 133 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitiphi

private

Definition at line 122 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitiphi_05

private

Definition at line 128 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitiphi_10

private

Definition at line 129 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitiphi_100

private

Definition at line 131 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitiphi_25

private

Definition at line 130 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

MonitorElement* HcalRecHitMonitor::h_rechitiphi_thresh

private

Definition at line 134 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

MonitorElement* HcalRecHitMonitor::h_TriggeredEvents

private

Definition at line 215 of file HcalRecHitMonitor.h.

Referenced by processEvent(), and setup().

double HcalRecHitMonitor::HB_occupancy_[260]

private

Definition at line 87 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HB_occupancy_thresh_[260]

private

Definition at line 88 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HBenergyThreshold_

private

Definition at line 54 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

double HcalRecHitMonitor::HBETThreshold_

private

Definition at line 60 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

int HcalRecHitMonitor::HBflagcounter_[32]

private

Definition at line 115 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

edm::InputTag HcalRecHitMonitor::hbheRechitLabel_

private

Definition at line 232 of file HcalRecHitMonitor.h.

Referenced by analyze(), and HcalRecHitMonitor().

bool HcalRecHitMonitor::HBpresent_

private

Definition at line 229 of file HcalRecHitMonitor.h.

double HcalRecHitMonitor::HBtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 85 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HBtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 86 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

std::vector<std::string> HcalRecHitMonitor::HcalHLTBits_

private

Definition at line 235 of file HcalRecHitMonitor.h.

Referenced by beginRun(), HcalRecHitMonitor(), and processEvent().

double HcalRecHitMonitor::HE_occupancy_[260]

private

Definition at line 91 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HE_occupancy_thresh_[260]

private

Definition at line 92 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HEenergyThreshold_

private

Definition at line 55 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

double HcalRecHitMonitor::HEETThreshold_

private

Definition at line 61 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

int HcalRecHitMonitor::HEflagcounter_[32]

private

Definition at line 114 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

bool HcalRecHitMonitor::HEpresent_

private

Definition at line 229 of file HcalRecHitMonitor.h.

double HcalRecHitMonitor::HEtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 89 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HEtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 90 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HF_occupancy_[174]

private

Definition at line 107 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HF_occupancy_thresh_[174]

private

Definition at line 108 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HFenergyLong_[200]

private

Definition at line 99 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFenergyLong_thresh_[200]

private

Definition at line 100 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFenergyShort_[200]

private

Definition at line 103 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFenergyShort_thresh_[200]

private

Definition at line 104 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFenergyThreshold_

private

Definition at line 57 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

double HcalRecHitMonitor::HFETThreshold_

private

Definition at line 63 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

int HcalRecHitMonitor::HFflagcounter_[32]

private

Definition at line 117 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HFlong_occupancy_[865]

private

Definition at line 109 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFlong_occupancy_thresh_[865]

private

Definition at line 110 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

MonitorElement* HcalRecHitMonitor::HFP_HFM_Energy

private

Definition at line 227 of file HcalRecHitMonitor.h.

Referenced by processEvent_rechit(), and setup().

bool HcalRecHitMonitor::HFpresent_

private

Definition at line 229 of file HcalRecHitMonitor.h.

edm::InputTag HcalRecHitMonitor::hfRechitLabel_

private

Definition at line 232 of file HcalRecHitMonitor.h.

Referenced by analyze(), and HcalRecHitMonitor().

double HcalRecHitMonitor::HFshort_occupancy_[865]

private

Definition at line 111 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFshort_occupancy_thresh_[865]

private

Definition at line 112 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 97 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HFtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 98 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HFtimeLong_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 101 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFtimeLong_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 102 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFtimeShort_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 105 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

double HcalRecHitMonitor::HFtimeShort_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 106 of file HcalRecHitMonitor.h.

Referenced by zeroCounters().

edm::InputTag HcalRecHitMonitor::hltresultsLabel_

private

Definition at line 234 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), and processEvent().

double HcalRecHitMonitor::HO_occupancy_[218]

private

Definition at line 95 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HO_occupancy_thresh_[218]

private

Definition at line 96 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HOenergyThreshold_

private

Definition at line 56 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

double HcalRecHitMonitor::HOETThreshold_

private

Definition at line 62 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().

int HcalRecHitMonitor::HOflagcounter_[32]

private

Definition at line 116 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

bool HcalRecHitMonitor::HOpresent_

private

Definition at line 229 of file HcalRecHitMonitor.h.

edm::InputTag HcalRecHitMonitor::hoRechitLabel_

private

Definition at line 232 of file HcalRecHitMonitor.h.

Referenced by analyze(), and HcalRecHitMonitor().

double HcalRecHitMonitor::HOtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 93 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::HOtime_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN]

private

Definition at line 94 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

edm::InputTag HcalRecHitMonitor::l1gtLabel_

private

Definition at line 233 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor().

std::vector<std::string> HcalRecHitMonitor::MinBiasHLTBits_

private

Definition at line 236 of file HcalRecHitMonitor.h.

Referenced by beginRun(), HcalRecHitMonitor(), and processEvent().

unsigned int HcalRecHitMonitor::occupancy_[85][72][4]

private

Definition at line 76 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

unsigned int HcalRecHitMonitor::occupancy_thresh_[85][72][4]

private

Definition at line 77 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

EtaPhiHists HcalRecHitMonitor::OccupancyByDepth

private

Definition at line 66 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

EtaPhiHists HcalRecHitMonitor::OccupancyThreshByDepth

private

Definition at line 67 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

bool HcalRecHitMonitor::setupDone_

private

Definition at line 230 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), and setup().

EtaPhiHists HcalRecHitMonitor::SqrtSumEnergy2ByDepth

private

Definition at line 70 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

EtaPhiHists HcalRecHitMonitor::SqrtSumEnergy2ThreshByDepth

private

Definition at line 72 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

EtaPhiHists HcalRecHitMonitor::SumEnergyByDepth

private

Definition at line 69 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

EtaPhiHists HcalRecHitMonitor::SumEnergyThreshByDepth

private

Definition at line 71 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

EtaPhiHists HcalRecHitMonitor::SumTimeByDepth

private

Definition at line 73 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and setup().

EtaPhiHists HcalRecHitMonitor::SumTimeThreshByDepth

private

Definition at line 74 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), and setup().

double HcalRecHitMonitor::time_[85][72][4]

private

Definition at line 82 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::time_thresh_[85][72][4]

private

Definition at line 83 of file HcalRecHitMonitor.h.

Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().

double HcalRecHitMonitor::timediffThresh_

private

Definition at line 182 of file HcalRecHitMonitor.h.

Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().