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

std::vector< double >	dEdX_weights_Nose_

double	fCperMIP_

std::vector< double >	fCperMIPnose_

double	LSB_scintillator_MIP_

double	LSB_silicon_fC_

double	thickCorr_

std::vector< double >	thickCorrNose_

HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 15 of file HGCalTriggerCellCalibration.h.

Constructor & Destructor Documentation

◆ HGCalTriggerCellCalibration()

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

Definition at line 6 of file HGCalTriggerCellCalibration.cc.

     : LSB_silicon_fC_(conf.getParameter<double>("siliconCellLSB_fC")),
       LSB_scintillator_MIP_(conf.getParameter<double>("scintillatorCellLSB_MIP")),
       fCperMIP_(conf.getParameter<double>("fCperMIP")),
       fCperMIPnose_(conf.getParameter<std::vector<double>>("fCperMIPnose")),
       thickCorr_(conf.getParameter<double>("thickCorr")),
       thickCorrNose_(conf.getParameter<std::vector<double>>("thickCorrNose")),
       dEdX_weights_(conf.getParameter<std::vector<double>>("dEdXweights")),
       dEdX_weights_Nose_(conf.getParameter<std::vector<double>>("dEdXweightsNose")) {
   if (fCperMIP_ <= 0 or fCperMIPnose_.at(0) <= 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 or thickCorrNose_.at(0) <= 0) {
     edm::LogWarning("DivisionByZero") << "WARNING: the cell-thickness correction factor is zero or negative. It won't "
                                          "be applied to correct trigger cell energies.";
   }
 }

References fCperMIP_, fCperMIPnose_, or, thickCorr_, and thickCorrNose_.

Member Function Documentation

◆ calibrateInGeV()

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

Definition at line 104 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 25 of file HGCalTriggerCellCalibration.cc.

                                                                                     {
   DetId trgdetid(trgCell.detId());
   bool isSilicon = triggerTools_.isSilicon(trgdetid);
   bool isNose = triggerTools_.isNose(trgdetid);
   int thickness = isSilicon ? triggerTools_.thicknessIndex(trgdetid, false) : 0;
  
   /* 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 (isNose and fCperMIPnose_[thickness] > 0) {
     trgCellMipP /= fCperMIPnose_[thickness];
   } else if (isSilicon && fCperMIP_ > 0) {
     trgCellMipP /= fCperMIP_;
   }
  
   /* compute the transverse-mip */
   double trgCellMipPt = trgCellMipP / cosh(trgCell.eta());
  
   /* setting pT [mip] */
   trgCell.setMipPt(trgCellMipPt);
 }

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

Referenced by calibrateInGeV().

◆ calibrateMipTinGeV()

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

Definition at line 56 of file HGCalTriggerCellCalibration.cc.

                                                                                        {
   const double MevToGeV(0.001);
   double trgCellEt(0.);
  
   DetId trgdetid(trgCell.detId());
   bool isNose = triggerTools_.isNose(trgdetid);
   unsigned trgCellLayer = triggerTools_.layerWithOffset(trgdetid);
   bool isSilicon = triggerTools_.isSilicon(trgdetid);
   int thickness = isSilicon ? triggerTools_.thicknessIndex(trgdetid, false) : 0;
  
   if (isNose) {
     if (dEdX_weights_Nose_.at(trgCellLayer) == 0.) {
       throw cms::Exception("BadConfiguration - HFNose")
           << "Trigger cell energy forced to 0 by calibration coefficients.\n"
           << "The configuration should be changed. \n";
     }
   } else {
     if (dEdX_weights_.at(trgCellLayer) == 0.) {
       throw cms::Exception("BadConfiguration - HGCAL ")
           << "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 */
   trgCellEt = trgCell.mipPt() * MevToGeV;
   trgCellEt *= (isNose) ? dEdX_weights_Nose_.at(trgCellLayer) : dEdX_weights_.at(trgCellLayer);
  
   /* correct for the cell-thickness */
   if (triggerTools_.isSilicon(trgdetid)) {
     if (isNose && thickCorrNose_[thickness] > 0)
       trgCellEt /= thickCorrNose_[thickness];
     else if (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);
 }