Return the energy deosite in the gas.
136{
139
140 if (mode == 0) return e3 * xi1 / (xi1 + xi2);
141
142 static TDatabasePDG* ptrTDatabasePDG = TDatabasePDG::Instance();
143 const double mass = ptrTDatabasePDG->GetParticle(pdg)->Mass();
144 const double chrg = fabs(ptrTDatabasePDG->GetParticle(pdg)->Charge() / 3);
145
147 if (mode == 2) {
149 return e3 * me1 / me0;
150 }
152 if (mode == 1) return e3 * me1 / (me1 + me2);
153
154
155
156
157 const double mom2 = mom * mom;
158 const double beta2 = mom2 / (mass * mass + mom2);
159 const double buf = 0.1535e-3 * dx * chrg * chrg / beta2;
160 xi1 *= buf;
161 xi2 *= buf;
162 const double xi1i = 1 / xi1;
163 const double xi2i = 1 / xi2;
164 const double lm1min = -me1 * xi1i;
165 const double lm1max = (e3 - me1) * xi1i;
166 double lm1 = (e3 * me1 / (me1 + me2) - me1) * xi1i;
167 lm1 = min(lm1, lm1max);
168 lm1 = max(lm1, lm1min);
169
170 int itry = 0;
171
172
173 const double lm2b = (e3 - me1 - me2) * xi2i;
174 const double r = xi1 * xi2i;
175 double deltal = 0;
176 double lm1old = lm1;
177 while (itry <= 2 || fabs(lm1 - lm1old) > 0.001) {
178 ++itry;
179 lm1old = lm1;
180 if (itry > 50) {
181
182 break;
183 }
184 lm1 -= deltal;
185 lm1 = max(lm1min, lm1);
186 lm1 = min(lm1max, lm1);
187 double lm2 = lm2b - r * lm1;
188 double emlm1 = exp(-lm1);
189 double emlm2 = exp(-lm2);
190 double dldlm1 = 1 - emlm1 - r + r * emlm2;
191 double dl2dlm12 = emlm1 + r * r * emlm2;
192 deltal = dldlm1 / dl2dlm12;
193
194 }
195
196 if (mode == 3) return xi1 * lm1 + me1;
197
198 double lm2 = lm2b - r * lm1;
199 double lmin = 0.5 * (lm1 + lm2 + exp(-lm1) + exp(-lm2));
200 const double ymax = exp(-lmin);
201
202
203
204 const double lm1w = lm1max - lm1min;
205 double y = 10;
206 double p = 1;
207 int jtry = 0;
208 double urndm[2];
209 while (y > p) {
210 ++jtry;
211 if (jtry > 500) {
212
213 lm1 = (e3 * me1 / (me1 + me2) - me1) * xi1i;
214 break;
215 }
216 gRandom->RndmArray(2, urndm);
217 lm1 = lm1min + lm1w * urndm[0];
218 lm2 = lm2b - r * lm1;
219 double l = 0.5 * (lm1 + lm2 + exp(-lm1) + exp(-lm2));
220 p = exp(-l);
221 y = ymax * urndm[1];
222 }
223
224 return xi1 * lm1 + me1;
225}
double m_i1
I of He-C2H6 gas.
double m_i0
I of wire+gas.
double m_z1
Z of He-C2H6 gas.
double m_i2
I of "wire gas".
double m_rho0
rho of wire+gas
double m_rho1
rho of He-C2H6 gas
double m_rho2
rho of "wire gas"
double m_a1
A of He-C2H6 gas.
double m_z2
Z of "wire gas".
double m_z0
Z of wire+gas.
double m_a0
A of wire+gas.
double getMostProbabEDep(double p, double mass, double zi, double dx, double z, double a, double i, double rho) const
Return the energy deosite in the material.
double m_a2
A of "wire gas".