|
|
|
#include <HcalNoiseInfoProducer.h>
Definition at line 38 of file HcalNoiseInfoProducer.h.
| HcalNoiseInfoProducer::HcalNoiseInfoProducer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 33 of file HcalNoiseInfoProducer.cc.
References adc2fC, calibdigiHBHEthreshold_, calibdigiHBHEtimeslices_, calibdigiHFthreshold_, calibdigiHFtimeslices_, caloTowerCollName_, digiCollName_, edm::ParameterSet::exists(), edm::ParameterSet::existsAs(), fillCaloTowers_, fillDigis_, fillRecHits_, fillTracks_, edm::ParameterSet::getParameter(), HcalAcceptSeverityLevel_, HcalRecHitFlagsToBeExcluded_, i, maxCaloTowerIEta_, maxProblemRBXs_, maxTrackEta_, minHighHitE_, minLowHitE_, minRecHitE_, minTrackPt_, recHitCollName_, python::multivaluedict::sort(), AlCaHLTBitMon_QueryRunRegistry::string, trackCollName_, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, TS4TS5UpperCut_, and useCalibDigi_.
: algo_(iConfig) { // set the parameters fillDigis_ = iConfig.getParameter<bool>("fillDigis"); fillRecHits_ = iConfig.getParameter<bool>("fillRecHits"); fillCaloTowers_ = iConfig.getParameter<bool>("fillCaloTowers"); fillTracks_ = iConfig.getParameter<bool>("fillTracks"); maxProblemRBXs_ = iConfig.getParameter<int>("maxProblemRBXs"); maxCaloTowerIEta_ = iConfig.getParameter<int>("maxCaloTowerIEta"); maxTrackEta_ = iConfig.getParameter<double>("maxTrackEta"); minTrackPt_ = iConfig.getParameter<double>("minTrackPt"); digiCollName_ = iConfig.getParameter<std::string>("digiCollName"); recHitCollName_ = iConfig.getParameter<std::string>("recHitCollName"); caloTowerCollName_ = iConfig.getParameter<std::string>("caloTowerCollName"); trackCollName_ = iConfig.getParameter<std::string>("trackCollName"); minRecHitE_ = iConfig.getParameter<double>("minRecHitE"); minLowHitE_ = iConfig.getParameter<double>("minLowHitE"); minHighHitE_ = iConfig.getParameter<double>("minHighHitE"); HcalAcceptSeverityLevel_ = iConfig.getParameter<uint32_t>("HcalAcceptSeverityLevel"); if (iConfig.exists("HcalRecHitFlagsToBeExcluded")) HcalRecHitFlagsToBeExcluded_ = iConfig.getParameter<std::vector<int> >("HcalRecHitFlagsToBeExcluded"); else{ edm::LogWarning("MisConfiguration")<<"the module is missing the parameter HcalAcceptSeverityLevel. created empty."; HcalRecHitFlagsToBeExcluded_.resize(0); } // Digi threshold and time slices to use for HBHE and HF calibration digis useCalibDigi_ = true; if(iConfig.existsAs<double>("calibdigiHBHEthreshold") == false) useCalibDigi_ = false; if(iConfig.existsAs<double>("calibdigiHFthreshold") == false) useCalibDigi_ = false; if(iConfig.existsAs<std::vector<int> >("calibdigiHBHEtimeslices") == false) useCalibDigi_ = false; if(iConfig.existsAs<std::vector<int> >("calibdigiHFtimeslices") == false) useCalibDigi_ = false; if(useCalibDigi_ == true) { calibdigiHBHEthreshold_ = iConfig.getParameter<double>("calibdigiHBHEthreshold"); calibdigiHBHEtimeslices_ = iConfig.getParameter<std::vector<int> >("calibdigiHBHEtimeslices"); calibdigiHFthreshold_ = iConfig.getParameter<double>("calibdigiHFthreshold"); calibdigiHFtimeslices_ = iConfig.getParameter<std::vector<int> >("calibdigiHFtimeslices"); } else { calibdigiHBHEthreshold_ = 0; calibdigiHBHEtimeslices_ = std::vector<int>(); calibdigiHFthreshold_ = 0; calibdigiHFtimeslices_ = std::vector<int>(); } TS4TS5EnergyThreshold_ = iConfig.getParameter<double>("TS4TS5EnergyThreshold"); std::vector<double> TS4TS5UpperThresholdTemp = iConfig.getParameter<std::vector<double> >("TS4TS5UpperThreshold"); std::vector<double> TS4TS5UpperCutTemp = iConfig.getParameter<std::vector<double> >("TS4TS5UpperCut"); std::vector<double> TS4TS5LowerThresholdTemp = iConfig.getParameter<std::vector<double> >("TS4TS5LowerThreshold"); std::vector<double> TS4TS5LowerCutTemp = iConfig.getParameter<std::vector<double> >("TS4TS5LowerCut"); for(int i = 0; i < (int)TS4TS5UpperThresholdTemp.size() && i < (int)TS4TS5UpperCutTemp.size(); i++) TS4TS5UpperCut_.push_back(std::pair<double, double>(TS4TS5UpperThresholdTemp[i], TS4TS5UpperCutTemp[i])); sort(TS4TS5UpperCut_.begin(), TS4TS5UpperCut_.end()); for(int i = 0; i < (int)TS4TS5LowerThresholdTemp.size() && i < (int)TS4TS5LowerCutTemp.size(); i++) TS4TS5LowerCut_.push_back(std::pair<double, double>(TS4TS5LowerThresholdTemp[i], TS4TS5LowerCutTemp[i])); sort(TS4TS5LowerCut_.begin(), TS4TS5LowerCut_.end()); // if digis are filled, then rechits must also be filled if(fillDigis_ && !fillRecHits_) { fillRecHits_=true; edm::LogWarning("HCalNoiseInfoProducer") << " forcing fillRecHits to be true if fillDigis is true.\n"; } const float adc2fCTemp[128]={-0.5,0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5,12.5, 13.5,15.,17.,19.,21.,23.,25.,27.,29.5,32.5,35.5,38.5,42.,46.,50.,54.5,59.5, 64.5,59.5,64.5,69.5,74.5,79.5,84.5,89.5,94.5,99.5,104.5,109.5,114.5,119.5, 124.5,129.5,137.,147.,157.,167.,177.,187.,197.,209.5,224.5,239.5,254.5,272., 292.,312.,334.5,359.5,384.5,359.5,384.5,409.5,434.5,459.5,484.5,509.5,534.5, 559.5,584.5,609.5,634.5,659.5,684.5,709.5,747.,797.,847.,897.,947.,997., 1047.,1109.5,1184.5,1259.5,1334.5,1422.,1522.,1622.,1734.5,1859.5,1984.5, 1859.5,1984.5,2109.5,2234.5,2359.5,2484.5,2609.5,2734.5,2859.5,2984.5, 3109.5,3234.5,3359.5,3484.5,3609.5,3797.,4047.,4297.,4547.,4797.,5047., 5297.,5609.5,5984.5,6359.5,6734.5,7172.,7672.,8172.,8734.5,9359.5,9984.5}; for(int i = 0; i < 128; i++) adc2fC[i] = adc2fCTemp[i]; // we produce a vector of HcalNoiseRBXs produces<HcalNoiseRBXCollection>(); // we also produce a noise summary produces<HcalNoiseSummary>(); }
| HcalNoiseInfoProducer::~HcalNoiseInfoProducer | ( | ) |
Definition at line 127 of file HcalNoiseInfoProducer.cc.
{
}
| void HcalNoiseInfoProducer::fillcalotwrs | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HcalNoiseRBXArray & | array, | ||
| HcalNoiseSummary & | summary | ||
| ) | const [private] |
Definition at line 643 of file HcalNoiseInfoProducer.cc.
References caloTowerCollName_, CaloTower::emEnergy(), HcalNoiseSummary::emenergy_, Exception, reco::HcalNoiseRBXArray::findHPD(), edm::Event::getByLabel(), CaloTower::hadEnergy(), HcalNoiseSummary::hadenergy_, patZpeak::handle, CaloTower::ietaAbs(), edm::HandleBase::isValid(), maxCaloTowerIEta_, and edm::errors::ProductNotFound.
Referenced by produce().
{
// get the calotowers
edm::Handle<CaloTowerCollection> handle;
iEvent.getByLabel(caloTowerCollName_, handle);
if(!handle.isValid()) {
throw edm::Exception(edm::errors::ProductNotFound)
<< " could not find CaloTowerCollection named " << caloTowerCollName_ << "\n.";
return;
}
summary.emenergy_ = summary.hadenergy_ = 0.0;
// loop over all of the calotower information
for(CaloTowerCollection::const_iterator it = handle->begin(); it!=handle->end(); ++it) {
const CaloTower& twr=(*it);
// create a persistent reference to the tower
edm::Ref<CaloTowerCollection> myRef(handle, it-handle->begin());
// get all of the hpd's that are pointed to by the calotower
std::vector<std::vector<HcalNoiseHPD>::iterator> hpditervec;
array.findHPD(twr, hpditervec);
// loop over the hpd's and add the reference to the RefVectors
for(std::vector<std::vector<HcalNoiseHPD>::iterator>::iterator it=hpditervec.begin();
it!=hpditervec.end(); ++it)
(*it)->calotowers_.push_back(myRef);
// skip over anything with |ieta|>maxCaloTowerIEta
if(twr.ietaAbs()>maxCaloTowerIEta_) {
summary.emenergy_ += twr.emEnergy();
summary.hadenergy_ += twr.hadEnergy();
}
}
return;
}
| void HcalNoiseInfoProducer::filldigis | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HcalNoiseRBXArray & | array, | ||
| HcalNoiseSummary & | summary | ||
| ) | [private] |
Definition at line 300 of file HcalNoiseInfoProducer.cc.
References ecalMGPA::adc(), HcalCoderDb::adc2fC(), adc2fC, reco::HcalNoiseRBX::allCharge_, edm::RefVector< C, T, F >::begin(), reco::HcalNoiseHPD::big5Charge_, reco::HcalNoiseHPD::bigCharge_, HcalNoiseSummary::calibCharge_, HcalNoiseSummary::calibChargegt15TS45_, HcalNoiseSummary::calibChargeTS45_, HcalNoiseSummary::calibCountgt15TS45_, HcalNoiseSummary::calibCountHF_, HcalNoiseSummary::calibCountTS45_, calibdigiHBHEthreshold_, calibdigiHBHEtimeslices_, calibdigiHFthreshold_, calibdigiHFtimeslices_, HcalCalibDetId::calibFlavor(), digiCollName_, HcalSeverityLevelComputer::dropChannel(), edm::RefVector< C, T, F >::end(), Exception, reco::HcalNoiseRBXArray::findHPD(), reco::HcalNoiseRBXArray::findRBX(), edm::EventSetup::get(), edm::Event::getByLabel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), patZpeak::handle, HcalBarrel, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary, HcalEndcap, HcalForward, HcalCalibDetId::HOCrosstalk, i, HBHEDataFrame::id(), edm::HandleBase::isValid(), reco::HcalNoiseHPD::maxZeros_, HcalCalibrations::pedestal(), edm::ESHandle< T >::product(), edm::errors::ProductNotFound, DetId::rawId(), HI_PhotonSkim_cff::rechits, reco::HcalNoiseHPD::rechits_, CaloSamples::size(), TotalCalibCharge, reco::HcalNoiseHPD::totalZeros_, and HBHEDataFrame::zsMarkAndPass().
Referenced by produce().
{
// Some initialization
TotalCalibCharge = 0;
// Starting with this version (updated by Jeff Temple on Dec. 6, 2012), the "TS45" names in the variables are mis-nomers. The actual time slices used are determined from the digiTimeSlices_ variable, which may not be limited to only time slices 4 and 5. For now, "TS45" name kept, because that is what is used in HcalNoiseSummary object (in GetCalibCountTS45, etc.). Likewise, the charge value in 'gt15' is now configurable, though the name remains the same. For HBHE, we track both the number of calibration channels (NcalibTS45) and the number of calibration channels above threshold (NcalibTS45gt15). For HF, we track only the number of channels above the given threshold in the given time window (NcalibHFgtX). Default for HF in 2012 is to use the full time sample with effectively no threshold (threshold=-999)
int NcalibTS45=0;
int NcalibTS45gt15=0;
int NcalibHFgtX=0;
double chargecalibTS45=0;
double chargecalibgt15TS45=0;
// get the conditions and channel quality
edm::ESHandle<HcalDbService> conditions;
iSetup.get<HcalDbRecord>().get(conditions);
edm::ESHandle<HcalChannelQuality> qualhandle;
iSetup.get<HcalChannelQualityRcd>().get(qualhandle);
const HcalChannelQuality* myqual = qualhandle.product();
edm::ESHandle<HcalSeverityLevelComputer> mycomputer;
iSetup.get<HcalSeverityLevelComputerRcd>().get(mycomputer);
const HcalSeverityLevelComputer* mySeverity = mycomputer.product();
// get the digis
edm::Handle<HBHEDigiCollection> handle;
iEvent.getByLabel(digiCollName_, handle);
if(!handle.isValid()) {
throw edm::Exception(edm::errors::ProductNotFound) << " could not find HBHEDigiCollection named " << digiCollName_ << "\n.";
return;
}
// loop over all of the digi information
for(HBHEDigiCollection::const_iterator it=handle->begin(); it!=handle->end(); ++it) {
const HBHEDataFrame &digi=(*it);
HcalDetId cell = digi.id();
DetId detcell=(DetId)cell;
// check on cells to be ignored and dropped
const HcalChannelStatus* mydigistatus=myqual->getValues(detcell.rawId());
if(mySeverity->dropChannel(mydigistatus->getValue())) continue;
if(digi.zsMarkAndPass()) continue;
// Drop if exclude bit set
if ((mydigistatus->getValue() & (1 <<HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary))==1) continue;
// get the calibrations and coder
const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
HcalCoderDb coder (*channelCoder, *shape);
// match the digi to an rbx and hpd
HcalNoiseRBX &rbx=(*array.findRBX(digi));
HcalNoiseHPD &hpd=(*array.findHPD(digi));
// determine if the digi is one the highest energy hits in the HPD
// this works because the rechits are sorted by energy (see fillrechits() below)
bool isBig=false, isBig5=false, isRBX=false;
int counter=0;
edm::RefVector<HBHERecHitCollection> &rechits=hpd.rechits_;
for(edm::RefVector<HBHERecHitCollection>::const_iterator rit=rechits.begin();
rit!=rechits.end(); ++rit, ++counter) {
if((*rit)->id() == digi.id()) {
if(counter==0) isBig=isBig5=true; // digi is also the highest energy rechit
if(counter<5) isBig5=true; // digi is one of 5 highest energy rechits
isRBX=true;
}
}
// loop over each of the digi's time slices
int totalzeros=0;
CaloSamples tool;
coder.adc2fC(digi,tool);
for(int ts=0; ts<tool.size(); ++ts) {
// count zero's
if(digi[ts].adc()==0) {
++hpd.totalZeros_;
++totalzeros;
}
// get the fC's
double corrfc = tool[ts]-calibrations.pedestal(digi[ts].capid());
// fill the relevant digi arrays
if(isBig) hpd.bigCharge_[ts]+=corrfc;
if(isBig5) hpd.big5Charge_[ts]+=corrfc;
if(isRBX) rbx.allCharge_[ts]+=corrfc;
}
// record the maximum number of zero's found
if(totalzeros>hpd.maxZeros_)
hpd.maxZeros_=totalzeros;
}
// get the calibration digis
edm::Handle<HcalCalibDigiCollection> hCalib;
iEvent.getByLabel("hcalDigis", hCalib);
// get total charge in calibration channels
if(hCalib.isValid() == true)
{
for(HcalCalibDigiCollection::const_iterator digi = hCalib->begin(); digi != hCalib->end(); digi++)
{
if(digi->id().hcalSubdet() == 0)
continue;
/*
HcalCalibDetId cell = digi->id();
DetId detcell = (DetId)cell;
const HcalChannelStatus* mydigistatus=myqual->getValues(detcell.rawId());
if(mySeverity->dropChannel(mydigistatus->getValue()))
continue;
if(digi->zsMarkAndPass())
continue;
const HcalQIECoder *channelCoder = conditions->getHcalCoder(cell);
const HcalQIEShape* shape = conditions->getHcalShape(channelCoder);
HcalCoderDb coder(*channelCoder, *shape);
CaloSamples tool;
coder.adc2fC(*digi, tool);
for(int i = 0; i < (int)digi->size(); i++)
TotalCalibCharge = TotalCalibCharge + tool[i];
*/
// Original code computes total calib charge over all digis. While I think it would be more useful to skip
// zs mark-and-pass channels, I keep this computation as is. Individual HBHE and HF variables do skip
// the m-p channels. -- Jeff Temple, 6 December 2012
for(int i = 0; i < (int)digi->size(); i++)
TotalCalibCharge = TotalCalibCharge + adc2fC[digi->sample(i).adc()&0xff];
HcalCalibDetId myid=(HcalCalibDetId)digi->id();
if ( myid.calibFlavor()==HcalCalibDetId::HOCrosstalk)
continue; // ignore HOCrosstalk channels
if(digi->zsMarkAndPass()) continue; // skip "mark-and-pass" channels when computing charge in calib channels
if (digi->id().hcalSubdet()==HcalForward) // check HF
{
double sumChargeHF=0;
for (unsigned int i=0;i<calibdigiHFtimeslices_.size();++i)
{
// skip unphysical time slices
if (calibdigiHFtimeslices_[i]<0 || calibdigiHFtimeslices_[i]>digi->size())
continue;
sumChargeHF+=adc2fC[digi->sample(calibdigiHFtimeslices_[i]).adc()&0xff];
}
if (sumChargeHF>calibdigiHFthreshold_) ++NcalibHFgtX;
} // end of HF check
else if (digi->id().hcalSubdet()==HcalBarrel || digi->id().hcalSubdet()==HcalEndcap) // now check HBHE
{
double sumChargeHBHE=0;
for (unsigned int i=0;i<calibdigiHBHEtimeslices_.size();++i)
{
// skip unphysical time slices
if (calibdigiHBHEtimeslices_[i]<0 || calibdigiHBHEtimeslices_[i]>digi->size())
continue;
sumChargeHBHE+=adc2fC[digi->sample(calibdigiHBHEtimeslices_[i]).adc()&0xff];
}
++NcalibTS45;
chargecalibTS45+=sumChargeHBHE;
if (sumChargeHBHE>calibdigiHBHEthreshold_)
{
++NcalibTS45gt15;
chargecalibgt15TS45+=sumChargeHBHE;
}
} // end of HBHE check
} // loop on HcalCalibDigiCollection
} // if (hCalib.isValid()==true)
summary.calibCharge_ = TotalCalibCharge;
summary.calibCountTS45_=NcalibTS45;
summary.calibCountgt15TS45_=NcalibTS45gt15;
summary.calibChargeTS45_=chargecalibTS45;
summary.calibChargegt15TS45_=chargecalibgt15TS45;
summary.calibCountHF_=NcalibHFgtX;
return;
}
| void HcalNoiseInfoProducer::fillOtherSummaryVariables | ( | HcalNoiseSummary & | summary, |
| const CommonHcalNoiseRBXData & | data | ||
| ) | const [private] |
Definition at line 209 of file HcalNoiseInfoProducer.cc.
References algo_, HcalNoiseSummary::cnthit10_, HcalNoiseSummary::cnthit25_, CommonHcalNoiseRBXData::e10ts(), CommonHcalNoiseRBXData::e2ts(), HcalNoiseSummary::filterstatus_, HcalNoiseSummary::hasBadRBXTS4TS5_, CommonHcalNoiseRBXData::highEHitTimeSqrd(), HcalNoiseSummary::hlnoisetwrs_, CommonHcalNoiseRBXData::HPDEMF(), join(), HcalNoiseSummary::loosenoisetwrs_, CommonHcalNoiseRBXData::lowEHitTimeSqrd(), HcalNoiseSummary::max10_, HcalNoiseSummary::max10GeVHitTime(), HcalNoiseSummary::max25_, HcalNoiseSummary::max25GeVHitTime(), HcalNoiseSummary::maxe10ts_, HcalNoiseSummary::maxE2Over10TS(), HcalNoiseSummary::maxe2ts_, CommonHcalNoiseRBXData::maxHighEHitTime(), HcalNoiseSummary::maxHPDHits(), HcalNoiseSummary::maxhpdhits_, HcalNoiseSummary::maxhpdhitsnoother_, HcalNoiseSummary::maxHPDNoOtherHits(), CommonHcalNoiseRBXData::maxLowEHitTime(), HcalNoiseSummary::maxRBXHits(), HcalNoiseSummary::maxrbxhits_, HcalNoiseSummary::maxZeros(), HcalNoiseSummary::maxzeros_, HcalNoiseSummary::min10_, HcalNoiseSummary::min10GeVHitTime(), HcalNoiseSummary::min25_, HcalNoiseSummary::min25GeVHitTime(), HcalNoiseSummary::mine10ts_, HcalNoiseSummary::minE2Over10TS(), HcalNoiseSummary::mine2ts_, CommonHcalNoiseRBXData::minHighEHitTime(), HcalNoiseSummary::minHPDEMF(), HcalNoiseSummary::minhpdemf_, CommonHcalNoiseRBXData::minLowEHitTime(), HcalNoiseSummary::minRBXEMF(), HcalNoiseSummary::minrbxemf_, CommonHcalNoiseRBXData::numHighEHits(), CommonHcalNoiseRBXData::numHPDHits(), CommonHcalNoiseRBXData::numHPDNoOtherHits(), CommonHcalNoiseRBXData::numLowEHits(), CommonHcalNoiseRBXData::numRBXHits(), CommonHcalNoiseRBXData::numZeros(), HcalNoiseAlgo::passEMFThreshold(), HcalNoiseAlgo::passHighLevelNoiseFilter(), HcalNoiseAlgo::passLooseHits(), HcalNoiseAlgo::passLooseNoiseFilter(), HcalNoiseAlgo::passLooseRatio(), HcalNoiseAlgo::passLooseTiming(), HcalNoiseAlgo::passLooseZeros(), HcalNoiseAlgo::passRatioThreshold(), HcalNoiseAlgo::passTightHits(), HcalNoiseAlgo::passTightNoiseFilter(), HcalNoiseAlgo::passTightRatio(), HcalNoiseAlgo::passTightTiming(), HcalNoiseAlgo::passTightZeros(), CommonHcalNoiseRBXData::PassTS4TS5(), HcalNoiseAlgo::passZerosThreshold(), CommonHcalNoiseRBXData::ratio(), CommonHcalNoiseRBXData::RBXEMF(), CommonHcalNoiseRBXData::rbxTowers(), HcalNoiseSummary::rms10_, HcalNoiseSummary::rms25_, HcalNoiseSummary::tightnoisetwrs_, and CommonHcalNoiseRBXData::validRatio().
Referenced by produce().
{
// charge ratio
if(algo_.passRatioThreshold(data) && data.validRatio()) {
if(data.ratio()<summary.minE2Over10TS()) {
summary.mine2ts_ = data.e2ts();
summary.mine10ts_ = data.e10ts(); }
if(data.ratio()>summary.maxE2Over10TS()) {
summary.maxe2ts_ = data.e2ts();
summary.maxe10ts_ = data.e10ts();
}
}
// ADC zeros
if(algo_.passZerosThreshold(data)) {
if(data.numZeros()>summary.maxZeros()) {
summary.maxzeros_ = data.numZeros();
}
}
// hits count
if(data.numHPDHits() > summary.maxHPDHits()) {
summary.maxhpdhits_ = data.numHPDHits();
}
if(data.numRBXHits() > summary.maxRBXHits()) {
summary.maxrbxhits_ = data.numRBXHits();
}
if(data.numHPDNoOtherHits() > summary.maxHPDNoOtherHits()) {
summary.maxhpdhitsnoother_ = data.numHPDNoOtherHits();
}
// TS4TS5
if(data.PassTS4TS5() == false)
summary.hasBadRBXTS4TS5_ = true;
// hit timing
if(data.minLowEHitTime()<summary.min10GeVHitTime()) {
summary.min10_ = data.minLowEHitTime();
}
if(data.maxLowEHitTime()>summary.max10GeVHitTime()) {
summary.max10_ = data.maxLowEHitTime();
}
summary.rms10_ += data.lowEHitTimeSqrd();
summary.cnthit10_ += data.numLowEHits();
if(data.minHighEHitTime()<summary.min25GeVHitTime()) {
summary.min25_ = data.minHighEHitTime();
}
if(data.maxHighEHitTime()>summary.max25GeVHitTime()) {
summary.max25_ = data.maxHighEHitTime();
}
summary.rms25_ += data.highEHitTimeSqrd();
summary.cnthit25_ += data.numHighEHits();
// EMF
if(algo_.passEMFThreshold(data)) {
if(summary.minHPDEMF() > data.HPDEMF()) {
summary.minhpdemf_ = data.HPDEMF();
}
if(summary.minRBXEMF() > data.RBXEMF()) {
summary.minrbxemf_ = data.RBXEMF();
}
}
// summary flag
if(!algo_.passLooseRatio(data)) summary.filterstatus_ |= 0x1;
if(!algo_.passLooseHits(data)) summary.filterstatus_ |= 0x2;
if(!algo_.passLooseZeros(data)) summary.filterstatus_ |= 0x4;
if(!algo_.passLooseTiming(data)) summary.filterstatus_ |= 0x8;
if(!algo_.passTightRatio(data)) summary.filterstatus_ |= 0x100;
if(!algo_.passTightHits(data)) summary.filterstatus_ |= 0x200;
if(!algo_.passTightZeros(data)) summary.filterstatus_ |= 0x400;
if(!algo_.passTightTiming(data)) summary.filterstatus_ |= 0x800;
if(!algo_.passHighLevelNoiseFilter(data)) summary.filterstatus_ |= 0x10000;
// summary refvectors
JoinCaloTowerRefVectorsWithoutDuplicates join;
if(!algo_.passLooseNoiseFilter(data))
join(summary.loosenoisetwrs_, data.rbxTowers());
if(!algo_.passTightNoiseFilter(data))
join(summary.tightnoisetwrs_, data.rbxTowers());
if(!algo_.passHighLevelNoiseFilter(data))
join(summary.hlnoisetwrs_, data.rbxTowers());
return;
}
| void HcalNoiseInfoProducer::fillrechits | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HcalNoiseRBXArray & | array, | ||
| HcalNoiseSummary & | summary | ||
| ) | const [private] |
Definition at line 487 of file HcalNoiseInfoProducer.cc.
References CaloRecHit::detid(), CaloRecHit::energy(), HcalNoiseSummary::energyInLaserRegion_, HcalNoiseSummary::energyInNonLaserRegion_, Exception, reco::HcalNoiseRBXArray::findHPD(), CaloRecHit::flags(), HcalNoiseSummary::flatnoisee_, HcalNoiseSummary::flatnoiseet_, edm::EventSetup::get(), edm::Event::getByLabel(), CaloGeometry::getPosition(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), patZpeak::handle, HcalCaloFlagLabels::HBHEFlatNoise, HcalCaloFlagLabels::HBHEIsolatedNoise, HcalCaloFlagLabels::HBHESpikeNoise, HcalCaloFlagLabels::HBHETriangleNoise, HcalCaloFlagLabels::HBHETS4TS5Noise, HcalAcceptSeverityLevel_, HcalChannelStatus::HcalBadLaserSignal, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummaryR45, HcalRecHitFlagsToBeExcluded_, HcalNoiseSummary::hitsInLaserRegion_, HcalNoiseSummary::hitsInNonLaserRegion_, i, HBHERecHit::id(), HcalNoiseSummary::isolnoisee_, HcalNoiseSummary::isolnoiseet_, edm::HandleBase::isValid(), HcalNoiseSummary::nflatnoise_, HcalNoiseSummary::nisolnoise_, HcalNoiseSummary::nspikenoise_, HcalNoiseSummary::ntrianglenoise_, HcalNoiseSummary::nts4ts5noise_, edm::ESHandle< T >::product(), edm::errors::ProductNotFound, recHitCollName_, HcalNoiseSummary::rechitCount15_, HcalNoiseSummary::rechitCount_, HcalNoiseSummary::rechitEnergy15_, HcalNoiseSummary::rechitEnergy_, HcalSeverityLevelComputer::recoveredRecHit(), reco::HcalNoiseHPD::refrechitset_, HcalNoiseSummary::spikenoisee_, HcalNoiseSummary::spikenoiseet_, HcalNoiseSummary::trianglenoisee_, HcalNoiseSummary::trianglenoiseet_, HcalNoiseSummary::ts4ts5noisee_, and HcalNoiseSummary::ts4ts5noiseet_.
Referenced by produce().
{
// get the HCAL channel status map
edm::ESHandle<HcalChannelQuality> hcalChStatus;
iSetup.get<HcalChannelQualityRcd>().get( hcalChStatus );
const HcalChannelQuality* dbHcalChStatus = hcalChStatus.product();
// get the severity level computer
edm::ESHandle<HcalSeverityLevelComputer> hcalSevLvlComputerHndl;
iSetup.get<HcalSeverityLevelComputerRcd>().get(hcalSevLvlComputerHndl);
const HcalSeverityLevelComputer* hcalSevLvlComputer = hcalSevLvlComputerHndl.product();
// get the calo geometry
edm::ESHandle<CaloGeometry> pG;
iSetup.get<CaloGeometryRecord>().get(pG);
const CaloGeometry* geo = pG.product();
// get the rechits
edm::Handle<HBHERecHitCollection> handle;
iEvent.getByLabel(recHitCollName_, handle);
if(!handle.isValid()) {
throw edm::Exception(edm::errors::ProductNotFound)
<< " could not find HBHERecHitCollection named " << recHitCollName_ << "\n.";
return;
}
summary.rechitCount_ = 0;
summary.rechitCount15_ = 0;
summary.rechitEnergy_ = 0;
summary.rechitEnergy15_ = 0;
summary.hitsInLaserRegion_=0;
summary.hitsInNonLaserRegion_=0;
summary.energyInLaserRegion_=0;
summary.energyInNonLaserRegion_=0;
// loop over all of the rechit information
for(HBHERecHitCollection::const_iterator it=handle->begin(); it!=handle->end(); ++it) {
const HBHERecHit &rechit=(*it);
// skip bad rechits (other than those flagged by the isolated noise, triangle, flat, and spike algorithms)
const DetId id = rechit.detid();
uint32_t recHitFlag = rechit.flags();
uint32_t isolbitset = (1 << HcalCaloFlagLabels::HBHEIsolatedNoise);
uint32_t flatbitset = (1 << HcalCaloFlagLabels::HBHEFlatNoise);
uint32_t spikebitset = (1 << HcalCaloFlagLabels::HBHESpikeNoise);
uint32_t trianglebitset = (1 << HcalCaloFlagLabels::HBHETriangleNoise);
uint32_t ts4ts5bitset = (1 << HcalCaloFlagLabels::HBHETS4TS5Noise);
for(unsigned int i=0; i<HcalRecHitFlagsToBeExcluded_.size(); i++) {
uint32_t bitset = (1 << HcalRecHitFlagsToBeExcluded_[i]);
recHitFlag = (recHitFlag & bitset) ? recHitFlag-bitset : recHitFlag;
}
const uint32_t dbStatusFlag = dbHcalChStatus->getValues(id)->getValue();
// Ignore rechit if exclude bit set, regardless of severity of other bits
if ((dbStatusFlag & (1 <<HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary))==1) continue;
int severityLevel = hcalSevLvlComputer->getSeverityLevel(id, recHitFlag, dbStatusFlag);
bool isRecovered = hcalSevLvlComputer->recoveredRecHit(id, recHitFlag);
if(severityLevel!=0 && !isRecovered && severityLevel>static_cast<int>(HcalAcceptSeverityLevel_)) continue;
// do some rechit counting and energies
summary.rechitCount_ = summary.rechitCount_ + 1;
summary.rechitEnergy_ = summary.rechitEnergy_ + rechit.energy();
if ((dbStatusFlag & (1 <<HcalChannelStatus::HcalBadLaserSignal))==1) // hit comes from a region where no laser calibration pulse is normally seen
{
++summary.hitsInNonLaserRegion_;
summary.energyInNonLaserRegion_+=rechit.energy();
}
else // hit comes from region where laser calibration pulse is seen
{
++summary.hitsInLaserRegion_;
summary.energyInLaserRegion_+=rechit.energy();
}
if(rechit.energy() > 1.5)
{
summary.rechitCount15_ = summary.rechitCount15_ + 1;
summary.rechitEnergy15_ = summary.rechitEnergy15_ + rechit.energy();
}
// if it was ID'd as isolated noise, update the summary object
if(rechit.flags() & isolbitset) {
summary.nisolnoise_++;
summary.isolnoisee_ += rechit.energy();
GlobalPoint gp = geo->getPosition(rechit.id());
double et = rechit.energy()*gp.perp()/gp.mag();
summary.isolnoiseet_ += et;
}
if(rechit.flags() & flatbitset) {
summary.nflatnoise_++;
summary.flatnoisee_ += rechit.energy();
GlobalPoint gp = geo->getPosition(rechit.id());
double et = rechit.energy()*gp.perp()/gp.mag();
summary.flatnoiseet_ += et;
}
if(rechit.flags() & spikebitset) {
summary.nspikenoise_++;
summary.spikenoisee_ += rechit.energy();
GlobalPoint gp = geo->getPosition(rechit.id());
double et = rechit.energy()*gp.perp()/gp.mag();
summary.spikenoiseet_ += et;
}
if(rechit.flags() & trianglebitset) {
summary.ntrianglenoise_++;
summary.trianglenoisee_ += rechit.energy();
GlobalPoint gp = geo->getPosition(rechit.id());
double et = rechit.energy()*gp.perp()/gp.mag();
summary.trianglenoiseet_ += et;
}
if(rechit.flags() & ts4ts5bitset) {
if ((dbStatusFlag & (1 <<HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummaryR45))==0) // only add to TS4TS5 if the bit is not marked as "HcalCellExcludeFromHBHENoiseSummaryR45"
{
summary.nts4ts5noise_++;
summary.ts4ts5noisee_ += rechit.energy();
GlobalPoint gp = geo->getPosition(rechit.id());
double et = rechit.energy()*gp.perp()/gp.mag();
summary.ts4ts5noiseet_ += et;
}
}
// find the hpd that the rechit is in
HcalNoiseHPD& hpd=(*array.findHPD(rechit));
// create a persistent reference to the rechit
edm::Ref<HBHERecHitCollection> myRef(handle, it-handle->begin());
// store it in a place so that it remains sorted by energy
hpd.refrechitset_.insert(myRef);
} // end loop over rechits
// loop over all HPDs and transfer the information from refrechitset_ to rechits_;
for(HcalNoiseRBXArray::iterator rbxit=array.begin(); rbxit!=array.end(); ++rbxit) {
for(std::vector<HcalNoiseHPD>::iterator hpdit=rbxit->hpds_.begin(); hpdit!=rbxit->hpds_.end(); ++hpdit) {
// loop over all of the entries in the set and add them to rechits_
for(std::set<edm::Ref<HBHERecHitCollection>, RefHBHERecHitEnergyComparison>::const_iterator
it=hpdit->refrechitset_.begin(); it!=hpdit->refrechitset_.end(); ++it) {
hpdit->rechits_.push_back(*it);
}
}
}
// now the rechits in all the HPDs are sorted by energy!
return;
}
| void HcalNoiseInfoProducer::filltracks | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup, | ||
| HcalNoiseSummary & | summary | ||
| ) | const [private] |
Definition at line 684 of file HcalNoiseInfoProducer.cc.
References reco::TrackBase::eta(), edm::Event::getByLabel(), patZpeak::handle, edm::HandleBase::isValid(), maxTrackEta_, minTrackPt_, reco::TrackBase::p(), reco::TrackBase::pt(), trackCollName_, and HcalNoiseSummary::trackenergy_.
Referenced by produce().
{
edm::Handle<reco::TrackCollection> handle;
iEvent.getByLabel(trackCollName_, handle);
// don't throw exception, just return quietly
if(!handle.isValid()) {
// throw edm::Exception(edm::errors::ProductNotFound)
// << " could not find trackCollection named " << trackCollName_ << "\n.";
return;
}
summary.trackenergy_=0.0;
for(reco::TrackCollection::const_iterator iTrack = handle->begin(); iTrack!=handle->end(); ++iTrack) {
reco::Track trk=*iTrack;
if(trk.pt()<minTrackPt_ || fabs(trk.eta())>maxTrackEta_) continue;
summary.trackenergy_ += trk.p();
}
return;
}
| void HcalNoiseInfoProducer::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [override, private, virtual] |
Implements edm::EDProducer.
Definition at line 138 of file HcalNoiseInfoProducer.cc.
References algo_, data, CommonHcalNoiseRBXData::energy(), fillCaloTowers_, fillcalotwrs(), filldigis(), fillDigis_, fillOtherSummaryVariables(), fillrechits(), fillRecHits_, filltracks(), fillTracks_, HcalNoiseAlgo::isProblematic(), maxProblemRBXs_, minHighHitE_, minLowHitE_, minRecHitE_, HcalNoiseSummary::nproblemRBXs_, HcalNoiseAlgo::passHighLevelNoiseFilter(), HcalNoiseAlgo::passLooseNoiseFilter(), HcalNoiseAlgo::passTightNoiseFilter(), edm::Event::put(), edmLumisInFiles::summary, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, and TS4TS5UpperCut_.
{
// this is what we're going to actually write to the EDM
std::auto_ptr<HcalNoiseRBXCollection> result1(new HcalNoiseRBXCollection);
std::auto_ptr<HcalNoiseSummary> result2(new HcalNoiseSummary);
// define an empty HcalNoiseRBXArray that we're going to fill
HcalNoiseRBXArray rbxarray;
HcalNoiseSummary &summary=*result2;
// fill them with the various components
// digi assumes that rechit information is available
if(fillRecHits_) fillrechits(iEvent, iSetup, rbxarray, summary);
if(fillDigis_) filldigis(iEvent, iSetup, rbxarray, summary);
if(fillCaloTowers_) fillcalotwrs(iEvent, iSetup, rbxarray, summary);
if(fillTracks_) filltracks(iEvent, iSetup, summary);
// Why is this here? Shouldn't it have been in the filldigis method? Any reason for TotalCalibCharge to be defined outside filldigis(...) ?-- Jeff, 7/2/12
//if(fillDigis_) summary.calibCharge_ = TotalCalibCharge;
// select those RBXs which are interesting
// also look for the highest energy RBX
HcalNoiseRBXArray::iterator maxit=rbxarray.begin();
double maxenergy=-999;
bool maxwritten=false;
for(HcalNoiseRBXArray::iterator rit = rbxarray.begin(); rit!=rbxarray.end(); ++rit) {
HcalNoiseRBX &rbx=(*rit);
CommonHcalNoiseRBXData data(rbx, minRecHitE_, minLowHitE_, minHighHitE_, TS4TS5EnergyThreshold_,
TS4TS5UpperCut_, TS4TS5LowerCut_);
// find the highest energy rbx
if(data.energy()>maxenergy) {
maxenergy=data.energy();
maxit=rit;
maxwritten=false;
}
// find out if the rbx is problematic/noisy/etc.
bool writerbx = algo_.isProblematic(data) || !algo_.passLooseNoiseFilter(data) ||
!algo_.passTightNoiseFilter(data) || !algo_.passHighLevelNoiseFilter(data);
// fill variables in the summary object not filled elsewhere
fillOtherSummaryVariables(summary, data);
if(writerbx) {
summary.nproblemRBXs_++;
if(summary.nproblemRBXs_<=maxProblemRBXs_) {
result1->push_back(rbx);
if(maxit==rit) maxwritten=true;
}
}
} // end loop over rbxs
// if we still haven't written the maximum energy rbx, write it now
if(!maxwritten) {
HcalNoiseRBX &rbx=(*maxit);
// add the RBX to the event
result1->push_back(rbx);
}
// put the rbxcollection and summary into the EDM
iEvent.put(result1);
iEvent.put(result2);
return;
}
float reco::HcalNoiseInfoProducer::adc2fC[128] [private] |
Definition at line 111 of file HcalNoiseInfoProducer.h.
Referenced by filldigis(), and HcalNoiseInfoProducer().
Definition at line 92 of file HcalNoiseInfoProducer.h.
Referenced by fillOtherSummaryVariables(), and produce().
double reco::HcalNoiseInfoProducer::calibdigiHBHEthreshold_ [private] |
Definition at line 97 of file HcalNoiseInfoProducer.h.
Referenced by filldigis(), and HcalNoiseInfoProducer().
std::vector<int> reco::HcalNoiseInfoProducer::calibdigiHBHEtimeslices_ [private] |
Definition at line 98 of file HcalNoiseInfoProducer.h.
Referenced by filldigis(), and HcalNoiseInfoProducer().
double reco::HcalNoiseInfoProducer::calibdigiHFthreshold_ [private] |
Definition at line 100 of file HcalNoiseInfoProducer.h.
Referenced by filldigis(), and HcalNoiseInfoProducer().
std::vector<int> reco::HcalNoiseInfoProducer::calibdigiHFtimeslices_ [private] |
Definition at line 101 of file HcalNoiseInfoProducer.h.
Referenced by filldigis(), and HcalNoiseInfoProducer().
std::string reco::HcalNoiseInfoProducer::caloTowerCollName_ [private] |
Definition at line 86 of file HcalNoiseInfoProducer.h.
Referenced by fillcalotwrs(), and HcalNoiseInfoProducer().
std::string reco::HcalNoiseInfoProducer::digiCollName_ [private] |
Definition at line 84 of file HcalNoiseInfoProducer.h.
Referenced by filldigis(), and HcalNoiseInfoProducer().
bool reco::HcalNoiseInfoProducer::fillCaloTowers_ [private] |
Definition at line 73 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
bool reco::HcalNoiseInfoProducer::fillDigis_ [private] |
Definition at line 71 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
bool reco::HcalNoiseInfoProducer::fillRecHits_ [private] |
Definition at line 72 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
bool reco::HcalNoiseInfoProducer::fillTracks_ [private] |
Definition at line 74 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
uint32_t reco::HcalNoiseInfoProducer::HcalAcceptSeverityLevel_ [private] |
Definition at line 108 of file HcalNoiseInfoProducer.h.
Referenced by fillrechits(), and HcalNoiseInfoProducer().
std::vector<int> reco::HcalNoiseInfoProducer::HcalRecHitFlagsToBeExcluded_ [private] |
Definition at line 109 of file HcalNoiseInfoProducer.h.
Referenced by fillrechits(), and HcalNoiseInfoProducer().
int reco::HcalNoiseInfoProducer::maxCaloTowerIEta_ [private] |
Definition at line 80 of file HcalNoiseInfoProducer.h.
Referenced by fillcalotwrs(), and HcalNoiseInfoProducer().
int reco::HcalNoiseInfoProducer::maxProblemRBXs_ [private] |
Definition at line 77 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
double reco::HcalNoiseInfoProducer::maxTrackEta_ [private] |
Definition at line 81 of file HcalNoiseInfoProducer.h.
Referenced by filltracks(), and HcalNoiseInfoProducer().
double reco::HcalNoiseInfoProducer::minHighHitE_ [private] |
Definition at line 91 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
double reco::HcalNoiseInfoProducer::minLowHitE_ [private] |
Definition at line 91 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
double reco::HcalNoiseInfoProducer::minRecHitE_ [private] |
Definition at line 91 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
double reco::HcalNoiseInfoProducer::minTrackPt_ [private] |
Definition at line 82 of file HcalNoiseInfoProducer.h.
Referenced by filltracks(), and HcalNoiseInfoProducer().
std::string reco::HcalNoiseInfoProducer::recHitCollName_ [private] |
Definition at line 85 of file HcalNoiseInfoProducer.h.
Referenced by fillrechits(), and HcalNoiseInfoProducer().
double reco::HcalNoiseInfoProducer::TotalCalibCharge [private] |
Definition at line 89 of file HcalNoiseInfoProducer.h.
Referenced by filldigis().
std::string reco::HcalNoiseInfoProducer::trackCollName_ [private] |
Definition at line 87 of file HcalNoiseInfoProducer.h.
Referenced by filltracks(), and HcalNoiseInfoProducer().
double reco::HcalNoiseInfoProducer::TS4TS5EnergyThreshold_ [private] |
Definition at line 104 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
std::vector<std::pair<double, double> > reco::HcalNoiseInfoProducer::TS4TS5LowerCut_ [private] |
Definition at line 106 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
std::vector<std::pair<double, double> > reco::HcalNoiseInfoProducer::TS4TS5UpperCut_ [private] |
Definition at line 105 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer(), and produce().
bool reco::HcalNoiseInfoProducer::useCalibDigi_ [private] |
Definition at line 94 of file HcalNoiseInfoProducer.h.
Referenced by HcalNoiseInfoProducer().
1.7.3