Possibility to switch off heavy flavor production added.

[u/mrichter/AliRoot.git] / TMEVSIM / TMevSimPartTypeParams.cxx

diff --git a/TMEVSIM/TMevSimPartTypeParams.cxx b/TMEVSIM/TMevSimPartTypeParams.cxx
index 8ce8fb1d5a6a8e191076ad134572d75876d8fc77..0ae4775735491f851dafbcb24a4bf4fc12f8e1f4 100644 (file)
--- a/TMEVSIM/TMevSimPartTypeParams.cxx
+++ b/TMEVSIM/TMevSimPartTypeParams.cxx
@@ -34,8 +34,8 @@ ClassImp(TMevSimPartTypeParams)
 
 //__________________________________________________________________________
 TMevSimPartTypeParams::TMevSimPartTypeParams(Int_t agpid, Int_t amultmean, Int_t amultvc, 
-                                            Float_t atempmean, Float_t atempstdev, Float_t asigmamean,
-                                            Float_t asigmastdev, Float_t aexpvelmean, Float_t aexpvelstdev) 
+                                             Float_t atempmean, Float_t atempstdev, Float_t asigmamean,
+                                             Float_t asigmastdev, Float_t aexpvelmean, Float_t aexpvelstdev) 
 {
 // Construct the particle type parametrs class. Use the values provide
 // by the user to initialize the internal variables. For the meaning of
@@ -59,10 +59,13 @@ TMevSimPartTypeParams::TMevSimPartTypeParams(Int_t agpid, Int_t amultmean, Int_t
 }
 
 //______________________________________________________________________________
-TMevSimPartTypeParams::~TMevSimPartTypeParams() {}
+TMevSimPartTypeParams::~TMevSimPartTypeParams() 
+{
+  //dtor
+}
 
 //______________________________________________________________________________
-TMevSimPartTypeParams::TMevSimPartTypeParams (const TMevSimPartTypeParams& pars) {
+TMevSimPartTypeParams::TMevSimPartTypeParams (const TMevSimPartTypeParams& pars) : TObject(pars) {
 // The copy constructor
   
    this->fGPid = pars.GetGPid();
@@ -104,72 +107,88 @@ TMevSimPartTypeParams& TMevSimPartTypeParams::operator=(const TMevSimPartTypePar
 
 //______________________________________________________________________________
 void    TMevSimPartTypeParams::SetGPid(Int_t gpid) {
+//Sets pid
    fGPid = gpid;
 }
 //______________________________________________________________________________
 Int_t   TMevSimPartTypeParams::GetGPid() const {
+//returns pid
    return fGPid;
 }
 //______________________________________________________________________________
 void    TMevSimPartTypeParams::SetMultMean(Int_t multmean) {
+//sets multiplicity mean
    fMultMean = multmean;
 }
 //______________________________________________________________________________
 Int_t TMevSimPartTypeParams::GetMultMean() const {
+//returns multiplicity mean
    return fMultMean;
 }
 //______________________________________________________________________________
 void    TMevSimPartTypeParams::SetMultVarianceControl(Int_t multvc) {
+//sets multiplicity variance
    fMultVarianceControl = multvc;
 }
 //______________________________________________________________________________
 Int_t   TMevSimPartTypeParams::GetMultVarianceControl() const {
+//returns multiplicity variance
    return fMultVarianceControl;
 }
 //______________________________________________________________________________
 void    TMevSimPartTypeParams::SetTempParams(Float_t tempmean, Float_t tempsdev) { 
+//Sets mean and variance of temperature
    fTempMean = tempmean;
    fTempStDev = tempsdev;
 }
 //______________________________________________________________________________
 Float_t TMevSimPartTypeParams::GetTempMean() const { 
+//returns mean temperature
    return fTempMean;
 }
 //______________________________________________________________________________
 Float_t TMevSimPartTypeParams::GetTempStDev() const { 
+//returns temperature variance
    return fTempStDev;
 }
 //______________________________________________________________________________
 void    TMevSimPartTypeParams::SetSigmaPrams(Float_t sigmamean, Float_t sigmastdev) {  
+//Sets mean and variance of rapidity
    fSigmaMean = sigmamean;
    fSigmaStDev = sigmastdev;
 }
 //______________________________________________________________________________
 Float_t TMevSimPartTypeParams::GetSigmaMean() const { 
+//returns mean rapidity
    return fSigmaMean;
 }
 //______________________________________________________________________________
 Float_t TMevSimPartTypeParams::GetSigmaStDev() const { 
+//returns rapidity variance
    return fSigmaStDev;
 }
 //______________________________________________________________________________
 void    TMevSimPartTypeParams::SetExpVelParams(Float_t expvelmean, Float_t expvelstdev) {  
+//Sets mean and variance of  Expansion velocity ala Scott Pratt (in units of c)
    fExpVelMean = expvelmean;
    fExpVelStDev = expvelstdev;
 }
 //______________________________________________________________________________
 Float_t TMevSimPartTypeParams::GetExpVelMean() const { 
+//Sets mean of Expansion velocity ala Scott Pratt (in units of c)
    return fExpVelMean;
 }
 //______________________________________________________________________________
 Float_t TMevSimPartTypeParams::GetExpVelStDev() const { 
+//Sets variance of  Expansion velocity ala Scott Pratt (in units of c)
    return fExpVelStDev;
 }
 //______________________________________________________________________________
 
 void TMevSimPartTypeParams::SetVnMeanComponent
 (Int_t nComponent, Float_t mean1, Float_t mean2,  Float_t mean3, Float_t mean4)  {  
-
+  //Sets VnMean: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6
+  //                 mean include all 6 sets of parameters, set them to 0 if not used                                
 
   if (nComponent < 0 || nComponent > NFLOWTERMS )
     Error("SetVnMeanComponent", "Wrong Vn component n = %d (must be [%d , %d])", 
@@ -184,6 +203,8 @@ void TMevSimPartTypeParams::SetVnMeanComponent
 
 void TMevSimPartTypeParams::SetVnStDevComponent
 (Int_t nComponent, Float_t stdev1, Float_t stdev2,Float_t stdev3, Float_t stdev4) {  
+  //Sets VnStDev: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6
+  //              standard deviation include all 6 sets of parameters, set them to 0 if not used                                
 
   if (nComponent < 0 || nComponent > NFLOWTERMS )
     Error("SetVnStDevComponent", "Wrong Vn component n = %d (must be [%d , %d])",
@@ -198,13 +219,14 @@ void TMevSimPartTypeParams::SetVnStDevComponent
 //______________________________________________________________________________
 
 void TMevSimPartTypeParams::SetVnMeanComponent (Int_t nComponent, Float_t mean[4]) {
-  
+  //Sets VnMean: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6
+  //                 mean include all 6 sets of parameters, set them to 0 if not used                                
   if (nComponent < 0 || nComponent > NFLOWTERMS )
     Error("SetVnMeanComponent", "Wrong Vn component (%d) (must be [%d,%d])",
        nComponent, 0, NFLOWTERMS);
   
-  Int_t N = 4;
-  for (Int_t i=0; i<N; i++) 
+  Int_t n = 4;
+  for (Int_t i=0; i<n; i++) 
     fVnMean[nComponent][i] = mean[i];
 
 }
@@ -212,25 +234,31 @@ void TMevSimPartTypeParams::SetVnMeanComponent (Int_t nComponent, Float_t mean[4
 //______________________________________________________________________________
 
 void TMevSimPartTypeParams::SetVnStDevComponent (Int_t nComponent, Float_t stdev[4]) {
+  //Sets VnStDev: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6
+  //              standard deviation include all 6 sets of parameters, set them to 0 if not used                                
 
   if (nComponent < 0 || nComponent > NFLOWTERMS )
     Error("SetVnStDevComponent", "Wrong Vn component n = %d (must be [%d , %d])",
          nComponent, 0, NFLOWTERMS);
   
-  Int_t N = 4;
-  for (Int_t i=0; i<N; i++) 
+  Int_t n = 4;
+  for (Int_t i=0; i<n; i++) 
     fVnStDev[nComponent][i] = stdev[i];  
 }
 
 //______________________________________________________________________________
 
 Float_t     TMevSimPartTypeParams::GetVnMeanComponent(Int_t nComponent, Int_t nMean) const {
+  //Returns VnMean: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6
+  //                 mean include all 6 sets of parameters, set them to 0 if not used                                
 
   if ((nComponent < 0) || (nComponent>NFLOWTERMS) || (nMean < 0) || (nMean > 3)) return 0.0;
   return fVnMean[nComponent][nMean];
 }
 //______________________________________________________________________________
 Float_t     TMevSimPartTypeParams::GetVnStDevComponent(Int_t nComponent, Int_t nStDev) const {
+  //Returns VnStDev: Anisotropic flow parameters for Fourier components NFLOWTERMS=1,6
+  //              standard deviation include all 6 sets of parameters, set them to 0 if not used                                
 
   if ((nComponent < 0) || (nComponent>NFLOWTERMS) || (nStDev < 0) || (nStDev > 3)) return 0.0;
   return fVnMean[nComponent][nStDev];