FemtoWeightGeneratorLednicky2.h

/*
 * FemtoWeightGeneratorLednicky2.h
 *
 *  Created on: 30 kwi 2018
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
#ifndef HALFEMTOWEIGHTGENERATORLEDNICKY2_H_
#define HALFEMTOWEIGHTGENERATORLEDNICKY2_H_
 
 
#include "FemtoWeightGenerator.h"
 
#include "FemtoConst.h"
namespace Hal {
  class FemtoFsiParsed;
 
  class FemtoWeightGeneratorLednicky2 : public FemtoWeightGenerator {
  protected:
    FemtoFsiParsed* fParser;
    double fWei;
    double fWein;
    double fWeif;
    double fWeightDen;
    int fItest;  //
 
    // int mNs;
    int fIch;
    int fIqs;
    int fIsi;
    int fI3c;
    double fNuclMass;
    double fNuclCharge;
    bool fSphereApp;
    bool fT0App;
 
    // Pair identification
    int fLL;
    short fNuclChargeSign;
    bool fSwap;
    int const fLLMax;
    char** fLLName;  
    int* fNumProcessPair;  
    // WTF? Streamer?
    int fNumbNonId;
 
    // Interface to the fortran functions
    void FsiInit();
    void FsiSetLL();
    void FsiNucl();
    bool SetPid(const int aPid1, const int aPid2);
    virtual void SetPairTypeFromPairAndConfigureFSI(FemtoPair* pair);
 
  public:
    FemtoWeightGeneratorLednicky2();  // call SetDefaultCalcPar
    FemtoWeightGeneratorLednicky2(const FemtoWeightGeneratorLednicky2& aWeight);  // call SetDefaultCalcPar
    // --- Destructor : nothing to explicitly delete
    FemtoWeightGeneratorLednicky2& operator=(const FemtoWeightGeneratorLednicky2& aWeight);
    virtual FemtoWeightGeneratorLednicky2* MakeClone() const { return new FemtoWeightGeneratorLednicky2(*this); };
    virtual ~FemtoWeightGeneratorLednicky2();
 
    virtual Double_t GenerateWeight(FemtoPair* pair);
 
    virtual void SetPairType(Femto::EPairType aPairType);
    virtual void SetPairTypeFromPair(FemtoPair* pair);
    virtual Femto::EPairType GetPairType() const;
 
    virtual Double_t GetKStar() const;
    virtual Double_t GetKStarOut() const;
    virtual Double_t GetKStarSide() const;
    virtual Double_t GetKStarLong() const;
    virtual Double_t GetRStar() const;
    virtual Double_t GetRStarOut() const;
    virtual Double_t GetRStarSide() const;
    virtual Double_t GetRStarLong() const;
 
    virtual FemtoWeightGenerator* Clone(const char* newname = "") const;
 
    // --- Setting
    Bool_t Init(Int_t taks_id, FemtoPair* pair);
    virtual Bool_t IsPairSupported(Femto::EPairType type) const;
    // >>> Calculation mode
    void SetDefaultCalcPar();  // Default is CoulOn, QuantumOn, StrongOn,
                               // 3BodyOff, Square, T0ApproxOff
    void SetCoulOn();
    void SetCoulOff();
    void SetQuantumOn();
    void SetQuantumOff();
    void SetStrongOn();
    void SetStrongOff();
    void Set3BodyOn();
    void Set3BodyOff();
    void SetSphere();  // use Spherical wave approximation
    void SetSquare();
    void SetT0ApproxOff();
    void SetT0ApproxOn();
    void PrintLambdas() { ; }
    void SetNuclCharge(const double aNuclCharge);  // for 3-body calculation
    void SetNuclMass(const double aNuclMass);
    ClassDef(FemtoWeightGeneratorLednicky2, 1)
  };
}  // namespace Hal
#endif /* HALFEMTOWEIGHTGENERATORLEDNICKY2_H_ */