TTherminator/Therminator/Particle.cxx

   1 /******************************************************************************
   2  *                      T H E R M I N A T O R                                 *
   3  *                   THERMal heavy-IoN generATOR                              *
   4  *                           version 1.0                                      *
   5  *                                                                            *
   6  * Authors of the model: Wojciech Broniowski, Wojciech.Broniowski@ifj.edu.pl, *
   7  *                       Wojciech Florkowski, Wojciech.Florkowski@ifj.edu.pl  *
   8  * Authors of the code:  Adam Kisiel, kisiel@if.pw.edu.pl                     *
   9  *                       Tomasz Taluc, ttaluc@if.pw.edu.pl                    *
  10  * Code designers: Adam Kisiel, Tomasz Taluc, Wojciech Broniowski,            *
  11  *                 Wojciech Florkowski                                        *
  12  *                                                                            *
  13  * For the detailed description of the program and furhter references         *
  14  * to the description of the model plesase refer to: nucl-th/0504047,         *
  15  * accessibile at: http://www.arxiv.org/nucl-th/0504047                       *
  16  *                                                                            *
  17  * Homepage: http://hirg.if.pw.edu.pl/en/therminator/                         *
  18  *                                                                            *
  19  * This code can be freely used and redistributed. However if you decide to   *
  20  * make modifications to the code, please contact the authors, especially     *
  21  * if you plan to publish the results obtained with such modified code.       *
  22  * Any publication of results obtained using this code must include the       *
  23  * reference to nucl-th/0504047 and the published version of it, when         *
  24  * available.                                                                 *
  25  *                                                                            *
  26  *****************************************************************************/
  27 #include "Particle.h"
  28 #include "ParticleType.h"
  29 #include <TMath.h>
  30 #define FIELDWIDTH 15
  31
  32 Particle::Particle()
  33 {
  34   mPartType  = 0;
  35   px = 0; py = 0; pz = 0;
  36   rx = 0; ry = 0; rz = 0; rt = 0;
  37   mDecayed = 0;
  38   mHasFather = -1;
  39 }
  40
  41 Particle::Particle(double aRapidity, double aPt, double aPhip,
  42                    double aAlfam, double aRho, double aPhis, double aTau,
  43                    ParticleType *aType)
  44 {
  45   mPartType = aType;
  46   mDecayed = 0;
  47   mHasFather = -1;
  48
  49   px = aPt*TMath::Cos(aPhip);
  50   py = aPt*TMath::Sin(aPhip);
  51   double tMt = TMath::Hypot(GetMass(),aPt);
  52   pz = tMt*TMath::SinH(aRapidity);
  53
  54   rx = aRho*TMath::Cos(aPhis);
  55   ry = aRho*TMath::Sin(aPhis);
  56
  57 // New method of calculating rz, rt
  58   rz = aTau*TMath::SinH(aAlfam);
  59   rt = aTau*TMath::CosH(aAlfam);
  60 }
  61
  62 Particle::Particle(ParticleType *aType,
  63                    double aPx, double aPy, double aPz,
  64                    double aRx, double aRy, double aRz,
  65                    double aTime)
  66 {
  67   mPartType = aType;
  68   mDecayed = 0;
  69   mHasFather = -1;
  70
  71   px = aPx;
  72   py = aPy;
  73   pz = aPz;
  74
  75   rx = aRx;
  76   ry = aRy;
  77   rz = aRz;
  78   rt = aTime;
  79 }
  80
  81
  82 Particle::Particle(const Particle& aParticle)
  83 {
  84   mPartType = aParticle.GetParticleType();
  85   mDecayed = 0;
  86   mHasFather = aParticle.GetFather();
  87
  88   px = aParticle.px;
  89   py = aParticle.py;
  90   pz = aParticle.pz;
  91
  92   rx = aParticle.rx;
  93   ry = aParticle.ry;
  94   rz = aParticle.rz;
  95   rt = aParticle.rt;
  96 }
  97
  98 Particle::~Particle()
  99 {
 100   /* no-op */
 101 }
 102
 103 void
 104 Particle::WriteParticle(ostream *aOuts)
 105 {
 106   aOuts->flags(ios::left | ios::scientific);
 107   aOuts->width(10);
 108   (*aOuts) << mPartType->GetPDGCode();
 109   aOuts->width(FIELDWIDTH);
 110   (*aOuts) << px;
 111   aOuts->width(FIELDWIDTH);
 112   (*aOuts) << py;
 113   aOuts->width(FIELDWIDTH);
 114   (*aOuts) << pz;
 115   aOuts->width(FIELDWIDTH);
 116   (*aOuts) << GetEnergy();
 117   aOuts->width(FIELDWIDTH);
 118   (*aOuts) << GetMass();
 119   aOuts->width(FIELDWIDTH);
 120   (*aOuts) << rx;
 121   aOuts->width(FIELDWIDTH);
 122   (*aOuts) << ry;
 123   aOuts->width(FIELDWIDTH);
 124   (*aOuts) << rz;
 125   aOuts->width(FIELDWIDTH);
 126   (*aOuts) << rt;
 127 //  (*aOuts) << mPartType->GetCharge() << "   ";
 128 //  if (mHasFather>-1) (*aOuts) << mHasFather << "   ";
 129   aOuts->width(6);
 130   (*aOuts) << (mHasFather);
 131   aOuts->width(2);
 132   (*aOuts) << (mDecayed);
 133 }
 134
 135
 136 double
 137 Particle::Pt()
 138 {
 139   return TMath::Hypot(px, py);
 140 }
 141
 142 double
 143 Particle::Rapidity()
 144 {
 145   double tE = TMath::Sqrt(px*px+py*py+pz*pz+GetMass()*GetMass());
 146   return 0.5*TMath::Log((tE-pz)/(tE+pz));
 147 }
 148
 149 ParticleType*
 150 Particle::GetParticleType() const
 151 {
 152   return mPartType;
 153 }
 154
 155 int
 156 Particle::HadDecayed()
 157 {
 158   return mDecayed;
 159 }
 160
 161 double
 162 Particle::GetMass()
 163 {
 164   return mPartType->GetMass();
 165 }
 166
 167 double
 168 Particle::GetEnergy()
 169 {
 170   return TMath::Sqrt(GetMass()*GetMass()+px*px+py*py+pz*pz);
 171 }
 172
 173
 174 double
 175 Particle::GetI3()
 176 {
 177   return mPartType->GetI3();
 178 }
 179
 180 double
 181 Particle::GetBarionN()
 182 {
 183   return mPartType->GetBarionN();
 184 }
 185
 186 double
 187 Particle::GetStrangeness()
 188 {
 189   return mPartType->GetStrangeness();
 190 }
 191
 192 void
 193 Particle::SetDecayed()
 194 {
 195   mDecayed = 1;
 196 }
 197
 198 void
 199 Particle::SetFather(int aFather)
 200 {
 201   mHasFather = aFather;
 202 }
 203
 204 int
 205 Particle::GetFather() const
 206 {
 207   return mHasFather;
 208 }
 209