]>
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: | |
45d54d9a | 23 | // $Rev:: 176 $: revision of last commit |
24 | // $Author:: jseger $: author of last commit | |
25 | // $Date:: 2014-06-20 22:15:20 +0200 #$: date of last commit | |
da32329d AM |
26 | // |
27 | // Description: | |
28 | // | |
29 | // | |
30 | // | |
31 | /////////////////////////////////////////////////////////////////////////// | |
32 | ||
33 | ||
34 | #ifndef STARLIGHTCONSTANTS_H_INCLUDE | |
35 | #define STARLIGHTCONSTANTS_H_INCLUDE | |
36 | ||
37 | ||
38 | /* | |
39 | * Constants are set here | |
40 | */ | |
41 | namespace starlightConstants | |
42 | { | |
43 | ||
44 | ||
45 | // constants | |
46 | static const double hbarc = 0.197327053; | |
47 | static const double hbarcmev = hbarc*1000.; | |
48 | static const double pi = 3.141592654; | |
49 | static const double twoPi = 2 * pi; | |
50 | static const double alpha = 1/137.0359895; | |
51 | ||
45d54d9a | 52 | // deuteron slope parameter |
53 | static const double deuteronSlopePar = 9.5; // [(GeV/c)^{-2}] | |
da32329d AM |
54 | // particle masses |
55 | static const double protonMass = 0.93827231; // [GeV/c^2] | |
56 | static const double pionChargedMass = 0.13956995; // [GeV/c^2] | |
57 | static const double pionNeutralMass = 0.1349766; // [GeV/c^2] | |
58 | static const double kaonChargedMass = 0.493677; // [GeV/c^2] | |
59 | static const double mel = 0.00051099907; // [GeV/c^2] | |
60 | static const double muonMass = 0.105658389; // [GeV/c^2] | |
61 | static const double tauMass = 1.777; // [GeV/c^2] | |
62 | ||
63 | static const double f0Mass = 0.980; // [GeV/c^2] | |
64 | static const double etaMass = 0.54730; // [GeV/c^2] | |
65 | static const double etaPrimeMass = 0.95766; // [GeV/c^2] | |
66 | static const double etaCMass = 2.9803; // [GeV/c^2] | |
67 | static const double f2Mass = 1.2751; // [GeV/c^2] | |
68 | static const double a2Mass = 1.3183; // [GeV/c^2] | |
69 | static const double f2PrimeMass = 1.525; // [GeV/c^2] | |
70 | ||
71 | enum particleTypeEnum { | |
72 | UNKNOWN = 0, | |
73 | ELECTRON = 11, | |
74 | MUON = 13, | |
75 | TAUON = 15, | |
45d54d9a | 76 | TAUONDECAY = 10015, |
da32329d AM |
77 | PROTON = 212, |
78 | PION = 211, | |
79 | KAONCHARGE = 321, | |
80 | KAONNEUTRAL = 310, | |
81 | A2 = 115, | |
82 | ETA = 221, | |
83 | F2 = 225, | |
84 | ETAPRIME = 331, | |
85 | F2PRIME = 335, | |
86 | ETAC = 441, | |
87 | F0 = 9010221, | |
88 | ZOVERZ03 = 33, | |
89 | RHO = 113, | |
90 | RHOZEUS = 913, | |
91 | FOURPRONG = 999, | |
92 | OMEGA = 223, | |
93 | PHI = 333, | |
94 | JPSI = 443, | |
95 | JPSI_ee = 443011, | |
96 | JPSI_mumu = 443013, | |
97 | JPSI2S = 444, | |
98 | JPSI2S_ee = 444011, | |
99 | JPSI2S_mumu = 444013, | |
100 | UPSILON = 553, | |
101 | UPSILON_ee = 553011, | |
102 | UPSILON_mumu = 553013, | |
103 | UPSILON2S = 554, | |
104 | UPSILON2S_ee = 554011, | |
105 | UPSILON2S_mumu = 554013, | |
106 | UPSILON3S = 555, | |
107 | UPSILON3S_ee = 555011, | |
108 | UPSILON3S_mumu = 555013 | |
109 | }; | |
110 | ||
111 | enum decayTypeEnum { | |
112 | NOTKNOWN = 0, | |
113 | NARROWVMDEFAULT = 1, | |
114 | WIDEVMDEFAULT = 2, | |
115 | PSIFAMILY = 3, | |
116 | LEPTONPAIR = 4, | |
117 | SINGLEMESON = 5 | |
118 | }; | |
119 | ||
120 | enum interactionTypeEnum { | |
121 | UNSPECIFIED = 0, | |
122 | PHOTONPHOTON = 1, | |
123 | PHOTONPOMERONNARROW = 2, | |
124 | PHOTONPOMERONWIDE = 3, | |
125 | PHOTONPOMERONINCOHERENT = 4, | |
126 | PHOTONUCLEARSINGLE = 5, | |
127 | PHOTONUCLEARDOUBLE = 6, | |
128 | PHOTONUCLEARSINGLEPA = 7, | |
129 | PHOTONUCLEARSINGLEPAPY = 8 | |
130 | ||
131 | }; | |
132 | //Structure for each event's set of tracks. | |
133 | struct event{ | |
134 | ||
135 | public: | |
136 | ||
137 | int _numberOfTracks; | |
138 | //Right now this is set up for a maximum of 4 tracks,if we want more, just increase the arrays | |
139 | //Moved it to 30, this way when pythia returns, it wont complain too much...hedging bets...not sure how many | |
140 | double px[30],py[30],pz[30]; | |
141 | //starlightConstants::particleTypeEnum _fsParticle[30]; | |
142 | int _fsParticle[30]; | |
143 | int _charge[30]; | |
144 | //To help track mothers and daughters produced through pythia. | |
145 | int _mother1[30]; | |
146 | int _mother2[30]; | |
147 | int _daughter1[30]; | |
148 | int _daughter2[30]; | |
149 | //Normally we just set vertices to 0 | |
150 | //But for pythia, we decay additional states | |
151 | int _numberOfVertices; | |
152 | double _vertx[10],_verty[10],_vertz[10]; | |
153 | }; | |
154 | ||
155 | ||
156 | } // starlightConstants | |
157 | ||
158 | ||
159 | #endif // STARLIGHTCONSTANTS_H_INCLUDE | |
160 |