-
-/*if(ncoll==1){
- printf(" Ncoll %d ", ncoll);
- printf(" Ngp %1.0f Nbp %1.0f", ngp, nbp);
- printf(" Ngn %1.0f Nbn %1.0f \n", ngn, nbn);
-}*/
-/*if(ngn+nbn<0.005){
- printf(" Ncoll %d ", ncoll);
- printf(" <Ngp> %f <Nbp> %f \t Ngp %1.0f Nbp %1.0f \t LCP %f\n", nGrayp, nBlackp, ngp, nbp,lcp);
- printf(" <Nslow,n> %f Ngn %1.0f Nbn %1.0f \n", nSlow, ngn, nbn);
-}*/
-
-}
-
-
-//--------------------------------------------------------------------------------------------------
-/* void AliCentralityGlauberFit::MakeSlowNucleonsAM(Int_t ncoll, Double_t alpha, Double_t k, Double_t bog,
- Double_t gamma, Double_t &nbn, Double_t &ngn, Double_t &nbp, Double_t &ngp)
-{
-// from AliGenSlowNucleonModelExp (Chiara Oppedisano)
-//
-// Return the number of black and gray nucleons
-//
-// Number of collisions
-
- // based on E910 model ================================================================
-
- Int_t fP = 82;
- Int_t fN = 126;
-
- Float_t nu = (Float_t) ncoll;
- Float_t sigmap = 0.25;
-
- // gray protons
- Float_t poverpd = 0.843;
- Float_t zAu2zPb = 82./79.;
- Float_t grayp = (-0.27 + 0.63 * nu - 0.0008 *nu *nu)*poverpd*zAu2zPb;
- Float_t nGrayp = gRandom->Gaus(grayp, sigmap);
- //if(nGrayp<0.) nGrayp=0.;
-
- Double_t p=0.;
- p = nGrayp/fP;
- ngp = gRandom->Binomial((Int_t) fP, p);
- //ngp = gRandom->Gaus(nGrayp, TMath::Sqrt(fP*p*(1-p)));
-
- // black protons
- //Float_t blackovergray = 3./7.;// from spallation
- //Float_t blackovergray = 0.65; // from COSY
- Float_t blackovergray = bog;
- //
- //Float_t blackp = blackovergray*nGrayp;
- //Float_t nBlackp = gRandom->Gaus(blackp, sigmap);
- Float_t nBlackp = blackovergray*nGrayp;
- //if(nBlackp<0.) nBlackp=0.;
-
- p = nBlackp/fP;
- nbp = gRandom->Binomial((Int_t) fP, p);
-