Modified parameters (Chiara)
authorsvnbits <svnbits@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 22 Jan 2013 15:44:56 +0000 (15:44 +0000)
committersvnbits <svnbits@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 22 Jan 2013 15:44:56 +0000 (15:44 +0000)
EVGEN/AliGenSlowNucleons.cxx
EVGEN/AliSlowNucleonModelExp.cxx
EVGEN/AliSlowNucleonModelExp.h

index a1f69f2..7b9a52b 100644 (file)
@@ -195,7 +195,7 @@ void AliGenSlowNucleons::Generate()
        GenerateSlow(fCharge, fTemperatureG, fBetaSourceG, p, theta);
        if (fDebug) fCosThetaGrayHist->Fill(TMath::Cos(theta));
        PushTrack(fTrackIt, -1, kf, p, origin, polar,
-                time, kPNoProcess, nt, 1.,-1);
+                time, kPNoProcess, nt, 1.,1);
        KeepTrack(nt);
     }
 //
@@ -207,7 +207,7 @@ void AliGenSlowNucleons::Generate()
        GenerateSlow(fCharge, fTemperatureG, fBetaSourceG, p, theta);
        if (fDebug) fCosThetaGrayHist->Fill(TMath::Cos(theta));
        PushTrack(fTrackIt, -1, kf, p, origin, polar,
-                time, kPNoProcess, nt, 1.,-1);
+                time, kPNoProcess, nt, 1.,1);
        KeepTrack(nt);
     }
 //
@@ -218,7 +218,7 @@ void AliGenSlowNucleons::Generate()
     for(i = 0; i < fNbp; i++) {
        GenerateSlow(fCharge, fTemperatureB, fBetaSourceB, p, theta);
        PushTrack(fTrackIt, -1, kf, p, origin, polar,
-                time, kPNoProcess, nt, 1.,-1);
+                time, kPNoProcess, nt, 1.,1);
        KeepTrack(nt);
     }
 //
@@ -229,7 +229,7 @@ void AliGenSlowNucleons::Generate()
     for(i = 0; i < fNbn; i++) {
        GenerateSlow(fCharge, fTemperatureB, fBetaSourceB, p, theta);
        PushTrack(fTrackIt, -1, kf, p, origin, polar,
-                time, kPNoProcess, nt, 1.,-1);
+                time, kPNoProcess, nt, 1.,1);
        KeepTrack(nt);
     }
 }
index 0fc562d..9598cb2 100644 (file)
@@ -33,8 +33,8 @@ ClassImp(AliSlowNucleonModelExp)
 AliSlowNucleonModelExp::AliSlowNucleonModelExp():
     fP(82),
     fN (208 - 82),
-    fAlphaGray(2.),
-    fAlphaBlack(4.),
+    fAlphaGray(2.3),
+    fAlphaBlack(3.6),
     fApplySaturation(kTRUE),
     fnGraySaturation(15),
     fnBlackSaturation(28)
@@ -42,6 +42,10 @@ AliSlowNucleonModelExp::AliSlowNucleonModelExp():
   //
   // Default constructor
   //
+  //
+  printf("\n\nInitializing slow nucleon model with parameters:\n");
+  printf(" \t alpha_{gray} %1.2f  alpha_{black} %1.2f\n",fAlphaGray, fAlphaBlack);
+  printf(" \t SATURATION %d w. %d (gray) %d (black) \n\n",fApplySaturation,fnGraySaturation,fnBlackSaturation);
 }
 
 
@@ -65,7 +69,7 @@ void AliSlowNucleonModelExp::GetNumberOfSlowNucleons(AliCollisionGeometry* geo,
     Float_t nBlack  = 0.;
     if(!fApplySaturation || (fApplySaturation && nGray<fnGraySaturation)) nBlack = fAlphaBlack * nu;
     else if(fApplySaturation && nGray>=fnGraySaturation) nBlack = fnBlackSaturation;
-    Float_t nBlackNeutrons = nBlack * fN / (fN + fP);
+    Float_t nBlackNeutrons = nBlack * 0.84;
     Float_t nBlackProtons  = nBlack - nBlackNeutrons;
 
 // Actual number (including fluctuations) from binomial distribution
index d49e0ae..83a240d 100644 (file)
@@ -24,7 +24,7 @@ class AliSlowNucleonModelExp : public AliSlowNucleonModel
                                         Int_t& ngp, Int_t& ngn, Int_t& nbp, Int_t& nbn) const;
     virtual void SetParameters(Float_t alpha1, Float_t alpha2);
     virtual void SetSaturation(Bool_t saturation) {fApplySaturation = saturation;}
-    virtual void SetSaturationParams(Float_t ngray=15, Float_t nblack=28) 
+    virtual void SetSaturationParams(Int_t ngray=15, Int_t nblack=28) 
                {fnGraySaturation=ngray; fnBlackSaturation=nblack;}
     
  protected: