STARLIGHT/starlight/include/.svn/text-base/gammaavm.h.svn-base

   1 ///////////////////////////////////////////////////////////////////////////
   2 //
   3 //    Copyright 2010
   4 //
   5 //    This file is part of starlight.
   6 //
   7 //    starlight is free software: you can redistribute it and/or modify
   8 //    it under the terms of the GNU General Public License as published by
   9 //    the Free Software Foundation, either version 3 of the License, or
  10 //    (at your option) any later version.
  11 //
  12 //    starlight is distributed in the hope that it will be useful,
  13 //    but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 //    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
  15 //    GNU General Public License for more details.
  16 //
  17 //    You should have received a copy of the GNU General Public License
  18 //    along with starlight. If not, see <http://www.gnu.org/licenses/>.
  19 //
  20 ///////////////////////////////////////////////////////////////////////////
  21 //
  22 // File and Version Information:
  23 // $Rev::                             $: revision of last commit
  24 // $Author::                          $: author of last commit
  25 // $Date::                            $: date of last commit
  26 //
  27 // Description:
  28 //
  29 //
  30 //
  31 ///////////////////////////////////////////////////////////////////////////
  32
  33
  34 #ifndef GAMMAAVM_H
  35 #define GAMMAAVM_H
  36
  37
  38 #include <vector>
  39
  40 #include "starlightconstants.h"
  41 #include "readinluminosity.h"
  42 #include "beambeamsystem.h"
  43 #include "randomgenerator.h"
  44 #include "eventchannel.h"
  45 #include "upcevent.h"
  46 #include "nBodyPhaseSpaceGen.h"
  47
  48
  49 class Gammaavectormeson : public eventChannel
  50 {
  51
  52  public:
  53   Gammaavectormeson(beamBeamSystem& bbsystem);
  54   virtual ~Gammaavectormeson();
  55   starlightConstants::event produceEvent(int &ievent);
  56
  57    upcEvent produceEvent();
  58
  59   void pickwy(double &W, double &Y);
  60   void momenta(double W,double Y,double &E,double &px,double &py,double &pz,int &tcheck);
  61   double pTgamma(double E);
  62   void vmpt(double W,double Y,double &E,double &px,double &py, double &pz,int &tcheck);
  63   void twoBodyDecay(starlightConstants::particleTypeEnum &ipid,double E,double W,double px0,double py0,double pz0,double &px1,double &py1,double&pz1,double &px2,double &py2,double &pz2,int &iFbadevent);
  64   bool fourBodyDecay(starlightConstants::particleTypeEnum& ipid, const double E, const double W, const double* p, lorentzVector* decayMoms, int& iFbadevent);
  65   double getMass();
  66   double getWidth();
  67   virtual double getTheta(starlightConstants::particleTypeEnum ipid);
  68   double getSpin();
  69   double _VMbslope;
  70   virtual double getDaughterMass(starlightConstants::particleTypeEnum &ipid);
  71   double pseudoRapidity(double px, double py, double pz);
  72
  73  private:
  74   starlightConstants::particleTypeEnum _VMpidtest;
  75   int _VMnumw;
  76   int _VMnumy;
  77   int _VMinterferencemode;
  78   int _VMCoherence;
  79   int _ProductionMode;
  80   double _VMCoherenceFactor;
  81   double _VMgamma_em;
  82   double _VMNPT;
  83   double _VMWmax;
  84   double _VMWmin;
  85   double _VMYmax;
  86   double _VMYmin;
  87   double _mass;
  88   double _width;
  89   double _VMptmax;
  90   double _VMdpt;
  91   nBodyPhaseSpaceGen* _phaseSpaceGen;
  92
  93 };
  94
  95 class Gammaanarrowvm : public Gammaavectormeson
  96 {
  97  public:
  98   Gammaanarrowvm(beamBeamSystem& bbsystem);
  99   virtual ~Gammaanarrowvm();
 100 };
 101
 102 class Gammaawidevm : public Gammaavectormeson
 103 {
 104  public:
 105   Gammaawidevm(beamBeamSystem& bbsystem);
 106   virtual ~Gammaawidevm();
 107 };
 108
 109 class Gammaaincoherentvm : public Gammaavectormeson
 110 {
 111  public:
 112   Gammaaincoherentvm(beamBeamSystem& bbsystem);
 113   virtual ~Gammaaincoherentvm();
 114 };
 115
 116 #endif  // GAMMAAVM_H