--- /dev/null
+///////////////////////////////////////////////////////////////////////////
+//
+// Copyright 2010
+//
+// This file is part of starlight.
+//
+// starlight is free software: you can redistribute it and/or modify
+// it under the terms of the GNU General Public License as published by
+// the Free Software Foundation, either version 3 of the License, or
+// (at your option) any later version.
+//
+// starlight is distributed in the hope that it will be useful,
+// but WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with starlight. If not, see <http://www.gnu.org/licenses/>.
+//
+///////////////////////////////////////////////////////////////////////////
+//
+// File and Version Information:
+// $Rev:: $: revision of last commit
+// $Author:: $: author of last commit
+// $Date:: $: date of last commit
+//
+// Description:
+//
+//
+//
+///////////////////////////////////////////////////////////////////////////
+
+
+#ifndef GAMMAGAMMASINGLE_H
+#define GAMMAGAMMASINGLE_H
+
+
+#include <vector>
+
+#include "starlightconstants.h"
+#include "readinluminosity.h"
+#include "beambeamsystem.h"
+#include "randomgenerator.h"
+#include "eventchannel.h"
+#include "starlightconfig.h"
+#ifdef ENABLE_PYTHIA
+#include "pythiadecayer.h"
+#endif
+
+class Gammagammasingle : public eventChannel
+{
+ public:
+ Gammagammasingle(beamBeamSystem& bbsystem);
+ ~Gammagammasingle();
+
+ void singleCrossSection();
+ starlightConstants::event produceEvent(int &ievent);
+
+ upcEvent produceEvent();
+
+ private:
+ double _sigmax[starlightLimits::MAXWBINS][starlightLimits::MAXYBINS];//=new double[500][500]; //decreased from 1000*1000; too big! causes fault!
+ double _sigmaSum;
+ double _wdelta; //Added 7/26/07 for passing sigmadelta to pickw
+ double _remainwd;// "
+ int _ivalwd; // "
+
+ void pickw(double &w);
+ void picky(double &y);
+
+ void parentMomentum(double w,double y,double &E,double &px,double &py,double&pz);
+ double pp(double E);
+ 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);
+ // void transform(double betax,double betay,double betaz,double &E,double &px,double &py,double &pz,int &iFbadevent);
+ void thephi(double W,double px,double py,double pz,double E,double &theta,double &phi);
+
+ double getMass();
+ double getWidth();
+ double getSpin();
+
+ starlightConstants::particleTypeEnum _GGsingInputpidtest;
+ int _GGsingInputnumw;
+ int _GGsingInputnumy;
+ double _GGsingInputGamma_em;
+#ifdef ENABLE_PYTHIA
+ pythiaDecayer _pyDecayer;
+#endif
+
+};
+
+
+#endif // GAMMAGAMMASINGLE_H