TWIST experiment GEANT Monte Carlo (UNITs/DEFAULTs are in square brackets) C LIST C C **** Geant FFKEYs: see GEANT manual for more details **** C C ========== RUNG: IDRUN IDEVT ========== C == IDRUN == User run number [1] C == IDEVT == User event number [0] C RUNG 1 C C ========== TIME: TIMINT TIMEND ITIME ========== C == TIMINT == Time used for initialization C NOTE: FFCARD input for TIMINT is ignored/overwritten C == TIMEND == Time required for termination [100 seconds] C == ITIME == Test every ITIME events C NOTE: User must optimize TIMEND/ITIME so that ITIME is C as large as is save! - Program termination is C initiated as soon as the time left on a particular C queue is smaller than TIMEND. TIME 0.0 100. -1 C C ========== TRIG: NEVENT ========== C C == NEVENT == Number of events to be processed C TRIG 999999 C TRIG 4938 C C ========== NETT: NETT ========== C C == NETT == Number of Events To Tape [0] C The GEANT run is terminated with IEORUN when C either: C NETT > 0 .and. Number of Events To Tape = NETT C or: C IEVENT = NEVENT (set by TRIG above) C which ever happens first. C NETT 0 C C ========== RLUX: LUX ISEED ========== C C The RANLUX random number generator C For more information, see: C http://wwwinfo.cern.ch/asdoc/shortwrupsdir/v115/top.html C C LUX: RANLUX Luxury Level C Higher means C ISEED: RANLUX initial seed (when seeded with only one INTEGER) C Note that *any* initial seed should be good. C ISEED=0 means "use RUNG as the seed". C RLUX 3 0 C C ========== RNDM: NRNDM(1) NRNDM(2) ========== C C NRNDM(1) and NRNDM(2) are used for restarting the generator C from a break point. RANLUX will then skip over C NRNDM(1) + 10**9*NRNDM(2) numbers to reach the break point C again (slow). C C Initial value of RANLUX extra seeds NRNDM(1), NRNDM(2) C in call to RLUXGO(Lux,ISEED,NRNDM(1),NRNDM(2)) C C To reproduce an event from a "blind analysis" run made with C micheld server, set IDECAY in the SRVU ffcard in addition to the C random seeds. RNDM 0 0 C C ========== RLXn: IVEC(25) ========== C C You can also use these to restart the generator from a break C point, using 25 32-bit integers. More efficient than RNDM, C but less convenient. C C Initial state of RANLUX random number generator C C RLX1 3581912 5574638 6409493 1188367 13993016 C RLX2 10082780 8942804 10259692 6705902 13910431 C RLX3 5182377 8365843 981434 10328133 5335947 C RLX4 12452237 5251499 8980566 6033477 15107937 C RLX5 9138526 1822325 15433747 1746609 3071024 C C ========= SORD: ISTORD ========= C C Stack ordering is required for mu->e g nu nu in order for the decay C gamma to be tracked after all other particles in the simulation; C i.e. daughters of muons are tracked in the order they are generated. C C 0 = Stack ordering is not activated C 1 = Stack ordering is activated (E614GEANT default) C C See GEANT manual BASE040-2 and TRAK310-1 C SORD 1 C C ========= AUTO: IGAUTO ========= C C 1 = Automatic computation of STMIN,STEMAX,DEEMAX,TMAXFD (default) C (except EPSIL) C 0 = Tracking media parameters taken from the argument list of GSTMED C (unless parameters are < 0 ) C AUTO 0 C C ========== PHYS: IPHYS ========== C C 0 = all physics processed off C 1 = physics processes switched on/off via GEANT cards (default) C C PHYS 0 C C ========== ANNI: IANNI ========== C C 0 = no positron annihilation effect C 1 = positron annihilation with generation of secondaries (default) C 2 = same without generation of secondaries C C ANNI 0 C C ========== BREM: IBREM ========== C C 0 = no Bremsstrahlung effect C 1 = Bremsstrahlung with generation of secondaries (default) C 2 = same without generation of secondaries C C BREM 0 C C ========== COMP: ICOMP ========== C C 0 = no Compton scattering C 1 = Compton scattering with generation of secondaries (default) C 2 = same without generation of secondaries C C COMP 0 C C ========== DCAY: IDCAY ========== C C 0 = no particle decays C 1 = unstable particle decay with generation of secondaries (default) C 2 = same without generation of secondaries C C DCAY 0 C C ========== DRAY: IDRAY ========== C C 0 = no delta rays effect C 1 = delta rays with generation of secondaries (default) C 2 = same without generation of secondaries C Note: DRAY 1 is only possible for reduced Landau fluctuations C (LOSS 1). When full Landau fluctuations (LOSS 2) then C IDRAY = 0 regardless of the setting here. C C DRAY 0 C C ========== HADR: IHADR ========== C C 0 = no hadron interactions effect C 1,4,5 = hadron interactions with generation of secondaries C 2 = same without generation of secondaries C C GHEISHA hadronic shower code if IHADR = 1 C FLUKA hadronic shower code if IHADR = 4 (default) C FLUKA/MICAP had. shower code if IHADR = 5 C C HADR 0 C C ========== LOSS: ILOSS ========== C C 0 = no energy loss effect C 1 = delta ray and reduced Landau fluctuations C 2 = full Landau fluctuations and no delta rays (default) C 3 = same as 1 C 4 = average Energy loss and no fluctuations C LOSS 1 C LOSS 2 C LOSS 0 C C ========== STRA: ISTRA ========== C C 0 = Urban model for energy loss for thin layer (default) C 1 = PAI model " " " " " " C 2 = ASHO model for 1e+e- pair production C 1 = pair production with generation of secondaries (default) C 2 = same without generation of secondaries C C PAIR 0 C C ========== PHOT: IPHOT ========== C C 0 = no photo-electric effect C 1 = photo-electric effect with generation of secondaries (default) C 2 = same without generation of secondaries C C PHOT 0 C C *** The ENERGY RANGE of the cross section and energy loss tables can C be fixed by the user with the data card : C 'ERAN' EKMIN EKMAX NEKBIN C which defines nkbin bins from Ekmin to Ekmax in a logarithmic scale C The default is, 196 bins from 10 KeV to 100 MeV kinetic energy but C in logarithmic scale. NEKBIN must be 50 ISWIT(3) are active only when DEBU is on C C == ISWIT(1) = 2: the content of the temporary stack for secondaries in C the common /GCKING/ is printed; C == ISWIT(2) = 1: the current point of the track is stored in the JDXYZ C bank via the routine GSXYZ; C = 2: the current information on the track is printed via C the routine GPCXYZ; C = 3: the current step is drawn via the routine GDCXYZ; C = 4: the current point of the track is stored in the JDXYZ C bank via the routine GSXYZ. When the particle stops C the track is drawn via the routine GDTRAK and the C space occupied by the track in the structure JDXYZ C released; C = 5: print GEANT vertex information via GPVERT at the end C of the event (in GUOUT) C == ISWIT(3) = 1: the current point of the track is stored in the JDXYZ C bank via the routine GSXYZ; C == ISWIT(4) = 0: no input RAYFILE C 1: read input RAYFILE C > 1: start reading at ISWIT(4)th input ray C == ISWIT(5) = 1: RAYFILE is M13GEANT format C = 2: RAYFILE is REVMOC format C = 3: RAYFILE is TECTRCK format C == ISWIT(8) = 0: Batch version C = 1: Interactive version C == ISWIT(10)= 0: no digitization performed (for fast execution) C = 1: writes simulated TDC data to disk file C 2: writes ASCII data to disk file C 3: writes simulated TDC data and ASCII data to disk files C C SWIT 0 2 0 1 1 0 0 0 0 3 SWIT 0 0 0 0 0 0 0 0 0 1 C C ========== HSTA: LHSTA C == LHSTA == NHSTA names of required standard histograms C C 'TIME' provides histogram of CPU time per event C C HSTA 'TIME' 'SIZE' 'MULT' 'NTRA' 'STAK' C C ========== PRIN: LPRIN C == LPRIN == NPRIN names of GEANT data structures to be printed C C PRIN 'PART' 'MATE' 'TMED' 'VOLU' 'SETS' C C ========== RGET: LRGET C == LRGET == NRGET names of GEANT data structures to fetch from RZ C C RGET 'INIT' C C ========== RSAV: LRSAV C == LRSAV == NRSAVE names of GEANT data structures to fetch from RZ C C RSAV 'INIT' C C ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ C C ***** E614GEANT Run directives ***** C C ========== MAXS: MAX_STEP MAX_LENGTH C C max_step: maximum number of steps for a track [1000] C max_length: maximum path length [4000 cm] C MAXS 60000 6000. C C ========== VCIN: IN_CUT(10) ========== C C If the first argument is 'ALL ', the card applies to all particles; C if not, it applies to only secondary (including decay) particles. C If such a particle enters one of these volumes, stop tracking it. C C in_cut: volume name where the track is cut when entering C example: 'HRNG' - a list of up to nine choices C VCIN 'SEC ' 'HRNG' 'GRFS' 'HSCV' 'HFNS' 'HEBG' 'EVOL' C C ========== VCOT: OUT_CUT(10) ========== C C If the first argument is 'ALL ', the card applies to all particles; C if not, it applies to only secondary (including decay) particles. C If such a particle leaves one of these volumes, stop tracking it. C C out_cut: volume name where the track is cut when exiting C example: 'HCTR' - a list of up to nine choices C C VCOT 'ALL ' 'HCTR' C C ========== STPL: ISTEP STEPLIM(9) ==> Step length limiting parameters C C ISTEP = -1 will cause the target tracking to be the same as in all C other materials, i.e. as given for the first 5 steplim C C STEPLIM C 1st-5th: tmaxfd, stemax, deemax, epsil and stmin for ALL materials C [-1.] [-0.5] [-0.25] [.0001] [-.1] C 6th-9th: stemax, deemax, epsil and stmin for 'TARGET' material C [.001] [-0.25] [.0001] [-.1] C STPL 0 -1.0 -0.5 -0.25 0.0001 -0.1 -0.1 -0.0001 0.0001 0.0001 C C ========== BFLD: JFIELD FIELD_MAX B_THETA B_PHI C DETECT_POS(1) DETECT_POS(2) DETECT_POS(3) ========== C C Properties of the magnetic field. C C jfield: 0 -> no magnetic field tracking (default) C 1 -> tracking performed with GRKUTA C 2 -> tracking performed with GHELIX C 3 -> tracking performed with GHELX3 C C From the GEANT code/manual: C C GRKUTA for inhomogeneous fields C GHELIX for quasi-homogeneous fields tilted w.r.t. the reference C GHELX3 for one-component fields along the z axis C C GRKUTA and GHELIX call the user subroutine GUFLD where the C components of the field at the given point are computed. C GHELX3 takes the value of the field in the tracking medium C parameter FIELDM. C C field_max: maximum field [20.0 kGauss] C C For jfield == 1 C C b_theta: polar angle of detector wrt B-field [deg] C b_phi: azimuthal angle of detector wrt B-field about z-axis [deg] C C For jfield == 2 C C b_theta: polar angle of B-field relative to z-axis [deg] C b_phi: azimuthal angle of B-field about z-axis [deg] C C These angles are used to establish the B-field components as -- C C Bz = field_max*cos(degrad*b_theta) C Bx = field_max*sin(degrad*b_theta)*cos(degrad*b_phi) C By = field_max*sin(degrad*b_theta)*sin(degrad*b_phi) C C in the case of jfield .eq. 2 C C For a uniform field along the z-axis use: BFLD 3 20.0 C Note: b_theta/b_phi are ignored for jfield .eq. 3 !!! C C detect_pos: x/y/z position of the detector wrt. B-field map center C BFLD 1 20.0 0.0 0.0 0.0 0.0 0.0 C BFLD 2 20.0 10.0 5.0 0.0 0.0 0.0 C setting for uniform 2 Tesla field: C BFLD 3 20.0 0.0 0.0 0.0 0.0 0.0 C C ========== FSTR: NSTRA M_STRA V_STRA C C NSTRA (max.12!) Number of media with specified non-default thin C layer option STRA C M_STRA(12) - List of tracking media numbers for thin layer option C V_STRA(12) - Value of thin layer option (must be REAL!); C C i.e. v_stra(i) goes with m_stra(i) C in the CALL GSTPAR(m_stra(i),'STRA',v_stra(i)) C C NOTE: What's required is the TRACKING MEDIA number (not the MATERIAL C number). So, for example, HE-GAS (material) is 30, but C 'HE IN MAGNET' is tracking medium 16. It is 16 that needs to C go into M_STRA. HOWEVER (!!!! NOTE !!!!) if you choose a C TRACKING MEDIA which refers to the same MATERIAL as another C TRACKING MEDIA, you MUST list all of them below, or else the C program will bomb!!! C C FSTR 7 30 31 33 36 15 16 19 0 0 0 0 0 1. 1. 1. 1. 1. 1. 1. 0. 0. 0. 0. 0. C C ======== ANCT: IANCT CTMIN CTMAX ======== C C ANgle CuT on cos(theta) of positron decay spectrum; especially C for eliminating tracks outside our acceptance that take a long C time to simulate. C C CTMIN and CTMAX are the limits on cos(theta) C [defaults: -0.2 and 0.2] C C IANCT: 0 or negative -> No cut. [default] C 1 -> Produce no positrons between CTMIN and CTMAX. C 2 -> Produce no positrons outside CTMIN and CTMAX. C ANCT 0 -0.2 0.2 C C ========== PMLT: NPRT(1) NPRT(2) ========== C C This card specifies the exact number of initial particles C (including the trigger particle): C >0 produces a fixed multiplicity C =0 statistical multiplicities according to the values C in the RATE card C C NOTE: nprt(1) < 0 is not allowed (that is to say, you don't even get C a trigger particle and hence no events) C C PMLT gives the initial particle multiplicity [1 -1] C when nprt(1) == 1 and nprt(2) < 0 no pile up C when nprt(1) == 0 and nprt(2) == 0 statistical multiplicity C when nprt(1) >= 0 and nprt(2) < 0 only muon pile up C when nprt(1) == 1 and nprt(2) >= 0 only positron pile up C PMLT 0 0 C C ========== RATE: RATE_MUON RATE_POSITRON TCOIN_PLUS TCOIN_MINUS ========== C C The 'rate' is the mean number of particles per second C The 'coincidence window' defines the trigger position relative C to the beginning and end of the event gate C C rate_muon: surface beam muon rate [MHz] C rate_positron: beam positron rate [MHz] C tcoin_plus: coincidence window after trigger [us] C tcoin_minus: coincidence window before trigger [us] C RATE 0.0022 0.021 10. -7. C C ========== MDCT: DECAY_TIME ========== C C The minimum trigger muon decay time in the simulation (in us) C C decay_time: minimum muon decay time [0 us] C C MDCT 1.0 C C ========== PGUN: IPGUN PGUN(13) ========== C C Particle GUN -- basically a particle beam, from any source C position, at any angle. C PGUN must be used in conjunction with one of the following cards: C MSOR, PNCL, SLGT, CONE C PGUN controls elements common to all these. C C Note: The element PGUN(11) controls the momentum bite; C PGUN(11)=1.007 means a 0.7% momentum bite. C C You can now have two different particle types for pile-up. C You specify the particle types on the PGUN line, then list C the PGUN properties on the following two lines. All the C various particle source cards (MSOR, etc) now have two lines, C in the same manor, corresponding to the particles specified C here. You can turn on different sources for each particle. C C FYI, particle type 65 is a mu+ with spin; type 2 is a e+. C C IPGUN = GEANT Particle type [65] (default) C PGUN( 1) = Particle gun origin x [0 cm] C PGUN( 2) = Particle gun origin y [0 cm] C PGUN( 3) = Particle gun origin z [-142.5 cm] C PGUN( 4) = Particle gun aim - rotation around x-axis [0 deg] C PGUN( 5) = Particle gun aim - rotation around y-axis [0 deg] C PGUN( 6) = Particle gun aim - rotation around z-axis [0 deg] C PGUN( 7) = sigma of x position [1 cm] C PGUN( 8) = sigma of y position [1 cm] C PGUN( 9) = sigma of z position [0 cm] C PGUN(10) = flat x width [0 cm] which is then smeared with sigma_x, only for IMSOR==3 C PGUN(11) = flat y width [0 cm] which is then smeared with sigma_y, only for IMSOR==3 C PGUN(12) = Particle momentum [29.79 MeV/c] C PGUN(13) = 1 +- deltaP/P [1.0 -> deltaP/P = 0] C C September 2002, info from Glen C FWHM_x = 27mm => sigma_x = 27mm/2.35 = 11.5mm C FWHM_y = 18mm => sigma_y = 18mm/2.35 = 7.7mm C PGUN 65 2 0. 0. -142.5 0. 0. 0. 1.15 0.77 0. 0. 0. 29.6 1.007 0. 0. -142.5 0. 0. 0. 1.15 0.77 0. 0. 0. 29.6 1.007 C PGUN (65) 0. 0. -153.0 0. 0. 0. 0.45 0.5 0. 0.9 0.9 29.79 1.0 C PGUN ( 2) 0. 0. 0. 0. 0. 0. 0. 0. 0. 0.9 0.9 52.83 1.0 C C ========== MSOR: IMSOR PSOR(9) ========== C C surface Muon SOuRce generator C C It is possible to specify some properties of the muon beam at C the focal point (with PGUN): i.e where the focal point is, the C overall aim of the beam, the beam cross section, particle momentum C and flat momentum spread, but launch these rays at the PSOR(1) C location, z0, with some flat distribution PSOR(2) around z0. The C rectangular angle spread in the horizontal and vertical direction, C in the absence of a magnetic field, is defined by PSOR(3)(4), and C a flag exists which eliminates all muon momenta above the surface C muon edge, PSOR(9). C C IMSOR = Use-Flag: any positive integer to use [1] C =1 : uniform direction into rectangular aperture, C gaussian beam spot from PGUN. C =2 : gaussian distributed direction, C gaussian beam spot from PGUN. C =3 : direction in combo with vertex(not at focus), C position and direction are correlated. C PSOR( 1) = Beam origin position in z [-300 cm] C PSOR( 2) = flat half-width of track start pos. [0 cm] C If IMSOR = 1 C PSOR( 3) = horizontal divergence [15 mrad] C PSOR( 4) = vertical divergence [10 mrad] C If IMSOR = 2 C PSOR( 3) = horizontal opening angle sigma [15 mrad] C PSOR( 4) = vertical opening angle sigma [10 mrad] C If IMSOR = 3 C PSOR( 3) = horizontal opening angle sigma [15 mrad] C PSOR( 4) = vertical opening angle sigma [10 mrad] C PSOR( 5) = error function width factor in dx [0 cm] C PSOR( 6) = error function width factor in dy [0 cm] C PSOR( 7) = x to dx correlation factor(not at focus) [0] C PSOR( 8) = y to dy correlation factor(not at focus) [0] C C All cases C PSOR( 9) = Surface muons: no = 0, yes > 0 (default) C This used to be PSOR( 5) C C September 2002, info from Glen: C From measurements with pion beam: FWHM_x = 15mrad, FWHM_y=18mrad C MSOR 1 1 -300. 0. 15. 18. 0. 0. 0. 0. 1. -300. 0. 15. 18. 0. 0. 0. 0. 0. C MSOR -182. 0. 13. 15. -0.009 -0.009 1. C MSOR -174.450836 0. 13. 15. -0.009 -0.009 1. C C ========== PNCL: IPNCL PNCL(2) ========== C C PeNCiL ray generator C C Origin is defined by PGUN. Direction specified here is relative C to that in PGUN; all particles produced are in the same direction. C C IPNCL = Use-Flag: any positive integer to use [OFF] C PNCL( 1) = RAYTRACE theta [mrad] C PNCL( 2) = RAYTRACE phi [mrad] C PNCL -1 -1 15. 10. 15. 10. C C ========== SLGT: ISLGT SLGT(1) ========== C C Search LiGhT beam generator C C A bunch of rays inside a cone (NOT a double-sided cone). C Cone's origin and central axis are defined by PGUN. C C ISLGT = Use-Flag: any positive integer to use [OFF] C SLGT( 1) = circular emittance [deg] C SLGT -1 -1 20.0 0.1 C C ========== CONE: ICONE CONE(2) ========== C C Rays on the surface of a CONE with opening angle theta, in a C range of +-delta_theta/2. The cos of the angle is sampled flat C within the specified limits. C C ICONE = Use-Flag: any positive integer to use [OFF] C CONE( 1) = theta [deg] C CONE( 2) = delta_theta [deg] C CONE -1 -1 -60. 0 2. 1. C C ========== CON2: ICON2 CON2(2) ========== C C Rays on the surface of a CONE with opening angle theta, in a C range of +-delta_theta/2. The *angle* (not the cosine) C is sampled flat within the specified limits. C C ICON2 = Use-Flag: any positive integer to use [OFF] C CON2( 1) = theta [deg] C CON2( 2) = delta_theta [deg] C CON2 -1 -1 90. 180. 90. 180. C C ========== COSM: CPART CYWIND CZWIND1 CZWIND2 CXWIND1 CXWIND2 ========== C C COSMic Ray generator C C mu+ and mu- (in ratio 1.25:1) uniformly distributed in x and z C within the specified window. See source file "cosang.F" for C details on momentum distribution. C C cpart: any positive integer to use cosmic generator [OFF] C cywind: 'window' height above detector (y) in cm C czwind: the z coord of ends of the rectangular 'window' C cxwind: the x coord of ends of the rectangular 'window' C COSM -1 130.8 -146.5 146.5 -132.0 132.0 C C ========== ROOM: P_ATM TEMP_C ========== C C Atmospheric pressure == P_atm [760 Torr] C Room temperature == temp_C [20 deg C] C ROOM 760.0 20.0 C C ========== GABS: IGAS GAS_MIXTURE ========== C C Properties of the gas degrader ("absorber") C C Identity of Gas == igas = 20 -> He/CO2 (default) C = 21 -> He/Ar C = 22 -> He/Xe C Gas mixture ratio == gas_mixture C C Note: This is the fraction by VOLUME of the heavier gas [0.5]; C (converted internally to fraction by weight) C GABS 20 0.42 C C ========== PABR: IPLASTIC PLAS_LNGTH PLAS_POSITN ========== C C Properties of the plastic degrader ("absorber") C C iplastic == 0 -> no plastic degrader C == 31 -> Mylar (default) C == 26 -> Scintillator C plas_thkns == [0.004 cm] (full length) C plas_side == [16 cm] (full length) C PABR 0 0.004 8.0 C PABR 31 0.004 8.0 C C ========== HEBG: IBAG_GAS ========== C C Gas in Helium bag between end of beampipe and detector volume C C ibag_gas == 13 --> magnetic air (default) C 16 --> magnetic helium C HEBG 13 C C ========== TECM: ITECM P_TECM NCEL_TEC TEC_CWDT TEC_GAP TEC_WDST C TEC_WDTH TEC_DIST TEC_WALL TEC_WNDW TEC_RIM C TEC_POSITN ======== C C Time Expansion Chamber properties C C itecm == 0/>0 means do not/do insert a TEC [0] C P_tecm == pressure in TEC [20 Torr] C ncel_tec == number of sense wires in each TEC section C tec_cwdt == separation of sense wires [cm] C tec_gap == gap between TEC sections [cm] C tec_wdst == distance between TEC sections and entry/exit window [cm] C tec_wdth == ion drift full width [cm] C tec_dist == distance: sense wires, field wires and TEC [cm] C tec_wall == aluminized mylar wall thickness [cm] C tec_wndw == mylar entrance/exit window thichness [cm] C tec_rim == dist.: TEC drift volume & end of unit [cm] C tec_positn == position of TEC (global z) [cm] C C TECM 1 20.0 48 0.1 6.0 2.5 2.0 0.3 0.01 0.000635 3.5 -143.2 C TECM 1 20.0 48 0.1 6.0 2.5 2.0 0.3 0.01 0.000635 3.5 -174.5 C C ========== TGRW: ITEC_GRDW NGRDW_TEC TEC_GRDW ========== C C Time Expansion Chamber Foil Grid C C itec_grdw == 0 -> no TEC Foil Grid (default) C 1 -> Aluminized Mylar Foil C 2 -> actual Foil Grid wires C ngrdw_tec == Number of foil grid wires C tec_grdw == TEC Foil Grid thickness/wire diameter C C TGRW 2 20 0.001 C C ========== TRSL: TEC_SIGMA0 TEC_SIGMA1 TEC_EFF ========== C C Time Expansion Chamber drift resolution C C tec_sigma0 == resolution next to field wire [cm] (default - 0.00626cm) C tec_sigma1 == dependence of resolution on sqrt(distance from field wire) C (default: tec_sigma0 = 0.00626cm, tec_sigma1 = 0.00701sqrt(cm)) C tec_eff: == efficiency of TEC [1.0] C C TRSL 0.00626 0.00701 1.0 C C ========== GCAP: N_ADMIX ========== C C n_admix == Nitrogen admixture (by VOLUME) in Helium gas C between chamber planes [0.02 -> 2%] C C September 2002, info from R.Openshaw: C There is about 2.7% N2 and 0.3% air in the He volume C GCAP 0.03 C C ========== YOKE: IYOKE ========== C C iyoke == 0 0/>0 means do not/do insert a YOKE [0] C C YOKE 1 C C ========== DSDK: IDSM DOWNSTREAM(1) DOWNSTREAM(2) DOWNSTREAM(3) ========== C C idsm == 0 -> no plastic disk downstream (default) C == 31 -> Mylar C == 26 -> Scintillator etc. C C downstream(1) == Plastic Disk Radius [cm] C downstream(2) == Plastic Disk Thickness [cm] C downstream(3) == Plastic Disk Global Position [cm] C C DSDK 26 18.5 0.5 100.0 C C ========== PIBR: IPIBR ========== C C Pion decay branch C C ipibr == 1 : pi+ -> mu+ nu mu has no spin (default) C == 2 : pi+ -> e+ nu C == 3 : pi+ -> mu+ n mu has spin C PIBR 1 C ========== MUBR: MDCAY_MOD ========== C C Muon decay branch C C mdcay_mod == 1 : mu+ -> e+ nu nub (Michel Decay - default) C == 2 : mu+ -> e+ gamma nu nub (rad. muon decay) C C MUBR 2 C C ========== MCHL: MICHEL_RHO MICHEL_DELTA MICHEL_XSI MICHEL_ETA ========== C C Set Michel parameters C C 1st through 4th entries: rho [0.75], delta [0.75], xi [1.0], eta [0.0] C MCHL 0.75 0.75 1.00 0.00 C C ========== SRVU: use micheld server ============= C ISRVUSE,ISPECTRUM,IBRUN,IDECAY C C ISRVUSE>0 to get a sample of Michel decays from a micheld server [0]. C C ISPECTRUM request the given spectrum number [0]. C C IBRUN: Base run number [0]. The sample number used is (RUNG-IBRUN). C Samples are numbered from zero within each spectrum. C C IDECAY is the initial offset in the sample buffer [0]. C To reproduce an event, set IDECAY as printed out by C photo in addition to setting random seeds in the RNDM card. C If the original run had default IDECAY==0, you can also C calculate IDECAY value to reproduce an event of that run C as (IEVENT-1), where IEVENT is printed out by GTRIGI. C SRVU 1 9 9000 0 C C ========== SRVC: micheld connection parameters ============= C ip1 ip2 ip3 ip4 port ntries delay_min delay_max C C Server address is: ip1.ip2.ip3.ip4:port C The default is 142.90.100.68:3469, that is, lin00.triumf.ca:3469 C C Up to ntries [1] connection attempts are made. The time interval C between two connection attempts starts at (randomized) delay_min [5s] C and doubles (with randomization) till it gets to delay_max [1000s], C then it stays around delay_max. C SRVC 142 90 100 68 3469 8 30 1000 C C ========== PMU0: polarization of muon at the origin ============= C C PMU0 sets polarization of muon with respect to its momentum in C initial event kinematics. The card has no effect on other muons, C e.g. those coming from pion decay. Default PMU0 is -1. C C PMU0 0.0 C C ========== DPOL: IF_DPOL TLIFE_DPOL ========== C C Set muon depolarization parameters C C if_dpol: depolarization > 0 = on, 0 or < 0 = off [0] C tlife_dpol: polarization mean life [0.0000022 sec] C DPOL 0 0.0000022 C C ========== DCIC: NCL_DC(2) STEPCL_DC(2) SLGTCL_DC(2) C C DC ion cluster simulation: C C ACL_DC: Number of clusters needed for hit [1] C STEPCL_DC: The ion cluster separation [0.03cm] C SLGTCL_DC: The ion cluster time length [10ns] C C (1) is for muon; (2) is for positron/electron C DCIC 1.0 1.6 0.0180 0.0140 100.0 100.0 C C ========== PCIC: ACL_PC(2) STEPCL_PC(2) SLGTCL_PC(2) C C PC ion cluster simulation: C C ACL_PC: Number of clusters needed for hit [1] C STEPCL_PC: The ion cluster separation [0.03cm] C SLGTCL_PC: The ion cluster time length [10ns] C C (1) is for muon; (2) is for positron/electron C PCIC 1.0 1.6 0.0180 0.0140 250.0 60.0 C C ========== DCEF: EFFS_DC_ON ========== C C Set flag to include DC chamber efficiency C C effs_dc_on == 0/>0 means do not/do invoke DC chamber efficiency [0] C DCEF 0 C C ========== DCDZ: DEAD_ZONE_DC_ON DEAD_ZONE_DC_FUZZ ========== C C dead_zone_dc_on == 0/>0 means do not/do invoke dead-zone C due to a through going muon in DC efficiency [0] C dead_zone_dc_fuzz == gives the extended dead region beyond the C muon track [0.05 cm] C DCDZ 0 0.05 C C ========== PCEF: EFFS_PC_ON ========== C C Set flag to include PC chamber efficiency C C effs_pc_on == 0/>0 means do not/do invoke PC chamber efficiency [0] C PCEF 0 C C ========== PCDZ: DEAD_ZONE_PC_ON DEAD_ZONE_PC_FUZZ ========== C C dead_zone_pc_on == 0/>0 means do not/do invoke dead-zone C due to a through going muon in PC efficiency [0] C dead_zone_pc_fuzz == gives the extended dead region beyond the C muon track [0.05 cm] C PCDZ 0 0.05 C C ========== DCRS: RESO_ON S0 DSDT ========== C C Set the drift chamber resolution funciton: C C sigma = sigma0 + d(sigma)/dT * T - where T is the drift time C C reso_on == 0/>0 means do not/do invoke DC resolution [0] C s0 == sigma0 [5.0 ns] C dsdT == d(sigma)/dT [0.01] C DCRS 1 3.0 0.0 C C ========== BOUT: M_FBU M_MCTR M_MCT2 M_MCSP M_MCS2 M_MCS3 M_MCWI M_MUSR ========== C C Set which event banks are written to YBOS file C C Digitized data (simulates Fastbus output): C M_FBU > 0 encodes the FBU1 and FBU2 banks in 'store_fbu' [1] C Monte Carlo "Truth" banks: C M_MCTR > 0 stores the MCTR bank in 'store_mctr' [1] C (Contains particle starting parameters) C M_MCT2 > 0 stores the MCT2 bank in 'store_mct2' [0] C (Contains both starting and stopping parameters) C Information about detector plane hits: C M_MCSP > 0 store sp_data into the MCSP bank in 'store_mcsp'[0] C M_MCS2 > 0 store sp_data into the MCS2 bank in 'store_mcs2'[0] C (MCS2 also contains momentum at each point.) C M_MCS3 > 0 store sp_data into the MCS3 bank in 'store_mcs3'[1] C (MCS3 also contains GEANT step point) C Information about wire hits (not used?). C M_MCWI > 0 encodes the MCWI bank in 'store_mcwi' [0] C M_MUSR > 0 produces the MUSR users MC bank in store_musr [0] C C Disable the space points bank BOUT 1 0 1 0 0 0 0 0 C Enable the space points bank C BOUT 1 0 1 0 0 1 0 0 C C ========== UCUT: UCUT(1) UCUT(2) ========== C C Set the minimum kinetic energy required of secondary to enter C into the JVERT and JKINE GEANT bank structure and MCTR bank C C UCUT(1) == Emin for gamma [MeV] (default - 1MeV) C UCUT(2) == Emin for electron [MeV] (default - 1MeV) C UCUT 1.0 1.0 C STOP