HGCalUncalibRecHitRecWeightsAlgo< C > Class Template Reference

#include <HGCalUncalibRecHitRecWeightsAlgo.h>

Public Member Functions
bool	isSiFESim () const
virtual HGCUncalibratedRecHit	makeRecHit (const C &dataFrame)
	Compute parameters. 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 22 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Constructor & Destructor Documentation

template<class C>

virtual HGCalUncalibRecHitRecWeightsAlgo< C >::~HGCalUncalibRecHitRecWeightsAlgo ( )

inlinevirtual

Definition at line 26 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ };

Member Function Documentation

template<class C>

bool HGCalUncalibRecHitRecWeightsAlgo< C >::isSiFESim ( ) const

inline

Definition at line 29 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ return isSiFESim_; }

template<class C>

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

inlinevirtual

Compute parameters.

Definition at line 53 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);
    
    // are we using the SiFE Simulation?
    if( isSiFESim_ ) {
      // mode == true means TDC readout was activated
      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_;
    jitter_    = double(sample.toa()) * toaLSBToNS_;
    LogDebug("HGCUncalibratedRecHit") << "TDC+: set the charge to: " << amplitude_ << ' ' << sample.data() 
                                          << ' ' << tdcLSB_ << std::endl
                                          << "TDC+: set the ToA to: " << jitter_ << ' ' 
                                          << sample.toa() << ' ' << toaLSBToNS_ << ' '
                                          << " flag=" << flag << std::endl;
      } else {
    amplitude_ = double(sample.data()) * adcLSB_;
    LogDebug("HGCUncalibratedRecHit") << "ADC+: set the charge to: " << amplitude_ << ' ' << sample.data() 
                                          << ' ' << adcLSB_ << ' ' << std::endl;
      }
    } else {
      amplitude_ = double(sample.data()) * adcLSB_;
      LogDebug("HGCUncalibratedRecHit") << "ADC+: set the charge to: " << amplitude_ << ' ' << sample.data() 
                         << ' ' << adcLSB_ << ' ' << std::endl;
    }
    
    int thickness = 1;
    if( ddd_ != nullptr ) {
      HGCalDetId hid(dataFrame.id());
      thickness = ddd_->waferTypeL(hid.wafer());
    }
    amplitude_ = amplitude_/fCPerMIP_[thickness-1];

    LogDebug("HGCUncalibratedRecHit") << "Final uncalibrated amplitude : " << amplitude_ << std::endl;
    
    return HGCUncalibratedRecHit( dataFrame.id(), amplitude_, pedestal_, jitter_, chi2_, flag);
   }

template<class C>

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

inline

Definition at line 31 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ adcLSB_ = adclsb; }

template<class C>

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

inline

Definition at line 38 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; 
  }

template<class C>

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

inline

Definition at line 28 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ isSiFESim_ = isSiFE; }

template<class C>

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

inline

Definition at line 32 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ tdcLSB_ = tdclsb; }

template<class C>

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

inline

Definition at line 36 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ tdcOnsetfC_ = tdcOnset; }

template<class C>

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

inline

Definition at line 34 of file HGCalUncalibRecHitRecWeightsAlgo.h.

{ toaLSBToNS_ = lsb2ns; }

template<class C>

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

inline

Definition at line 47 of file HGCalUncalibRecHitRecWeightsAlgo.h.

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

Member Data Documentation

template<class C>

double HGCalUncalibRecHitRecWeightsAlgo< C >::adcLSB_

private

Definition at line 102 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::set_ADCLSB().

template<class C>

const HGCalDDDConstants* HGCalUncalibRecHitRecWeightsAlgo< C >::ddd_

private

Definition at line 105 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry().

template<class C>

std::vector<double> HGCalUncalibRecHitRecWeightsAlgo< C >::fCPerMIP_

private

Definition at line 104 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::set_fCPerMIP().

template<class C>

bool HGCalUncalibRecHitRecWeightsAlgo< C >::isSiFESim_

private

Definition at line 103 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::isSiFESim(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::set_isSiFESim().

template<class C>

double HGCalUncalibRecHitRecWeightsAlgo< C >::tdcLSB_

private

Definition at line 102 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::set_TDCLSB().

template<class C>

double HGCalUncalibRecHitRecWeightsAlgo< C >::tdcOnsetfC_

private

Definition at line 102 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::set_tdcOnsetfC().

template<class C>

double HGCalUncalibRecHitRecWeightsAlgo< C >::toaLSBToNS_

private

Definition at line 102 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::set_toaLSBToNS().