//////////////////////////////////////////////
// Global info
//////////////////////////////////////////////

// constants
const Double_t HbarC = 0.197; // GeV fm
const Double_t AlphaEM = 1.0/137.0;

// masses
const Double_t RhoMass = 0.7755; // GeV
const Double_t ProtonMass = 0.938; // GeV
const Double_t JPsiMass = 3.0969; // GeV
const Double_t Psi2sMass = 3.686; // GeV
const Double_t Upsilon1sMass = 9.460; // GeV

// Parameters
//const Double_t pBeamEnergy = 6500; // GeV
// gold: 79, 197; Pb: 82, 208; Xe:54, 129
const Double_t IonChargePb = 82;
const Double_t IonChargeAu = 79;
const Double_t IonChargeXe = 54; 
const Double_t IonAtomicNumberPb = 208;
const Double_t IonAtomicNumberAu = 197;
const Double_t IonAtomicNumberXe = 129;
Double_t BeamEnergy = 0;
Double_t LorentzBoost = 0;
Double_t IonRadius = 0;
Double_t IonAtomicNumber = 0;
Double_t pBeamEnergy = 0; // GeV
Double_t IonCharge = 0;
Double_t mVM = 0;


// global constants: form factor a la SL 
const Double_t a_FormFactor = 0.7; // fermi
const Double_t rho_FormFactor = 0.13815; //  fixed by normalization to lead
const Double_t R_A_Pb = 1.2*TMath::Power(IonAtomicNumber,1.0/3.0); // fm

// global constants: Woods-Saxon for lead
const Double_t WSc_Pb = 6.624; // fm
const Double_t WSz_Pb = 0.549; // fm
const Double_t WSrho0_Pb = 0.16; // fm-3
const Double_t WSw_Pb = 0; 

// global constants: Woods-Saxon for gold
const Double_t WSc_Au = 6.38; // fm
const Double_t WSz_Au = 0.535; // fm
const Double_t WSrho0_Au = 0.1693; // fm-3
const Double_t WSw_Au = 0; 

// nuclear density for Xe (129,54)
const Double_t Xe_R = 5.36; // fm
const Double_t Xe_a = 0.59; // fm
const Double_t Xe_b2 = 0.18;
const Double_t Xe_rho0 = 0.19868; // fm^-3