#include <HGCalTriggerCellCalibration.h>

Public Member Functions
void	calibrateInGeV (l1t::HGCalTriggerCell &) const

void	calibrateInMipT (l1t::HGCalTriggerCell &) const

void	calibrateMipTinGeV (l1t::HGCalTriggerCell &) const

void	eventSetup (const edm::EventSetup &es)

	HGCalTriggerCellCalibration (const edm::ParameterSet &conf)

void	print ()

Private Attributes
std::vector< double >	dEdX_weights_

double	fCperMIP_

double	LSB_scintillator_MIP_

double	LSB_silicon_fC_

double	thickCorr_

HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 15 of file HGCalTriggerCellCalibration.h.

Constructor & Destructor Documentation

HGCalTriggerCellCalibration::HGCalTriggerCellCalibration ( const edm::ParameterSet & conf )

Definition at line 4 of file HGCalTriggerCellCalibration.cc.

References fCperMIP_, and thickCorr_.

     : LSB_silicon_fC_(conf.getParameter<double>("siliconCellLSB_fC")),
       LSB_scintillator_MIP_(conf.getParameter<double>("scintillatorCellLSB_MIP")),
       fCperMIP_(conf.getParameter<double>("fCperMIP")),
       thickCorr_(conf.getParameter<double>("thickCorr")),
       dEdX_weights_(conf.getParameter<std::vector<double>>("dEdXweights")) {
   if (fCperMIP_ <= 0) {
     edm::LogWarning("DivisionByZero") << "WARNING: the MIP->fC correction factor is zero or negative. It won't be "
                                          "applied to correct trigger cell energies.";
   }
   if (thickCorr_ <= 0) {
     edm::LogWarning("DivisionByZero") << "WARNING: the cell-thickness correction factor is zero or negative. It won't "
                                          "be applied to correct trigger cell energies.";
   }
 }

Member Function Documentation

void HGCalTriggerCellCalibration::calibrateInGeV ( l1t::HGCalTriggerCell & trgCell ) const

Definition at line 73 of file HGCalTriggerCellCalibration.cc.

References calibrateInMipT(), and calibrateMipTinGeV().

Referenced by HGCalConcentratorCoarsenerImpl::assignCoarseTriggerCellEnergy(), HGCalConcentratorSuperTriggerCellImpl::assignSuperTriggerCellEnergyAndPosition(), and eventSetup().

                                                                                    {
   /* calibrate from ADC count to transverse mip */
   calibrateInMipT(trgCell);
 
   /* calibrate from mip count to GeV */
   calibrateMipTinGeV(trgCell);
 }

void HGCalTriggerCellCalibration::calibrateInMipT ( l1t::HGCalTriggerCell & trgCell ) const

Definition at line 20 of file HGCalTriggerCellCalibration.cc.

References CustomPhysics_cfi::amplitude, l1t::HGCalTriggerCell::detId(), reco::LeafCandidate::eta(), fCperMIP_, l1t::L1Candidate::hwPt(), HGCalTriggerTools::isSilicon(), LSB_scintillator_MIP_, LSB_silicon_fC_, l1t::HGCalTriggerCell::setMipPt(), and triggerTools_.

Referenced by calibrateInGeV(), and eventSetup().

                                                                                     {
   bool isSilicon = triggerTools_.isSilicon(trgCell.detId());
 
   /* get the hardware pT in ADC counts: */
   int hwPt = trgCell.hwPt();
 
   // Convert ADC to charge in fC (in EE+FH) or in MIPs (in BH)
   double amplitude = hwPt * (!isSilicon ? LSB_scintillator_MIP_ : LSB_silicon_fC_);
 
   // The responses of the different cell thicknesses have been equalized
   // to the 200um response in the front-end. So there is only one global
   // fCperMIP and thickCorr here
   /* convert the charge amplitude in MIP: */
   double trgCellMipP = amplitude;
   if (isSilicon && fCperMIP_ > 0) {
     trgCellMipP /= fCperMIP_;
   }
 
   /* compute the transverse-mip */
   double trgCellMipPt = trgCellMipP / cosh(trgCell.eta());
 
   /* setting pT [mip] */
   trgCell.setMipPt(trgCellMipPt);
 }

void HGCalTriggerCellCalibration::calibrateMipTinGeV ( l1t::HGCalTriggerCell & trgCell ) const

Definition at line 45 of file HGCalTriggerCellCalibration.cc.

References dEdX_weights_, l1t::HGCalTriggerCell::detId(), reco::LeafCandidate::eta(), Exception, HGCalTriggerTools::isSilicon(), HGCalTriggerTools::layerWithOffset(), l1t::HGCalTriggerCell::mipPt(), reco::LeafCandidate::phi(), reco::LeafCandidate::setP4(), thickCorr_, and triggerTools_.

Referenced by calibrateInGeV(), and eventSetup().

                                                                                        {
   const double MevToGeV(0.001);
 
   DetId trgdetid(trgCell.detId());
   unsigned trgCellLayer = triggerTools_.layerWithOffset(trgdetid);
 
   if (dEdX_weights_.at(trgCellLayer) == 0.) {
     throw cms::Exception("BadConfiguration")
         << "Trigger cell energy forced to 0 by calibration coefficients.\n"
         << "The configuration should be changed. "
         << "Discarded layers should be defined in hgcalTriggerGeometryESProducer.TriggerGeometry.DisconnectedLayers "
            "and not with calibration coefficients = 0\n";
   }
 
   /* weight the amplitude by the absorber coefficient in MeV/mip + bring it in GeV */
   double trgCellEt = trgCell.mipPt() * dEdX_weights_.at(trgCellLayer) * MevToGeV;
 
   /* correct for the cell-thickness */
   if (triggerTools_.isSilicon(trgdetid) && thickCorr_ > 0) {
     trgCellEt /= thickCorr_;
   }
   /* assign the new energy to the four-vector of the trigger cell */
   math::PtEtaPhiMLorentzVector calibP4(trgCellEt, trgCell.eta(), trgCell.phi(), 0.);
 
   /* overwriting the 4p with the calibrated 4p */
   trgCell.setP4(calibP4);
 }