FemtoWeightGeneratorLednicky.h

/*
 * FemtoWeightGeneratorLednicky.h
 *
 *   derived from ALICE AliFemtoWeightGeneratorLednicky
 *  Created on: 21-08-2013
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
 
#ifndef HALFEMTOWEIGHTGENERATORLEDNICKY_H_
#define HALFEMTOWEIGHTGENERATORLEDNICKY_H_
 
#include "FemtoWeightGenerator.h"
 
namespace Hal {
  class FemtoWeightGeneratorLednicky : public FemtoWeightGenerator {
    friend class FemtoWeightGeneratorResidual;
 
  protected:
    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;
    int fNS;
    int* fNumProcessPair;  
    // WTF? Streamer?
    int fNumbNonId;
 
    int fKpKmModel;  // ij (i=1..4, j=1..4; see AliFemtoFsiWeightLednicky.F)
    int fPhi_OffOn;  // 0->Phi Off,1->Phi On
    int fNS_4;       // set NS is equal to 4
 
    // Interface to the fortran functions
    void FsiInit();
    void FsiSetLL();
    void FsiNucl();
    bool SetPid(const int aPid1, const int aPid2);
    void FsiSetKpKmModelType();
    virtual void SetPairTypeFromPairAndConfigureFSI(FemtoPair* pair);
    void PairIDToLL(Femto::EPairType pair_pid, Int_t& LL, Int_t& tCharge) const;
    Femto::EPairType LLToPairID(Int_t ll) const;
    inline void SetFlags(Int_t qs, Int_t s, Int_t c) {
      fIqs = qs;
      fIch = c;
      fIsi = s;
    };
    void RawCalc(FemtoPair* pair);
 
  public:
    FemtoWeightGeneratorLednicky();  // call SetDefaultCalcPar
    FemtoWeightGeneratorLednicky(const FemtoWeightGeneratorLednicky& aWeight);  // call SetDefaultCalcPar
    // --- Destructor : nothing to explicitly delete
    FemtoWeightGeneratorLednicky& operator=(const FemtoWeightGeneratorLednicky& aWeight);
    virtual ~FemtoWeightGeneratorLednicky();
 
    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 FemtoWeightGeneratorLednicky* MakeCopy() const { return new FemtoWeightGeneratorLednicky(*this); };
 
    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);
    void SetNS(int nNs) { fNS = nNs; };
 
    void FsiSetKpKmModelType(const int aModelType, const int aPhi_OffOn);
    virtual void Print(Option_t* option = "") const;
    virtual Package* Report() const;
    ClassDef(FemtoWeightGeneratorLednicky, 1)
  };
}  // namespace Hal
#endif /* HALFEMTOWEIGHTGENERATORLEDNICKY_H_ */