#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)

Private Attributes
std::vector< double >	chargeCollectionEfficiency_

std::vector< double >	dEdX_weights_

std::vector< double >	fCperMIP_

double	lsb_

std::vector< double >	thicknessCorrection_

HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 9 of file HGCalTriggerCellCalibration.h.

Constructor & Destructor Documentation

◆ HGCalTriggerCellCalibration()

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

Definition at line 5 of file HGCalTriggerCellCalibration.cc.

     : lsb_(conf.getParameter<double>("lsb")),
       fCperMIP_(conf.getParameter<std::vector<double>>("fCperMIP")),
       chargeCollectionEfficiency_(conf.getParameter<edm::ParameterSet>("chargeCollectionEfficiency")
                                       .getParameter<std::vector<double>>("values")),
       thicknessCorrection_(conf.getParameter<std::vector<double>>("thicknessCorrection")),
       dEdX_weights_(conf.getParameter<std::vector<double>>("dEdXweights")) {
   for (const auto& fCperMIP : fCperMIP_) {
     if (fCperMIP <= 0) {
       edm::LogWarning("DivisionByZero") << "WARNING: zero or negative MIP->fC correction factor. It won't be "
                                            "applied to correct trigger cell energies.";
     }
   }
   for (const auto& cce : chargeCollectionEfficiency_) {
     if (cce <= 0) {
       edm::LogWarning("DivisionByZero") << "WARNING: zero or negative cell-thickness correction factor. It won't be "
                                            "applied to correct trigger cell energies.";
     }
   }
   for (const auto& thickCorr : thicknessCorrection_) {
     if (thickCorr <= 0) {
       edm::LogWarning("DivisionByZero") << "WARNING: zero or negative cell-thickness correction factor. It won't be "
                                            "applied to correct trigger cell energies.";
     }
   }
 }

References chargeCollectionEfficiency_, hgcalVFEProducer_cfi::fCperMIP, fCperMIP_, and thicknessCorrection_.

Member Function Documentation

◆ calibrateInGeV()

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

Definition at line 99 of file HGCalTriggerCellCalibration.cc.

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

References calibrateInMipT(), and calibrateMipTinGeV().

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

◆ calibrateInMipT()

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

Definition at line 32 of file HGCalTriggerCellCalibration.cc.

                                                                                     {
   DetId trgdetid(trgCell.detId());
   bool isSilicon = triggerTools_.isSilicon(trgdetid);
   constexpr int kScintillatorIndex = 0;
   unsigned thickness = isSilicon ? triggerTools_.thicknessIndex(trgdetid, true) : kScintillatorIndex;
   if (thickness >= fCperMIP_.size()) {
     throw cms::Exception("OutOfBound") << "Trying to access thickness index " << thickness
                                        << " in fCperMIP, which is of size " << fCperMIP_.size();
   }
   if (thickness >= chargeCollectionEfficiency_.size()) {
     throw cms::Exception("OutOfBound") << "Trying to access thickness index " << thickness
                                        << " in chargeCollectionEfficiency, which is of size "
                                        << chargeCollectionEfficiency_.size();
   }
  
   /* 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 * lsb_;
  
   if (chargeCollectionEfficiency_[thickness] > 0) {
     amplitude /= chargeCollectionEfficiency_[thickness];
   }
  
   /* convert the charge amplitude in MIP: */
   double trgCellMipP = amplitude;
  
   if (fCperMIP_[thickness] > 0) {
     trgCellMipP /= fCperMIP_[thickness];
   }
  
   /* compute the transverse-mip */
   double trgCellMipPt = trgCellMipP / std::cosh(trgCell.eta());
  
   /* setting pT [mip] */
   trgCell.setMipPt(trgCellMipPt);
 }

References CustomPhysics_cfi::amplitude, chargeCollectionEfficiency_, l1t::HGCalTriggerCell::detId(), reco::LeafCandidate::eta(), Exception, fCperMIP_, l1t::L1Candidate::hwPt(), HGCalTriggerTools::isSilicon(), lsb_, l1t::HGCalTriggerCell::setMipPt(), Calorimetry_cff::thickness, HGCalTriggerTools::thicknessIndex(), and triggerTools_.

Referenced by calibrateInGeV().

◆ calibrateMipTinGeV()

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

Definition at line 71 of file HGCalTriggerCellCalibration.cc.

                                                                                        {
   constexpr double MevToGeV(0.001);
   double trgCellEt(0.);
  
   DetId trgdetid(trgCell.detId());
   unsigned trgCellLayer = triggerTools_.layerWithOffset(trgdetid);
   bool isSilicon = triggerTools_.isSilicon(trgdetid);
   constexpr int kScintillatorIndex = 0;
   unsigned thickness = isSilicon ? triggerTools_.thicknessIndex(trgdetid, true) : kScintillatorIndex;
   if (thickness >= thicknessCorrection_.size()) {
     throw cms::Exception("OutOfBound") << "Trying to access thickness index " << thickness
                                        << " in thicknessCorrection, which is of size " << thicknessCorrection_.size();
   }
  
   /* weight the amplitude by the absorber coefficient in MeV/mip + bring it in
    * GeV */
   trgCellEt = trgCell.mipPt() * MevToGeV;
   trgCellEt *= dEdX_weights_.at(trgCellLayer);
  
   /* correct for the cell-thickness */
   if (thicknessCorrection_[thickness] > 0) {
     trgCellEt /= thicknessCorrection_[thickness];
   }
  
   math::PtEtaPhiMLorentzVector calibP4(trgCellEt, trgCell.eta(), trgCell.phi(), 0.);
   trgCell.setP4(calibP4);
 }