PndSdsFE.h

Go to the documentation of this file.

/*
 * SdsFEAmpModel.h
 *
 *  Created on: Apr 4, 2011
 *      Author: esch
 *
 *      Class to calculate ToT, TimeWalk, TimeStamp and Time till the signal is back to baseline from a TF1 amplifier model.
 *
 *
 */

#ifndef PndSdsFE_H_
#define PndSdsFE_H_
#include "TObject.h"
#include "PndSdsFEAmpModelSimple.h"
#include "TMath.h"
#include "TF1.h"
#include "TRandom2.h"
#include "Math/Interpolator.h"
#include "Math/InterpolationTypes.h"
#include "PndSdsTotChargeConversion.h"
#include "TVectorT.h"

class PndSdsFE: public TObject {
public:
        PndSdsFE(): fFrontEndModel(0), fNumberOfSupportPoints(0), fNumberOfMaxElectons(0), fFunction(0), fMaximumAmplitude(0),fThreshold(0), fBaselineEpsilon(1), fToF(0), fEventTime(0),       fTimeOffSet(0), fTimeStep(0),fFunctionRange(0), fRand(0), i(0), stepsize(0), fCharge_list(0), fTot_list(0), inter(0) 
        {
                CalcTimeOffSet();
        };
  
  PndSdsFE(PndSdsFE& other):TObject(other),
      fFrontEndModel(other.fFrontEndModel),
      fNumberOfSupportPoints(other.fNumberOfSupportPoints),
      fNumberOfMaxElectons(other.fNumberOfMaxElectons),
      fFunction(other.fFunction),
      fMaximumAmplitude(other.fMaximumAmplitude),
      fThreshold(other.fThreshold),
      fBaselineEpsilon(other.fBaselineEpsilon),
      fToF(other.fToF),
      fEventTime(other.fEventTime),
          fTimeOffSet(other.fTimeOffSet),
      fTimeStep(other.fTimeStep),
      fFunctionRange(other.fFunctionRange),
      fRand(other.fRand),
      i(other.i),
      stepsize(other.stepsize),
      fCharge_list(other.fCharge_list),
      fTot_list(other.fTot_list),
      inter(other.inter)
  {};
//      PndSdsFE(double charingtime, double constcurrent, double threshold, double frequency, int verbose);
        virtual ~PndSdsFE();

        virtual double GetTotFromCharge(Double_t charge);                       // Calculates the ToT value from the charge
        virtual double GetChargeFromTot(double tot);                            // Calculates the charge from the tot value
        virtual double GetTimeWalkFromCharge(double charge);            // Calculates the TimeWalk from the charge
        virtual double GetTimeWalkFromTot(double tot);                          // Calculates the TimeWalk from the TOT value
        virtual double GetTimeBackToBaseline(double charge);            // Calculates the time from start of the signal till the signal is fBaselineEpsilon close to baseline
        virtual double GetTimeStamp(double eventtime, double tof, double charge);  // Calculates the TimeStamp of the signal
        virtual double GetTimeStep(){return fTimeStep;}
        virtual void SetParameter(TString parName, double parValue);
        virtual void SetThreshold(double threshold){fThreshold = threshold;}

  PndSdsFE& operator=(PndSdsFE& other)
  {
    if (this != &other) // protect against invalid self-assignment
    {   
      fFrontEndModel=other.fFrontEndModel;
      fNumberOfSupportPoints=other.fNumberOfSupportPoints;
      fNumberOfMaxElectons=other.fNumberOfMaxElectons;
      fFunction=other.fFunction;
      fMaximumAmplitude=other.fMaximumAmplitude;
      fThreshold=other.fThreshold;
      fBaselineEpsilon=other.fBaselineEpsilon;
      fToF=other.fToF;
      fEventTime=other.fEventTime;
          fTimeOffSet=other.fTimeOffSet;
      fTimeStep=other.fTimeStep;
      fFunctionRange=other.fFunctionRange;
      fRand=other.fRand;
      i=other.i;
      stepsize=other.stepsize;
      fCharge_list=other.fCharge_list;
      fTot_list=other.fTot_list;
      inter=other.inter;
    }
    return *this;
  };
protected:

  double DigitizeTime(double time);
        void CalcTimeOffSet();
  void CreateInterpolatorList();  // Creates the Interpolatorlist for GetTotFromCharge calculation
        void SaveInterpolatorList(std::vector<double> charge, std::vector<double> tot);
        void LoadInterpolatorList();
        void GetInterpolatorList();


        PndSdsFEAmpModelSimple *fFrontEndModel;

        int fNumberOfSupportPoints;
        int fNumberOfMaxElectons ;

        TF1 *fFunction;
        double fMaximumAmplitude;      // Saves position of the maximum amplitude of the signal
        double fThreshold;                              // Threshold for Signal from Parameter database
        double fBaselineEpsilon;                //
        double fToF;                                    // Time of Flight from MC
        double fEventTime;                              // EventTime from MC
        double fTimeOffSet;                             // Random offset
        double fTimeStep;
        double fFunctionRange;                  //Range of the function. max 100 fC => 624 150.9 e => max.tot 10503 ns (with thr.= 1100 e)
        TRandom2 fRand;
        int i;
        double stepsize;

        std::vector<double> fCharge_list;
        std::vector<double> fTot_list;
        ROOT::Math::Interpolator *inter;

private:




        ClassDef(PndSdsFE,1);
};

#endif /* PndSdsFE_H_ */