#include <HGCalRecHitWorkerSimple.h>

Inheritance diagram for HGCalRecHitWorkerSimple:

Public Member Functions
	HGCalRecHitWorkerSimple (const edm::ParameterSet &)

bool	run (const edm::Event &evt, const HGCUncalibratedRecHit &uncalibRH, HGCRecHitCollection &result) override

void	set (const edm::EventSetup &es) override

	~HGCalRecHitWorkerSimple () override

Public Member Functions inherited from HGCalRecHitWorkerBaseClass
	HGCalRecHitWorkerBaseClass (const edm::ParameterSet &)

virtual	~HGCalRecHitWorkerBaseClass ()

Protected Attributes
std::array< const HGCalDDDConstants *, 3 >	ddds_

std::vector< double >	hgcEE_cce_

std::vector< double >	hgcEE_fCPerMIP_

bool	hgcEE_isSiFE_

double	hgcEE_keV2DIGI_

std::vector< double >	hgcEE_noise_fC_

double	hgceeUncalib2GeV_

bool	hgcHEB_isSiFE_

double	hgcHEB_keV2DIGI_

double	hgcHEB_noise_MIP_

double	hgchebUncalib2GeV_

std::vector< double >	hgcHEF_cce_

std::vector< double >	hgcHEF_fCPerMIP_

bool	hgcHEF_isSiFE_

double	hgcHEF_keV2DIGI_

std::vector< double >	hgcHEF_noise_fC_

double	hgchefUncalib2GeV_

bool	killDeadChannels_

uint32_t	rangeMask_

uint32_t	rangeMatch_

std::vector< double >	rcorr_

std::unique_ptr< HGCalRecHitSimpleAlgo >	rechitMaker_

std::unique_ptr< hgcal::RecHitTools >	tools_

std::vector< int >	v_chstatus_

std::vector< int >	v_DB_reco_flags_

std::vector< float >	weights_

Detailed Description

Definition at line 19 of file HGCalRecHitWorkerSimple.h.

Constructor & Destructor Documentation

HGCalRecHitWorkerSimple::HGCalRecHitWorkerSimple ( const edm::ParameterSet & ps )

Definition at line 11 of file HGCalRecHitWorkerSimple.cc.

References constexpr, corr, MillePedeFileConverter_cfg::e, f, edm::ParameterSet::getParameter(), hgcEE_cce_, hgcEE_fCPerMIP_, hgcEE_isSiFE_, hgcEE_keV2DIGI_, hgcEE_noise_fC_, hgceeUncalib2GeV_, hgcHEB_isSiFE_, hgcHEB_keV2DIGI_, hgcHEB_noise_MIP_, hgchebUncalib2GeV_, hgcHEF_cce_, hgcHEF_fCPerMIP_, hgcHEF_isSiFE_, hgcHEF_keV2DIGI_, hgcHEF_noise_fC_, hgchefUncalib2GeV_, rangeMask_, rangeMatch_, rcorr_, rechitMaker_, tools_, and weights_.

                                                                          :
         HGCalRecHitWorkerBaseClass(ps)
 {
     rechitMaker_.reset(new HGCalRecHitSimpleAlgo());
     tools_.reset(new hgcal::RecHitTools());
     constexpr float keV2GeV = 1e-6;
     // HGCee constants
     hgcEE_keV2DIGI_ = ps.getParameter<double>("HGCEE_keV2DIGI");
     hgcEE_fCPerMIP_ = ps.getParameter < std::vector<double> > ("HGCEE_fCPerMIP");
     hgcEE_isSiFE_ = ps.getParameter<bool>("HGCEE_isSiFE");
     hgceeUncalib2GeV_ = keV2GeV / hgcEE_keV2DIGI_;
 
     // HGChef constants
     hgcHEF_keV2DIGI_ = ps.getParameter<double>("HGCHEF_keV2DIGI");
     hgcHEF_fCPerMIP_ = ps.getParameter < std::vector<double> > ("HGCHEF_fCPerMIP");
     hgcHEF_isSiFE_ = ps.getParameter<bool>("HGCHEF_isSiFE");
     hgchefUncalib2GeV_ = keV2GeV / hgcHEF_keV2DIGI_;
 
     // HGCheb constants
     hgcHEB_keV2DIGI_ = ps.getParameter<double>("HGCHEB_keV2DIGI");
     hgcHEB_isSiFE_ = ps.getParameter<bool>("HGCHEB_isSiFE");
     hgchebUncalib2GeV_ = keV2GeV / hgcHEB_keV2DIGI_;
 
     // layer weights (from Valeri/Arabella)
     const auto& dweights = ps.getParameter < std::vector<double> > ("layerWeights");
     for (auto weight : dweights)
     {
         weights_.push_back(weight);
     }
 
     rechitMaker_->setLayerWeights(weights_);
 
     // residual correction for cell thickness
     const auto& rcorr = ps.getParameter < std::vector<double> > ("thicknessCorrection");
     rcorr_.clear();
     rcorr_.push_back(1.f);
     for (auto corr : rcorr)
     {
         rcorr_.push_back(1.0 / corr);
     }
 
 
     hgcEE_noise_fC_ = ps.getParameter < std::vector<double> > ("HGCEE_noise_fC");
     hgcEE_cce_ = ps.getParameter< std::vector<double> > ("HGCEE_cce");
     hgcHEF_noise_fC_ = ps.getParameter < std::vector<double> > ("HGCHEF_noise_fC");
     hgcHEF_cce_ = ps.getParameter< std::vector<double> > ("HGCHEF_cce");
     hgcHEB_noise_MIP_ = ps.getParameter<double>("HGCHEB_noise_MIP");
 
     // don't produce rechit if detid is a ghost one
     rangeMatch_ = ps.getParameter<uint32_t>("rangeMatch");
     rangeMask_  = ps.getParameter<uint32_t>("rangeMask");
 }

