HGCalTriggerCellCalibration.cc

Go to the documentation of this file.

#include "L1Trigger/L1THGCal/interface/HGCalTriggerCellCalibration.h"
#include "DataFormats/ForwardDetId/interface/HFNoseTriggerDetId.h"
#include "DataFormats/ForwardDetId/interface/HGCalTriggerDetId.h"
 
// class constructor
HGCalTriggerCellCalibration::HGCalTriggerCellCalibration(const edm::ParameterSet& conf)
    : 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.";
  }
}
 
void HGCalTriggerCellCalibration::calibrateInMipT(l1t::HGCalTriggerCell& trgCell) const {
  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);
}
 
void HGCalTriggerCellCalibration::calibrateMipTinGeV(l1t::HGCalTriggerCell& trgCell) const {
  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);
}
 
void HGCalTriggerCellCalibration::calibrateInGeV(l1t::HGCalTriggerCell& trgCell) const {
  /* calibrate from ADC count to transverse mip */
  calibrateInMipT(trgCell);
 
  /* calibrate from mip count to GeV */
  calibrateMipTinGeV(trgCell);
}