CorrFit1DCFCumac.h

/*
 * CorrFIt1DCFCumac.h
 *
 *  Created on: 13 mar 2018
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
#ifndef HALCORRFIT1DCFCUMAC_H_
#define HALCORRFIT1DCFCUMAC_H_
 
#include "CorrFit1DCF.h"
#include "Object.h"
 
namespace Hal {
  class CorrFit1DCFCumac : public CorrFit1DCF {
  protected:
    const Double_t fPis;
    const Double_t fHc;
    mutable Double_t fZ;
    virtual Double_t Ff1(Double_t x) const;
    virtual Double_t Integr() const;
    Double_t Sign(Double_t a, Double_t b) const;
    Double_t
    Simps(Double_t a1, Double_t b1, Double_t h1, Double_t reps1, Double_t aeps1, Double_t x, Double_t aih, Double_t aiabs) const;
    Double_t F1(Double_t Z) const;
    Double_t F2(Double_t Z) const;
    CorrFit1DCFCumac(Int_t params);
 
  public:
    CorrFit1DCFCumac();
    virtual Double_t Get(Double_t q, Double_t r);
    virtual ~CorrFit1DCFCumac();
    ClassDef(CorrFit1DCFCumac, 1)
  };
 
  class CorrFIt1DCFCumacLamLam : public CorrFit1DCFCumac {
  protected:
    virtual Double_t CalculateCF(const Double_t* x, const Double_t* params) const;
 
  public:
    CorrFIt1DCFCumacLamLam();
    Int_t ScatteringLengthID() const { return 3; };
    Int_t EffectiveRadiusID() const { return 4; };
    Int_t LambdaPolarizationID() const { return 5; };
    Int_t ResidualAmplitudeID() const { return 6; };
    Int_t ResidualGaussWidhtID() const { return 7; };
    virtual ~CorrFIt1DCFCumacLamLam();
    ClassDef(CorrFIt1DCFCumacLamLam, 1)
  };
 
  class CorrFit1DCFCumacPLam : public CorrFit1DCFCumac {
  protected:
    virtual Double_t CalculateCF(const Double_t* x, const Double_t* params) const;
 
  public:
    CorrFit1DCFCumacPLam();
    Int_t SingletScatteringLengthID() const { return 3; }
    Int_t SingletEffectiveRadiusID() const { return 4; }
    Int_t TripletScatteringLenghtID() const { return 5; };
    Int_t TripletEffectiveRadiusID() const { return 6; };
    Int_t LambdaPolarizationID() const { return 7; };
    virtual ~CorrFit1DCFCumacPLam();
    ClassDef(CorrFit1DCFCumacPLam, 1)
  };
 
  class CorrFit1DCFCumacDLam : public CorrFit1DCFCumac {
  protected:
    Double_t F(Double_t d, Double_t r0) const;
    virtual Double_t CalculateCF(const Double_t* x, const Double_t* params) const;
 
  public:
    CorrFit1DCFCumacDLam();
    Int_t DoubletScatteringLengthID() const { return 3; }
    Int_t DoubletEffectiveRadiusID() const { return 4; }
    Int_t QuartetScatteringLenghtID() const { return 5; };
    Int_t QuartetEffectiveRadiusID() const { return 6; };
    virtual ~CorrFit1DCFCumacDLam() {};
    ClassDef(CorrFit1DCFCumacDLam, 1)
  };
 
  class CorrFit1DCFCumacK0K0 : public CorrFit1DCFCumac {
  protected:
    virtual Double_t CalculateCF(const Double_t* x, const Double_t* params) const;
 
  public:
    CorrFit1DCFCumacK0K0();
    Int_t Mf0ID() const { return 3; };
    Int_t Gamma_f0KKID() const { return 4; }
    Int_t Ma0ID() const { return 5; };
    Int_t Gamma_f0pipiID() const { return 6; };
    Int_t Gamma_a0PiEtaID() const { return 7; }
    Int_t Gamma_a0KKID() const { return 8; }
    Int_t AssymetryID() const { return 9; }
 
    void SetDefParams(Int_t opt);
    void UseLCMS() { fKinematics = Femto::EKinematics::kLCMS; }
    virtual ~CorrFit1DCFCumacK0K0() {};
    ClassDef(CorrFit1DCFCumacK0K0, 1)
  };
  class CorrFit1DCFCumacK0Kch : public CorrFit1DCFCumac {
  protected:
    Hal::Femto::CorrFitGammaCalc fGammaCalc;
 
    virtual Double_t CalculateCF(const Double_t* x, const Double_t* params) const;
 
  public:
    CorrFit1DCFCumacK0Kch();
    Int_t Ma0ID() const { return 0; };
    Int_t Gamma_a0KKID() const { return 1; }
    Int_t Gamma_a0PiEtaID() const { return 2; }
    void SetDefParams(Int_t opt);
    void UseLCMS() { fKinematics = Femto::EKinematics::kLCMS; }
    virtual ~CorrFit1DCFCumacK0Kch() {};
    ClassDef(CorrFit1DCFCumacK0Kch, 1)
  };
}  // namespace Hal
#endif /* HALCORRFIT1DCFCUMAC_H_ */