MC Routines
From magneticfields
To generate synthetic data, we adopt the fiducial model that the diffuse gamma ray background is isotropic.
Tashiro's MC routines in C
Here are Tashiro's MC routines in C:
Ferrer's MC Python routines
Here are Ferrer's MC routines in Python:
The script Media:mc.py generates 10000 synthetic samples with the same number of events in the real data. The events are uniformly distributed in the spherical caps and do not fall close to a Fermi source. The samples are saved to disk in numpy array format. Their statistics (average and standard deviation is computed using this script Media:qstat.py. For b>80deg we obtain the following MC samples Media:mc80.tar.gz.