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

void	set (const edm::EventSetup &es)

virtual	~HGCalRecHitWorkerSimple ()

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

virtual	~HGCalRecHitWorkerBaseClass ()

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

std::vector< double >	HGCEE_fCPerMIP_

bool	HGCEE_isSiFE_

double	HGCEE_keV2DIGI_

double	hgceeUncalib2GeV_

bool	HGCHEB_isSiFE_

double	HGCHEB_keV2DIGI_

double	hgchebUncalib2GeV_

std::vector< double >	HGCHEF_fCPerMIP_

bool	HGCHEF_isSiFE_

double	HGCHEF_keV2DIGI_

double	hgchefUncalib2GeV_

bool	killDeadChannels_

std::vector< float >	rcorr_

std::unique_ptr< HGCalRecHitSimpleAlgo >	rechitMaker_

std::vector< int >	v_chstatus_

std::vector< int >	v_DB_reco_flags_

Detailed Description

Definition at line 18 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_fCPerMIP_, HGCEE_isSiFE_, HGCEE_keV2DIGI_, hgceeUncalib2GeV_, HGCHEB_isSiFE_, HGCHEB_keV2DIGI_, hgchebUncalib2GeV_, HGCHEF_fCPerMIP_, HGCHEF_isSiFE_, HGCHEF_keV2DIGI_, hgchefUncalib2GeV_, rcorr_, rechitMaker_, and create_public_lumi_plots::weights.

                                                                          :
   HGCalRecHitWorkerBaseClass(ps) {
   rechitMaker_.reset( new HGCalRecHitSimpleAlgo() );
   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)
   std::vector<float> weights;
   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);
   }
   
 }

HGCalRecHitWorkerSimple::~HGCalRecHitWorkerSimple ( )

virtual

Definition at line 113 of file HGCalRecHitWorkerSimple.cc.

References DEFINE_EDM_PLUGIN.

113 {

114 }

Member Function Documentation

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

virtual

Implements HGCalRecHitWorkerBaseClass.

Definition at line 70 of file HGCalRecHitWorkerSimple.cc.

References ddds_, Exception, HcalEndcap, HGCEE, hgceeUncalib2GeV_, HGCHEB, hgchebUncalib2GeV_, HGCHEF, hgchefUncalib2GeV_, HGCUncalibratedRecHit::id(), edm::SortedCollection< T, SORT >::push_back(), DetId::rawId(), rcorr_, rechitMaker_, DetId::subdetId(), and HGCalDetId::wafer().

                                                              {
   DetId detid=uncalibRH.id();  
   uint32_t recoFlag = 0;
   //const std::vector<double>* fCPerMIP = nullptr;
     
   switch( detid.subdetId() ) {
   case HGCEE:
     rechitMaker_->setADCToGeVConstant(float(hgceeUncalib2GeV_) );
     //fCPerMIP = &HGCEE_fCPerMIP_;
     break;
   case HGCHEF:
     rechitMaker_->setADCToGeVConstant(float(hgchefUncalib2GeV_) );
     //fCPerMIP = &HGCHEF_fCPerMIP_;
     break;
   case HcalEndcap:
   case HGCHEB:
     rechitMaker_->setADCToGeVConstant(float(hgchebUncalib2GeV_) );
     break;  
   default:
     throw cms::Exception("NonHGCRecHit")
       << "Rechit with detid = " << detid.rawId() << " is not HGC!";
   }
   
   // make the rechit and put in the output collection
   if (recoFlag == 0) {    
     HGCRecHit myrechit( rechitMaker_->makeRecHit(uncalibRH, 0) );    
     int thk = -1;
     if( detid.subdetId() != HcalEndcap ) {
       HGCalDetId hid(detid);
       thk = ddds_[hid.subdetId()-3]->waferTypeL(hid.wafer());
       // units out of rechit maker are MIP * (GeV/fC)
       // so multiple
     }
     const double new_E = myrechit.energy()*(thk == -1 ? 1.0 : rcorr_[thk]);
     myrechit.setEnergy(new_E);
     result.push_back(myrechit);
   }
 
   return true;
 }

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

virtual

Implements HGCalRecHitWorkerBaseClass.

Definition at line 50 of file HGCalRecHitWorkerSimple.cc.

References HGCalTopology::dddConstants(), ddds_, edm::EventSetup::get(), HGCEE_isSiFE_, HGCHEF_isSiFE_, and HGCalGeometry::topology().

                                                          {
   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;  
 }