\START MUCOOL_RCP
\SIZE 166 81
!----------------------------------------------------------------------
! Name: MUCOOL.RCP
! Package:
! Library:
! Purpose: DEFINES MUCOOL FOR Muon COoling Experiment
! Version
! Created: 26-Jan-01 Rajendran Raja and Zafar Usubov
!----------------------------------------------------------------------
!----------------------------------------------------------------------
! Bend dipole magnetic field definitions
!----------------------------------------------------------------------
\ARRAY DIPOLE_BY_NUMERATORS
! 2 4 6 8 POWERS OF EXPANSION
 -1.0 25.0 -45.0 7605.0
\END
\ARRAY DIPOLE_BY_DENOMINATORS
! 2 4 6 8 POWERS OF EXPANSION
 8.0 384.0 1024.0 229376.0
\END
\ARRAY DIPOLE_BR_NUMERATORS
 ! 2 4 6 8 POWERS OF EXPANSION
 -5.0 15.0 -585. 14365.
\END
\ARRAY DIPOLE_BR_DENOMINATORS
! 2 4 6 8 POWERS OF EXPANSION
 24.0 128.0 14336.0 458752.0
\END
DIPOLE_BEND_RADIUS 52.0 !IN CENTIMETERS PHYSICAL R OF DIPOLE.
DIPOLE_BZERO -14.5485 !B0 FIELD OF DIPOLE IN KILOGAUSS !! -14.53 Old f, 14.5485 best
!DIPOLE_BZERO -14.53 !B0 FIELD OF DIPOLE IN KILOGAUSS !! -14.53 Old f, 14.5485 best
CONSTANT_DIPOLE_FIELD F !IF TRUE WILL HAVE CONSTANT FIELD IN DIPOLE
LONG_SECTION_SCALE_FACTOR 1.0 !WILL SCALE CANONICAL FIELD BY THIS FACTOR
CONSTANT_LONG_FIELD 0.0 !WILL CREATE CONSTANT FIELD IF NON ZERO
SHORT_SECTION_SCALE_FACTOR 1.0 !WILL SCALE CANONICAL FIELD BY THIS FACTOR
CONSTANT_SHORT_FIELD 0.0 !WILL CREATE CONSTANT FIELD IF NON ZERO
!----------------------------------------------------------------------
! rf PARAMETERS
!----------------------------------------------------------------------
RF_FREQUENCY 201.0E6 !IN Hz
! RF_GRADIENT 15.0E-5 !IN GIGA_VOLTS/CM
 RF_GRADIENT 12.704E-5 !IN GIGA_VOLTS/CM NEw value for new 
! RF_GRADIENT 12.702E-5 !IN GIGA_VOLTS/CM NEw value for new good !!
!RF_GRADIENT 0.0 !IN GIGA_VOLTS/CM
! rf_grad 12.702(31,29 for pt), 12.703(30,32) 
\ARRAY RF_TIMES
! times in NANOSECONDS of entry into rf in the 64 rf
! modules starting with RFCV copy 1-16 in RBY4 copy 1-4
!TARGET_TIME 0.11700E+02 !time of target in nsec 
 0.00000E+00 0.11962E+01 0.23923E+01 0.35885E+01
 0.47847E+01 0.59808E+01 0.71770E+01 0.83731E+01
 0.14998E+02 0.16194E+02 0.17390E+02 0.18587E+02
 0.19783E+02 0.20979E+02 0.22175E+02 0.23371E+02
 
 0.33974E+02 0.35170E+02 0.36366E+02 0.37562E+02
 0.38759E+02 0.39955E+02 0.41151E+02 0.42347E+02
 0.48972E+02 0.50168E+02 0.51364E+02 0.52560E+02
 0.53757E+02 0.54953E+02 0.56149E+02 0.57345E+02
 
 0.67948E+02 0.69144E+02 0.70340E+02 0.71536E+02
 0.72732E+02 0.73929E+02 0.75125E+02 0.76321E+02
 0.82946E+02 0.84142E+02 0.85338E+02 0.86534E+02
 0.87731E+02 0.88927E+02 0.90123E+02 0.91319E+02
 
 0.10192E+03 0.10312E+03 0.10431E+03 0.10551E+03
 0.10671E+03 0.10790E+03 0.10910E+03 0.11029E+03
 0.11692E+03 0.11812E+03 0.11931E+03 0.12051E+03
 0.12170E+03 0.12290E+03 0.12410E+03 0.12529E+03
