#include <HGCalTriggerCellCalibration.h>

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

void	calibrateInMipT (l1t::HGCalTriggerCell &)

void	calibrateMipTinGeV (l1t::HGCalTriggerCell &)

	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_

Static Private Attributes
static const unsigned	kLayers_ = kLayersEE_+kLayersFH_+kLayersBH_

static const unsigned	kLayersBH_ = 12

static const unsigned	kLayersEE_ = 28

static const unsigned	kLayersFH_ = 12

Detailed Description

Definition at line 17 of file HGCalTriggerCellCalibration.h.

Constructor & Destructor Documentation

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

Definition at line 4 of file HGCalTriggerCellCalibration.cc.

References dEdX_weights_, fCperMIP_, edm::ParameterSet::getParameter(), LSB_scintillator_MIP_, LSB_silicon_fC_, and thickCorr_.

                                                                                             {
     
     LSB_silicon_fC_ = beCodecConfig.getParameter<double>("siliconCellLSB_fC"); 
     LSB_scintillator_MIP_ = beCodecConfig.getParameter<double>("scintillatorCellLSB_MIP"); 
     fCperMIP_ = beCodecConfig.getParameter<double>("fCperMIP");
     dEdX_weights_ = beCodecConfig.getParameter<std::vector<double>>("dEdXweights");
     thickCorr_ = beCodecConfig.getParameter<double>("thickCorr");
     
 
 
     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 )

Definition at line 87 of file HGCalTriggerCellCalibration.cc.

References calibrateInMipT(), and calibrateMipTinGeV().

Referenced by HGCalTriggerBackend::HGCalTriggerSimCluster< FECODEC, DATA >::run().

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

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

Definition at line 24 of file HGCalTriggerCellCalibration.cc.

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

Referenced by calibrateInGeV().

 {
     
     HGCalDetId trgdetid( 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 * (trgdetid.subdetId()==HGCHEB ? 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( trgdetid.subdetId()!=HGCHEB && 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 )

Definition at line 52 of file HGCalTriggerCellCalibration.cc.

References dEdX_weights_, l1t::HGCalTriggerCell::detId(), reco::LeafCandidate::eta(), HGCHEB, HGCHEF, kLayersEE_, kLayersFH_, HGCalDetId::layer(), l1t::HGCalTriggerCell::mipPt(), reco::LeafCandidate::phi(), reco::LeafCandidate::setP4(), and thickCorr_.

Referenced by calibrateInGeV().

 {
     const double MevToGeV(0.001);
 
     HGCalDetId trgdetid( trgCell.detId() );
     int trgCellLayer = trgdetid.layer();
     int subdet = trgdetid.subdetId();
 
 
     if( subdet == HGCHEF ){
             trgCellLayer = trgCellLayer + kLayersEE_;
     }
     else if( subdet == HGCHEB ){
             trgCellLayer = trgCellLayer + kLayersEE_ + kLayersFH_;
     }
    
     /* 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( subdet!=HGCHEB && 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 );
 
 }