]>
Commit | Line | Data |
---|---|---|
da32329d AM |
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:: 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 | |
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); | |
75ce6a3a | 73 | void twoBodyDecay(starlightConstants::particleTypeEnum &ipid,double E,double W, |
74 | double px0,double py0,double pz0, | |
75 | double &px1,double &py1,double &pz1, double &pt1, | |
76 | double &px2,double &py2,double &pz2, double &pt2, | |
77 | double &mdec, | |
78 | int &iFbadevent); | |
da32329d AM |
79 | // void transform(double betax,double betay,double betaz,double &E,double &px,double &py,double &pz,int &iFbadevent); |
80 | void thephi(double W,double px,double py,double pz,double E,double &theta,double &phi); | |
81 | ||
82 | double getMass(); | |
83 | double getWidth(); | |
84 | double getSpin(); | |
85 | ||
86 | starlightConstants::particleTypeEnum _GGsingInputpidtest; | |
87 | int _GGsingInputnumw; | |
88 | int _GGsingInputnumy; | |
89 | double _GGsingInputGamma_em; | |
90 | #ifdef ENABLE_PYTHIA | |
91 | pythiaDecayer _pyDecayer; | |
92 | #endif | |
93 | ||
94 | }; | |
95 | ||
96 | ||
97 | #endif // GAMMAGAMMASINGLE_H |