#include <HGCalUncalibRecHitRecWeightsAlgo.h>

Public Member Functions
virtual HGCUncalibratedRecHit	makeRecHit (const C &dataFrame)
	Compute parameters. More...

virtual void	set_ADCLSB (const double adclsb)

virtual void	set_isSiFESim (const bool isSiFE)

virtual void	set_TDCLSB (const double tdclsb)

virtual void	set_tdcOnsetfC (const double tdcOnset)

virtual void	set_toaLSBToNS (const double lsb2ns)

virtual	~HGCalUncalibRecHitRecWeightsAlgo ()

Private Attributes
double	adcLSB_

double	fCToMIP_

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 20 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Constructor & Destructor Documentation

template<class C>

virtual HGCalUncalibRecHitRecWeightsAlgo< C >::~HGCalUncalibRecHitRecWeightsAlgo ( )

inlinevirtual

Definition at line 24 of file HGCalUncalibRecHitRecWeightsAlgo.h.

24 { };

Member Function Documentation

template<class C>

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

inlinevirtual

Compute parameters.

Definition at line 36 of file HGCalUncalibRecHitRecWeightsAlgo.h.

Referenced by HGCalUncalibRecHitWorkerWeights::run1(), HGCalUncalibRecHitWorkerWeights::run2(), and HGCalUncalibRecHitWorkerWeights::run3().

                                                                  {
     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
     amplitude_ = ( std::floor(tdcOnsetfC_/adcLSB_) + 1.0 )* adcLSB_ + double(sample.data()) * tdcLSB_;
     jitter_    = double(sample.toa()) * toaLSBToNS_;
     LogDebug("HGCUncalibratedRecHit") << "TDC+: set the energy to: " << amplitude_ << ' ' << sample.data() 
                                           << ' ' << tdcLSB_ << std::endl
                                           << "TDC+: set the jitter to: " << jitter_ << ' ' 
                                           << sample.toa() << ' ' << toaLSBToNS_ << ' '
                                           << " flag=" << flag << std::endl;
       } 
       else {
     amplitude_ = double(sample.data()) * adcLSB_;
     LogDebug("HGCUncalibratedRecHit") << "ADC+: set the energy to: " << amplitude_ << ' ' << sample.data() 
                                           << ' ' << adcLSB_ << ' ' << std::endl;
       }
     }
     else {
       amplitude_ = double(sample.data()) * adcLSB_;
       LogDebug("HGCUncalibratedRecHit") << "ADC+: set the energy to: " << amplitude_ << ' ' << sample.data() 
                                         << ' ' << adcLSB_ << ' ' << std::endl;
     }
     
     LogDebug("HGCUncalibratedRecHit") << "Final uncalibrated amplitude : " << amplitude_ << std::endl;
     return HGCUncalibratedRecHit( dataFrame.id(), amplitude_, pedestal_, jitter_, chi2_, flag);
    }