#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 Types
enum	detectortype { hgcee = 1, hgcfh = 2, hgcbh = 3, hgchfnose = 4 }

Protected Attributes
std::array< const HGCalDDDConstants *, 4 >	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_

std::vector< double >	hgcHFNose_cce_

std::vector< double >	hgcHFNose_fCPerMIP_

bool	hgcHFNose_isSiFE_

double	hgcHFNose_keV2DIGI_

std::vector< double >	hgcHFNose_noise_fC_

double	hgchfnoseUncalib2GeV_

bool	killDeadChannels_

uint32_t	rangeMask_

uint32_t	rangeMatch_

std::vector< double >	rcorr_

std::vector< double >	rcorrNose_

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_

std::vector< float >	weightsNose_

Detailed Description

Definition at line 19 of file HGCalRecHitWorkerSimple.h.

Member Enumeration Documentation

enum HGCalRecHitWorkerSimple::detectortype

protected

Enumerator
hgcee
hgcfh
hgcbh
hgchfnose

Definition at line 28 of file HGCalRecHitWorkerSimple.h.

28 { hgcee = 1, hgcfh = 2, hgcbh = 3, hgchfnose = 4 };

HGCalRecHitWorkerSimple::hgcfh

Definition: HGCalRecHitWorkerSimple.h:28

HGCalRecHitWorkerSimple::hgcbh

Definition: HGCalRecHitWorkerSimple.h:28

HGCalRecHitWorkerSimple::hgchfnose

Definition: HGCalRecHitWorkerSimple.h:28

HGCalRecHitWorkerSimple::hgcee

Definition: HGCalRecHitWorkerSimple.h:28

Constructor & Destructor Documentation

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

Definition at line 10 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_, hgcHFNose_cce_, hgcHFNose_fCPerMIP_, hgcHFNose_isSiFE_, hgcHFNose_keV2DIGI_, hgcHFNose_noise_fC_, hgchfnoseUncalib2GeV_, rangeMask_, rangeMatch_, rcorr_, rcorrNose_, rechitMaker_, tools_, weights_, and weightsNose_.

                                                                           : 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_;
 
   // HGChfnose constants
   hgcHFNose_keV2DIGI_ = ps.getParameter<double>("HGCHFNose_keV2DIGI");
   hgcHFNose_fCPerMIP_ = ps.getParameter<std::vector<double> >("HGCHFNose_fCPerMIP");
   hgcHFNose_isSiFE_ = ps.getParameter<bool>("HGCHFNose_isSiFE");
   hgchfnoseUncalib2GeV_ = keV2GeV / hgcHFNose_keV2DIGI_;
 
   // layer weights (from Valeri/Arabella)
   const auto& dweights = ps.getParameter<std::vector<double> >("layerWeights");
   for (auto weight : dweights) {
     weights_.push_back(weight);
   }
   const auto& weightnose = ps.getParameter<std::vector<double> >("layerNoseWeights");
   for (auto const& weight : weightnose)
     weightsNose_.emplace_back(weight);
 
   rechitMaker_->setLayerWeights(weights_);
   rechitMaker_->setNoseLayerWeights(weightsNose_);
 
   // 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);
   }
   const auto& rcorrnose = ps.getParameter<std::vector<double> >("thicknessNoseCorrection");
   rcorrNose_.clear();
   rcorrNose_.push_back(1.f);
   for (auto corr : rcorrnose) {
     rcorrNose_.push_back(1.0 / corr);
   }
 
   hgcEE_noise_fC_ = ps.getParameter<edm::ParameterSet>("HGCEE_noise_fC").getParameter<std::vector<double> >("values");
   hgcEE_cce_ = ps.getParameter<edm::ParameterSet>("HGCEE_cce").getParameter<std::vector<double> >("values");
   hgcHEF_noise_fC_ = ps.getParameter<edm::ParameterSet>("HGCHEF_noise_fC").getParameter<std::vector<double> >("values");
   hgcHEF_cce_ = ps.getParameter<edm::ParameterSet>("HGCHEF_cce").getParameter<std::vector<double> >("values");
   hgcHEB_noise_MIP_ = ps.getParameter<edm::ParameterSet>("HGCHEB_noise_MIP").getParameter<double>("noise_MIP");
   hgcHFNose_noise_fC_ =
       ps.getParameter<edm::ParameterSet>("HGCHFNose_noise_fC").getParameter<std::vector<double> >("values");
   hgcHFNose_cce_ = ps.getParameter<edm::ParameterSet>("HGCHFNose_cce").getParameter<std::vector<double> >("values");
 
   // 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 200 of file HGCalRecHitWorkerSimple.cc.

References DEFINE_EDM_PLUGIN.

200 {}

Member Function Documentation

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

overridevirtual

Implements HGCalRecHitWorkerBaseClass.

Definition at line 104 of file HGCalRecHitWorkerSimple.cc.

