#include <HcalNoiseInfoProducer.h>

Inheritance diagram for reco::HcalNoiseInfoProducer:

Public Member Functions
	HcalNoiseInfoProducer (const edm::ParameterSet &)

	~HcalNoiseInfoProducer ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Member Functions
virtual void	beginJob ()

virtual void	beginRun (edm::Run &, const edm::EventSetup &)

virtual void	endJob ()

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

void	fillcalotwrs (edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &) const

void	filldigis (edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &)

void	fillOtherSummaryVariables (HcalNoiseSummary &summary, const CommonHcalNoiseRBXData &data) const

void	fillrechits (edm::Event &, const edm::EventSetup &, HcalNoiseRBXArray &, HcalNoiseSummary &) const

void	filltracks (edm::Event &, const edm::EventSetup &, HcalNoiseSummary &) const

virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
float	adc2fC [128]

HcalNoiseAlgo	algo_

std::string	caloTowerCollName_

std::string	digiCollName_

bool	fillCaloTowers_

bool	fillDigis_

bool	fillRecHits_

bool	fillTracks_

uint32_t	HcalAcceptSeverityLevel_

std::vector< int >	HcalRecHitFlagsToBeExcluded_

int	maxCaloTowerIEta_

int	maxProblemRBXs_

double	maxTrackEta_

double	minHighHitE_

double	minLowHitE_

double	minRecHitE_

double	minTrackPt_

std::string	recHitCollName_

double	TotalCalibCharge

std::string	trackCollName_

double	TS4TS5EnergyThreshold_

std::vector< std::pair< double, double > >	TS4TS5LowerCut_

std::vector< std::pair< double, double > >	TS4TS5UpperCut_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 38 of file HcalNoiseInfoProducer.h.

Constructor & Destructor Documentation

HcalNoiseInfoProducer::HcalNoiseInfoProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 33 of file HcalNoiseInfoProducer.cc.

References adc2fC, caloTowerCollName_, digiCollName_, edm::ParameterSet::exists(), fillCaloTowers_, fillDigis_, fillRecHits_, fillTracks_, edm::ParameterSet::getParameter(), HcalAcceptSeverityLevel_, HcalRecHitFlagsToBeExcluded_, i, maxCaloTowerIEta_, maxProblemRBXs_, maxTrackEta_, minHighHitE_, minLowHitE_, minRecHitE_, minTrackPt_, recHitCollName_, python.multivaluedict::sort(), trackCollName_, TS4TS5EnergyThreshold_, TS4TS5LowerCut_, and TS4TS5UpperCut_.

                                                                            : 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);
   }
 
   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 105 of file HcalNoiseInfoProducer.cc.

106 {

107 }

Member Function Documentation

void HcalNoiseInfoProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 185 of file HcalNoiseInfoProducer.cc.

 {
   return;
 }

void HcalNoiseInfoProducer::beginRun	(	edm::Run &	,
		const edm::EventSetup &
	)

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 201 of file HcalNoiseInfoProducer.cc.

 {
   return;
 }

void HcalNoiseInfoProducer::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 192 of file HcalNoiseInfoProducer.cc.

 {
   return;
 }

void HcalNoiseInfoProducer::endRun	(	edm::Run &	,
		const edm::EventSetup &
	)

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 208 of file HcalNoiseInfoProducer.cc.

 {
   return;
 }

void HcalNoiseInfoProducer::fillcalotwrs	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		HcalNoiseRBXArray &	array,
		HcalNoiseSummary &	summary
	)		const

private

Definition at line 564 of file HcalNoiseInfoProducer.cc.

References caloTowerCollName_, CaloTower::emEnergy(), HcalNoiseSummary::emenergy_, edm::hlt::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
	)

private

Definition at line 306 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_, digiCollName_, HcalSeverityLevelComputer::dropChannel(), edm::RefVector< C, T, F >::end(), edm::hlt::Exception, reco::HcalNoiseRBXArray::findHPD(), reco::HcalNoiseRBXArray::findRBX(), edm::EventSetup::get(), edm::Event::getByLabel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), patZpeak::handle, i, HBHEDataFrame::id(), edm::HandleBase::isValid(), reco::HcalNoiseHPD::maxZeros_, HcalCalibrations::pedestal(), edm::ESHandle< class >::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;
 
   // get the conditions and channel quality
   edm::ESHandle<HcalDbService> conditions;
   iSetup.get<HcalDbRecord>().get(conditions);
   const HcalQIEShape* shape = conditions->getHcalShape();
   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;
 
     // get the calibrations and coder
     const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
     const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
     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);
         HcalCoderDb coder(*channelCoder, *shape);
 
         CaloSamples tool;
         coder.adc2fC(*digi, tool);
 
         for(int i = 0; i < (int)digi->size(); i++)
            TotalCalibCharge = TotalCalibCharge + tool[i];
         */
 
         for(int i = 0; i < (int)digi->size(); i++)
            TotalCalibCharge = TotalCalibCharge + adc2fC[digi->sample(i).adc()&0xff];
      }
   }
 
   return;
 }

void HcalNoiseInfoProducer::fillOtherSummaryVariables	(	HcalNoiseSummary &	summary,
		const CommonHcalNoiseRBXData &	data
	)		const

private

Definition at line 215 of file HcalNoiseInfoProducer.cc.

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 433 of file HcalNoiseInfoProducer.cc.

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;
 
   // 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();
     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(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) {
       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 605 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
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 116 of file HcalNoiseInfoProducer.cc.

References algo_, HcalNoiseSummary::calibCharge_, 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, TotalCalibCharge, 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);
   if(fillCaloTowers_) fillcalotwrs(iEvent, iSetup, rbxarray, summary);
   if(fillTracks_)     filltracks(iEvent, iSetup, summary);
 
   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;
 }