xs.h

#ifndef CRXS__XS_H
#define CRXS__XS_H

#include "linAlg_tools.h"

namespace CRXS {
    
    enum parametrization{
        KORSMEIER_I    =  1,
        KORSMEIER_II   =  2,
        WINKLER        =  3,
        DI_MAURO_I     =  4,
        DI_MAURO_II    =  5,
        ANDERSON       =  6,
        WINKLER_SELF   =  7,
        DI_MAURO_SELF  =  8,
        APPROX_1_OVER_T=  9,
        WINKLER_II     = 10,
        KORSMEIER_III  = 11
    };
    enum product{
        P_BAR    = 1,
        D_BAR    = 2,
        HE3_BAR  = 3,
        HE4_BAR  = 4
    };
    
    enum coalescence{
        FIXED_P0                   = 1,
        ENERGY_DEP__VAN_DOETINCHEM = 2,
        PT_DEP                     = 3
    };
    
    class XS{
        
    public:
        
        

        static bool convert_LAB_to_CM( const double T_p_LAB, const double T_prod_LAB, const double eta_LAB, double &s, double &E_prod, double &pT_prod, double &x_F, int product=P_BAR );
        

        static double inv_AA_pbar_CM( double s, double xF, double pT_pbar, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization);


        static double inv_AA_pbar_LAB( double Tn_proj_LAB, double T_pbar_LAB, double eta_LAB, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization );
        
        

        static double dE_AA_pbar_LAB( double Tn_proj_LAB, double T_pbar_LAB, int A_projectile=1, int N_projectile=0, int A_target=1, int N_target=0, int parametrization=KORSMEIER_II);
        static double integrand__dE_AA_pbar_LAB (double eta_LAB, void* parameters  );
        
        
        

        static double dE_AA_pbar_LAB_incNbarAndHyperon(double Tn_proj_LAB, double T_pbar_LAB, int A_projectile=1, int N_projectile=0, int A_target=1, int N_target=0, int parametrization=KORSMEIER_II);
        
        
        

        static double inv_AA_p_CM( double s, double xF, double pT_p, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization);
        

        static double inv_AA_p_LAB( double Tn_proj_LAB, double T_p_LAB, double eta_LAB, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization );
        
        

        static double dE_AA_p_LAB( double Tn_proj_LAB, double T_p_LAB, int A_projectile=1, int N_projectile=0, int A_target=1, int N_target=0, int parametrization=ANDERSON);
        static double integrand__dE_AA_p_LAB (double eta_LAB, void* parameters  );
        
        
        

        static double p_coal__VonDoetinchen( double s );
        
        

       static double p_coal__pTdep( double pToverA, double p0_val=0.160 );
        
        

        static double inv_AA_Dbar_CM( double s, double xF_Dbar, double pT_Dbar, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization, int coalescence, double p0_val=0.160 );
        

        static double inv_AA_Dbar_LAB( double Tn_proj_LAB, double Tn_Dbar_LAB, double eta_LAB, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization, int coalescence, double p0_val=0.160 );
        static double integrand__dE_AA_Dbar_LAB (double eta_LAB, void* parameters  );
        
        

        static double dEn_AA_Dbar_LAB( double Tn_proj_LAB, double Tn_Dbar_LAB, int A_projectile=1, int N_projectile=0, int A_target=1, int N_target=0, int parametrization=KORSMEIER_II, int coalescence=ENERGY_DEP__VAN_DOETINCHEM, double p0_val=0.160 );
        
        

        static double dEn_DbarA_Dbar_LAB( double Tn_Dbar_proj_LAB, double Tn_Dbar_prod_LAB, int A_target=1, int N_target=0, int parametrization=ANDERSON );
        
        

        static double inv_AA_He3bar_CM( double s, double xF_Hebar, double pT_Hebar, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization, int coalescence, double p0_val=0.160 );
        

        static double inv_AA_He3bar_LAB( double Tn_proj_LAB, double Tn_Hebar_LAB, double eta_LAB, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization, int coalescence, double p0_val=0.160 );
        static double integrand__dE_AA_He3bar_LAB (double eta_LAB, void* parameters  );
        
        

        static double dEn_AA_He3bar_LAB( double Tn_proj_LAB, double Tn_Hebar_LAB, int A_projectile=1, int N_projectile=0, int A_target=1, int N_target=0, int parametrization=KORSMEIER_II, int coalescence=ENERGY_DEP__VAN_DOETINCHEM, double p0_val=0.160 );
        
        

        static double dEn_He3barA_He3bar_LAB( double Tn_Hebar_proj_LAB, double Tn_Hebar_prod_LAB, int A_target=1, int N_target=0, int parametrization=ANDERSON );
        
        

       static double inv_AA_He4bar_CM( double s, double xF_Hebar, double pT_Hebar, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization, int coalescence, double p0_val=0.160 );
       

       static double inv_AA_He4bar_LAB( double Tn_proj_LAB, double Tn_Hebar_LAB, double eta_LAB, int A_projectile, int N_projectile, int A_target, int N_target, int parametrization, int coalescence, double p0_val=0.160 );
       static double integrand__dE_AA_He4bar_LAB (double eta_LAB, void* parameters  );
       
       

       static double dEn_AA_He4bar_LAB( double Tn_proj_LAB, double Tn_Hebar_LAB, int A_projectile=1, int N_projectile=0, int A_target=1, int N_target=0, int parametrization=KORSMEIER_II, int coalescence=ENERGY_DEP__VAN_DOETINCHEM, double p0_val=0.160 );
       
       

       static double dEn_He4barA_He4bar_LAB( double Tn_Hebar_proj_LAB, double Tn_Hebar_prod_LAB, int A_target=1, int N_target=0, int parametrization=ANDERSON );

        static void Set_SELF_C_parameters_diMauro(double* C);
        
        

        static void Set_SELF_C_parameters_Winkler(double* C);
        
        

        static void Set_SELF_D_parameters_diMauro(double* D);
        
        

        static void Set_SELF_D_parameters_Winkler(double* D);
        
        
        
        static bool     fIsRestricted_pp;
        static void     SetRestricted_pp     ( bool is_pp ){fIsRestricted_pp=is_pp;};
        
        static bool     isInRestricted_CM                  ( double s, double xf, double pT );
        static void     SetRestrictedParameterSpace_CM     ( double s, double xf, double pT );
        static void     RemoveRestrictedParameterSpace_CM  (  );
        static int      fRestrictedParameterSpace_CM;
        static double   fRestrictedParameterSpace_CM__s    [103];
        static double   fRestrictedParameterSpace_CM__xf   [103];
        static double   fRestrictedParameterSpace_CM__pT   [103];
        
        
        
        static bool     isInRestricted_LAB                  ( double Tp, double Tpbar, double eta );
        static void     SetRestrictedParameterSpace_LAB     ( double Tp, double Tpbar, double eta );
        static void     RemoveRestrictedParameterSpace_LAB  (  );
        static int      fRestrictedParameterSpace_LAB;
        static double   fRestrictedParameterSpace_LAB__Tp   [103];
        static double   fRestrictedParameterSpace_LAB__Tpbar[103];
        static double   fRestrictedParameterSpace_LAB__eta  [103];
        
        //        static double fRestrictedParameterSpace_LAB__Tp_Tpbar_eta[100][3];
        
        

    private:
        
        


        
    };
}

#endif