| 1 | #  SKYMIX  NbComponents  NSide_HealPix | 
|---|
| 2 | @SKYMIX  1 16 | 
|---|
| 3 |  | 
|---|
| 4 | #  READMAP  FitsName | 
|---|
| 5 | #  To Read the already prepared (mixed) Output MAP from FITS file | 
|---|
| 6 | #  FitsName should be the complete path to the FITS file | 
|---|
| 7 | READMAP xxx.fits | 
|---|
| 8 |  | 
|---|
| 9 | #  Defining the Detection filter pass - band | 
|---|
| 10 | #  GAUSSFILTER or FILTERFITSFILE  card | 
|---|
| 11 | #  GAUSSFILTER -> Gaussian filter | 
|---|
| 12 | #  FILTERFITSFILE -> Filter (nu, T(nu)) from FITS file | 
|---|
| 13 | #  GAUSSFILTER  Nu0 Sigma_Nu Tmax NuMin NuMax | 
|---|
| 14 | #  FILTERFITSFILE FileName NuMin NuMax | 
|---|
| 15 | #       FileName should be a complete path | 
|---|
| 16 | @GAUSSFILTER 143 12 0.84 50. 300. | 
|---|
| 17 |  | 
|---|
| 18 | #  MAPPATH  PathForFITS | 
|---|
| 19 | #  path for input Sky Component FITS files , and EmissionSpectra files | 
|---|
| 20 | MAPPATH /data/opus/Rz/SOPHYA/ | 
|---|
| 21 |  | 
|---|
| 22 | #  ----        Sky components  --------- | 
|---|
| 23 | #      ----     HealPix maps   ---- | 
|---|
| 24 | #  MAPFITSFILEi FITSfilename Normalisation | 
|---|
| 25 | @MAPFITSFILE4  gauss.fits  1.0 | 
|---|
| 26 |  | 
|---|
| 27 | #      ---   Emission Spectra definition ----- | 
|---|
| 28 | #   SPECTRAFITSFILEi  FITSfilename  Fmin Fmax | 
|---|
| 29 | @SPECTRAFITSFILE4  pl.fits 100 600 | 
|---|
| 30 |  | 
|---|
| 31 | #   Other possible definition of emission spectra | 
|---|
| 32 | #   BLACKBODYi  temperature | 
|---|
| 33 | #  The blackbody object returns flux in units of watt/m2/sr/hz | 
|---|
| 34 | #   We thus scale the map (x 1.e17) | 
|---|
| 35 | @MAPFITSFILE1  cossin.fits  1.e17 | 
|---|
| 36 | @DIPOLE1 2.726 | 
|---|
| 37 |  | 
|---|
| 38 | #   For power-law spectra  f = a ((nu-nu0)/dnu)^b) | 
|---|
| 39 | #   POWERLAWSPECTRAi  a nu0 dnu b Fmin Fmax | 
|---|
| 40 | @MAPFITSFILE2  noise.fits    0.1 | 
|---|
| 41 | @POWERLAWSPECTRA2 1. 150. 50. -0.5 100. 500. | 
|---|
| 42 |  | 
|---|
| 43 | #  For an emission spectrum which depends on the position | 
|---|
| 44 | #  in the map -> synchrotron for instance nu^(-beta(theta,phi)) | 
|---|
| 45 | #  MAPFITSFILE will still correspond to the brightness temperature | 
|---|
| 46 | @MAPFITSFILE3 /exp/planck/Sync/Template408MHz.fits 0.2 | 
|---|
| 47 | #@MAPFITSFILE1 gauss.fits 0.2 | 
|---|
| 48 | # and the BETAFITSFILE will correspond to a map with the spectral | 
|---|
| 49 | # index beta(theta,phi) | 
|---|
| 50 | # 408 is the normalisation frequency of the map | 
|---|
| 51 | # 32 is the nside of the map on which the integration over the | 
|---|
| 52 | #       spectralResponse of the detector will be performed | 
|---|
| 53 | #       for the full resolution, write the nside of the map | 
|---|
| 54 | #       or a negative value. | 
|---|
| 55 | @BETAFITSFILE3 /exp/planck/Sync/SpectralIndex.fits  408e-3 16 | 
|---|
| 56 | #@BETAFITSFILE1 gauss.fits  408e-3 32 | 
|---|
| 57 |  | 
|---|
| 58 | #  Define the Debug level | 
|---|
| 59 | @DEBUGLEVEL  0 | 
|---|
| 60 | #  Define the Print level | 
|---|
| 61 | @PRINTLEVEL  5 | 
|---|