FemtoConst.h

/*
 * FemtoConst.h
 *
 *  Created on: 26-11-2013
 *      Author: Daniel Wielanek
 *      E-mail: daniel.wielanek@gmail.com
 *      Warsaw University of Technology, Faculty of Physics
 */
 
#ifndef HALFEMTOCONST_H_
#define HALFEMTOCONST_H_
 
#include "Array.h"
#include "Object.h"
#include "TString.h"
class TVector3;
namespace Hal {
  namespace Femto {
 
    inline Double_t FmToGeV() { return 5.06842372; }
    inline Double_t FmToGeV(Double_t fm) { return fm * FmToGeV(); }
 
    enum class EKinematics {
      kPRF,   
      kLCMS,  
      kSH_LCMS,
      kSH_PRF,
      kPHIETA,
      kPRFL
    };
    enum class ECFType { kOneDim, kThreeDim, kSpherical, kPhiEta, kUnkown };
    enum class EPairType {
      kUnknown           = 100,
      kPionPlusPionPlus  = 101,
      kPionPlusPionMinus = 102,
      kKaonPlusKaonPlus  = 103,
      kKaonPlusKaonMinus = 4,
 
      kProtonProton      = 105,
      kProtonAntiproton  = 6,
      kPionPlusKaonPlus  = 107,
      kPionPlusKaonMinus = 108,
 
      kPionPlusProton     = 109,
      kPionPlusAntiproton = 110,
      kKaonPlusProton     = 111,
      kKaonPlusAntiproton = 112,
 
      kProtonLambda        = 113,
      kLambdaLambda        = 114,
      kKaonZeroKaonZero    = 15,
      kKaonZeroKaonZeroBar = 16,
 
      kNeutronNeutron   = 117,
      kNeutronProton    = 118,
      kPionZeroPionZero = 19,
      kNeutronLambda    = 20,
 
      // TPI weights only
      kProtonSigmaPlus     = 201,
      kProtonAntiSigmaPlus = 202,
      kProtonAntiLambda    = 203,
      kSigmaPlusSigmaPlus  = 204,
 
      kSigmaPlusAntiSigmaPlus = 205,
      kProtonXiZero           = 206,
      kNeutronXiMinus         = 207,
      kProtonXIMInus          = 208,
 
      kNeutronXiZero      = 209,
      kProtonSigmaZero    = 210,
      kSigmaZeroSigmaZero = 211,
      kLambdaSigmaZero    = 212,
 
      kLambdaAntiLambda = 213
 
    };
 
 
  }  // namespace Femto
 
  class FemtoPair;
  class DividedHisto1D;
  class XMLNode;
  class FemtoWeightGenerator;
  namespace Femto {
    Bool_t IsPairIdentical(EPairType pt);
    std::pair<Int_t, Int_t> PairTypeToPid(EPairType pt);
    EPairType PidToPairType(Int_t pid1, Int_t pid2);
    TString PairTypeToString(EPairType type);
    EKinematics LabelToKinematics(TString label);
    EKinematics CodeLabelToKinematics(TString label);
    FemtoPair* MakePair(EKinematics kin, Bool_t use_fake);
    DividedHisto1D* GetHistoFromXML(XMLNode* nod);
    FemtoWeightGenerator* GetWeightGeneratorFromXLM(XMLNode* nod);
    TString KinematicsToLabel(EKinematics kin);
    TString KinematicsToAxisLabel(EKinematics kin, Int_t ax, Int_t ndim);
    void FillRandomPair(FemtoPair& p, Int_t pid1, Int_t pid2, Double_t sigmaq = 0.1, Double_t sigmar = 1);
    void FillRandomKinematics(FemtoPair& p, const TVector3& sum, const TVector3& diff, EKinematics kin);
    ECFType GetCFType(TObject* obj);
    class CorrFitGammaCalc : public Object {
      Double_t fA2    = {0};
      Double_t fB2    = {0};
      Double_t fConsA = {0};
      Double_t fConsB = {0};
 
    public:
      CorrFitGammaCalc(Double_t a = 0, Double_t b = 0, Double_t c = 0, Double_t d = 0) { ReInit(a, b, c, d); };
      void ReInit(Double_t a, Double_t b, Double_t c, Double_t d);
      Double_t Calculate(Double_t kstar) const;
      virtual ~CorrFitGammaCalc() {};
      ClassDef(CorrFitGammaCalc, 0)
    };
  }  // namespace Femto
}  // namespace Hal
 
#endif /* HALFEMTOCONST_H_ */