FitterFuncs::PulseShapeFunctor Class Reference

#include <PulseShapeFitOOTPileupCorrection.h>

Public Member Functions
double	doublePulseShapeFunc (const double *x)
double	EvalPulse (const std::vector< double > &pars)
int	getcntNANinfit ()
	PulseShapeFunctor (const HcalPulseShapes::Shape &pulse, bool iPedestalConstraint, bool iTimeConstraint, bool iAddPulseJitter, bool iAddTimeSlew, double iPulseJitter, double iTimeMean, double iTimeSig, double iPedMean, double iPedSig, double iNoise)
void	setDefaultcntNANinfit ()
void	setpsFiterry (double *erry)
void	setpsFiterry2 (double *erry2)
void	setpsFitslew (double *slew)
void	setpsFitx (double *x)
void	setpsFity (double *y)
double	sigma (double ifC)
double	singlePulseShapeFunc (const double *x)
double	triplePulseShapeFunc (const double *x)
	~PulseShapeFunctor ()

Private Member Functions
void	funcHPDShape (std::array< float, HcalConst::maxSamples > &ntmpbin, const double &pulseTime, const double &pulseHeight, const double &slew)

Private Attributes
std::vector< float >	acc25nsVec
std::vector< float >	accVarLenIdxMinusOneVec
std::vector< float >	accVarLenIdxZEROVec
bool	addPulseJitter_
int	cntNANinfit
std::vector< float >	diff25nsItvlVec
std::vector< float >	diffVarItvlIdxMinusOneVec
std::vector< float >	diffVarItvlIdxZEROVec
double	invertpedSig2_
double	invertpedSig_
double	inverttimeSig2_
double	inverttimeSig_
double	noise_
bool	pedestalConstraint_
double	pedMean_
double	pedSig_
double	psFit_erry [HcalConst::maxSamples]
double	psFit_erry2 [HcalConst::maxSamples]
double	psFit_slew [HcalConst::maxSamples]
double	psFit_x [HcalConst::maxSamples]
double	psFit_y [HcalConst::maxSamples]
std::array< float, HcalConst::maxPSshapeBin >	pulse_hist
double	pulseJitter_
bool	timeConstraint_
double	timeMean_
double	timeShift_
double	timeSig_
bool	unConstrainedFit_

Detailed Description

Definition at line 33 of file PulseShapeFitOOTPileupCorrection.h.

Constructor & Destructor Documentation

FitterFuncs::PulseShapeFunctor::PulseShapeFunctor	(	const HcalPulseShapes::Shape &	pulse,
		bool	iPedestalConstraint,
		bool	iTimeConstraint,
		bool	iAddPulseJitter,
		bool	iAddTimeSlew,
		double	iPulseJitter,
		double	iTimeMean,
		double	iTimeSig,
		double	iPedMean,
		double	iPedSig,
		double	iNoise
	)

Definition at line 10 of file PulseShapeFitOOTPileupCorrection.cc.