!
\END
!!TARGET_TIME 11.782 !time of target IN NANOSECONDS 
!TARGET_TIME 11.6979 !time of target IN NANOSECONDS 
TARGET_TIME 11.7000 !time of target IN NANOSECONDS 
! (currently MABS middle) wrt to 1st rf.
!RF_PHASE 210.0 !phase of rf at which to accelerate in degrees
RF_PHASE 30.0 !phase of rf at which to accelerate in degrees
!----------------------------------------------------------------------
! rotation matrices
!----------------------------------------------------------------------
\ARRAY ROTATION_MATRICES
! 2001 90.000 0.000 90.000 90.000 0.000 0.000
! 2002 90.000 180.000 90.000 90.000 180.000 0.000
\END
!----------------------------------------------------------------------
! material and mixture definitions
!----------------------------------------------------------------------
!
\ARRAY MATERIAL_LIST
!
!----------------------------------------------------------------------
! ADD NEW MATERIALS HERE. CALL TO GSMATE
! FORMAT
!CHNAME = 'MATERIAL NAME',
!IMATE = MATERIAL NUMBER
! A = ATOMIC WT
! Z = ATOMIC NUMBER
! DENS = DENSITY
! RADL = RADIATION LENGTH
! ABSL = ABSORPTION LENGTH (IGNORED BY GEANT)
! NWBUF = NUMBER OF EXTRA USER WORDS
! UWBUF = NWBUF EXTRA USER WORDS
!----------------------------------------------------------------------
!
! CHNAME IMATE A Z DENS RADL ABSL NWBUF UWBUF
\END
\ARRAY MIXTURE_LIST
!----------------------------------------------------------------------
! MIXTURE DEFINITIONS . CALL TO GSMIXT
! FORMAT IS AS FOLLOWS.
! MIXTURE NAME (CAN BE UP TO 32 CHARS), MIXTURE NUMBER (OBSERVE CONVENTION HERE)
! NUMBER OF MATERIALS IN MIXTURE (NLMAT),
! NLMAT MATERIAL NUMBERS,
! WMAT(NLMAT) WEIGHTS.
! IF NLMAT>0, WMAT CONTAINS THICKNESSES OF MATERIALS IN PROPORTION.
! density is calculated.
! IF NLMAT<0, WMAT CONTAINS PROPORTION BY ATOMS OF MATERIALS IN MIXTURE
! followed by density of mixture.
!----------------------------------------------------------------------
!'IRON_COPPER_TITANIUM' 201 3 10 11 19 0.5 0.25 0.25
\END
!
\ARRAY MEDIA_LIST
! Media name, media number, material number,ISVOL IFIELD FIELDM TMAXFD STEMAX DEEMAX EPSIL STMIN
! THE REMAINING PARAMETERS ARE THE ARGUMENTS TO THE CALL TO GSTMED
! ISVOL 1 !1= sensitive volume. 0=not sensitive
! IFIELD 3 !0=no field
 !1=tracking with Runge Kutta
! !2=inhomegeneous magnetic field. Field map from GUFLD
 !3=Uniform magnetic field along z axis of strenght FIELDM
