#include <HcalNoiseInfoProducer.h>

Inheritance diagram for reco::HcalNoiseInfoProducer:

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

	~HcalNoiseInfoProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

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 &) override

Private Attributes
float	adc2fC [128]

HcalNoiseAlgo	algo_

double	calibdigiHBHEthreshold_

std::vector< int >	calibdigiHBHEtimeslices_

double	calibdigiHFthreshold_

std::vector< int >	calibdigiHFtimeslices_

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_

bool	useCalibDigi_

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

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::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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, 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.

128 {

129 }

Member Function Documentation

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_, 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,
		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(), counter, 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, HcalBarrel, HcalChannelStatus::HcalCellExcludeFromHBHENoiseSummary, HcalEndcap, HcalForward, HcalCalibDetId::HOCrosstalk, 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;
 
   // 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.

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_, edm::hlt::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< class >::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
	)

overrideprivatevirtual

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;
 }