Modifications

July 7, 2010
Began adding documentation.


CRT Input

This page will provide you with basic information you need to efficiently construct an input file for use with CRT. Full options are listed in the documentation included with the source code. If you have any questions please contact us.


Input

Overview

A certain format is defined for reading a file containing a list of source types and magnetic field configurations. Comments in the input file are defined as lines starting with a #. Since the particle only 'feels' the total magnetic field, the program will add together any individual field components that the user specifies. A wide variety of sources can be generated, including multi-component sources (proton plus iron, for example). The order in which fields and sources does not matter, along as each declaration is complete. The delimiter between lines is a CARRIAGE RETURN on the system in which CRT was compiled. Example input format follows below.


Sources

Templates

  • isotropic:
    I NPSRC CRSHIFT CRINDEX EMIN EMAX PMASS PCHARGE
  • isotropic-mono-plaw-auger:
    X NPSRC CRSHIFT CRINDEX EMIN EMAX PMASS PCHARGE
  • point:
    P NPSRC CRSHIFT CRINDEX EMIN EMAX D LONG LAT BEAMANGLE PMASS PCHARGE
  • long:
    L NPSRC CRSHIFT ENERGY STARTANGLE PMASS PCHARGE
  • long-mono-plaw-auger:
    Y NPSRC CRSHIFT ENERGY STARTANGLE PMASS PCHARGE
  • lat:
    B NPSRC CRSHIFT ENERGY STARTANGLE PMASS PCHARGE
  • lat-mono-plaw-auger:
    Z NPSRC CRSHIFT ENERGY STARTANGLE PMASS PCHARGE
  • CR:
    C ENERGY LONG LAT PMASS PCHARGE
  • ArbTraj:
    A ENERGY PX PY PZ VX VY VZ PMASS PCHARGE
  • EGsrc:
    E NPSRC CRSHIFT CRINDEX EMIN EMAX D LONG LAT BRAND LCORR LEG PMASS PCHARGE
  • isotropic-broken-plaw:
    BPI NPSRC CRSHIFT INDEX0 INDEX1 EBREAK EMIN EMAX PMASS PCHARGE
  • isotropic-broken-plaw-auger (backtrk and Auger exposure):
    BPX NPSRC CRSHIFT INDEX0 INDEX1 EBREAK EMIN EMAX PMASS PCHARGE

Units and descriptions

  • NPSRC - number of particles from this source
  • CRSHIFT - [kpc], size of disk within which to generate particles
  • CRINDEX - spectral index for plaw spectrum (over entire energy range), DO NOT SET = (-1)
  • EMIN - [EeV], minimum energy of spectrum
  • EMAX - [EeV], maximum energy of spectrum
  • EBREAK - [EeV], break energy of broken-plaw spectrum
  • INDEX0 - spectral index below EBREAK
  • INDEX1 - spectral index above EBREAK
  • PMASS - mass of particle, set as multiples of proton mass (PMASS=1 -> m=1.67*10^-27 kg, PMASS=26 -> 26*1.67*10^-27 kg)
  • PCHARGE - charge of particle, set as multiples of proton charge (same as PMASS)
  • D - [kpc], distance to source. if D=(-999), then source is assumed EXG and distance automatically set to 50kpc
  • LONG - [degrees], Galactic longitude coordinate of source
  • LAT - [degrees], Galactic latitude coordinate of source
  • STARTANGLE - [degrees], coordinate of strip for Long/Lat sources
  • BEAMANGLE - [degrees], opening angle of 'jet' at the source, uniform over solid angle. All particle velocity vectors will lie within directions encompassed by the 'jet.'
  • ENERGY - [EeV], energy of particle
  • P{xyz} - [kpc], the {xyz}-component of the source location for the ArbTraj source class in System-A
  • V{XYZ} - the unnormalized {xyz}-component of the particle velocity
  • BRAND - [microGauss], magnitude of turbulent magnetic field outside Galaxy
  • LCORR - [kpc], correlation length of turbulent magnetic field outside Galaxy
  • LEG - [kpc], distance to source


Magnetic Fields

Templates

  • ASS_A:
    F ass_a model NORM PITCH SCALE1 SCALE2 SCALE3 SCALE4 ZCUT
  • ASS_S:
    F ass_s model NORM PITCH SCALE1 SCALE2 SCALE3 SCALE4 ZCUT
  • BSS_A:
    F bss_a model NORM PITCH SCALE1 SCALE2 SCALE3 SCALE4 ZCUT
  • BSS_S:
    F bss_s model NORM PITCH SCALE1 SCALE2 SCALE3 SCALE4 ZCUT
  • UNIFORM:
    F uniform BX BY BZ
  • DIPOLE:
    F dipole NORM
  • SIMRNDM:
    F simprand NORM corrlen sigmalen sigmanorm
  • RING:
    F ring NORM INNEREDGE OUTEREDGE SCALE5
  • TOROID:
    F toroidal TORRAD SCALE6 LORWIDTH BMAX

Units and descriptions

  • 'model' - select either 'stanev' or 'hmr'. Descriptions of these fields can be found in: Stanev. Astropart. Journ. 479:290, 1997 (arXiv:astro-ph/9607086) Harari, Mollerach, Roulet. JHEP 9908:022, 1999 (arXiv:astro-ph/9906309)
  • NORM - [microGauss], the total magnitude of the spiral field at the position of the Sun for *SS_* models. While value of (1) will produce a dipole field with components Bx=By=0, Bz=(1 microGauss).
  • PITCH - [degrees], looking from the Galactic North Pole, the pitch angle is positive if the clockwise tangent to the spiral is outside the circle with radius R
  • SCALE1 - [kpc], Galactocentric distance of the location of maximum field strength in the direction (0,0)
  • SCALE2 - [kpc], scale length for tanh attenuation factor in the Galactic disk (not used for stanev field, but an entry here is required anyway)
  • SCALE3 - [kpc], scale length for exponential decay along z-direction
  • SCALE4 - [kpc], scale length for exponential decay along z-direction
  • ZCUT - [kpc], distance in z-direction where exponential decay scale length changes from SCALE3 to SCALE4 (not used for hmr field, but an entry here is required anyway)
  • B{X,Y,Z} - [microGauss], same as for NORM
  • INNEREDGE - [kpc], inner edge of annulus of nonzero field magnitude
  • OUTEREDGE - [kpc], outer edge of annulus
  • SCALE5 - [kpc], scale length for exponential decay along z-direction
  • TORRAD - [kpc], radius of circle for constant field magnitude
  • SCALE6 - [kpc], scale height above Galactic plane
  • LORWIDTH - [kpc], half-width of Lorentzian distribution
  • BMAX - [kpc], maximum value of field magnitude
  • corrlen - [kpc], nominal cell size.
  • sigmalen - [kpc], sigma of a gaussian defining the variation of cell to cell size.
  • sigmanorm - [microGauss], sigma of a gaussian defining the variation of cell to cell field strength.

 Note:SCALE2 is not used for 'stanev' fields. ZCUT is not used for 'hmr' fields. Simply set these values to (1) in these cases.


Example Input File

I 10 0.01 -2.7 3 200 1 1
P 10 0 -2.7 3 200 1 -5 10 0 1 1
BPI 10 0.01 -2.7 -3 50 3 200 1 1
F jf2012 0 0