Select Git revision
galplotdef_55_reltest19~
galplotdef_55_reltest19~ 40.60 KiB
1234567890123456789012
======================value
Title = galplot test case
E_gamma_min = 1.e0 min gamma-ray energy (MeV)
E_gamma_max = 1.e8 max gamma-ray energy (MeV)
E_gamma_factor = 2.0 gamma-ray energy factor
nu_synch_min = 1.0e6 min synchrotron frequency (Hz)
nu_synch_max = 1.0e10 max synchrotron frequency (Hz)
nu_synch_factor = 2.0 synchrotron frequency factor
NOT Used for profiles:
long_min =330.5 gamma-ray intensity skymap longitude minimum (deg)
long_max = 30.5 gamma-ray intensity skymap longitude maximum (deg)
lat_min =-89.75 gamma-ray intensity skymap latitude minimum (deg)
lat_max =+89.75 gamma-ray intensity skymap latitude maximum (deg)
d_long = 10. gamma-ray intensity skymap longitude binsize (deg) << used for Xco fitting only
d_lat = 2. gamma-ray intensity skymap latitude binsize (deg) << used for Xco fitting only
verbose = 0 -4000 -1002 verbosity: 0=min,10=max <0: selected debugs -1002= no curve labelling -4000=population synthesis only
test_suite = 0 run test suite instead of normal run
screen_output = 0 0=no screen output 1=screen output
psfile_tag = P7.7_37m string appended to plot files
output_format = 3 0=none 1=ps 2=gif 3=both
gamma_spectrum = 1 2 1 1=plot model & data, 2=model only 3=isotropic only 4=data only
gamma_long_profile = 1 0 1
gamma_lat_profile = 1 0 1
gamma_IC_selectcomp = 1 for spectra: 0=plot total only, 1=plot components and total ; for profiles: 0=total, 1,2,3=use one component only (for healpix:always total)
Xlong_min1 = 0.5 region X
Xlong_max1 = 10.5
Xlong_min2 = 350.5
Xlong_max2 = 359.5
Xlat_min1 = -05.5
Xlat_max1 = -00.5
Xlat_min2 = +00.5
Xlat_max2 = +05.5
Ylong_min1 = 10.5 region Y
Ylong_max1 = 30.5
Ylong_min2 = 330.5
Ylong_max2 = 350.5
Ylat_min1 = -05.5
Ylat_max1 = -00.5
Ylat_min2 = +00.5
Ylat_max2 = +05.5
test for Olaf
xlong_min1 = 0.25
xlong_max1 = 19.75
xlong_min2 = 340.25
xlong_max2 = 359.75
xlat_min1 = -05.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +05.75
test deBoer small region
xlong_min1 = 12.25
xlong_max1 = 14.75
xlong_min2 = 15.25
xlong_max2 = 19.75
xlat_min1 = -80.75
xlat_max1 = -20.25
xlat_min2 = +20.25
xlat_max2 = +80.75
xAlong_min1 = 0.5
xAlong_max1 = 30.5
xAlong_min2 = 330.5
xAlong_max2 = 359.5
xAlat_min1 = -05.5
xAlat_max1 = -00.5
xAlat_min2 = +00.5
xAlat_max2 = +05.5
region A
xlong_min1 = 0.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -04.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +04.75
region A rationalized
xlong_min1 = 0.00
xlong_max1 = 30.00
xlong_min2 = 330.00
xlong_max2 = 360.00
xlat_min1 = -05.00
xlat_max1 = -00.00
xlat_min2 = +00.00
xlat_max2 = +05.00
Blong_min1 = 30.5
Blong_max1 = 179.5
Blong_min2 = 180.5
Blong_max2 = 330.5
Blat_min1 = -05.5
Blat_max1 = -00.5
Blat_min2 = +00.5
Blat_max2 = +05.5
region C
Clong_min1 = 90.5
Clong_max1 = 179.5
Clong_min2 = 180.5
Clong_max2 = 270.5
Clat_min1 = -10.5
Clat_max1 = -00.5
Clat_min2 = +00.5
Clat_max2 = +10.5
outer Galaxy avoiding Crab, Geminga, Velaxlong_min1 = 90.25
xlong_max1 = 175.25
xlong_min2 = 205.25
xlong_max2 = 250.25
xlat_min1 = -10.25
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +10.25
region D all longitudes 10<|b|<20
Dlong_min1 = 0.5
Dlong_max1 = 179.5
Dlong_min2 = 180.5
Dlong_max2 = 359.5
Dlat_min1 = -20.5
Dlat_max1 = -10.5
Dlat_min2 = +10.5
Dlat_max2 = +20.5
region D more precisely defined
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -19.75
xlat_max1 = -10.25
xlat_min2 = +10.25
xlat_max2 = +19.75
all longitudes -20<b<-10
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -19.75
xlat_max1 = -10.25
xlat_min2 = -19.75
xlat_max2 = -10.25
all longitudes +10<b<+20
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = +10.25
xlat_max1 = +19.75
xlat_min2 = +10.25
xlat_max2 = +19.75
all longitudes 15<|b|<30
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -29.75
xlat_max1 = -15.25
xlat_min2 = +15.25
xlat_max2 = +29.75
all longitudes -30<b<-15
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -29.75
xlat_max1 = -15.25
xlat_min2 = -29.75
xlat_max2 = -15.25
all longitudes +15<b<+30
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = +15.25
xlat_max1 = +29.75
xlat_min2 = +15.25
xlat_max2 = +29.75
Elong_min1 = 0.5
Elong_max1 = 179.5
Elong_min2 = 180.5
Elong_max2 = 359.5
Elat_min1 = -60.5
Elat_max1 = -20.5
Elat_min2 = +20.5
Elat_max2 = +60.5
Region F
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -89.75
xlat_max1 = -60.25
xlat_min2 = +60.25
xlat_max2 = +89.75
Region F rationalized
xlong_min1 = 0.00
xlong_max1 = 180.00
xlong_min2 = 180.00
xlong_max2 = 360.00
xlat_min1 = -90.00
xlat_max1 = -60.00
xlat_min2 = +60.00
xlat_max2 = +90.00
Region F South
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -89.75
xlat_max1 = -60.25
xlat_min2 = -89.75
xlat_max2 = -60.25
Region F North
xlong_min1 = 0.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = +60.25
xlat_max1 = +89.75
xlat_min2 = +60.25
xlat_max2 = +89.75
xlong_min1 = 00.5
xlong_max1 = 30.5
xlong_min2 = 330.5
xlong_max2 = 359.5
xlat_min1 = -89.5
xlat_max1 = -60.5
xlat_min2 = +60.5
xlat_max2 = +89.5
Region H
xlong_min1 = 00.25
xlong_max1 = 59.75
xlong_min2 = 300.25
xlong_max2 = 359.75
xlat_min1 = -9.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +9.75
Region H rationalized ========================================
long_min1 = 00.00
long_max1 = 60.00
long_min2 = 300.00
long_max2 = 360.00
lat_min1 = -10.00
lat_max1 = -0.00
lat_min2 = +0.00
lat_max2 = +10.00
H but 5< |b| <10
xlong_min1 = 00.25
xlong_max1 = 59.75
xlong_min2 = 300.25
xlong_max2 = 359.75
xlat_min1 = -9.75
xlat_max1 = -5.25
xlat_min2 = +5.25
xlat_max2 = +9.75
H but +10<b<+20
xlong_min1 = 00.25
xlong_max1 = 59.75
xlong_min2 = 300.25
xlong_max2 = 359.75
xlat_min1 = +10.25
xlat_max1 = +19.75
xlat_min2 = +10.25
xlat_max2 = +19.75
H but -20<b<-10
xlong_min1 = 00.25
xlong_max1 = 59.75
xlong_min2 = 300.25
xlong_max2 = 359.75
xlat_min1 = -19.75
xlat_max1 = -10.25
xlat_min2 = -19.75
xlat_max2 = -10.75
all longitudes -60<b<-20
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -59.75
xlat_max1 = -20.25
xlat_min2 = -59.75
xlat_max2 = -20.75
all longitudes +20<b<+60
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = +20.25
xlat_max1 = +59.75
xlat_min2 = +20.25
xlat_max2 = +59.75
longitudes 30-90 |b|<5
xlong_min1 = 29.75
xlong_max1 = 89.75
xlong_min2 = 29.25
xlong_max2 = 89.75
xlat_min1 = -04.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +04.75
longitudes 270-330 |b|<5
xlong_min1 = 270.25
xlong_max1 = 329.75
xlong_min2 = 270.25
xlong_max2 = 329.75
xlat_min1 = -04.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +04.75
270<l<90 -30<b<-15 for Troy
xlong_min1 = 270.25
xlong_max1 = 359.75
xlong_min2 = 00.25
xlong_max2 = 89.75
xlat_min1 = -29.75
xlat_max1 = -15.25
xlat_min2 = -29.75
xlat_max2 = -15.75
270<l<90 +15<b<+30 for Troy
xlong_min1 = 270.25
xlong_max1 = 359.75
xlong_min2 = 00.25
xlong_max2 = 89.75
xlat_min1 = +15.25
xlat_max1 = +29.75
xlat_min2 = +15.25
xlat_max2 = +29.75
90<l<270 30<|b|<15 for Troy
xlong_min1 = 90.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 269.75
xlat_min1 = -29.75
xlat_max1 = -15.25
xlat_min2 = +15.25
xlat_max2 = +29.75
90<l<270 -30<b<-15 for Troy
xlong_min1 = 90.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 269.75
xlat_min1 = -29.75
xlat_max1 = -15.25
xlat_min2 = -29.75
xlat_max2 = -15.75
90<l<270 +15<b<+30 for Troy
xlong_min1 = 90.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 269.75
xlat_min1 = +15.25
xlat_max1 = +29.75
xlat_min2 = +15.25
xlat_max2 = +29.75
xlong_min1 = 00.5 FIT REGION HIGH LATITUDES ==========
xlong_max1 = 179.5
xlong_min2 = 180.5
xlong_max2 = 359.5
xlat_min1 = -80.5
xlat_max1 = -10.5
xlat_min2 = +10.5xlat_max2 = +80.5
xlong_min1 = 00.5 FIT REGION FULL SKY
xlong_max1 = 179.5
xlong_min2 = 180.5
xlong_max2 = 359.5
xlat_min1 = -80.5
xlat_max1 = -00.5
xlat_min2 = +00.5
xlat_max2 = +80.5
xlong_min1 = 00.5 FIT REGION FULL SKY |b|<30
xlong_max1 = 179.5
xlong_min2 = 180.5
xlong_max2 = 359.5
xlat_min1 = -29.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +29.75
xlong_min1 = 00.5 FIT REGION FULL SKY |b|<10
xlong_max1 = 179.5
xlong_min2 = 180.5
xlong_max2 = 359.5
xlat_min1 = -09.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +09.75
xlong_min1 = 00.5 FIT REGION FULL SKY |b|<5
xlong_max1 = 179.5
xlong_min2 = 180.5
xlong_max2 = 359.5
xlat_min1 = -04.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +04.75
xlong_min1 = 160.25 FIT REGION SMALL ANTICENTRE
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 199.75
xlat_min1 = -09.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +09.75
xlong_min1 = 183.00 REGION WHERE SUN CUTS ANTICENTRE around l=186
xlong_max1 = 186.00
xlong_min2 = 186.00
xlong_max2 = 189.00
xlat_min1 = -03.00
xlat_max1 = -00.00
xlat_min2 = +00.00
xlat_max2 = +03.00
xlong_min1 = 197.00 REGION OUT OF PLANE WHERE SUN PASSES
xlong_max1 = 200.00
xlong_min2 = 200.00
xlong_max2 = 203.00
xlat_min1 = +19.00xlat_max1 = +22.00
xlat_min2 = +22.00
xlat_max2 = +25.00
SPI regions
xlong_min1 = 00.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -9.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +9.75
xlong_min1 = 00.25
xlong_max1 = 09.75
xlong_min2 = 350.25
xlong_max2 = 359.75
xlat_min1 = -9.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +9.75
xlong_min1 = 09.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 350.75
xlat_min1 = -9.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +9.75
Bouchet 2008
xlong_min1 = 00.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -14.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +14.75
xlong_min1 = 00.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -09.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +09.75
rationalized version
xlong_min1 = 00.00
xlong_max1 = 30.00
xlong_min2 = 330.00
xlong_max2 = 360.00
xlat_min1 = -10.00
xlat_max1 = -0.00
xlat_min2 = +0.00
xlat_max2 = +10.00
rationalized version =============================================================
xlong_min1 = 00.00xlong_max1 = 30.00
xlong_min2 = 330.00
xlong_max2 = 360.00
xlat_min1 = -15.00
xlat_max1 = -0.00
xlat_min2 = +0.00
xlat_max2 = +15.00
xlong_min1 = 00.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -04.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +04.75
AO6-KP proposal 1
xlong_min1 = 00.25
xlong_max1 = 29.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -80.75
xlat_max1 = -44.25
xlat_min2 = +44.25
xlat_max2 = +80.75
AO6-KP proposal 2
xlong_min1 = 60.25
xlong_max1 = 89.75
xlong_min2 = 270.25
xlong_max2 = 299.75
xlat_min1 = -09.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +09.75
reference region for Petry spimodfit flux
xlong_min1 = 00.25
xlong_max1 = 89.75
xlong_min2 = 270.25
xlong_max2 = 359.75
xlat_min1 = -25.25
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +25.25
MILAGRO PRL paper:
xlong_min1 = 40.25
xlong_max1 = 99.75
xlong_min2 = 40.25
xlong_max2 = 99.75
xlat_min1 = -04.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +04.75
MILAGRO Cygnus:
xlong_min1 = 65.25
xlong_max1 = 75.75
xlong_min2 = 76.25
xlong_max2 = 84.75
xlat_min1 = -02.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +02.75
MILAGRO inner Galaxy
xlong_min1 = 30.25
xlong_max1 = 49.75
xlong_min2 = 50.25
xlong_max2 = 64.75
xlat_min1 = -02.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +02.75
Cygnus:
xlong_min1 = 70.00
xlong_max1 = 75.25
xlong_min2 = 75.75
xlong_max2 = 80.00
xlat_min1 = -05.00
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +05.00
GC 20 by 10 deg
xlong_min1 = 00.25
xlong_max1 = 09.75
xlong_min2 = 350.25
xlong_max2 = 359.75
xlat_min1 = -4.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +4.75
GC 10 by 5 deg
xlong_min1 = 00.25
xlong_max1 = 04.75
xlong_min2 = 355.25
xlong_max2 = 359.75
xlat_min1 = -2.75
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +2.75
GC for Dogiel
xlong_min1 = 00.25
xlong_max1 = 02.25
xlong_min2 = 357.75
xlong_max2 = 359.75
xlat_min1 = -2.25
xlat_max1 = -0.25
xlat_min2 = +0.25
xlat_max2 = +2.25
Chandra
xlong_min1 = 27.75
xlong_max1 = 28.25
xlong_min2 = 28.25
xlong_max2 = 28.75
xlat_min1 = -1.25
xlat_max1 = +0.25
xlat_min2 = +0.25
xlat_max2 = +1.25
inner South
xlong_min1 = 330.25
xlong_max1 = 359.75
xlong_min2 = 330.25
xlong_max2 = 359.75
xlat_min1 = -1.25
xlat_max1 = +0.25
xlat_min2 = +0.25
xlat_max2 = +1.25
synchrotron whole sky
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -80.25
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +80.25
synchrotron inner low lats
xlong_min1 = 00.25
xlong_max1 = 59.75
xlong_min2 = 300.25
xlong_max2 = 359.75
xlat_min1 = -00.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +00.75
synchrotron
xlong_min1 = 09.25
xlong_max1 = 49.75
xlong_min2 = 310.25
xlong_max2 = 349.75
xlat_min1 = -29.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +29.75
synchrotron North
xlong_min1 = 19.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 180.75
xlat_min1 = +00.25
xlat_max1 = +59.75
xlat_min2 = +00.25
xlat_max2 = +59.75
synchrotron North low lats
xlong_min1 = 19.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 240.75
xlat_min1 = -02.75
xlat_max1 = -00.25xlat_min2 = +00.25
xlat_max2 = +02.75
synchrotron North low lats wide
xlong_min1 = 19.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 240.75
xlat_min1 = -19.75
xlat_max1 = -00.25
xlat_min2 = +00.25
xlat_max2 = +19.75
synchrotron North high lats
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = +60.25
xlat_max1 = +80.75
xlat_min2 = +60.25
xlat_max2 = +80.75
synchrotron North lats 20-60 long 0-240
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 240.75
xlat_min1 = +20.25
xlat_max1 = +60.75
xlat_min2 = +20.25
xlat_max2 = +60.75
synchrotron North lats 20-60 long 0-360
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = +20.25
xlat_max1 = +60.75
xlat_min2 = +20.25
xlat_max2 = +60.75
synchrotron South lats -20/-60 long 0-360
xlong_min1 = 00.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 359.75
xlat_min1 = -59.75
xlat_max1 = -19.75
xlat_min2 = +19.75
xlat_max2 = +59.75
synchrotron inner low lats
xlong_min1 = 00.25
xlong_max1 = 59.75
xlong_min2 = 300.25
xlong_max2 = 359.75
xlat_min1 = -19.75
xlat_max1 = -00.25xlat_min2 = +00.25
xlat_max2 = +19.75
synchrotron
xlong_min1 = 00.25
xlong_max1 = 04.75
xlong_min2 = 355.25
xlong_max2 = 359.75
xlat_min1 = -30.75
xlat_max1 = -20.25
xlat_min2 = +20.25
xlat_max2 = +30.75
synchrotron
xlong_min1 = 170.25
xlong_max1 = 179.75
xlong_min2 = 180.25
xlong_max2 = 189.75
xlat_min1 = -30.75
xlat_max1 = -20.25
xlat_min2 = +20.25
xlat_max2 = +30.75
long_binsize = 2.0 binsize for longitude profiles (degrees)
long_log_scale = 0 longitude profile intensity scale: 0=linear 1=log
lat_binsize = 1.0 binsize for latitude profiles (degrees)
lat_log_scale = 1 latitude profile intensity scale: 0=linear 1=log
# for gammas
gamma_spectrum_Emin = 1 1.0e-2 1.0e1 1.0e-2 1.0e+2 minimum energy on spectrum plot 0.0-> all energies
gamma_spectrum_Emax = 5.0e5 1.0e6 40 1.0e+5 40 1.0e+5 1.0e6 1.0e+7 0.6e+5 2.0e+5 maximum energy on spectrum plot 0.0-> all energies
gamma_spectrum_Imin = 1.0e-4 1.0e-6 1.0e-4 1.0e-6 1.0e-4 1.0e-5 minimum intensity*E^2 on spectrum plot
gamma_spectrum_Imax = 1.0 2.0e-2 1.0e-0 2.0e-2 maximum intensity*E^2 on spectrum plot
# for MW
#gamma_spectrum_Emin = 1.0e-15 1.0e+2 1.0e-3 1.0e+2 1.0e+2 1.0e-2 minimum energy on spectrum plot 0.0-> all energies
#gamma_spectrum_Emax = 1.0e+6 1.0e5 1.0e+7 0.6e+5 2.0e+5 maximum energy on spectrum plot 0.0-> all energies
#gamma_spectrum_Imin = 1.0e-6 1.0e-5 minimum intensity*E^2 on spectrum plot
#gamma_spectrum_Imax = 1.0e+4 2.0e-1 2.0e-2 maximum intensity*E^2 on spectrum plot
sources_EGRET = 0 200 0=sources removed 1=sources in 28,29 special data 100,101=GLAST DC2 simulations 200=GLAST SC2
energies_EGRET = 13 13 number of EGRET energies (10=30-10000 MeV,13= >10000 MeV)
convolve_EGRET = 0 2 0=no EGRET convolution, 1/2/3=convolution (without/with/with for ratios>1) renormalization
convolve_EGRET_iEmax = 7 maximum energy range to convolve EGRET (1-10)
convolve_EGRET_thmax = 20.0 maximum angle (deg) for EGRET convolution
error_bar_EGRET = 0.15 EGRET error bar, as fraction
spectrum_style_EGRET = 2 1=segments 2=at mean energy 3=both
spectrum_50GeV_EGRET = 0 1=plot EGRET spectra in 10-50 GeV from Hunter et al. 1997 (300<l<60, |b|<10 only)
spectrum_index_EGRET = 2.0 0=automatic index computation, otherwise use this index to plot EGRET points
spectrum_cut_total = 0 1=cut total prediction outside EGRET energy range
isotropic_use = 1 0=ignore isotropic, 1=add to prediction
isotropic_type = 2 1=power law, 2=explicit EGRET ranges only 1 will work for GLAST at present. in fact 2=read from isotropic file defined below
isotropic_const = 4.0e-3 1.55e-3 gardian fit 2.74e-3 4.33e-3 5.496e-3 2.91e-3 1.0e-3 sree=2.74e-3 isotropic background = isotropic_const * E^(-isotropic_g)
isotropic_g = 2.2 2.095 gardian fit 2.247 2.300 2.1925 2.1 cm^-2 sr^-1 MeV^-1 s^-1 fitted_isotropic_spectrum_500030
isotropic_EGRET = 1.68e-5,1.06e-5,6.66e-6,4.48e-6,3.92e-6,1.20e-6,7.6e-7,3.20e-7,2.4e-7,1.17e-7,3.4e-8,1.1e-8,0.0 isotropic background in EGRET ranges, cm^-2 s^-1
isotropic_sree_plot = 0 (2) 1=plot Sreekumar et al. isotropic background separately 2=plot isotropic_EGRET 3=both
fit_EGRET = 0 1=fit EGRET data components+EGB, 2=Xco
fit_nrebin_EGRET = 1 20 2880 number of 0.5*0.5 bins for EGRET fitting rebinning eg 2880
fit_counts_min_EGRET = 0 5 10 minimum counts for EGRET fitting after rebinningfit_function_EGRET = 3 1=background only, 2=scaling+background 3=both
fit_options_EGRET = 0 other options 1=write convolved EGRET skymaps to files
convolve_GLAST = 0 1 0 0=no convolution, 1=convolution
data_EGRET = 0 1=plot EGRET data
data_GLAST = 101 >=1:plot GLAST data 100,101=GLAST DC2 simulations 200=GLAST SC2 201=GLAST SC2 13 energies (absent)
# dataclean 10m
# for profiles healpix order from counts, need same order as for exposure hence o6 but has to be same as galprop maps so this should work gives nans in Exposure convert. does not work with galprop o7 maps
xGLAST_counts_file = GLAST/dataclean/counts_10m_dataclean_z100_t65_cl3_healpix_o6_EBOUNDS.fits
xGLAST_exposure_file = GLAST/dataclean/expcube_10m_dataclean_z100_t65_cl3_p6_v3_dataclean_healpix_o6.fits
xGLAST_psf_file = GLAST/dataclean/psf_10m_dataclean_z100_t65_cl3_p6_v3_dataclean.fits
# this was normally used:
xLAST_counts_file = GLAST/dataclean/16m/maps/counts_16m_dataclean_z100_t65_cl3_healpix_o7_gardian_30bins.fits
# but since healpix order is based on counts for profiles and expcube o7 is not available, use o6 for profiles: but the 30bin gardian counts are not available
xLAST_exposure_file = GLAST/dataclean/16m/irfs/expcube_16m_datacl_z100_t65_cl3_p6_v3_dataclean_healpix_o6.fits
xLAST_psf_file = GLAST/dataclean/16m/irfs/psf_16m_datacl_z100_t65_cl3_p6_v3_dataclean.fits
# pass7.4
#GLAST_counts_file = GLAST/24m_pass7.4/24m_pass7.4_clean_z90_t72_cl0/maps/counts_24m_pass7.4_clean_z90_t72_cl0_healpix_o7_EBOUNDS.fits
#GLAST_exposure_file = GLAST/24m_pass7.4/24m_pass7.4_clean_z90_t72_cl0/irfs/expcube_24m_pass7.4_clean_z90_t72_cl0_p7_v4_clean_healpix_o6.fits
#GLAST_psf_file = GLAST/24m_pass7.4/24m_pass7.4_clean_z90_t72_cl0/irfs/psf_24m_pass7.4_clean_z90_t72_cl0_p7_v4_clean.fits
# 37m pass 7.7
xxGLAST_counts_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx/maps/counts_37m_pass7.7_clean_z90_t72_ec1xxx_healpix_o7_EBOUNDS.fits
xxGLAST_exposure_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx/irfs/expcube_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v7_healpix_o6.fits.gz
xxGLAST_psf_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx/irfs/psf_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v7.fits
# 37m pass 7.7 down to 10 MeV FRONT+BACK
GLAST_counts_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/maps/counts_37m_pass7.7_clean_z90_t72_ec1xxx_emin10_healpix_o7_columnnames.fits
GLAST_exposure_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/irfs/expcube_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v6_healpix_o6.fits.gz
GLAST_psf_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/irfs/psf_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v6.fits
# 37m pass 7.7 down to 10 MeV FRONT
xGLAST_counts_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/maps/counts_37m_pass7.7_clean_z90_t72_ec1xxx_emin10_front_healpix_o7_columnnames.fits
xGLAST_exposure_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/irfs/expcube_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v6_front_healpix_o6.fits
xGLAST_psf_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/irfs/psf_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v6_front.fits
# 37m pass 7.7 down to 10 MeV BACK last used 27 Aug 2013 OK
xGLAST_counts_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/maps/counts_37m_pass7.7_clean_z90_t72_ec1xxx_emin10_back_healpix_o7_columnnames.fits
xGLAST_exposure_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/irfs/expcube_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v6_back_healpix_o6.fits
xGLAST_psf_file = GLAST/37m_pass7.7_clean_z90_t72_ec1xxx_emin10/irfs/psf_37m_pass7.7_clean_z90_t72_ec1xxx_p7clean_v6_back.fits
# P8 1st year test. Generated by Johann and Gulli
xGLAST_counts_file = GLAST/P300x_10MeV_1TeV_zmax100_evclass4_mktime_healpix7_counts.fits
xGLAST_exposure_file = GLAST/expcube_P300x_10MeV_1TeV_zmax100_evclass4_mktime_healpix7.fits
xGLAST_psf_file = GLAST/psf_P300x_10MeV_1TeV_zmax100_evclass4_mktime.fits
# this one looks much to low, does not fit at high lats.
xisotropic_bgd_file = GLAST/isotropic_24month_P74_clean_jeanmarc.fits
# the new EGB from Markus looks quite similar to the old one so still use it:
# isotropic_dataclean_3.fits same as isotropic_dataclean_2.fits but copied 57.5 MeV to new 25 MeV
isotropic_bgd_file = GLAST/isotropic_dataclean_3.fits GLAST/isotropic_dataclean_2.fits GLAST/isotropic_dataclean_1.fits
Fermi_cat_file = GLAST/gll_psc_v08_NoNullFluxDensity.fit official 2FGL, 2 NULLs elimnated by AWS (blue points)
Fermi_cat_file2 = GLAST/gll_psc_v08_NoNullFluxDensity_NoExtragalactic.fit 2FGL with all extragalactic sources removed (red points)
Fermi_cat_file3 = GLAST/gll_psc_v08_NoNullFluxDensity_NoExtragalactic_NoUnassoc.fit 2FGL with all extragalactic and unassociated sources removed (black points)
data_COMPTEL = 2 0 2 >=1=plot COMPTEL data (1,2,3 different styles)
data_OSSE = 0 1=plot OSSE data
data_RXTE = 0 1=plot RXTE data
data_INTEGRAL = -2012 -2011 -100 -224 233 232 224 N=plot INTEGRAL SPI data list N (>0) -100=Bouchet -2011=Bouchet new |b|<15, -2012=|b|<10
data_INTEGRAL_comp = 01,09,-1,150 INTEGRAL SPI components range e.g. 1,3,4,95 11,102 -ve suppresses
data_INTEGRAL_mode = 0 INTEGRAL SPI 1=data 2=fit 3=both 0=none (but plot other data)
data_INTEGRAL_syserr = 0.150 INTEGRAL SPI systematic error (fraction of total) 0.25
data_INTEGRAL_sysemn = 0.100 INTEGRAL SPI systematic error maximum! energy (MeV) to apply 0.100
data_IBIS = 0 1=plot IBIS data
data_MILAGRO = 6 1-6 =plot MILAGRO data new 2007: 5=Cyg 6=30<l<65
data_HEGRA = 0 1=plot HEGRA data
data_WHIPPLE = 0 1=plot WHIPPLE data
model_ridge = 0 >=1=plot ridge models (1,2)
spiskymax_profile = 0 1=plot SPI spiskmax profiles
spiskymax_image_ID = 2610 spiskymax_image.IDIDID 2200 2612 2611 2610
spiskymax_iteration = 40
# population synthesis parameters
sourcepop1_verbose = 0 1
sourcepop1_density0 = 10
#sourcepop1_oversampl = 1 for realistic catalogue
sourcepop1_oversampl = 10 100 10 1 .1 .001 .01 .1 100 1 10 100 1 10 100 100 use small values for dim numerous sources
sourcepop1_L_min = 6e34 6e33 6e34 6e33 6e32 6e31 6e30 6e31 6e34 6e33 6e32 6e33 6e34 1e35 1e36 1e37
sourcepop1_L_max = 6e37 6e36 6e37 6e36 6e35 6e34 6e33 6e34 6e37 6e36 6e35 6e36 6e37 1e38 1e39 1e39
sourcepop1_alpha_L = -1.01 -2.0 avoid -1.0 exactly
sourcepop1_fluxlimit = 5e-9 1e-8 5e-9 5e-9 1e-9 1e-8 1e-7 0 gives empty NS plot!
sourcepop1_alpha_R = 2.35 0.0 0.0 pulsars:2.35 rho= (R/Ro)^alpha_R * exp(- beta_R*(R-Ro)/Ro) * exp(-|z|/zscale)
sourcepop1_beta_R = 5.56283 1.70 1.0625
sourcepop1_zscale = 0.5 1.0 0.5 0.18 0.5 0.1
sourcepop1_alpha_z = 0.0
sourcepop1_specmodel = 3 1 2 3 2 1=2 breaks, fixed 2=2 breaks with uncorr dispersion, 3=cutoff powerlaw with dispersion
sourcepop1_spectrumg = 1.5 2.01 1.8 OBSOLETE
sourcepop1_spect_g_0 = 1.40 mean index for correlated dispersion, spectral model 3
sourcepop1_spect_br0 = 3.32 mean log10 Ecutoff for correlated dispersion, spectral model 3
sourcepop1_sp_g0_sig = 0.43 sigma index for correlated dispersion, spectral model 3
sourcepop1_spbr0_sig = 0.16 sigma log10 Ecutoff for correlated dispersion, spectral model 3
sourcepop1_spect_g_0 = 1.4 1.8 2.2 2.0 1.5 1.5 spectral index below break 0 spectral model 2
sourcepop1_spect_br0 = 2089 2e3 1.0 2e3 break 0 (MeV)
sourcepop1_spect_g_1 = 3.0 2.2 2.0 3.0 spectral index from break 0 to break 1
sourcepop1_spect_br1 = 2e4 break 1 (MeV)
sourcepop1_spect_g_2 = 3.5 2.2 2.0 3.0 3.5 spectral index above break 1
sourcepop1_sp_g0_sig = 0.34 sigma in g_0 for break with dispersion spectral model 2
sourcepop1_spbr0_sig = 0.29 sigma in log10(br0) for break with dispersion spectral model 2
sourcepop1_sp_g1_sig = 0.0 sigma in g_1 for break with dispersion
sourcepop1_spbr1_sig = 0 sigma in br1 for break with dispersion
sourcepop1_sp_g2_sig = 0 sigma in g_2 for break with dispersion
sourcepop1_E_ref_low = 1e3 50 1e2 1e3 1e2 1e3 100 .018 100 reference energy range low MeV
sourcepop1_E_ref_hig = 2e4 10000 .028 10000 reference energy range high MeV
sourcepop_SPI_cat = catalogue.spimodfit.2481.fits
sourcepop_total = 3 1 1 0 1=both methods 3=use healpix method 4=lb method ? NO 1=healpix method
sourcepop_sublimit = 3 1 0
sourcepop_soplimit = 3 1 0
gcr_spectra = 1
gcr_spectra_r_min = 7.5
gcr_spectra_r_max = 8.5
gcr_spectra_z_min = 0.0
gcr_spectra_z_max = +0.1
gcr_spectra_Z = 01,02,-1,+1 list of Z (separate by comma, no blanks)
gcr_spectra_A = 01,04,00,00 list of A (separate by comma, no blanks)
gcr_spectra_n_ratio = 3 number of sec/prim ratios
gcr_spectra_sec_Z_01 = 05,05 list of secondary Z for sec/prim ratio 1
gcr_spectra_sec_A_01 = 10,11 list of secondary A for sec/prim ratio 1
gcr_spectra_pri_Z_01 = 06,06 list of primary Z for sec/prim ratio 1
gcr_spectra_pri_A_01 = 12,13 list of primary A for sec/prim ratio 1
gcr_spectra_ramin_01 = 0.0 minimum value for sec/prim ratio 1
gcr_spectra_ramax_01 = 1.0 maximum value for sec/prim ratio 1
gcr_spectra_sec_Z_02 = 04 list of secondary Z for sec/prim ratio 1
gcr_spectra_sec_A_02 = 10 list of secondary A for sec/prim ratio 1
gcr_spectra_pri_Z_02 = 04 list of primary Z for sec/prim ratio 1
gcr_spectra_pri_A_02 = 09 list of primary A for sec/prim ratio 1
gcr_spectra_ramin_02 = 0.0 minimum value for sec/prim ratio 1
gcr_spectra_ramax_02 = 0.6 maximum value for sec/prim ratio 1
gcr_spectra_sec_Z_03 = 05 list of secondary Z for sec/prim ratio
gcr_spectra_sec_A_03 = -1 list of secondary A for sec/prim ratio -1 does not work for model
gcr_spectra_pri_Z_03 = 06 list of primary Z for sec/prim ratio
gcr_spectra_pri_A_03 = -1 list of primary A for sec/prim ratio -1 does not work for model
gcr_spectra_ramin_03 = 0.0 minimum value for sec/prim ratio
gcr_spectra_ramax_03 = 1.0 maximum value for sec/prim ratio
#gcr_spectra_sec_Z_01 = 05,04 list of secondary Z for sec/prim ratio 1
#gcr_spectra_sec_A_01 = -1,-1 list of secondary A for sec/prim ratio 1
#gcr_spectra_pri_Z_01 = 06,08 list of primary Z for sec/prim ratio 1
#gcr_spectra_pri_A_01 = -1,-1 list of primary A for sec/prim ratio 1
#gcr_spectra_ramin_01 = 0.0 minimum value for sec/prim ratio 1
#gcr_spectra_ramax_01 = 0.4 maximum value for sec/prim ratio 1
gcr_spectra_Emin = 1.0e2 1.0e3 0.0e+1 minimum energy on spectrum plot 0.0-> all energies
gcr_spectra_Emax = 1.0e7 1.0e6 1.0e+7 maximum energy on spectrum plot 0.0-> all energies
gcr_spectra_Imin = 1.0e1 1.0e4 1.0e2 1.0e5 minimum intensity*E^3 on spectrum plot (protons, electrons different ranges)
gcr_spectra_Imax = 1.0e8 1.0e8 1.0e5 1.0e7 maximum intensity*E^3 on spectrum plot
gcr_spectra_mod_phi = 0.,500.,600 modulation phi (MV)
gcr_database_file = GCR_data_11.dat GCR experiments database file
sync_index_nu = 38,408,408,1420,408,30000,20000,30000 frequencies (MHz) in pairs, for synchrotron index (=4 pairs at present)
sync_data____10MHz = 1 1=plot data at this frequency
sync_data____22MHz = 1 1=plot data at this frequency
sync_data____45MHz = 1 1=plot data at this frequency
sync_data___150MHz = 1 1=plot data at this frequency
sync_data___408MHz = 1 1=plot data at this frequency
sync_data___820MHz = 1 1=plot data at this frequency
sync_data__1420MHz = 1 1=plot data at this frequency
sync_data__2326MHz = 1 1=plot data at this frequency
sync_data_22800MHz = 1 1=plot data at this frequency
sync_data_33000MHz = 1 1=plot data at this frequency
sync_data_41000MHz = 1 1=plot data at this frequency
sync_data_61000MHz = 1 1=plot data at this frequency
sync_data_94000MHz = 1 1=plot data at this frequency
sync_data_WMAP = 4 0=WMAP team 3yr sync maps, 1-4= Miville-Deschenes beta-scaled 408MHz
luminosity = 1 only if ISRF_read=1
galdef_series = 1 for several galdef runs plot_luminosity_multiple_galdef 0=no plot
ISRF_luminosity_R = 30 R (kpc) for ISRF luminosity calculation
ISRF_luminosity_z = 30 z (kpc) for ISRF luminosity calculation
ISRF_read = 0 1