HGCalRecHitWorkerSimple::~HGCalRecHitWorkerSimple ( )

override

Definition at line 147 of file HGCalRecHitWorkerSimple.cc.

References DEFINE_EDM_PLUGIN.

148 {

149 }

Member Function Documentation

bool HGCalRecHitWorkerSimple::run	(	const edm::Event &	evt,
		const HGCUncalibratedRecHit &	uncalibRH,
		HGCRecHitCollection &	result
	)

overridevirtual

Implements HGCalRecHitWorkerBaseClass.

Definition at line 90 of file HGCalRecHitWorkerSimple.cc.

References ddds_, MillePedeFileConverter_cfg::e, CaloRecHit::energy(), Exception, HcalEndcap, HGCEE, hgcEE_cce_, hgcEE_fCPerMIP_, hgcEE_noise_fC_, hgceeUncalib2GeV_, HGCHEB, hgcHEB_noise_MIP_, hgchebUncalib2GeV_, HGCHEF, hgcHEF_cce_, hgcHEF_fCPerMIP_, hgcHEF_noise_fC_, hgchefUncalib2GeV_, HGCUncalibratedRecHit::id(), edm::SortedCollection< T, SORT >::push_back(), rangeMask_, rangeMatch_, DetId::rawId(), rcorr_, rechitMaker_, DetId::subdetId(), tools_, HGCalDetId::wafer(), and weights_.

 {
     DetId detid = uncalibRH.id();
 // don't produce rechit if detid is a ghost one
 
     if((detid & rangeMask_) == rangeMatch_)
         return false;
 
     int thickness = -1;
     float sigmaNoiseGeV = 0.f;
     unsigned int layer = tools_->getLayerWithOffset(detid);
     HGCalDetId hid;
     float cce_correction = 1.0;
 
     switch (detid.subdetId())
     {
     case HGCEE:
         rechitMaker_->setADCToGeVConstant(float(hgceeUncalib2GeV_));
         hid = detid;
         thickness = ddds_[hid.subdetId() - 3]->waferTypeL(hid.wafer());
         cce_correction = hgcEE_cce_[thickness - 1];
         sigmaNoiseGeV = 1e-3 * weights_[layer] * rcorr_[thickness]
                     * hgcEE_noise_fC_[thickness - 1] / hgcEE_fCPerMIP_[thickness - 1];
         break;
     case HGCHEF:
         rechitMaker_->setADCToGeVConstant(float(hgchefUncalib2GeV_));
         hid = detid;
         thickness = ddds_[hid.subdetId() - 3]->waferTypeL(hid.wafer());
         cce_correction = hgcHEF_cce_[thickness - 1];
         sigmaNoiseGeV = 1e-3 * weights_[layer] * rcorr_[thickness]
                     * hgcHEF_noise_fC_[thickness - 1] / hgcHEF_fCPerMIP_[thickness - 1];
         break;
     case HcalEndcap:
     case HGCHEB:
         rechitMaker_->setADCToGeVConstant(float(hgchebUncalib2GeV_));
         hid = detid;
         sigmaNoiseGeV = 1e-3 * hgcHEB_noise_MIP_ * weights_[layer];
         break;
     default:
         throw cms::Exception("NonHGCRecHit") << "Rechit with detid = " << detid.rawId()
                 << " is not HGC!";
     }
 
     // make the rechit and put in the output collection
 
     HGCRecHit myrechit(rechitMaker_->makeRecHit(uncalibRH, 0));
     const double new_E = myrechit.energy() * (thickness == -1 ? 1.0 : rcorr_[thickness])/cce_correction;
 
 
     myrechit.setEnergy(new_E); 
     myrechit.setSignalOverSigmaNoise(new_E/sigmaNoiseGeV);
     result.push_back(myrechit);
 
     return true;
 }

void HGCalRecHitWorkerSimple::set ( const edm::EventSetup & es )

overridevirtual

Implements HGCalRecHitWorkerBaseClass.

Definition at line 64 of file HGCalRecHitWorkerSimple.cc.

References HGCalTopology::dddConstants(), ddds_, edm::EventSetup::get(), hgcEE_isSiFE_, hgcHEF_isSiFE_, tools_, and HGCalGeometry::topology().

 {
     tools_->getEventSetup(es);
     if (hgcEE_isSiFE_)
     {
         edm::ESHandle < HGCalGeometry > hgceeGeoHandle;
         es.get<IdealGeometryRecord>().get("HGCalEESensitive", hgceeGeoHandle);
         ddds_[0] = &(hgceeGeoHandle->topology().dddConstants());
     }
     else
     {
         ddds_[0] = nullptr;
     }
     if (hgcHEF_isSiFE_)
     {
         edm::ESHandle < HGCalGeometry > hgchefGeoHandle;
         es.get<IdealGeometryRecord>().get("HGCalHESiliconSensitive", hgchefGeoHandle);
         ddds_[1] = &(hgchefGeoHandle->topology().dddConstants());
     }
     else
     {
         ddds_[1] = nullptr;
     }
     ddds_[2] = nullptr;
 }