References ddds_, DetId::det(), MillePedeFileConverter_cfg::e, CaloRecHit::energy(), Exception, DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, hgcbh, HGCEE, hgcee, hgcEE_cce_, hgcEE_fCPerMIP_, hgcEE_noise_fC_, hgceeUncalib2GeV_, hgcfh, HGCHEB, hgcHEB_noise_MIP_, hgchebUncalib2GeV_, HGCHEF, hgcHEF_cce_, hgcHEF_fCPerMIP_, hgcHEF_noise_fC_, hgchefUncalib2GeV_, hgchfnose, hgcHFNose_cce_, hgcHFNose_fCPerMIP_, hgcHFNose_noise_fC_, hgchfnoseUncalib2GeV_, HGCUncalibratedRecHit::id(), edm::SortedCollection< T, SORT >::push_back(), rangeMask_, rangeMatch_, DetId::rawId(), rcorr_, rcorrNose_, rechitMaker_, DetId::subdetId(), Calorimetry_cff::thickness, tools_, HFNoseDetId::type(), HGCSiliconDetId::type(), weights_, and weightsNose_.

                                                                {
   DetId detid = uncalibRH.id();
   // don't produce rechit if detid is a ghost one
 
   if (detid.det() == DetId::Forward || detid.det() == DetId::Hcal) {
     if ((detid & rangeMask_) == rangeMatch_)
       return false;
   }
 
   int thickness = -1;
   float sigmaNoiseGeV = 0.f;
   unsigned int layer = tools_->getLayerWithOffset(detid);
   float cce_correction = 1.0;
   int idtype(0);
 
   switch (detid.det()) {
     case DetId::HGCalEE:
       idtype = hgcee;
       thickness = 1 + HGCSiliconDetId(detid).type();
       break;
     case DetId::HGCalHSi:
       idtype = hgcfh;
       thickness = 1 + HGCSiliconDetId(detid).type();
       break;
     case DetId::HGCalHSc:
       idtype = hgcbh;
       break;
     default:
       switch (detid.subdetId()) {
         case HGCEE:
           idtype = hgcee;
           thickness = ddds_[detid.subdetId() - 3]->waferTypeL(HGCalDetId(detid).wafer());
           break;
         case HGCHEF:
           idtype = hgcfh;
           thickness = ddds_[detid.subdetId() - 3]->waferTypeL(HGCalDetId(detid).wafer());
           break;
         case HcalEndcap:
           [[fallthrough]];
         case HGCHEB:
           idtype = hgcbh;
           break;
         case HFNose:
           idtype = hgchfnose;
           thickness = 1 + HFNoseDetId(detid).type();
           break;
         default:
           break;
       }
   }
   switch (idtype) {
     case hgcee:
       rechitMaker_->setADCToGeVConstant(float(hgceeUncalib2GeV_));
       cce_correction = hgcEE_cce_[thickness - 1];
       sigmaNoiseGeV =
           1e-3 * weights_[layer] * rcorr_[thickness] * hgcEE_noise_fC_[thickness - 1] / hgcEE_fCPerMIP_[thickness - 1];
       break;
     case hgcfh:
       rechitMaker_->setADCToGeVConstant(float(hgchefUncalib2GeV_));
       cce_correction = hgcHEF_cce_[thickness - 1];
       sigmaNoiseGeV = 1e-3 * weights_[layer] * rcorr_[thickness] * hgcHEF_noise_fC_[thickness - 1] /
                       hgcHEF_fCPerMIP_[thickness - 1];
       break;
     case hgcbh:
       rechitMaker_->setADCToGeVConstant(float(hgchebUncalib2GeV_));
       sigmaNoiseGeV = 1e-3 * hgcHEB_noise_MIP_ * weights_[layer];
       break;
     case hgchfnose:
       rechitMaker_->setADCToGeVConstant(float(hgchfnoseUncalib2GeV_));
       cce_correction = hgcHFNose_cce_[thickness - 1];
       sigmaNoiseGeV = 1e-3 * weightsNose_[layer] * rcorrNose_[thickness] * hgcHFNose_noise_fC_[thickness - 1] /
                       hgcHFNose_fCPerMIP_[thickness - 1];
       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));
   double new_E = myrechit.energy();
   if (detid.det() == DetId::Forward && detid.subdetId() == ForwardSubdetector::HFNose) {
     new_E *= (thickness == -1 ? 1.0 : rcorrNose_[thickness]) / cce_correction;
   } else {
     new_E *= (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 77 of file HGCalRecHitWorkerSimple.cc.

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

                                                          {
   tools_->getEventSetup(es);
   rechitMaker_->set(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;
   if (hgcHFNose_isSiFE_) {
     edm::ESHandle<HGCalGeometry> hgchfnoseGeoHandle;
     es.get<IdealGeometryRecord>().get("HGCalHFNoseSensitive", hgchfnoseGeoHandle);
     ddds_[3] = &(hgchfnoseGeoHandle->topology().dddConstants());
   } else {
     ddds_[3] = nullptr;
   }
 }