1 ///////////////////////////////////////////////////////////////////////////
5 // This file is part of starlight.
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.
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.
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/>.
20 ///////////////////////////////////////////////////////////////////////////
22 // File and Version Information:
23 // $Rev:: 164 $: revision of last commit
24 // $Author:: odjuvsla $: author of last commit
25 // $Date:: 2013-10-06 16:18:08 +0200 #$: date of last commit
31 ///////////////////////////////////////////////////////////////////////////
34 #ifndef GAMMAGAMMASINGLE_H
35 #define GAMMAGAMMASINGLE_H
40 #include "starlightconstants.h"
41 #include "readinluminosity.h"
42 #include "beambeamsystem.h"
43 #include "randomgenerator.h"
44 #include "eventchannel.h"
45 #include "starlightconfig.h"
47 #include "pythiadecayer.h"
50 class Gammagammasingle : public eventChannel
53 Gammagammasingle(beamBeamSystem& bbsystem);
56 void singleCrossSection();
57 starlightConstants::event produceEvent(int &ievent);
59 upcEvent produceEvent();
62 double _sigmax[starlightLimits::MAXWBINS][starlightLimits::MAXYBINS];//=new double[500][500]; //decreased from 1000*1000; too big! causes fault!
64 double _wdelta; //Added 7/26/07 for passing sigmadelta to pickw
68 void pickw(double &w);
69 void picky(double &y);
71 void parentMomentum(double w,double y,double &E,double &px,double &py,double&pz);
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);
81 starlightConstants::particleTypeEnum _GGsingInputpidtest;
84 double _GGsingInputGamma_em;
86 pythiaDecayer _pyDecayer;
92 #endif // GAMMAGAMMASINGLE_H