References acc25nsVec, accVarLenIdxMinusOneVec, accVarLenIdxZEROVec, addPulseJitter_, diff25nsItvlVec, diffVarItvlIdxMinusOneVec, diffVarItvlIdxZEROVec, i, invertpedSig2_, invertpedSig_, inverttimeSig2_, inverttimeSig_, j, HcalConst::maxPSshapeBin, HcalConst::maxSamples, noise_, HcalConst::nsPerBX, pedestalConstraint_, pedMean_, pedSig_, psFit_erry, psFit_erry2, psFit_slew, psFit_x, psFit_y, pulse_hist, pulseJitter_, timeConstraint_, timeMean_, timeShift_, and timeSig_.

                                      : 
    cntNANinfit(0),
    acc25nsVec(HcalConst::maxPSshapeBin), diff25nsItvlVec(HcalConst::maxPSshapeBin),
    accVarLenIdxZEROVec(HcalConst::nsPerBX), diffVarItvlIdxZEROVec(HcalConst::nsPerBX), 
    accVarLenIdxMinusOneVec(HcalConst::nsPerBX), diffVarItvlIdxMinusOneVec(HcalConst::nsPerBX) {
    //The raw pulse
    for(int i=0;i<HcalConst::maxPSshapeBin;++i) pulse_hist[i] = pulse(i);
    // Accumulate 25ns for each starting point of 0, 1, 2, 3...
    for(int i=0; i<HcalConst::maxPSshapeBin; ++i){
      for(int j=i; j<i+HcalConst::nsPerBX; ++j){  //sum over HcalConst::nsPerBXns from point i
    acc25nsVec[i] += ( j < HcalConst::maxPSshapeBin? pulse_hist[j] : pulse_hist[HcalConst::maxPSshapeBin-1]);
      }
      diff25nsItvlVec[i] = ( i+HcalConst::nsPerBX < HcalConst::maxPSshapeBin? pulse_hist[i+HcalConst::nsPerBX] - pulse_hist[i] : pulse_hist[HcalConst::maxPSshapeBin-1] - pulse_hist[i]);
    }
    // Accumulate different ns for starting point of index either 0 or -1
    for(int i=0; i<HcalConst::nsPerBX; ++i){
      if( i==0 ){
    accVarLenIdxZEROVec[0] = pulse_hist[0];
    accVarLenIdxMinusOneVec[i] = pulse_hist[0];
      } else{
    accVarLenIdxZEROVec[i] = accVarLenIdxZEROVec[i-1] + pulse_hist[i];
    accVarLenIdxMinusOneVec[i] = accVarLenIdxMinusOneVec[i-1] + pulse_hist[i-1];
      }
      diffVarItvlIdxZEROVec[i] = pulse_hist[i+1] - pulse_hist[0];
      diffVarItvlIdxMinusOneVec[i] = pulse_hist[i] - pulse_hist[0];
    }
    for(int i = 0; i < HcalConst::maxSamples; i++) { 
      psFit_x[i]      = 0;
      psFit_y[i]      = 0;
      psFit_erry[i]   = 1.;
      psFit_erry2[i]  = 1.;
      psFit_slew [i]  = 0.;
    }
    //Constraints
    pedestalConstraint_ = iPedestalConstraint;
    timeConstraint_     = iTimeConstraint;
    addPulseJitter_     = iAddPulseJitter;
    pulseJitter_        = iPulseJitter;
    timeMean_           = iTimeMean;
    timeSig_            = iTimeSig;
    pedMean_            = iPedMean;
    pedSig_             = iPedSig;
    noise_              = iNoise;
    timeShift_          = 100.;
    timeShift_ += 12.5;//we are trying to get BX 

    inverttimeSig_ = 1./timeSig_;
    inverttimeSig2_ = inverttimeSig_*inverttimeSig_;
    invertpedSig_ = 1./pedSig_;
    invertpedSig2_ = invertpedSig_*invertpedSig_;
  }

FitterFuncs::PulseShapeFunctor::~PulseShapeFunctor

(

)

Definition at line 102 of file PulseShapeFitOOTPileupCorrection.cc.

                                        {
  }

Member Function Documentation

double FitterFuncs::PulseShapeFunctor::doublePulseShapeFunc

(

const double *

x

)

Definition at line 156 of file PulseShapeFitOOTPileupCorrection.cc.

References EvalPulse().

Referenced by PulseShapeFitOOTPileupCorrection::resetPulseShapeTemplate(), and PulseShapeFitOOTPileupCorrection::setPulseShapeTemplate().

                                                                   {
      std::vector<double> pars(x, x+5);
      return EvalPulse(pars);
   }

double FitterFuncs::PulseShapeFunctor::EvalPulse

(

const std::vector< double > &

pars

)

Definition at line 105 of file PulseShapeFitOOTPileupCorrection.cc.

References addPulseJitter_, cntNANinfit, constexpr, funcHPDShape(), i, HcalConst::invertnsPerBx, invertpedSig2_, inverttimeSig2_, edm::isNotFinite(), j, HcalConst::maxSamples, pileupCalc::nbins, HcalConst::nsPerBX, pedestalConstraint_, pedMean_, psFit_erry, psFit_erry2, psFit_slew, psFit_y, pulseJitter_, cond::rpcobgas::time, timeConstraint_, timeMean_, and timeShift_.

