const Float_t kParamSkinPb[2]={0.762408, 20.};
// Thickness of nuclear surface
- const Float_t kNuclearThick = 0.52;
+ //const Float_t kNuclearThick = 0.52;
// Maximum impact parameter for U [r0*A**(1/3)]
const Float_t kbMaxU = 14.87;
// Maximum impact parameter for Pb [r0*A**(1/3)]
//Float_t bCore = kbMaxPb-2*kNuclearThick;
if(fB>kbMaxPb){
fZbAverage = TMath::Exp(-kParamSkinPb[0]*(fB-kParamSkinPb[1]));
- printf(" b = %1.2f fm Z_bound %1.2f\n", fB, fZbAverage);
+ //printf(" b = %1.2f fm Z_bound %1.2f\n", fB, fZbAverage);
}
if(fZbAverage>kZProj) fZbAverage = kZProj;
Float_t zbNorm = fZbAverage/kZProj;