#include <HGCalTriggerCellCalibration.h>

Public Member Functions
void	calibrateInGeV (l1t::HGCalTriggerCell &, int cellThickness)

void	calibrateInMipT (l1t::HGCalTriggerCell &, int cellThickness)

void	calibrateMipTinGeV (l1t::HGCalTriggerCell &)

	HGCalTriggerCellCalibration (const edm::ParameterSet &conf)

void	print ()

Private Attributes
std::vector< double >	dEdX_weights_

std::vector< double >	fCperMIP_ee_

std::vector< double >	fCperMIP_fh_

double	LSB_

std::vector< double >	thickCorr_

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 corr, dEdX_weights_, fCperMIP_ee_, fCperMIP_fh_, edm::ParameterSet::getParameter(), LSB_, and thickCorr_.

                                                                                             {
     
     LSB_ = beCodecConfig.getParameter<double>("cellLSB"); 
     fCperMIP_ee_ = beCodecConfig.getParameter<std::vector<double>>("fCperMIPee");
     fCperMIP_fh_ = beCodecConfig.getParameter<std::vector<double>>("fCperMIPfh");
     dEdX_weights_ = beCodecConfig.getParameter<std::vector<double>>("dEdXweights");
     thickCorr_ = beCodecConfig.getParameter<std::vector<double>>("thickCorr");
     
 
     for(auto corr : thickCorr_){
         if(corr <= 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,
		int	cellThickness
	)

Definition at line 85 of file HGCalTriggerCellCalibration.cc.

References calibrateInMipT(), and calibrateMipTinGeV().

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

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

void HGCalTriggerCellCalibration::calibrateInMipT	(	l1t::HGCalTriggerCell &	trgCell,
		int	cellThickness
	)

Definition at line 22 of file HGCalTriggerCellCalibration.cc.

References CustomPhysics_cfi::amplitude, l1t::HGCalTriggerCell::detId(), reco::LeafCandidate::eta(), fCperMIP_ee_, fCperMIP_fh_, HGCEE, HGCHEB, HGCHEF, l1t::L1Candidate::hwPt(), LSB_, l1t::HGCalTriggerCell::setMipPt(), DetId::subdetId(), and thickCorr_.

Referenced by calibrateInGeV().

 {
     
     HGCalDetId trgdetid( trgCell.detId() );
     int subdet = trgdetid.subdetId();
 
     /* get the hardware pT in ADC counts: */
     int hwPt = trgCell.hwPt();
 
     /* set the lowest signal bit and convert in charge amplitude: */
     double amplitude = hwPt * LSB_;  
 
     /* convert the charge amplitude in MIP: */
     if( subdet == HGCEE ){ 
         amplitude = amplitude / fCperMIP_ee_.at(cellThickness-1);
     }else if( subdet == HGCHEF ){
         amplitude = amplitude / fCperMIP_fh_.at(cellThickness-1);
     }else if( subdet == HGCHEB ){
         edm::LogWarning("DataNotFound") << "WARNING: the BH trgCells are not yet implemented";
     }
 
     /* correct the charge amplitude for the sensor thickness */
     double trgCellMipP = amplitude;
     if( thickCorr_.at( cellThickness-1 ) > 0 ){
         trgCellMipP = trgCellMipP / thickCorr_.at( cellThickness-1 ); 
     }
      
     double trgCellMipPt = trgCellMipP/cosh( trgCell.eta() ); 
 
     /* setting pT [mip] */
     trgCell.setMipPt( trgCellMipPt ) ;
 } 

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

Definition at line 56 of file HGCalTriggerCellCalibration.cc.

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

Referenced by calibrateInGeV().

 {
     const double MevToGeV(0.001);
 
     HGCalDetId trgdetid( trgCell.detId() );
     int trgCellLayer = trgdetid.layer();
     int subdet = trgdetid.subdetId();
 
     /* get the transverse momentum in mip units */
     double mipP = trgCell.mipPt() * cosh( trgCell.eta() );
 
     if( subdet == HGCHEF ){
             trgCellLayer = trgCellLayer + 28;
     }
    
     //weight the amplitude by the absorber coefficient in MeV/mip + bring it in GeV
     double trgCellE = mipP * dEdX_weights_.at(trgCellLayer) * MevToGeV;
 
     //assign the new energy to the four-vector of the trigger cell
     math::PtEtaPhiMLorentzVector calibP4(trgCellE/cosh( trgCell.eta() ), 
                                          trgCell.eta(), 
                                          trgCell.phi(), 
                                          trgCell.p4().M() );
     
     // overwriting the 4p with the calibrated 4p     
     trgCell.setP4( calibP4 );
 
 }