Referenced by doublePulseShapeFunc(), singlePulseShapeFunc(), and triplePulseShapeFunc().

                                                                   {
      constexpr unsigned nbins = (unsigned) HcalConst::maxSamples;
      unsigned i =0, j=0;
      //Stop crashes
      for(i =0; i < pars.size(); ++i ) if( edm::isNotFinite(pars[i]) ){ ++ cntNANinfit; return 1e10; }
      
      //calculate chisquare
      double chisq  = 0;
      double delta2 =0;
      std::array<float,HcalConst::maxSamples> pulse_shape_sum;
      for(i=0; i < (pars.size()-1)/2; ++i ){
         int time = (pars[i*2]+timeShift_-timeMean_)*HcalConst::invertnsPerBx;
         //Interpolate the fit (Quickly)
         std::array<float,HcalConst::maxSamples> pulse_shape;
         funcHPDShape(pulse_shape, pars[i*2],pars[i*2+1],psFit_slew[time]);
         // add an uncertainty from the pulse (currently noise * pulse height =>Ecal uses full cov)
         if(addPulseJitter_) {
            for (j=0; j<nbins; ++j) {
               if( i==0 ) psFit_erry2[j]  = psFit_erry[j]*psFit_erry[j];
               psFit_erry2[j] += pulse_shape[j]*pulse_shape[j]*pulseJitter_*pulseJitter_;
            }
         }
         for(j=0; j<nbins; ++j){
            if( i==0 ) pulse_shape_sum[j] = pulse_shape[j] + pars.back();
            else pulse_shape_sum[j] += pulse_shape[j];
         }
      }
      for (i=0;i<nbins; ++i) {
        delta2 = (psFit_y[i]- pulse_shape_sum[i])*(psFit_y[i]- pulse_shape_sum[i])/psFit_erry2[i];
        chisq += delta2;
      }
      //Add the pedestal Constraint to chi2
      if(pedestalConstraint_) {
    chisq += invertpedSig2_*(pars.back() - pedMean_)*(pars.back()- pedMean_);
      }
      //Add the time Constraint to chi2
      if(timeConstraint_) {
    for(j=0; j< (pars.size()-1)/2; ++j ){
      int time = (pars[j*2]+timeShift_-timeMean_)/(double)HcalConst::nsPerBX;
      double time1 = -100.+time*HcalConst::nsPerBX;
      chisq += inverttimeSig2_*(pars[j*2] - timeMean_ - time1)*(pars[j*2] - timeMean_ - time1);
    }
      }
      return chisq;
   }

void FitterFuncs::PulseShapeFunctor::funcHPDShape ( std::array< float, HcalConst::maxSamples > &  ntmpbin, const double &  pulseTime, const double &  pulseHeight, const double &  slew  )

private

Definition at line 65 of file PulseShapeFitOOTPileupCorrection.cc.

References funct::abs(), acc25nsVec, accVarLenIdxMinusOneVec, accVarLenIdxZEROVec, cntNANinfit, constexpr, diff25nsItvlVec, diffVarItvlIdxMinusOneVec, diffVarItvlIdxZEROVec, i, HcalConst::iniTimeShift, edm::isNotFinite(), HcalConst::maxSamples, and HcalConst::nsPerBX.

