StFcsDbPulse.h

/*
David Kapukchyan
April 21, 2021
First instance of pulse db class to contain all the basic variables for the pulse shape from DEP boards as well as the functional forms
 */

#ifndef STFCSDBPULSE_H
#define STFCSDBPULSE_H

#include <iostream>

#include "TDataSet.h"
#include "TMath.h"
#include "TGraphAsymmErrors.h"
#include "TF1.h"
#include "St_base/StMessMgr.h"

class StFcsDbPulse : public TDataSet {

 public:
  StFcsDbPulse(const char* name = "fcsPulse"); 
  virtual ~StFcsDbPulse();                     
  int Init();                                  

  //Basic contants
  static double sqrtpi(){return sqrt(TMath::Pi());}        
  static double sqrt2pi(){return sqrt(2.0*TMath::Pi());}   
  
  //Constants for simulating pulses
  void setTBPerRC(double v){mTBPerRC = v;}         
  double TBPerRC()const{return mTBPerRC;}               
  double nsecPerTB()const{return 107.0/TBPerRC();}      
  double BeamLengthSig()const{return 10.0/nsecPerTB();} 

  void setTail(int tail);

  void setGSigma(double v){mGSigma = v;} 
  void setA1(double v){mA1 = v;}         
  void setA2(double v){mA2 = v;}         
  void setXoff1(double v){mXoff1 = v;}   
  void setXoff2(double v){mXoff2 = v;}   
  void setTau1(double v){mTau1 = v;}     
  void setTau2(double v){mTau2 = v;}     
  void setP1(double v){mP1 = v;}         
  void setP2(double v){mP2 = v;}         

  //Variables related to the tail function
  double GSigma()const{return mGSigma;}  
  double A1()const{return mA1;}          
  double A2()const{return mA2;}          
  double Xoff1()const{return mXoff1;}    
  double Xoff2()const{return mXoff2;}    
  double Tau1()const{return mTau1;}      
  double Tau2()const{return mTau2;}      
  double P1()const{return mP1;}          
  double P2()const{return mP2;}          

  static void setTGraphAsymmErrors(TGraphAsymmErrors* gae, const int &i, const double &adc, double Yerr, double YerrSat);

  double pulseShape(double* x, double* p);
  
  double multiPulseShape(double* x, double* p);
  
  TF1* createPulse(double xlow=0, double xhigh=1,int npars=5 );

  static int GenericPadPos(int value, int Nvals, int PadNums );
  
  static int PadNum4x4(int det, int col, int row);

  static Int_t getYMinMax(TGraphAsymmErrors* gae, Double_t &Ymin, Double_t &Ymax, Double_t xmin=-5, Double_t xmax=2000);

  virtual void Print(Option_t* opt = "") const; 

 protected:
  double mTBPerRC;  
  double mGSigma;   
  //Pulse shape constants
  double mA1;       
  double mA2;       
  double mXoff1;    
  double mXoff2;    
  double mTau1;     
  double mTau2;     
  double mP1;       
  double mP2;       

 private:
  int mTail;  
  static const int mAdcSaturation = 4000; 

  ClassDef(StFcsDbPulse, 1);
};

#endif