Flushing the buffer for debug purposes (R. Preghenella)

[u/mrichter/AliRoot.git] / EVGEN / AliSlowNucleonModelExp.h

diff --git a/EVGEN/AliSlowNucleonModelExp.h b/EVGEN/AliSlowNucleonModelExp.h
index 2a14b7415d5ba64e1fd07060d6d8f0da12da3d95..69aef96ef1d994ab685bb3b28ddfc19a525fa6f0 100644 (file)
--- a/EVGEN/AliSlowNucleonModelExp.h
+++ b/EVGEN/AliSlowNucleonModelExp.h
@@ -5,6 +5,12 @@
 
 /* $Id$ */
 
+//
+// Experimental data inspired Gray Particle Model for p-Pb collisions
+// Fluctuations are calculated from a binomial distribution.
+// Author: A.Morsch
+//
+
 #include "AliSlowNucleonModel.h"
 
 class AliCollisionGeometry;
@@ -15,16 +21,39 @@ class AliSlowNucleonModelExp : public AliSlowNucleonModel
     AliSlowNucleonModelExp();
     virtual ~AliSlowNucleonModelExp(){;}
     virtual void GetNumberOfSlowNucleons(AliCollisionGeometry* geo,
-                                        Int_t& ngp, Int_t& ngn, Int_t& nbp, Int_t& nbn);
+                                        Int_t& ngp, Int_t& ngn, Int_t& nbp, Int_t& nbn) const;
+    virtual void GetNumberOfSlowNucleons2(AliCollisionGeometry* geo,
+                                        Int_t& ngp, Int_t& ngn, Int_t& nbp, Int_t& nbn) const;
+    virtual void GetNumberOfSlowNucleons2s(AliCollisionGeometry* geo,
+                                        Int_t& ngp, Int_t& ngn, Int_t& nbp, Int_t& nbn) const;
+    // 1st model
     virtual void SetParameters(Float_t alpha1, Float_t alpha2);
+    virtual void SetSaturation(Bool_t saturation) {fApplySaturation = saturation;}
+    virtual void SetSaturationParams(Int_t ngray=15, Int_t nblack=28) 
+               {fnGraySaturation=ngray; fnBlackSaturation=nblack;}
+    // 2nd model
+    virtual void SetLCPparam(Float_t al) {fLCPparam=al;}
+    virtual void SetNslowParams(Float_t a, Float_t b, Float_t c) 
+                {fSlownparam[0]=a; fSlownparam[1]=b; fSlownparam[2]=c;}
     
  protected:
     Float_t  fP;          // Number of protons  in the target 
     Float_t  fN;          // Number of neutrons in the target
     Float_t  fAlphaGray;  // Proportionality between gray   particles and number of collisions
     Float_t  fAlphaBlack; // Proportionality between black  particles and number of collisions
+    Bool_t   fApplySaturation;  // If true apply satoration to N_black vs. N_gray
+    Int_t    fnGraySaturation;  // N_gray value for N_black saturation
+    Int_t    fnBlackSaturation; // N_black saturation value
+    //
+    // Adding parameters for 2nd model that can be tuned during config
+    Float_t  fLCPparam;                // parameter to calculate LCP from <Nslow p>
+    Float_t  fSlownparam[3];   // parameters to calculate <Nslow n> from LCP
+    //
+    // Adding parameter to smear the number of slow nucleons
+    Float_t  fSigmaSmear;
+    
     
-  ClassDef(AliSlowNucleonModelExp,1) // Gray Particle Model (Experiment inspired)
+  ClassDef(AliSlowNucleonModelExp, 4) // Gray Particle Model (Experiment inspired)
 };
 #endif
 
@@ -33,3 +62,4 @@ class AliSlowNucleonModelExp : public AliSlowNucleonModel
 
 
 
+