Referenced by EvalPulse().

                                                                                                                                                           { 
    // pulse shape components over a range of time 0 ns to 255 ns in 1 ns steps
    constexpr int ns_per_bx = HcalConst::nsPerBX;
    constexpr int num_ns = HcalConst::nsPerBX*HcalConst::maxSamples;
    constexpr int num_bx = num_ns/ns_per_bx;
    // zeroing output binned pulse shape
    ntmpbin = { {0.0f,0.0f,0.0f,0.0f,0.0f,0.0f,0.0f,0.0f,0.0f,0.0f} };
    //Get the starting time
    int i_start         = ( -HcalConst::iniTimeShift - pulseTime - slew >0 ? 0 : (int)std::abs(-HcalConst::iniTimeShift-pulseTime-slew) + 1);
    double offset_start = i_start - HcalConst::iniTimeShift - pulseTime - slew; //-199-2*pars[0]-2.*slew (for pars[0] > 98.5) or just -98.5-pars[0]-slew;
    if( offset_start == 1.0 ){ offset_start = 0.; i_start-=1; } //Deal with boundary
    const int bin_start        = (int) offset_start; //bin off to integer
    const int bin_0_start      = ( offset_start < bin_start + 0.5 ? bin_start -1 : bin_start ); //Round it
    const int iTS_start        = i_start/ns_per_bx;         //Time Slice for time shift
    const int distTo25ns_start = HcalConst::nsPerBX - 1 - i_start%ns_per_bx;    //Delta ns 
    const double factor = offset_start - bin_0_start - 0.5; //Small correction?
    
    if( edm::isNotFinite(offset_start) ){ //Check for nan
      ++ cntNANinfit;
    }else{
      //Build the new pulse
      ntmpbin[iTS_start] = (bin_0_start == -1 ? // Initial bin (I'm assuming this is ok)
                  accVarLenIdxMinusOneVec[distTo25ns_start] + factor * diffVarItvlIdxMinusOneVec[distTo25ns_start]
                : accVarLenIdxZEROVec    [distTo25ns_start] + factor * diffVarItvlIdxZEROVec    [distTo25ns_start]);
      //Fill the rest of the bins
      for(int iTS = iTS_start+1; iTS < num_bx; ++iTS){
    int bin_idx = distTo25ns_start + 1 + (iTS-iTS_start-1)*ns_per_bx + bin_0_start;
    ntmpbin[iTS] = acc25nsVec[bin_idx] + factor * diff25nsItvlVec[bin_idx];
      }
    }
    //Scale the pulse 
    for(int i=0; i < num_bx; ++i) {
      ntmpbin[i]     *= pulseHeight;
    }
    return;
  }

int FitterFuncs::PulseShapeFunctor::getcntNANinfit ( )

inline

Definition at line 43 of file PulseShapeFitOOTPileupCorrection.h.

References cntNANinfit.

{ return cntNANinfit; }

void FitterFuncs::PulseShapeFunctor::setDefaultcntNANinfit ( )

inline

Definition at line 42 of file PulseShapeFitOOTPileupCorrection.h.

References cntNANinfit.

{ cntNANinfit =0; }

void FitterFuncs::PulseShapeFunctor::setpsFiterry ( double *  erry )

inline

Definition at line 47 of file PulseShapeFitOOTPileupCorrection.h.

References i, HcalConst::maxSamples, and psFit_erry.

{ for(int i=0; i<HcalConst::maxSamples; ++i) psFit_erry  [i] = erry [i]; }

void FitterFuncs::PulseShapeFunctor::setpsFiterry2 ( double *  erry2 )

inline

Definition at line 48 of file PulseShapeFitOOTPileupCorrection.h.

References i, HcalConst::maxSamples, and psFit_erry2.

{ for(int i=0; i<HcalConst::maxSamples; ++i) psFit_erry2 [i] = erry2[i]; }

void FitterFuncs::PulseShapeFunctor::setpsFitslew ( double *  slew )

inline

Definition at line 49 of file PulseShapeFitOOTPileupCorrection.h.

References i, HcalConst::maxSamples, and psFit_slew.

{ for(int i=0; i<HcalConst::maxSamples; ++i) {psFit_slew [i] = slew [i]; } }

void FitterFuncs::PulseShapeFunctor::setpsFitx ( double *  x )

inline

Definition at line 45 of file PulseShapeFitOOTPileupCorrection.h.

References i, HcalConst::maxSamples, and psFit_x.

{ for(int i=0; i<HcalConst::maxSamples; ++i) psFit_x[i] = x[i]; }

void FitterFuncs::PulseShapeFunctor::setpsFity ( double *  y )

inline

Definition at line 46 of file PulseShapeFitOOTPileupCorrection.h.

References i, HcalConst::maxSamples, and psFit_y.

{ for(int i=0; i<HcalConst::maxSamples; ++i) psFit_y[i] = y[i]; }

double FitterFuncs::PulseShapeFunctor::sigma

(

double

ifC

)

Definition at line 166 of file PulseShapeFitOOTPileupCorrection.cc.

References alignCSCRings::e.

                                             { 
      if(ifC < 75) return (0.577 + 0.0686*ifC)/3.; 
      return (2.75  + 0.0373*ifC + 3e-6*ifC*ifC)/3.; 
   }

double FitterFuncs::PulseShapeFunctor::singlePulseShapeFunc

(

const double *

x

)

Definition at line 151 of file PulseShapeFitOOTPileupCorrection.cc.

References EvalPulse().

