STARLIGHT/starlight/include/.svn/text-base/gammagammasingle.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 GAMMAGAMMASINGLE_H
  35 #define GAMMAGAMMASINGLE_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 "starlightconfig.h"
  46 #ifdef ENABLE_PYTHIA
  47 #include "pythiadecayer.h"
  48 #endif
  49
  50 class Gammagammasingle : public eventChannel
  51 {
  52  public:
  53   Gammagammasingle(beamBeamSystem& bbsystem);
  54   ~Gammagammasingle();
  55
  56   void singleCrossSection();
  57   starlightConstants::event produceEvent(int &ievent);
  58
  59   upcEvent produceEvent();
  60
  61  private:
  62   double _sigmax[starlightLimits::MAXWBINS][starlightLimits::MAXYBINS];//=new double[500][500];   //decreased from 1000*1000; too big! causes fault!
  63   double _sigmaSum;
  64   double _wdelta;  //Added 7/26/07 for passing sigmadelta to pickw
  65   double _remainwd;// "
  66   int _ivalwd;     // "
  67
  68   void pickw(double &w);
  69   void picky(double &y);
  70
  71   void parentMomentum(double w,double y,double &E,double &px,double &py,double&pz);
  72   double pp(double E);
  73   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 &py2,*/double &pz2,int &iFbadevent);
  74   // void transform(double betax,double betay,double betaz,double &E,double &px,double &py,double &pz,int &iFbadevent);
  75   void thephi(double W,double px,double py,double pz,double E,double &theta,double &phi);
  76
  77   double getMass();
  78   double getWidth();
  79   double getSpin();
  80
  81   starlightConstants::particleTypeEnum _GGsingInputpidtest;
  82   int _GGsingInputnumw;
  83   int _GGsingInputnumy;
  84   double _GGsingInputGamma_em;
  85 #ifdef ENABLE_PYTHIA
  86   pythiaDecayer _pyDecayer;
  87 #endif
  88
  89 };
  90
  91
  92 #endif  // GAMMAGAMMASINGLE_H