PndMvdCalcTot.h

Go to the documentation of this file.

#ifndef PNDMVDCALCTOT_H
#define PNDMVDCALCTOT_H

#include "TObject.h"
#include "TString.h"
#include "TRandom2.h"
#include <iostream>
#include <iomanip>


class PndMvdCalcTot
{
  public :
      PndMvdCalcTot(){
          ftr = 100.;
          fa = 60.;
          fthreshold = 3000;
          fclockfrequency = 50.;
          ftimestep = 1. / fclockfrequency * 1000.;
          SetStartOffset();
          std::cout<<"-I- <PndMvdCalcTot::PndMvdCalcTot()>: default constructor used!"<<std::endl;
      };

      PndMvdCalcTot(Double_t tr, Double_t a, Double_t threshold, Double_t clockfrequency){
          ftr = tr;
          fa = a;
          fthreshold = threshold;
          fclockfrequency = clockfrequency;
          ftimestep = 1. / fclockfrequency * 1000.;
          SetStartOffset();
      };

          ~PndMvdCalcTot(){};
                
            void SetStartOffset();

                Double_t GetTot(Double_t fcharge);
                void SetParameter(Double_t tr, Double_t a, Double_t threshold);
                Double_t GetPileUpTime(Double_t fcharge);               //returns the time the capacitor is loaded and therefore the time the hit MVD is blind for other events
                
  private : 
                Double_t fthreshold;            //in electrons
        Double_t ftr;                           //raising time [ns]
        Double_t fa;                            //falling ration [e/ns]
        Double_t Q, Qt;                         //max charge, threshold charge
                Double_t t1e, t2e;                      //exact time point when signal is over threshold/under threshold again
        Double_t fclockfrequency;       //the frequency the readout chip works with
        Double_t ftimestep;                     //one time step of the clock
        Double_t fstarttime;            //absolute point when TOT begins
        Double_t fstoptime;                     //absolute point when TOT ends
        Double_t ftimeoffset;           //difference between rising clock edge and when detector gets hit
        
        Double_t GetTotWC();            //assumes a clock which results in a quantized TOT
        
        TRandom2 fRand;
};

#endif /*MVDCALCTOT_H*/