]>
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:: $: 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 |