|
42 | 42 | generator = CobremsGenerator(E0, Epeak)
|
43 | 43 |
|
44 | 44 | def usage():
|
| 45 | + """ |
| 46 | + Print a brief summary of the internal configuration variables |
| 47 | + that regulate the behavior of the CobremsGenerator and exit. |
| 48 | + """ |
45 | 49 | print """
|
46 | 50 | Usage: cobrems.init(key=value, ...)
|
47 | 51 | where key should be one of the following, [defaults in backets]:
|
@@ -322,18 +326,48 @@ def plotPolarization(collimated=1):
|
322 | 326 | return polarH1
|
323 | 327 |
|
324 | 328 | def acceptance(vars):
|
| 329 | + """ |
| 330 | + TF1 user function that can be used to plot the collimator acceptance |
| 331 | + for photons emitted at lab polar angle vars[0] relative to the incident |
| 332 | + electron beam direction at the radiator. Both beam emittance and |
| 333 | + multiple-scattering in the target contribute to smearing of the angular |
| 334 | + acceptance at the the collimator edge. The scattering angle is contained |
| 335 | + in list argument vars[0], expressed in units of (me/fBeamEnergy). |
| 336 | + """ |
325 | 337 | return generator.Acceptance(vars[0] ** 2)
|
326 | 338 |
|
327 | 339 | def polarization(vars):
|
328 |
| - return generator.Polarization(vars[0]) |
| 340 | + """ |
| 341 | + TF1 user function that can be used to plot the linear polarization of |
| 342 | + the photon beam at energy k = vars[0] * fBeamEnergy, and production angle |
| 343 | + vars[1] expressed in units of (me/fBeamEnergy). |
| 344 | + """ |
| 345 | + return generator.Polarization(vars[0], vars[1] ** 2) |
329 | 346 |
|
330 | 347 | def dRtdx(vars):
|
| 348 | + """ |
| 349 | + TF1 user function that can be used to plot the total beam flux spectrum |
| 350 | + of the collimated coherent bremsstrahlung photon beam at photon energy |
| 351 | + k = vars[0] * fBeamEnergy, in units of (/GeV/s). |
| 352 | + """ |
331 | 353 | return generator.Rate_dNtdx(vars[0]) * cur / 1.6e-13
|
332 | 354 |
|
333 | 355 | def dRcdx(vars):
|
| 356 | + """ |
| 357 | + TF1 user function that can be used to plot the beam flux spectrum |
| 358 | + (coherent component only) of the collimated coherent bremsstrahlung |
| 359 | + photon beam at photon energy k = vars[0] * fBeamEnergy, in units |
| 360 | + of (/GeV/s). |
| 361 | + """ |
334 | 362 | return generator.Rate_dNcdx(vars[0]) * cur / 1.6e-13
|
335 | 363 |
|
336 | 364 | def dRidx(vars):
|
| 365 | + """ |
| 366 | + TF1 user function that can be used to plot the beam flux spectrum |
| 367 | + (incoherent component only) of the collimated coherent bremsstrahlung |
| 368 | + photon beam at photon energy k = vars[0] * fBeamEnergy, in units |
| 369 | + of (/GeV/s). |
| 370 | + """ |
337 | 371 | return generator.Rate_dNidx(vars[0]) * cur / 1.6e-13
|
338 | 372 |
|
339 | 373 | def plotTotal_rc(rchist, collimated=1):
|
|
0 commit comments