#include <HGCalUncalibRecHitRecWeightsAlgo.h>

Public Member Functions
bool	isSiFESim () const

virtual HGCUncalibratedRecHit	makeRecHit (const C &dataFrame)
	Compute HGCUncalibratedRecHit from DataFrame. More...

void	set_ADCLSB (const double adclsb)

void	set_fCPerMIP (const std::vector< double > &fCPerMIP)

void	set_isSiFESim (const bool isSiFE)

void	set_TDCLSB (const double tdclsb)

void	set_tdcOnsetfC (const double tdcOnset)

void	set_toaLSBToNS (const double lsb2ns)

void	setGeometry (const HGCalGeometry *geom)

virtual	~HGCalUncalibRecHitRecWeightsAlgo ()

Private Attributes
double	adcLSB_

const HGCalDDDConstants *	ddd_

std::vector< double >	fCPerMIP_

bool	isSiFESim_

double	tdcLSB_

double	tdcOnsetfC_

double	toaLSBToNS_

Detailed Description

template<class C>
class HGCalUncalibRecHitRecWeightsAlgo< C >

Template used to produce fast-track HGCAL Reco, weight=1

Author: Valeri Andreev

Definition at line 21 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Constructor & Destructor Documentation

◆ ~HGCalUncalibRecHitRecWeightsAlgo()

template<class C>

virtual HGCalUncalibRecHitRecWeightsAlgo< C >::~HGCalUncalibRecHitRecWeightsAlgo ( )

inlinevirtual

Definition at line 24 of file HGCalUncalibRecHitRecWeightsAlgo.h.

24 {};

Member Function Documentation

◆ isSiFESim()

template<class C>

bool HGCalUncalibRecHitRecWeightsAlgo< C >::isSiFESim ( ) const

inline

Definition at line 27 of file HGCalUncalibRecHitRecWeightsAlgo.h.

27 { return isSiFESim_; }

◆ makeRecHit()

template<class C>

virtual HGCUncalibratedRecHit HGCalUncalibRecHitRecWeightsAlgo< C >::makeRecHit ( const C & dataFrame )

inlinevirtual

Compute HGCUncalibratedRecHit from DataFrame.

Definition at line 51 of file HGCalUncalibRecHitRecWeightsAlgo.h.

                                                                {
     double amplitude_(-1.), pedestal_(-1.), jitter_(-1.), chi2_(-1.);
     uint32_t flag = 0;
  
     constexpr int iSample = 2;  //only in-time sample
     const auto& sample = dataFrame.sample(iSample);
  
     // Were digis done w/ the complete digitization and using signal shape?
     // (originally done only for the silicon, while for scitillator it was trivial. Fomr 11_ also scinti uses shape)
     if (isSiFESim_) {
       // mode == true: TDC readout was activated and amplitude comes from TimeOverThreshold
       if (sample.mode()) {
         flag = !sample.threshold();  // raise flag if busy cell
         // LG (23/06/2015):
         // to get a continuous energy spectrum we must add here the maximum value in fC ever
         // reported by the ADC. Namely: floor(tdcOnset/adcLSB_) * adcLSB_
         // need to increment by one so TDC doesn't overlap with ADC last bin
         // LG (11/04/2016):
         // offset the TDC upwards to reflect the bin center
         amplitude_ = (std::floor(tdcOnsetfC_ / adcLSB_) + 1.0) * adcLSB_ + (double(sample.data()) + 0.5) * tdcLSB_;
  
         if (sample.getToAValid()) {
           jitter_ = double(sample.toa()) * toaLSBToNS_;
         }
       } else {
         amplitude_ = double(sample.data()) * adcLSB_;  // why do we not have +0.5 here ?
         if (sample.getToAValid()) {
           jitter_ = double(sample.toa()) * toaLSBToNS_;
         }
       }  //isSiFESim_
     }    //mode()
  
     // trivial digitization, i.e. no signal shape
     else {
       amplitude_ = double(sample.data()) * adcLSB_;
     }
  
     int thickness = (ddd_ != nullptr) ? ddd_->waferType(dataFrame.id()) : 0;
     amplitude_ = amplitude_ / fCPerMIP_[thickness];
  
     LogDebug("HGCUncalibratedRecHit") << "isSiFESim_: " << isSiFESim_ << " ADC+: set the charge to: " << amplitude_
                                       << ' ' << sample.data() << ' ' << adcLSB_ << ' '
                                       << "               TDC+: set the ToA to: " << jitter_ << ' ' << sample.toa()
                                       << ' ' << toaLSBToNS_ << ' ' << tdcLSB_
                                       << "               getToAValid(): " << sample.getToAValid()
                                       << " mode(): " << sample.mode() << std::endl;
     LogDebug("HGCUncalibratedRecHit") << "Final uncalibrated amplitude : " << amplitude_ << std::endl;
  
     return HGCUncalibratedRecHit(dataFrame.id(), amplitude_, pedestal_, jitter_, chi2_, flag);
   }

◆ set_ADCLSB()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::set_ADCLSB ( const double adclsb )

inline

Definition at line 29 of file HGCalUncalibRecHitRecWeightsAlgo.h.

29 { adcLSB_ = adclsb; }

Referenced by configureIt().

◆ set_fCPerMIP()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::set_fCPerMIP ( const std::vector< double > & fCPerMIP )

inline

Definition at line 36 of file HGCalUncalibRecHitRecWeightsAlgo.h.

                                                        {
     if (std::any_of(fCPerMIP.cbegin(), fCPerMIP.cend(), [](double conv) { return conv <= 0.0; })) {
       throw cms::Exception("BadConversionFactor") << "At least one of fCPerMIP is zero!" << std::endl;
     }
     fCPerMIP_ = fCPerMIP;
   }

Referenced by configureIt().

◆ set_isSiFESim()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::set_isSiFESim ( const bool isSiFE )

inline

Definition at line 26 of file HGCalUncalibRecHitRecWeightsAlgo.h.

26 { isSiFESim_ = isSiFE; }

Referenced by configureIt().

◆ set_TDCLSB()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::set_TDCLSB ( const double tdclsb )

inline

Definition at line 30 of file HGCalUncalibRecHitRecWeightsAlgo.h.

30 { tdcLSB_ = tdclsb; }

Referenced by configureIt().

◆ set_tdcOnsetfC()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::set_tdcOnsetfC ( const double tdcOnset )

inline

Definition at line 34 of file HGCalUncalibRecHitRecWeightsAlgo.h.

34 { tdcOnsetfC_ = tdcOnset; }

Referenced by configureIt().

◆ set_toaLSBToNS()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::set_toaLSBToNS ( const double lsb2ns )

inline

Definition at line 32 of file HGCalUncalibRecHitRecWeightsAlgo.h.

32 { toaLSBToNS_ = lsb2ns; }

Referenced by configureIt().

◆ setGeometry()

template<class C>

void HGCalUncalibRecHitRecWeightsAlgo< C >::setGeometry ( const HGCalGeometry * geom )

inline

Definition at line 43 of file HGCalUncalibRecHitRecWeightsAlgo.h.

                                               {
     if (geom)
       ddd_ = &(geom->topology().dddConstants());
     else
       ddd_ = nullptr;
   }