c===============================================================================
c ERLAEUTERUNGEN : siehe Datei-Ende
c===============================================================================
c KOMMENTARZEILEN: (Zeilen beginnend mit '@')
c===============================================================================

 $parameter_liste

   TestRun    =   T             ! -> fileNummern zwischen 9900 und 9999  [false]

 !==============================================================================
 !			WEITERFUEHRUNG FRUEHERER SIMULATIONEN
 !==============================================================================

  fileName_ACCEL      = 'ac_9915'	! ACCEL-file
  fileName_MUTRACK    = 'mu_9915'	! prev. MUTRACK sim. up to TD-Foil

  previousSimulation =   1      ! fruehere Simulation fortfuehren?       [   0 ]
                                ! 0 : fromScratch
                                ! 1 : ACCEL
                                ! 2 : MUTRACK ('FoilFile')

 ! im Fall von 'priviousSimulation' = 1 oder = 2 (ACCEL oder MUTRACK):

  previousSettings =      T     ! uebernehme fruehere Schleifensettings  [false]
     ow_U_Tgt   =      f        ! - aber ueberschreibe 'U_Tgt_'          [false]
     ow_U_Gua   =      f        ! - aber ueberschreibe 'U_Gua_'          [false]
     ow_U_G1    =      f        ! - aber ueberschreibe 'U_G1_'           [false]
     ow_U_L1    =      f        ! - aber ueberschreibe 'U_L1_'           [false]
     ow_U_Sp    =      f        ! - aber ueberschreibe 'U_Sp_'           [false]
     ow_U_L2    =      f        ! - aber ueberschreibe 'U_L2_'           [false]
     ow_U_Folie =      f        ! - aber ueberschreibe 'U_Folie_'        [false]
    !ow_B_TD    =  NICHT IMPLEMENTIERT   ueberschreibe 'BTD_'            [false]
    !ow_B_Helm  =  NICHT IMPLEMENTIERT   ueberschreibe 'BHelm_'          [false]
     ow_masse   =      f        ! - aber ueberschreibe 'masse_'          [false]
     ow_ladung  =      f        ! - aber ueberschreibe 'ladung_'         [false]
     ow_E0      =         T     ! - aber ueberschreibe 'E0_'             [false]
     ow_y0      =      f        ! - aber ueberschreibe 'y0_'             [false]
     ow_z0      =      f        ! - aber ueberschreibe 'z0_'             [false]
     ow_theta0  =      f        ! - aber ueberschreibe 'theta0_'         [false]
     ow_phi0    =      f        ! - aber ueberschreibe 'phi0_'           [false]
     ow_alfaTgt =      f        ! - aber ueberschreibe 'alfaTgt_'        [false]
     ow_alfaSp  =      f        ! - aber ueberschreibe 'alfaSp_'         [false]
     ow_alfaTD  =      f        ! - aber ueberschreibe 'alfaTD_'         [false]
     ow_deltaL1 =      f        ! - aber ueberschreibe 'deltaL1_'        [false]
     ow_artList =      f        ! - aber ueberschreibe 'artlist'         [false]


 !==============================================================================
 !			KAMMER-SETTINGS UND TEILCHEN-PARAMETER
 !==============================================================================
 !   Einheiten: Spannungen in kV, Winkel in degree, Masse in keV/c**2,
 !              Ladung in pos.Elementarlad., Energie in keV, Laengen in mm,
 !		Zeiten in ns, Geschwindigkeiten in mm/ns
 !------------------------------------------------------------------------------

 ! - - - - - - - - - - Schleifen-Parameter - - - - - - - - - - - - - - - - - - -

   U_Tgt_     =   18                   ! U(Target)                       [ 0. ]
   U_Gua_     =   15                   ! U(Guard)                        [ 0. ]
   U_G1_      =   12                   ! U(Gitter)                       [ 0. ]
   U_L1_      =   5                    ! U(kalte Linse)                  [ 0. ]
   U_Sp_      =   13.5                 ! U(Spiegel)                      [ 0. ]
   U_L2_      =   4.5                  ! U(warme Linse)                  [ 0. ]
   U_L3_      =   4.5          	       ! U(warme Linse)                  [ 0. ]

   U_Folie_   =   -2.980               ! U(Folie)                        [ 0. ]
   U_Vorne_   =   -2.790               ! U(Vorne)                        [ 0. ]
   U_Hinten_  =   -3.230               ! U(Hinten)                       [ 0. ]
   U_MCP3_    =   -2.600               ! U(MCP3)                         [ 0. ]
   
   U_MCP2_    =   -1.87                ! U(MCP2)			 [ 0. ]

  !B_TD_      =   NICHT IMPLEMENTIERT  ! BFeld(Kompensationsspule)       [ 0. ]
  !B_Helm_    =   NICHT IMPLEMENTIERT  ! BFeld(Helmholtzpule)            [ 0. ]

   alfaTgtVertically = f               ! alfa(Target) = vertikale Verkippung
   alfaTgt_   =    0                   ! alfa(Target)  (x=0,y=90)	 [ 0. ]
   alfaSp_    =   45                   ! alfa(Spiegel) (x=0,y=90)        [45. ]
   alfaTD_    =    0                   ! alfa(Trigger) (x=0,y=90)	 [ 0. ]

   Masse_     =    0                   ! (*) Masse                      [105659]
   Ladung_    =    0                   ! (*) Ladung                      [ 1. ]

   E0_        =    0.01                ! (*) Startenergie (in keV!)      [ 0. ]
   y0_        =    0                   ! (*) y0 (Startposition)          [ 0. ]
   z0_        =    0                   ! (*) z0 (Startposition)          [ 0. ]
   theta0_    =    0                   ! (*) theta0 (Startwinkel)(x=0)   [ 0. ]
   phi0_      =    0                   ! (*) phi0 (Startwinkel)(y=0,z=90)[ 0. ]

   Thickness_ =   3                    ! Foliendicke in ug/cm**2         [ 3.0]
   mean_Eloss_=    15                  ! mittlerer Energieverlust in keV [ 0.0]

   DeltaL1_   = 0!  -3,3,3             ! Verlaengerung 1 in mm           [ 0.0]
   DeltaL2_   = 0! -5,+5,1             ! Verlaengerung 2 in mm           [ 0.0]

				       ! (*) <-> siehe Dateiende

 ! - - - - - - - - - - Zufallsverteilte Startparameter - - - - - - - - - - - - -

	randomStarts   = 5000        ! Anzahl zufallsverteilter Starts    [ 50 ]

 ! Start-Energie:
	random_energy   =    1      ! 0: gemaess 'E0_'                   [  0 ]
                                    ! 1: E0_ + gleichverteilter Offset
                                    ! 2: E0_ + gaussverteilter Offset
	lowerE0         = .000      ! min. Energie (in keV!)             [0.000]
	upperE0         = .001      ! max. Energie (in keV!)             [0.010]
	sigmaE0         = .005      ! Breite der Gaussverteilung         [0.010]

 ! Start-Position:
	random_position =    0      ! 0: gemaess 'y0_' und 'z0_'	 [  0 ]
				    ! 1: gleichverteilt auf Viereckflaeche
				    ! 2: gleichverteilt auf Kreisflaeche
				    ! 3: Gaussverteilt auf Viereckflaeche
				    ! 4: Gaussverteilt auf Kreisflaeche
	StartBreite     =   -1      !        			  [TargetBreite]
	StartHoehe      =   -1      !        			  [TargetHoehe ]
	StartRadius     =   -1      !       			         [ 40 ]
	sigmaPosition   =   15      ! Breite der Gaussverteilung	 [ 15 ]


 ! Start-Winkel:
	random_winkel   =    0      ! 0: gemaess 'theta0_' und 'phi0_'	 [  0 ]
				    ! 1: Lambertverteilt
				    ! 2: Gaussverteilt
	StartLambertOrd =    1      ! Ordnung der Lambert-Verteilung	 [  1. ]
	sigmaWinkel     =    5      ! Breite der Gaussverteilung	 [  1. ]

 ! - - - - - - - - - - Projektile (case-sensitive) - - - - - - - - - - - - - - -

        !   maximal 9 verschiedene aus:

        !   e+, e-, m+, m-, Mu, Mu-, H+, H, H-, H2+, H2, H2-, alfa      ![ ' ' ]
        !   H2O1, H2O2, H2O3, H2O4, H2O5,
        !   Hyd1, Hyd2, Hyd3, Hyd4, Hyd5, Hyd6, Hyd7, Hyd8, Hyd9,
 	!   A11+, A12+, A21+, A31+, A32+, N11+, N21+, K11+, K12+

        !   H2O1,H2O2,..,H2Oi,..: Wassercluster (H2O)i
        !   Hyd1,Hyd2,..,Hydi,..: protonierte Wassercluster [(H2O)i]H+
        !   A: Argon, N: Stickstoff, K: Krypton

        !   (alle Zeichen in artList nach einem '!' werden ignoriert)

        artList = 'm+'

	!   nur verwendet falls 'log_neutralize = t': neutrale Anteile in %
        !   (eine Angabe je Projektil, '-1' -> verwende Funktion von M. Gonin
	!   fuer Neutralisierung von Protonen in Graphit).

	neutral_fract	  =  0


 ! - - - - - - - - - - Trigger-Detektor- - - - - - - - - - - - - - - - - - - - -
                                                                                
	TriggerInBeam     =   t     ! TriggerDetektor im Strahlweg?      [false]
	transTDFoil       =  88.    ! Transmission der Folie in %	 [100.]

  	log_neutralize	  =   F     ! Neutralisierung in der Folie?	 [false]

	n_E_Verlust       =  5      ! 0: kein Energieverlust             [  0 ]
				    !<0: mittl. Verlust gemaess 'mean_Eloss_'
				    !>0: mittl. Verlust gemaess ICRU-Tabelle
				    !|1| keine Energieaufstreuung
				    !|2| Streuung gaussvert. gemaess 'sigmaE'
				    !|3| Streuung gleichvert. [lowerE,upperE]
				    !|4| Streuung gaussv. nach Lindhard/Scharff
				    !|5| Streuung gaussv. nach Yang (recommended
	calculate_each    =    F    ! bei ICRU: fuer jedes Teilchen neu? [false]
	graphitData       =    F    ! bei ICRU: Elossdaten fuer Graphit? [false]
				    ! (Alternative: amorpher Kohlenstoff)
	sigmaE            =  0.5    ! GaussBreite des Energieverlustes   [  .5]
	lowerE            =  0      ! untere Grenze fuer Energieverlust  [ 0. ]
	upperE            =  0.050  ! obere Grenze fuer Energieverlust   [  .1]

	n_Aufstreu        =    3    ! 0: keine Winkelaufstreuung	 [  0 ]
				    ! 1: gemaess 'sigmaAufstreu'
				    ! 2: gemaess Meyer (gaussverteilte Winkel)
				    ! 3: gemaess Meyer (nach f-Funktion)
	sigmaAufstreu     =   2.    ! Gaussbreite der Aufweitung         [ 5. ]

	generate_FE       =   F     ! FolienElektronen generieren?	 [false]


 !==============================================================================
 !			AUSGABE - STEUERUNG
 !==============================================================================

   writeTraj2File = f	! Um das rausschreiben der Trajektorie in ein separates
			! file 'OUT.TRAJ' zu ermoeglichen, muessen im SourceCode
			! alle mit 'cMBc' herauskommentierten Zeilen aktiviert
			! werden. (Dieses feature wurden fuer M.B. implementiert
			! aber wegen Rechenzeitargumenten herauskommentiert) 

 ! - - - - - - - - PHYSICA-Tabelle erstellen?- - - - - - - - - - - - - - - - - -

 	createPhysTab =  f          !                 			 [false]

 ! welche Statistiken?

	PHY_S1xM2     =   t         ! S1xM2     im PHYSICA-file?         [false]
	PHY_S1M2      =   t         ! S1M2      im PHYSICA-file?         [false]
	PHY_S1Fo      =   t         ! S1Fo      im PHYSICA-file?         [false]
	PHY_FoM2      =   t         ! FoM2      im PHYSICA-file?         [false]
	PHY_S1M3      =   t         ! S1M3      im PHYSICA-file?         [false]
	PHY_M3M2      =   t         ! M3M2      im PHYSICA-file?         [false]
	PHY_t_FE      =   t         ! tFE       im PHYSICA-file?         [false]
	PHY_y_Fo      =   t         ! y(Folie)  im PHYSICA-file?         [false]
	PHY_z_Fo      =   t         ! z(Folie)  im PHYSICA-file?         [false]
	PHY_r_Fo      =   t         ! r(Folie)  im PHYSICA-file?         [false]
	PHY_y_M2      =   t         ! y(MCP2)   im PHYSICA-file?         [false]
	PHY_z_M2      =   t         ! z(MCP2)   im PHYSICA-file?         [false]
	PHY_r_M2      =   t         ! r(MCP2)   im PHYSICA-file?         [false]
	PHY_y_xM2     =   t         ! y(xMCP2)  im PHYSICA-file?         [false]
	PHY_z_xM2     =   t         ! z(xMCP2)  im PHYSICA-file?         [false]
	PHY_r_xM2     =   t         ! r(xMCP2)  im PHYSICA-file?         [false]

 ! welche statistischen Groessen?

	PHY_mean      =  t          ! Mittelwerte?                       [false]
	PHY_variance  =  t          ! Varianzen?                         [false]
	PHY_minimum   =  t          ! Minimalwerte?                      [false]
	PHY_maximum   =  t          ! Maximalwerte?                      [false]
	PHY_percent   =  t          ! rel. Anteil beitragender Teilchen? [false]

 ! - - - - - - - - NTupel-relevante Vorgaben - - - - - - - - - - - - - - - - - -

	createNTP     =   t         ! NTP-File erzeugen			 [false]

        Fo_triggered  =   F         ! TD-Folientreffer verlangt          [false]
        xM2_triggered =   F         ! Erreichte MCP2-Ebene verlangt      [false]
        M2_triggered  =   F         ! MCP2-Treffer verlangt              [false]

	smearS1Fo     =  f	    ! Verschmierung von S1Fo in NTP (zur
	sigmaS1fo     = 1.6	    !    Beruecksichtigung der Zeitaufloesung)

	NTP_S1xM2     =   F         ! - Flugzeit bis MCP2-Ebene?	 [false]
	NTP_times     =   F         ! - sonstige Flugzeiten?		 [false]
	NTP_FoM2Only  =   F         ! - von 'times' nur FoM2 (rezessiv!) [false]
	NTP_lifetime  =   F         ! - Lebenszeiten?			 [false]
	NTP_start     =   F         ! - Startgroessen?  		 [false]
	NTP_stop      =   F         ! - Stopgroessen?			 [false]
	NTP_40mm      =   F         ! - Groessen bei x=40mm?		 [false]
  	NTP_Folie     =   t         ! - Everlust und Aufstreuung?	 [false]
	NTP_charge    =   F         ! - Ladungszustand nach Folie?	 [false]
	NTP_steps     =   F         ! - min. und max. Schrittweiten?	 [false]

 ! - - - - - - - - - - Vorgaben fuer das Summary - - - - - - - - - - - - - - - -
                                                                                
	n_outWhere   =  2           ! 0-> keine Ausgabe des Summarys   	 [  2 ]
                                    ! 1-> Ausgabe nur in Datei
                                    ! 2-> in Datei und auf Schirm
                                    ! 3-> nur auf Schirm
                                    !-1-> nur kleines Logfile
                                                                                
	SUM_S1xM2    =   t          ! S1xM2     im Summary? 		 [false]
	SUM_S1M2     =   t          ! S1M2      im Summary? 		 [false]
	SUM_S1Fo     =   t          ! S1Fo      im Summary? 		 [false]
	SUM_FoM2     =   t          ! FoM2      im Summary?		 [false]
	SUM_S1M3     =   t          ! S1M3      im Summary?		 [false]
	SUM_M3M2     =   t          ! M3M2      im Summary?		 [false]
	SUM_t_FE     =   t          ! tFE       im Summary?		 [false]
	SUM_y_Fo     =   t          ! y(Folie)  im Summary?		 [false]
	SUM_z_Fo     =   t          ! z(Folie)  im Summary?		 [false]
	SUM_r_Fo     =   t          ! r(Folie)  im Summary?		 [false]
	SUM_y_M2     =   t          ! y(MCP2)   im Summary?		 [false]
	SUM_z_M2     =   t          ! z(MCP2)   im Summary?		 [false]
	SUM_r_M2     =   t          ! r(MCP2)   im Summary?		 [false}SUM_FoM2SUM_FoM2
   	SUM_y_xM2    =   t          ! y(xMCP2)  im Summary?		 [false]
	SUM_z_xM2    =   t          ! z(xMCP2)  im Summary?		 [false]
	SUM_r_xM2    =   t          ! r(MxCP2)  im Summary?		 [false}SUM_FoM2SUM_FoM2
                                                                                
	log_out_FE        =    F    ! die FE-Statistik im Summary?	 [false]
	log_out_pfosten(1)=    F    ! P : Pfosten einzeln im Summary?    [false]
	log_out_pfosten(2)=    F    ! FE: Pfosten einzeln im Summary?    [false]

 ! - - - - - - - - File mit zugrundeliegender Geometrie erstellen? - - - - - - -

	write_geo     =        T    !					 [false]

 ! - - - - - - - - zu erstellende Tabellenfiles- - - - - - - - - - - - - - - - -

	Tab_S1xM2     =     F       ! Tabelle S1xM2			 [false]
	Tab_S1M2      =     F       ! Tabelle S1M2			 [false]
	Tab_S1Fo      =     F       ! Tabelle S1Fo			 [false]
	Tab_FoM2      =     F       ! Tabelle FoM2			 [false]
 	Tab_S1M3      =     F       ! Tabelle S1M3			 [false]
	Tab_M3M2      =     F       ! Tabelle M3M2			 [false]
	Tab_t_FE      =     F       ! Tabelle tFE			 [false]
	Tab_y_Fo      =     F       ! Tabelle y(Folie)			 [false]
	Tab_z_Fo      =     F       ! Tabelle z(Folie)			 [false]
	Tab_r_Fo      =     F       ! Tabelle r(Folie)			 [false]
	Tab_y_M2      =     F       ! Tabelle y(MCP2)			 [false]
	Tab_z_M2      =     F       ! Tabelle z(MCP2)			 [false]
	Tab_r_M2      =     F       ! Tabelle r(MCP2)			 [false]
	Tab_y_xM2     =     F       ! Tabelle y(xMCP2)			 [false]
	Tab_z_xM2     =     F       ! Tabelle z(xMCP2)			 [false]
	Tab_r_xM2     =     F       ! Tabelle r(xMCP2)			 [false]

 ! - - - - - - - - - - Graphikausgabe- - - - - - - - - - - - - - - - - - - - - -
                                                                                
	GRAPHICS          =   F     ! Graphikausgabe der Trajektorien?   [false]
	GRAPHICS_Anzahl	  =  40     ! fuer wieviele Starts je Schleife?  [  25 ]
  	vertical          =   T     ! Blick auf die Kammer VON OBEN?     [true!]
	log_marker	  =   F     ! Marker am Trajektorienende?	 [false]
	plot_FE           =   F     ! FE-Trajektorien plotten?           [false]

	n_postSkript      =     0   ! Postskript-files erstellen?        [  1 ]
	                            ! 0 : nicht erstellen
	                            ! 1 : Abfrage nach jeder Schleife
	                            ! 2 : immer erstellen

	iMonitor          =	3   ! Abtastfrequenz fuer Graphik        [  2 ]
	color             =     1   ! 0=schwarzweiss, 1=farbig           [  0 ]
	schnitt_p         =     2   ! Kammerteil: 1=horiz., 2=vertikal   [  2 ]
	schnitt_x         =    -1   ! x-Koord. der Schnittebene	         [MCP2]

 !==============================================================================
 !			GRUNDEINSSTELLUNGEN
 !==============================================================================

 ! - - - - - - - - - - Programmsteuerung - - - - - - - - - - - - - - - - - - - -

	TestOnWireHit     =    t    ! Drahttreffer beruecksichtigen?     [false]
	UseDecay          =    t    ! MYONEN-Zerfall beruecksichtigen?   [false]

	upToTDFoilOnly	  =  F      ! Integration nur bis zur TD-Folie   [false]

	createFoilFile    =  F      ! neues 'FoilFile' erstellen?	 [false]

	idealMirror	  =  F      ! Spiegel als ideal betrachten?	 [false]

	DEBUG             =  F      ! Debug-Version?			 [false]
	DEBUG_FE          =  F      ! FE in Debug-Info?			 [false]
	DEBUG_Anzahl	  =  3	    ! fuer wieviele Starts je Schleife?  [   2 ]

 ! - - - - - - - - - - Geometrie - - - - - - - - - - - - - - - - - - - - - - - -

  	Geo_fileName      = 'GEO_KAMMER_RUN9_NEW_MUONS'
 ! 	Geo_fileName      = 'GEO_KAMMER_RUN9_NEW_PROTONS'
	
 ! 	Geo_fileName      = '9nM'   ! RunNummer bzw. Name des Geo-files
 !	Geo_fileName      = '9'     ! RunNummer bzw. Name des Geo-files
				    ! '9n': mit laengerer Mappe fuer MCP2

	gridInFrontOfFoil =  f      ! Gitter vor der Triggerfolie?       [false]
                                    ! ('grindInFrontOfFoil' nur bei Run 9)
	xBlende		  = -1      ! Abstand der Blende zum MCP2   (-1 => none)
	radius_Blende	  =  0      ! Radius der Blende vor MCP2 

 ! - - - - - - - - - - Startflaeche- - - - - - - - - - - - - - - - - - - - - - -
                                                                                
	StartFlaeche    =  0        !-1: entsprechend x0_ und Kammerteil [  0 ]
	                            ! 0: Target
                                    ! 1: 1. Gitter
                                    ! 2: Folie (Startwinkel im Kammersystem)
                                    ! 3: Folie (Startwinkel im Triggersystem)

	x0_             =  0        ! Startort rel zu. Tgt bzw. Spiegel  [ 0. ]
	Kammerteil      =  1        ! 1: horizontal, 2: vertikal         [  1 ]

 ! - - - - - - - - - - Fehlerkontrolle - - - - - - - - - - - - - - - - - - - - -

	eps_x       =  1.e-6        ! max. Fehler im Ort                 [1.e-5]
	eps_v       =  1.e-6        ! max. Fehler in Geschw.             [1.e-5]
	log_relativ =    F          ! relative Fehlerbetrachtung?        [false]

	maxStep     =  7000         ! max. Anzahl Integrationsschritte   [ 1000]

	dtsmall	    =  0.00001	    ! kleinster Integrations-Zeitschritt [0.001]

	maxBelowDtSmall = 700 	    ! max. Anzahl von dtSmall Unterschr. [  50 ]

 ! - - - - - - - - - - seed- - - - - -  - - - - - - - - - - - - - - - - - - - -

  	seed_	    =    0          ! > 1E6, ungerade   (if .LE. 0 => automatic)

 $END

c********************************************************************************
c*
c*  In dieser Datei werden die Parameter fuer die Flugbahnberechnung von 
c*  Teilchen in der NEMU-Apparatur mittels des Fortranprogramms 'MUTRACK'
c*  festgelegt. 
c*
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Alle Zeilen vor der Zeile ' $parameter_liste', die nicht mit 'c ' beginnen
c*  und keine Leerzeilen sind, werden als Kommentarzeilen fuer das Mutrack-
c*  Programm betrachtet und im Kopf des Logfiles ausgegeben.
c*
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Bleibt die rechte Seite einer Zuweisung frei (variable =    ), so wird der
c*  am Zeilenende in eckigen Klammern angegebene Defaultwert verwendet.
c* 
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Jeder Parameter des Abschnitts 'SCHLEIFEN-PARAMETER' kann einen, zwei
c*  oder mehrere aequidistante Werte annehmen. Die Eingabe geschieht
c*  entsprechend in einer der Formen
c*
c*		_parameter = wert
c*		_parameter = wert1    , wert2
c*		_parameter = wert_min , wert_max , wert_step.
c*
c*  Die Zuweisungen in den mit (*) markierten Zeilen des gleichen Abschnitts
c*  kommen nur zum Tragen, wenn sie nicht durch prioritaere Zuweisungen
c*  ueberdeckt werden. Fuer die Massen- und die Ladungsschleife geschieht dies
c*  durch Spezifizieren von 'artList' im Abschnitt 'PROJEKTILE' (vgl. unten),
c*  fuer die Energie-, Positions- und Winkel-Schleifen durch die Wahl einer
c*  Zufallsverteilung fuer die entsprechende Groesse im Abschnitt 'ZUFALLS-
c*  VERTEILTE STARTPARAMETER' ('random_energy', 'random_position' bzw.
c*  'random_winkel' > 0).
c*
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Wird die Variable 'artList' im Abschnitt 'PROJEKTILE' mit einer Komma-
c*  getrennten Liste aus
c* 
c*    m+, m-, Mu, Mu-, e+, e-, H+, H, H-, H2+, H2, H2-, alfa,
c*    A11+, A12+, A21+, A31+, A32+,   N11+, N21+,   K11+, K12+,   X11+, X12+
c*
c*  belegt (als String: z.B. 'm+,H+' ), so werden die Parameter 'Ladung' und
c*  'Masse' der Reihe nach mit den zugehoerigen Werten belegt.
c*  (m=Myon, Mu=Myonium, e=Elektron, H=Wasserstoff, alfa=alpha-Teilchen).
c*  Die zweite der beiden Zeilen gibt Kuerzel fuer Ionen schwererer Atome an:
c*  'Znm+' steht fuer einen "Cluster" aus n Atomen der Art 'Z' im Ladungszustand
c*  m+.
c*  A: Argon, N: Stickstoff, K: Krypton, X: Xenon.
c*
c*  Soll diese Moeglichkeit nicht genutzt werden, so ist 'artList' undefiniert
c*  zu belassen.
c*
c*  Das Einlesen der Teilchenliste geschieht CASE-SENSITIVE. Ein '!' an
c*  beliebiger Stelle der Teilchenliste bewirkt, dass alle Zeichen nach
c*  dem ! ignoriert werden.
c*  
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Befindet sich der Trigger im Strahlweg, so kann mit n_E_Verlust die Art
c*  der Beruecksichtigung des Energieverlustes in der Triggerfolie gewaehlt
c*  werden:
c*
c*	n_E_Verlust=0: kein Energieverlust
c*	n_E_Verlust=1: mittlerer Energieverlust gemaess 'mean_E_Verlust',
c*		       Breite der Energieverlustverteilung gemaess 'sigmaE'.
c*	n_E_Verlust=2: Funktioniert nur mit einfach positiv geladenen
c*		       Teilchen und wenn die Startflaeche entweder der Moderator
c*		       oder das 1. Gitter ist. Der Energieverlust wird hier
c*		       ueber Geschwindigkeitsskalierung aus Protonendaten
c*		       gewonnen (ICRU-Tabelle). Im so erhaltenen Energieverlust
c*		       ist auch die Wegstreckenverlaengerung in der Folie
c*		       bei nicht senkrechtem Einfall beruecksichtigt (Einfalls-
c*		       richtung!). Mit 'sigmaE' ungleich Null kann die Energie-
c*		       verlustverteilung zusaetlich verbreitert werden.
c*		       Bei gewuerfelten Startenergien fuer die Projektile
c*		       wird fuer die Berechnung des mittleren Energieverlustes
c*		       der Mittelwert von 'lowerE0' und 'upperE0' verwendet, es
c*		       sei denn, dass 'calculate_each'=.true. gesetzt wurde.
c*		       In diesem Fall wird der Energieverlust fuer jedes
c*		       Teilchen neu berechnet. Dies soll es ermoeglichen,
c*		       groessere Bereiche fuer die Startenergie vorzugeben und
c*		       trotzdem sinnvolle Energieverluste zu erhalten.
c*		       Allerdings kann dies die Rechenzeit u. U. stark erhoehen.
c*
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Wird das logical 'log_marker' im Abschnitt 'Graphikausgabe' gesetzt, so
c*  wird gegebenenfalls bei der Graphikdarstellung der Projektiltrajektorien
c*  das Trajektorienende mit einem Marker versehen, wobei der Markertyp vom
c*  Teilchenschicksal abhengt. Dabei gibt es aber momentan (Stand: Feb.96) noch
c*  Probleme, da vor allem fuer Teilchen, die auf der aktiven Flaeche des
c*  MCP2 aufschlagen, falsche Markertypen gewaehlt werden.
c*  Scheint ein spezielles Problem der CERN-library zu sein und wird bis auf
c*  weiteres nicht wieter verfolgt. Dieses feature wurde speziell fuer DEBUG-
c*  Taetigkeiten entwickelt und sollte ansonsten nicht von besonderer Bedeutung
c*  sein.
c*  Die Teilchenschicksale und die zugehoerigen Markertypen sind:
c*
c*      auf Rand des MCP2           : gefuellter Kreis
c*      im TD abgebrochen           : Stern (3 gekreuzte Linien)
c*      Spiegel durchquert          : gefuellter Kreis
c*      MCP2 getroffen              : leere Raute
c*      zerfallen                   : Schweizerkreuz
c*      Element verfehlt            : leerer Kreis
c*      reflektiert                 : leeres Dreieck
c*      aufgeschlagen               : gefuellter Kreis
c*      verloren (steps>maxsteps)   : gefuellter fuenfstrahliger Stern
c*      Zeitschritt kleiner dtSmall : leerer fuenfstrahliger Stern
c*
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Der Name der Ausgabedateien lautet 'MU_nnnn', wobei 'nnnn' eine durch-
c*  laufende Nummer ist. Die Ausgabedateien erhalten je nach Inhalt eine der
c*  Extensionen
c* 
c*	.LOG	 ->  Summary-file
c*	.NTP	 ->  PAW-Ntupel
c*	._???    ->  Statistik-Tabellen
c*	.PHYSICA ->  Gestamttabelle fuer Graphische Darstellung mit PHYSICA
c*	.PS	 ->  Postskript-Files mit Graphik der Bahnkurven. Bei diesen
c*		     Dateien wird der Dateiname um die Nummer der jeweiligen
c*		     Schleife, zu der die Graphik gehoert, erweitert. (Insofern
c*		     mehr als eine Schleife durchlaufen wird).
c*	.GEO	 ->  file mit der verwendeten Kammergeometrie.
c*	.FOIL	 ->  file mit Trajektoriendaten am Ort der Triggerfolie. Kann
c*		     von spaeteren MUTRACK runs eingelesen werden.
c*	.INFO	 ->  Info-file zu .FOIL mit den verwendeten Parametersettings.
c*
c*  Die Dateien werden im Directory 'mutrack$OUTdirectory' abgelegt.
c* 
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*  Das Achsensystem ist mit dem Ursprung in die Mitte des Moderators
c*  gelegt. Die x-Achse verlaeuft in Strahlrichtung, die z-Achse nach oben,
c*  die y-Achse so, dass ein rechtshaendiges System besteht.
c*  Nach dem Spiegel wird das Koordinatendreibein so gedreht, dass die x-Achse
c*  wiederum in Strahlrichtung weist. Der Ursprung wird dabei in das Zentrum
c*  des Doppelkreuzes der Spiegelaufhaengung verschoben.
c* 
c*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
c*
c*	*****************************************************************
c*	*			     Die Zeilen				*
c*	* $PARAMETER_LISTE						*
c*	*				und				*
c*	* $END								*
c*	*	         MUESSEN IN SPALTE ZWEI BEGINNEN!		*
c*	*****************************************************************
c*	*  26.09.1995			>>  anselm.hofer@psi.ch		*
c*******************************************************************************