! FIELDM 60.0000 !Maximum field value in Kilo Gauss
! TMAXFD 1.0000 !Maximum angle due to field in degrees.
! STEMAX 0.1000 !Maximum displacement in one step for mult. scat
! DEEMAX 0.1000 !Maximum fractional energy loss per step
! EPSIL 2.0000E-02 !Tracking precision (cm)
! STMIN 2.0000E-02 !Minimum step due to energy loss or mult. scat.(cm)
! NU !NUMBER OF USER WORDS
! W(1:NU) !USER WORDS
!----------------------------------------------------------------------
! MEDIA DEFINITION. CALL TO GSTMED
!----------------------------------------------------------------------
! MEDNAME IMED IMAT ISVL IFLD FIELDM TMAXFD STEMAX DEEMAX EPSIL STMIN NU W
!
\END
!----------------------------------------------------------------------
! Geometry definitions
!----------------------------------------------------------------------
!
! EVERY VOLUME XXXX HAS AN ARRAY XXXX THAT DESCRIBES THE GSVOLU PARAMETERS
! AND AN ARRAY ACT_XXXX THAT DESCRIBES THE ACTION ON THE VOLUME .GSPOS,GSPOSP
! GSDVN,GSDVN2,GSDVT,GSDVX,OR GSDVT2.
!
! FORMAT OF ACTIONS:-
!**** ALL CHARACTER VARIABLES IN THESE FIELDS ARE 4 CHARACTERS LONG****
!THE FIRST 4 CHARACTER VARIABLE OF THE ACTION FIELD IS 'VOLU','POS ','POSP',
!'DVN ','DVN2','DVT ','DVT2','DVX ' WHICH IS EQUIVALENT TO THE CALL TO GEANT
!SUBROUTINES GSVOLU, GSPOS ETC.
!THE ORDER OF THE VARIABLES IS THE SAME AS IN THE GEANT SUBROUTINE CALL EXCEPT
!FOR THE FOLLOWING EXCEPTIONS. IF ONE IS PERFORMING THE ACTION ON A VOLUME NAME
!(SAY 'CEM1'),
!THAT INFORMATION IS CONTAINED IN THE ARRAY NAME (ACT_CEM1),
!SO IT IS NOT REPEATED IN THE ACTION FIELD. ALSO, IN THE ACSE OF VARIABLE
!PARAMETERS AS IN THE VOLU AND POSP ACTIONS, THE NUMBER OF PARAMETERS (NPARS)
!IS DECLARED BEFORE THA PARAMETERS THEMSELVES. THIS IS OPPOSITE TO THAT IN THE
!SUBROUTINE CALL. THE READER SHOULD READ THE GEANT GEOMETRY MANUAL FOR A
!COMPLETE SPECIFICATION OF THE OTHER VARIABLES.
!----------------------------------------------------------------------
!----------------------------------------------------------------------
! TRACKING MEDIA NUMBERS
!----------------------------------------------------------------------
 air 1 ! AIR 
 lead 2 ! LEAD 
 scin 3 ! SCINTILLATOR 
 ledd 4 ! LEADD 
 sci1 5 ! SCINTI  
 iron 6 ! IRON 
 copp 7 ! COOPER 
 argn 8 ! ARGON 
 plex 9 ! PLEX 
 null 10 ! 
 null 11 ! 
 null 12 ! 
 hydr 13 ! LIQUID HYDROGEN 
 vac0 14 ! VACUUM0- NO MAGNETIC FIELD 
 vacm 15 ! VACUUMM- WITH MAGNETIC FIELD 
 null 16 ! 
 null 17 ! 
!----------------------------------------------------------------------
! VOLME DEFINTIONS
!----------------------------------------------------------------------
\ARRAY VOLUME_LIST
 'ARTE' !Mothervolume 
 'RBY4' !4 fold symmetric trapezium mother 
 'BMAG' !Bending magnets  
 'BMA1' !Bending magnets 
 'LCOI' !Long section solenoid 
 'LCO1' !Long section solenoid
 'MABS' !Long section absorber 
 'RFCV' !rf Cavity
 'ELCO' !Long section solenoid end volume
 'SCOI' !Short Section solenoid 
 'SCO1' !Short section solenoid
 'SCO2' !Short section solenoid 
 'ESCO' !Short section solenoid end volume
 'WABS' !Wedge absorber 1
\END
 P001 0.7500000E+03
 P002 0.3338950E+03
 P003 0.9150000E+02
 P004 0.0000000E+00
 P005 0.4673037E+03
 P006 -0.4673037E+03
 P007 0.9000000E+02
 P008 0.3000000E+02
 P009 0.4500000E+02
 P010 -0.1009122E+02
 P011 0.3095000E+02
 P012 -0.4151054E+01
! P012 -0.2443480E+01
 P013 0.2895000E+02
 P014 0.4190878E+02
 P015 21.88495 ! P017 Absorber half width in z has to fit in SCOI
 P016 0.8950000E+02
 P017 -21.88495 !linked to P011/sqrt(2)
 P018 43.59934 !SCOI/2.0 Delta Z
 P019 43.76991 !twice P015
 P020 12.7 !half angle of wedge absorber in degrees.
\ARRAY ACT_ARTE
 'VOLU' 'BOX ' 'vac0' 3 'P001' 'P001' 'P001'
