|
|
|
#include <HcalRecHitMonitor.h>
| 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::NLumiBlocks_, HcalBaseDQMonitor::Online_, HcalBaseDQMonitor::prefixME_, HcalBaseDQMonitor::skipOutOfOrderLS_, 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.);
} //constructor
| HcalRecHitMonitor::~HcalRecHitMonitor | ( | ) |
Definition at line 61 of file HcalRecHitMonitor.cc.
{
} //destructor
| void HcalRecHitMonitor::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | s | ||
| ) | [virtual] |
Reimplemented from HcalBaseDQMonitor.
Definition at line 576 of file HcalRecHitMonitor.cc.
References edm::EventBase::bunchCrossing(), gather_cfg::cout, HcalBaseDQMonitor::currentLS, HcalBaseDQMonitor::debug_, alignCSCRings::e, MonitorElement::Fill(), edm::Event::getByLabel(), h_LumiPlot_BX_allevents, h_LumiPlot_LS_allevents, hbheRechitLabel_, hfRechitLabel_, hoRechitLabel_, HcalBaseDQMonitor::IsAllowedCalibType(), HcalBaseDQMonitor::LumiInOrder(), edm::EventBase::luminosityBlock(), and processEvent().
{
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 486 of file HcalRecHitMonitor.cc.
References DQMStore::bookString(), gather_cfg::cout, HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::debug_, HcalHLTBits_, i, HcalBaseDQMonitor::mergeRuns_, MinBiasHLTBits_, reset(), DQMStore::setCurrentFolder(), setup(), 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 549 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 539 of file HcalRecHitMonitor.cc.
References 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 1225 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 517 of file HcalRecHitMonitor.cc.
References gather_cfg::cout, HcalBaseDQMonitor::debug_, and HcalBaseDQMonitor::subdir_.
| void HcalRecHitMonitor::fill_Nevents | ( | void | ) | [private] |
Definition at line 1236 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 615 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 679 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(), 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 527 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 69 of file HcalRecHitMonitor.cc.
References HcalCaloFlagLabels::ADCSaturationBit, DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookFloat(), DQMStore::bookProfile(), gather_cfg::cout, HcalLogicalMapGenerator::createMap(), HcalBaseDQMonitor::dbe_, HcalBaseDQMonitor::debug_, MonitorElement::Fill(), HcalCaloFlagLabels::Fraction2TS, relval_steps::gen(), 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, logicalMap, HcalBaseDQMonitor::NLumiBlocks_, OccupancyByDepth, OccupancyThreshByDepth, RECHITMON_TIME_MAX, RECHITMON_TIME_MIN, reset(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), HcalBaseDQMonitor::SetupEtaPhiHists(), SqrtSumEnergy2ByDepth, SqrtSumEnergy2ThreshByDepth, HcalBaseDQMonitor::subdir_, SumEnergyByDepth, SumEnergyThreshByDepth, SumTimeByDepth, SumTimeThreshByDepth, timediffThresh_, HcalCaloFlagLabels::TimingAddedBit, HcalCaloFlagLabels::TimingErrorBit, HcalCaloFlagLabels::TimingSubtractedBit, and zeroCounters().
Referenced by beginRun().
{
HcalBaseDQMonitor::setup();
if (debug_>0)
std::cout <<"<HcalRecHitMonitor::setup> Setting up histograms"<<std::endl;
// Can we include this just in the setup, or do we need to get a new logical map with every run?
HcalLogicalMapGenerator gen;
logicalMap=new HcalLogicalMap(gen.createMap());
// 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);
// clear all counters, reset histograms
this->zeroCounters();
this->reset();
return;
} //void HcalRecHitMonitor::setup(...)
| void HcalRecHitMonitor::zeroCounters | ( | void | ) |
Definition at line 1395 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)
double HcalRecHitMonitor::energy2_[85][72][4] [private] |
Definition at line 82 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::energy2_thresh_[85][72][4] [private] |
Definition at line 84 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::energy_[85][72][4] [private] |
Definition at line 81 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::energy_thresh_[85][72][4] [private] |
Definition at line 83 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::energyThreshold_ [private] |
Definition at line 54 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor().
double HcalRecHitMonitor::ETThreshold_ [private] |
Definition at line 60 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor().
Definition at line 171 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 169 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 172 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 170 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 173 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 175 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 174 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HBflagcounter [private] |
Definition at line 164 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HBHE_FlagCorr [private] |
Definition at line 183 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 179 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 180 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 224 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HBOccupancy [private] |
Definition at line 146 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
Definition at line 223 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HBsizeVsLS [private] |
Definition at line 139 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 147 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_HBThreshTime [private] |
Definition at line 145 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HBTime [private] |
Definition at line 144 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 226 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 195 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 194 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 190 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HEflagcounter [private] |
Definition at line 165 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 222 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HEOccupancy [private] |
Definition at line 151 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
Definition at line 221 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HEsizeVsLS [private] |
Definition at line 140 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 152 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_HEThreshTime [private] |
Definition at line 150 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HETime [private] |
Definition at line 149 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 189 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 227 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HF_FlagCorr [private] |
Definition at line 182 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 193 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 192 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 178 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 188 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HFflagcounter [private] |
Definition at line 167 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 177 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 220 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HFOccupancy [private] |
Definition at line 161 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
Definition at line 219 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HFsizeVsLS [private] |
Definition at line 142 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 162 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_HFThreshTime [private] |
Definition at line 160 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HFTime [private] |
Definition at line 159 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 187 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 229 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HOflagcounter [private] |
Definition at line 166 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HOOccupancy [private] |
Definition at line 156 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_HOsizeVsLS [private] |
Definition at line 141 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 157 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_HOThreshTime [private] |
Definition at line 155 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_HOTime [private] |
Definition at line 154 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 228 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 208 of file HcalRecHitMonitor.h.
Definition at line 211 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 212 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 209 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 210 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 197 of file HcalRecHitMonitor.h.
Definition at line 200 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 201 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 198 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 199 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 214 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 206 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 203 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 204 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 215 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 216 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitieta [private] |
Definition at line 124 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitieta_05 [private] |
Definition at line 127 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitieta_10 [private] |
Definition at line 128 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 130 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitieta_25 [private] |
Definition at line 129 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 136 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitiphi [private] |
Definition at line 125 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitiphi_05 [private] |
Definition at line 131 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitiphi_10 [private] |
Definition at line 132 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 134 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
MonitorElement* HcalRecHitMonitor::h_rechitiphi_25 [private] |
Definition at line 133 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
Definition at line 137 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
Definition at line 218 of file HcalRecHitMonitor.h.
Referenced by processEvent(), and setup().
double HcalRecHitMonitor::HB_occupancy_[260] [private] |
Definition at line 90 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HB_occupancy_thresh_[260] [private] |
Definition at line 91 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HBenergyThreshold_ [private] |
Definition at line 55 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
double HcalRecHitMonitor::HBETThreshold_ [private] |
Definition at line 61 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
int HcalRecHitMonitor::HBflagcounter_[32] [private] |
Definition at line 118 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
Definition at line 234 of file HcalRecHitMonitor.h.
Referenced by analyze(), and HcalRecHitMonitor().
bool HcalRecHitMonitor::HBpresent_ [private] |
Reimplemented from HcalBaseDQMonitor.
Definition at line 232 of file HcalRecHitMonitor.h.
double HcalRecHitMonitor::HBtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 88 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 89 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
std::vector<std::string> HcalRecHitMonitor::HcalHLTBits_ [private] |
Definition at line 237 of file HcalRecHitMonitor.h.
Referenced by beginRun(), HcalRecHitMonitor(), and processEvent().
double HcalRecHitMonitor::HE_occupancy_[260] [private] |
Definition at line 94 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HE_occupancy_thresh_[260] [private] |
Definition at line 95 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HEenergyThreshold_ [private] |
Definition at line 56 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
double HcalRecHitMonitor::HEETThreshold_ [private] |
Definition at line 62 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
int HcalRecHitMonitor::HEflagcounter_[32] [private] |
Definition at line 117 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
bool HcalRecHitMonitor::HEpresent_ [private] |
Reimplemented from HcalBaseDQMonitor.
Definition at line 232 of file HcalRecHitMonitor.h.
double HcalRecHitMonitor::HEtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 92 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 93 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HF_occupancy_[174] [private] |
Definition at line 110 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HF_occupancy_thresh_[174] [private] |
Definition at line 111 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HFenergyLong_[200] [private] |
Definition at line 102 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFenergyLong_thresh_[200] [private] |
Definition at line 103 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFenergyShort_[200] [private] |
Definition at line 106 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFenergyShort_thresh_[200] [private] |
Definition at line 107 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFenergyThreshold_ [private] |
Definition at line 58 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
double HcalRecHitMonitor::HFETThreshold_ [private] |
Definition at line 64 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
int HcalRecHitMonitor::HFflagcounter_[32] [private] |
Definition at line 120 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HFlong_occupancy_[865] [private] |
Definition at line 112 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFlong_occupancy_thresh_[865] [private] |
Definition at line 113 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
MonitorElement* HcalRecHitMonitor::HFP_HFM_Energy [private] |
Definition at line 230 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
bool HcalRecHitMonitor::HFpresent_ [private] |
Reimplemented from HcalBaseDQMonitor.
Definition at line 232 of file HcalRecHitMonitor.h.
Definition at line 234 of file HcalRecHitMonitor.h.
Referenced by analyze(), and HcalRecHitMonitor().
double HcalRecHitMonitor::HFshort_occupancy_[865] [private] |
Definition at line 114 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFshort_occupancy_thresh_[865] [private] |
Definition at line 115 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 100 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 101 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 104 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFtimeLong_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 105 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFtimeShort_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 108 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
double HcalRecHitMonitor::HFtimeShort_thresh_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 109 of file HcalRecHitMonitor.h.
Referenced by zeroCounters().
Definition at line 236 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), and processEvent().
double HcalRecHitMonitor::HO_occupancy_[218] [private] |
Definition at line 98 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HO_occupancy_thresh_[218] [private] |
Definition at line 99 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::HOenergyThreshold_ [private] |
Definition at line 57 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
double HcalRecHitMonitor::HOETThreshold_ [private] |
Definition at line 63 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
int HcalRecHitMonitor::HOflagcounter_[32] [private] |
Definition at line 119 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
bool HcalRecHitMonitor::HOpresent_ [private] |
Reimplemented from HcalBaseDQMonitor.
Definition at line 232 of file HcalRecHitMonitor.h.
Definition at line 234 of file HcalRecHitMonitor.h.
Referenced by analyze(), and HcalRecHitMonitor().
double HcalRecHitMonitor::HOtime_[RECHITMON_TIME_MAX-RECHITMON_TIME_MIN] [private] |
Definition at line 96 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 97 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
edm::InputTag HcalRecHitMonitor::l1gtLabel_ [private] |
Definition at line 235 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor().
HcalLogicalMap* HcalRecHitMonitor::logicalMap [private] |
Definition at line 66 of file HcalRecHitMonitor.h.
Referenced by processEvent_rechit(), and setup().
std::vector<std::string> HcalRecHitMonitor::MinBiasHLTBits_ [private] |
Definition at line 238 of file HcalRecHitMonitor.h.
Referenced by beginRun(), HcalRecHitMonitor(), and processEvent().
unsigned int HcalRecHitMonitor::occupancy_[85][72][4] [private] |
Definition at line 79 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
unsigned int HcalRecHitMonitor::occupancy_thresh_[85][72][4] [private] |
Definition at line 80 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
Definition at line 69 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 70 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 73 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 75 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
Definition at line 72 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 74 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
EtaPhiHists HcalRecHitMonitor::SumTimeByDepth [private] |
Definition at line 76 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and setup().
Definition at line 77 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), and setup().
double HcalRecHitMonitor::time_[85][72][4] [private] |
Definition at line 85 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::time_thresh_[85][72][4] [private] |
Definition at line 86 of file HcalRecHitMonitor.h.
Referenced by fill_Nevents(), processEvent_rechit(), and zeroCounters().
double HcalRecHitMonitor::timediffThresh_ [private] |
Definition at line 185 of file HcalRecHitMonitor.h.
Referenced by HcalRecHitMonitor(), processEvent_rechit(), and setup().
1.7.3