CorrFit3DCFBowlerSinyukov.h

/*
 * CorrFit3DCFBowlerSinyukov.h
 *
 *  Created on: 22 sty 2018
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
#ifndef HALCORRFIT3DCFBOWLERSINYUKOV_H_
#define HALCORRFIT3DCFBOWLERSINYUKOV_H_
 
#include "CorrFit3DCF.h"
#include "Splines.h"
 
#include "TF1.h"
 
namespace Hal {
  class Femto3DCF;
  class CorrFit3DCFBowlerSinyukov : public CorrFit3DCF {
  protected:
    Spline2D* fCFs;
    Spline3D* fQinvMap;
    Bool_t fTStarZero;
    mutable Double_t fRinv;
    Double_t fAverageBeta;
    Double_t fGamma;
    Double_t fQSMode;
    const Double_t fSqrt3;
    Femto::EKinematics fFrame;
    //  Double_t Integral(Double_t *x, Double_t *params)const;
    virtual void Fit(TObject* histo);
    virtual void FitDummy(TObject* histo);
    virtual Double_t GetRoutPRF(Double_t r_out, Double_t tau) const;
    inline Double_t Gaus(Double_t r, Double_t R) const { return TMath::Exp(-r * r / (4.0 * R * R)) * r * r; }
    virtual Double_t CalculateCF(const Double_t* x, const Double_t* params) const;
    virtual void RecalculateFunction() const;
    CorrFit3DCFBowlerSinyukov(Int_t extra_params);
 
  public:
    CorrFit3DCFBowlerSinyukov();
    void AsFermions() { fQSMode = -0.5; };
    void AsBosons() { fQSMode = 1; }
    void AsNonIdentical() { fQSMode = 0.0; };
    void LoadMap(TH2D* h, Option_t* opt = "");
    void SetAveragePairBeta(Double_t beta);
    void SetAveragePairBeta(Double_t m1, Double_t m2, Double_t avkt);
    void SetTStarZero() { fTStarZero = kTRUE; };
    Spline2D* GetKqMap() const { return fCFs; };
    Femto3DCF* UnBowlerize() const;
    static Int_t Tau() { return 5; };
    virtual ~CorrFit3DCFBowlerSinyukov();
    ClassDef(CorrFit3DCFBowlerSinyukov, 1)
  };
}  // namespace Hal
#endif /* HALCORRFIT3DCFBOWLERSINYUKOV_H_ */