\END
\ARRAY ACT_RBY4
 'VOLU' 'TRD1' 'vacm' 4  'P002' 0.600712556E+03 'P003' 0.133408778E+03
 'POS ' 1 'ARTE' 'P004' 'P004' 'P005' 0 'ONLY'
 'POS ' 2 'ARTE' 'P005' 'P004' 'P004' 27 'ONLY'
 'POS ' 3 'ARTE'  'P004' 'P004' 'P006' 25 'ONLY'
 'POS ' 4 'ARTE' 'P006' 'P004' 'P004' 28 'ONLY'
\END
\ARRAY ACT_BMAG
 'VOLU' 'TUBS' 'vacm' 5 'P004' 104.000 'P008' 'P009'  'P007'
 'POS ' 1 'RBY4' 'P002' 'P004' 'P010' 24 'ONLY'
\END
\ARRAY ACT_BMA1
 'VOLU' 'TUBS' 'vacm' 5 'P004' 104.000 'P008' 'P007' 0.1350000E+03
 'POS ' 1 'RBY4' -0.3338950E+03 'P004' 'P010' 24 'ONLY'
\END
\ARRAY ACT_SCOI
 'VOLU' 'TUBE' 'vacm' 3 'P004' 'P011' 'P018'
 'POS ' 1 'RBY4' -0.40149394E+03 'P004' 'P012' 22 'ONLY'
 'POS ' 2 'RBY4' 0.40149394E+03 'P004' 'P012' 23 'ONLY'
\END
\ARRAY ACT_SCO1
 'VOLU' 'TUBE' 'iron' 3 'P013' 'P011' 0.1500000E+02
 'POS ' 1 'SCOI' 'P004' 'P004' 0.2140066E+02 0 'ONLY'
\END
\ARRAY ACT_SCO2
 'VOLU' 'TUBE' 'iron' 3 'P013' 'P011' 0.2500500E+02
 'POS ' 1 'SCOI' 'P004' 'P004' -0.1859934E+02 0 'ONLY'
\END
\ARRAY ACT_ESCO
 'VOLU' 'TUBE' 'vacm' 3 'P004' 'P011' 0.05
 'POS ' 1 'SCOI' 'P004' 'P004' -43.54934 0 'ONLY'
 'POS ' 2 'SCOI' 'P004' 'P004' 43.54934 0 'ONLY'
\END
\ARRAY ACT_MABS
 'VOLU' 'TUBE' 'hydr' 3 'P004' 'P011' 0.6650000E+02
 'POS ' 1 'LCOI' 'P004' 'P004' 'P004' 0 'ONLY'
\END
\ARRAY ACT_WABS
 'VOLU' 'TUBS' 'hydr' 5 'P004' 30.95000 'P015' 'P004' 'P020'
!'VOLU' 'TUBS' 'hydr' 5 'P004' 'P019' 'P015' 'P004' 'P020'
 'POS ' 1 'SCOI' 'P004' 'P004' 43.49900 31 'ONLY'
!'POS ' 1 'SCOI' 'P004' 'P017' 43.54934 31 'ONLY'
\END
\ARRAY ACT_LCOI
 'VOLU' 'TUBE' 'vacm' 3 'P004' 'P003' 'P002'
 'POS ' 1 'RBY4' 'P004' 'P004' 'P014' 27 'ONLY'
\END
\ARRAY ACT_LCO1
 'VOLU' 'TUBE' 'iron' 3 'P016' 'P003' 0.2200000E+02
 'POS ' 1 'LCOI' 'P004' 'P004' -0.3018950E+03 0 'ONLY'
 'POS ' 2 'LCOI' 'P004' 'P004' -0.2378950E+03 0 'ONLY'
 'POS ' 3 'LCOI' 'P004' 'P004' -0.1738950E+03 0 'ONLY'
 'POS ' 4 'LCOI' 'P004' 'P004' -0.1098950E+03 0 'ONLY'
 'POS ' 5 'LCOI' 'P004' 'P004' -0.4589500E+02 0 'ONLY'
 'POS ' 6 'LCOI' 'P004' 'P004' 0.4589500E+02 0 'ONLY'
 'POS ' 7 'LCOI' 'P004' 'P004' 0.1098950E+03 0 'ONLY'
 'POS ' 8 'LCOI' 'P004' 'P004' 0.1738950E+03 0 'ONLY'
 'POS ' 9 'LCOI' 'P004' 'P004' 0.2378950E+03  0 'ONLY'
 'POS ' 10 'LCOI' 'P004' 'P004' 0.3018950E+03 0 'ONLY'
