Line | |
---|
1 | #include "MCascade.h"
|
---|
2 |
|
---|
3 | void phys()
|
---|
4 | {
|
---|
5 | MCascade cascade;
|
---|
6 |
|
---|
7 | cascade.SetSourceZ(0.1); // Readshift of source
|
---|
8 | cascade.SetB(0/*1e-6*/); // [G] mean magnetic field
|
---|
9 | cascade.SetBradius(50); // [Mpc]
|
---|
10 | cascade.SetRuntime(8*60); // [min] maximum time to run the simulation
|
---|
11 | cascade.SetEnergyBins(18, 1e2, 1e5); // [GeV]
|
---|
12 | cascade.SetMaxInvCompton(512); // [#] maximum number of inv. Compton (-1 means infinite)
|
---|
13 | cascade.SetRatioInvCompton(0); // [%] allowed Emis of electron (0 means disabled)
|
---|
14 | cascade.SetSpectralIndex(-1); // -1 means a E^2 plot
|
---|
15 | cascade.SetDisplayIndex(1); // 1 means a E^-2 spectrum
|
---|
16 |
|
---|
17 | //
|
---|
18 | // Run the simulation: filename, eventdisply (on/off)
|
---|
19 | //
|
---|
20 | //cascade.Run("cascade_0.03_18_1e2_1e5_B1e-6_50Mpc_256_1.root", kTRUE);
|
---|
21 | cascade.Run("cascade_0.1_18_1e2_1e5_B0_512_02.root", kFALSE);
|
---|
22 | }
|
---|
23 |
|
---|
Note:
See
TracBrowser
for help on using the repository browser.