#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 }

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.

Member Enumeration Documentation

enum HGCalRecHitWorkerSimple::detectortype

protected

Enumerator
hgcee
hgcfh
hgcbh

Definition at line 29 of file HGCalRecHitWorkerSimple.h.

29 { hgcee=1, hgcfh=2, hgcbh=3};

HGCalRecHitWorkerSimple::hgcfh

Definition: HGCalRecHitWorkerSimple.h:29

HGCalRecHitWorkerSimple::hgcbh

Definition: HGCalRecHitWorkerSimple.h:29

HGCalRecHitWorkerSimple::hgcee

Definition: HGCalRecHitWorkerSimple.h:29

Constructor & Destructor Documentation

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

Definition at line 9 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<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>("value");
 
     // 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 169 of file HGCalRecHitWorkerSimple.cc.

References DEFINE_EDM_PLUGIN.

170 {

171 }

Member Function Documentation

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

overridevirtual

Implements HGCalRecHitWorkerBaseClass.

Definition at line 88 of file HGCalRecHitWorkerSimple.cc.

References ddds_, DetId::det(), MillePedeFileConverter_cfg::e, CaloRecHit::energy(), Exception, DetId::Forward, DetId::Hcal, HcalEndcap, 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_, HGCUncalibratedRecHit::id(), edm::SortedCollection< T, SORT >::push_back(), rangeMask_, rangeMatch_, DetId::rawId(), rcorr_, rechitMaker_, DetId::subdetId(), tools_, HGCSiliconDetId::type(), and weights_.

 {
     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:
       case HGCHEB:
     idtype = hgcbh; 
     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;
     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 62 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;
 }