\END
\ARRAY ACT_RFCV
 'VOLU' 'TUBE' 'vacm' 3 'P004' 60.0 0.15900E+02
 'POS ' 1 'LCOI' 'P004' 'P004' -317.895  0 'ONLY'
 'POS ' 2 'LCOI' 'P004' 'P004' -285.395 0 'ONLY'
 'POS ' 3 'LCOI' 'P004' 'P004' -252.895 0 'ONLY'
 'POS ' 4 'LCOI' 'P004' 'P004' -220.395 0 'ONLY'
 'POS ' 5 'LCOI' 'P004' 'P004' -187.895 0 'ONLY'
 'POS ' 6 'LCOI' 'P004' 'P004' -155.395 0 'ONLY'
 'POS ' 7 'LCOI' 'P004' 'P004' -122.895 0 'ONLY'
 'POS ' 8 'LCOI' 'P004' 'P004' -90.395 0 'ONLY'
!
 'POS ' 9 'LCOI' 'P004' 'P004' 90.395 0 'ONLY'
 'POS ' 10 'LCOI' 'P004' 'P004' 122.895 0 'ONLY'
 'POS ' 11 'LCOI' 'P004' 'P004' 155.395 0 'ONLY'
 'POS ' 12 'LCOI' 'P004' 'P004' 187.895 0 'ONLY'
 'POS ' 13 'LCOI' 'P004' 'P004' 220.395 0 'ONLY'
 'POS ' 14 'LCOI' 'P004' 'P004' 252.895 0 'ONLY'
 'POS ' 15 'LCOI' 'P004' 'P004' 285.395 0 'ONLY'
 'POS '  16 'LCOI' 'P004' 'P004' 317.895 0 'ONLY'
!
 \END
\ARRAY ACT_ELCO
 'VOLU' 'TUBE' 'vacm' 3 'P004' 'P003' 0.05
 'POS ' 1 'LCOI' 'P004' 'P004' -333.845 0 'ONLY'
 'POS ' 2 'LCOI' 'P004' 'P004' 333.845 0 'ONLY'
\END
!----------------------------------------------------------------------
! Detector definition
!----------------------------------------------------------------------
\ARRAY DETECTOR_LIST
 'DET_MABS' ! main absorbers
 'DET_WABS' ! wedge absorbers
 'DET_RFCV' ! RF Cavities
\END
\ARRAY DET_MABS
!
!Arguments for multi-calls to GSDET
! CHSET NWHI NWDI
 'MABS' 1024 1024
! Number of detectors belonging to this CHSET
 1
! IDTYPE CHDET NV CHNSMV(1:nv) NBITSV Hits definition pointer
 2001 'MABS' 4 'ARTE' 'RBY4' 'LCOI' 'MABS' 16 'CHIT'
\END 
\ARRAY DET_WABS
!
!Arguments for multi-calls to GSDET
!  CHSET NWHI NWDI
 'WABS' 1024 1024
! Number of detectors belonging to this CHSET
 1
! IDTYPE CHDET NV CHNSMV(1:nv) NBITSV Hits definition pointer
 2002 'WABS' 4 'ARTE' 'RBY4' 'SCOI' 'WABS' 16 'CHIT'
\END 
\ARRAY DET_RFCV
!
!Arguments for multi-calls to GSDET
! CHSET NWHI NWDI
 'RFCV' 1024 1024
! Number of detectors belonging to this CHSET
 1
! IDTYPE CHDET NV CHNSMV(1:nv) NBITSV Hits definition pointer
 2002 'RFCV' 4 'ARTE' 'RBY4' 'LCOI' 'RFCV' 16 'CHIT'
\END 
\ARRAY CHIT
!
!HITS DEFINITION. Arguments for GSDETH
! NH
 14
!CHNAMH(1:NH) followed by NH values of NBITSH, ORIG and FACT
'ITR' 'SLN' 'DE' 'TSLN' 'XE' 'YE' 'ZE' 'TFE' 'ENE' 'XX' 'YX' 'ZX' 'TFX' 'ENX' 
 14*32
 0.0 0.0 0.0 0.0 1.E4 1.E4 1.E4 0.0 0.0 1.E4 1.E4 1.E4 0.0 0.0
 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5 1.E5
\END
\STOP ! MUCOOL_RCP