Referenced by PulseShapeFitOOTPileupCorrection::resetPulseShapeTemplate(), and PulseShapeFitOOTPileupCorrection::setPulseShapeTemplate().

                                                                   {
      std::vector<double> pars(x, x+3);
      return EvalPulse(pars);
   }

double FitterFuncs::PulseShapeFunctor::triplePulseShapeFunc

(

const double *

x

)

Definition at line 161 of file PulseShapeFitOOTPileupCorrection.cc.

References EvalPulse().

Referenced by PulseShapeFitOOTPileupCorrection::resetPulseShapeTemplate(), and PulseShapeFitOOTPileupCorrection::setPulseShapeTemplate().

                                                                   {
      std::vector<double> pars(x, x+7);
      return EvalPulse(pars);
   }

Member Data Documentation

std::vector<float> FitterFuncs::PulseShapeFunctor::acc25nsVec

private

Definition at line 59 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by funcHPDShape(), and PulseShapeFunctor().

std::vector<float> FitterFuncs::PulseShapeFunctor::accVarLenIdxMinusOneVec

private

Definition at line 61 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by funcHPDShape(), and PulseShapeFunctor().

std::vector<float> FitterFuncs::PulseShapeFunctor::accVarLenIdxZEROVec

private

Definition at line 60 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by funcHPDShape(), and PulseShapeFunctor().

bool FitterFuncs::PulseShapeFunctor::addPulseJitter_

private

Definition at line 67 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

int FitterFuncs::PulseShapeFunctor::cntNANinfit

private

Definition at line 58 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), funcHPDShape(), getcntNANinfit(), and setDefaultcntNANinfit().

std::vector<float> FitterFuncs::PulseShapeFunctor::diff25nsItvlVec

private

Definition at line 59 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by funcHPDShape(), and PulseShapeFunctor().

std::vector<float> FitterFuncs::PulseShapeFunctor::diffVarItvlIdxMinusOneVec

private

Definition at line 61 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by funcHPDShape(), and PulseShapeFunctor().

std::vector<float> FitterFuncs::PulseShapeFunctor::diffVarItvlIdxZEROVec

private

Definition at line 60 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by funcHPDShape(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::invertpedSig2_

private

Definition at line 78 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::invertpedSig_

private

Definition at line 78 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::inverttimeSig2_

private

Definition at line 77 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::inverttimeSig_

private

Definition at line 77 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::noise_

private

Definition at line 74 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor().

bool FitterFuncs::PulseShapeFunctor::pedestalConstraint_

private

Definition at line 65 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::pedMean_

private

Definition at line 72 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::pedSig_

private

Definition at line 73 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::psFit_erry[HcalConst::maxSamples]

private

Definition at line 63 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), PulseShapeFunctor(), and setpsFiterry().

double FitterFuncs::PulseShapeFunctor::psFit_erry2[HcalConst::maxSamples]

private

Definition at line 63 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), PulseShapeFunctor(), and setpsFiterry2().

double FitterFuncs::PulseShapeFunctor::psFit_slew[HcalConst::maxSamples]

private

Definition at line 63 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), PulseShapeFunctor(), and setpsFitslew().

double FitterFuncs::PulseShapeFunctor::psFit_x[HcalConst::maxSamples]

private

Definition at line 63 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor(), and setpsFitx().

double FitterFuncs::PulseShapeFunctor::psFit_y[HcalConst::maxSamples]

private

Definition at line 63 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), PulseShapeFunctor(), and setpsFity().

std::array<float,HcalConst::maxPSshapeBin> FitterFuncs::PulseShapeFunctor::pulse_hist

private

Definition at line 56 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::pulseJitter_

private

Definition at line 69 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

bool FitterFuncs::PulseShapeFunctor::timeConstraint_

private

Definition at line 66 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::timeMean_

private

Definition at line 70 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::timeShift_

private

Definition at line 75 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by EvalPulse(), and PulseShapeFunctor().

double FitterFuncs::PulseShapeFunctor::timeSig_

private

Definition at line 71 of file PulseShapeFitOOTPileupCorrection.h.

Referenced by PulseShapeFunctor().

bool FitterFuncs::PulseShapeFunctor::unConstrainedFit_

private

Definition at line 68 of file PulseShapeFitOOTPileupCorrection.h.