physics QGSP_BIC #param histoFile=MC3_7M_P32.root # killshld used to kill tracks in shielding, significantly faster running # set to nonzero to kill all particles that enter param killshld=1 g4ui when=4 "/vis/viewer/set/background 1 1 1" ############################################################################## ######## Use this beam if you want to hit the target ###################### ########### Run 100000 events ############################################################################## #beam gaussian particle=pi+ nEvents=100 sigmaX=2.0 sigmaY=2.0 beamZ=-500 \ # meanMomentum=20000.0 rotation=y0 beamX=0.0 ############################################################################## virtualdetector DetBeam width=50.0 height=50.0 color=0,1,0 length=25.0 virtualdetector DetBeam2 width=50.0 height=50.0 color=1,1,0 length=25.0 virtualdetector DetBeam3 width=200.0 height=200.0 color=1,1,0 length=25.0 # place DetBeam z=-200 # place DetBeam z=0 # place DetBeam z=150.0 # place DetBeam z=200 ################################################################################################## ############### Target ################################################################ ################################################################################################## param scale=16.0/16.0 #31.75 mm is actually the width....the length is 320.802 mm #Flat surface is 211.328 mm extrusion Target length=31.75 vertexes='0,0;109.728,31.75;321.056,31.75;211.328,0;' material=Cu color=1,0,1 param targetPosition=350 place Target x=-(31.75/2) y=0 z=$targetPosition rotation=x0,y90,z-90 #corner z=10000 rotation=Y10 cornerarc z=$targetPosition-(200) angle=-16.0*$scale centerRadius=0.001 # place DetBeam z=215.0 # place DetBeam2 z=700.0 # place DetBeam2 z=800.0 ################################################################################################## ####################################################################################### ######### Use this beam command to display a beam of specific energy ########### ########## After the target thus enabling you to see where it travels ############# ####################################################################################### beam gaussian particle=pi+ nEvents=99 sigmaX=2.0 sigmaY=2.0 beamZ=250 \ meanMomentum=2050.0 sigmaP=100 rotation=y0 beamX=0.0 ####################################################################################### ####################################################################################### ###############--------- collimator -------------------################################ ####################################################################################### # The Fe length = 42.76 inches ~= 1086 mm box topbot height=180 length=1086 width=800 material=Fe color=0,1,1 kill=$killshld box side height=50 length=1086 width=400 material=Fe color=0,1,1 kill=$killshld * us end of coll is 2.83", ds is 5.65" length is 42.71/cos(16)=44.43 * so opening angle of collimator is 1.8176 deg * half width of collimator at the center is 2.12 inches = 53.85 mm #place topbot rename=bot y=-115 z=882 #place topbot rename=top y=+115 z=882 place side rename=right x=-254.0 z=882.0 rotation=y-1.81 place side rename=left x=+254.0 z=882.0 rotation=y1.81 ######################################################################################## param zPosition=((4092.2+4776.9)/2) param halfSteel=(1066.8/2) ############################################################################## ######### FTBF Drift Chambers 1 and 2 ################################## ##### FTBF MWPC has 128 wires 1 mm apart for X and Y planes ############## ############################################################################## virtualdetector detector1 width=128.0 height=128.0 color=0,0,1 length=25.0 virtualdetector detector2 width=128.0 height=128.0 color=0,0,1 length=25.0 place detector1 z=$zPosition-$halfSteel-100-1000 material=POLYSTYRENE format=ascii place detector2 z=$zPosition-$halfSteel-100 material=POLYSTYRENE format=ascii # place DetBeam3 z=$zPosition-$halfSteel-100-1000 #I put these in to be more visible # place DetBeam3 z=$zPosition-$halfSteel-100 ############################################################################## ######################################################################################## ########## 2 NDB magnets which ben 5 degrees each ################################ ##########------------ total 10 deg bend --------------------------------############# ######################################################################################## genericbend NDB fieldWidth=317.5 fieldHeight=142.24 fieldLength=591.0 kill=$killshld \ ironColor=1,0,0 ironWidth=533.4 ironHeight=474.675 ironLength=466.725 param Bscale=1.0 param B1=-0.35*$Bscale param B1z=4092.2 param B2=-0.35*$Bscale param B2z=4776.9 #place NDB rename=NDB22 By=$B1 z=$B1z rotation=Y2.5 # cornerarc z=$B1z-233.0 angle=5 centerRadius=3600.0/$Bscale #place NDB rename=NDB21 By=$B2 z=$B2z rotation=Y5 # cornerarc z=$B2z-233.0 angle=5 centerRadius=3600.0/$Bscale ######################################################################################### ######################################################################################### ######################################################################################### ################## Nova Test Beam Magnets Using ################################# ################# 1 42" Wideband Sweeper ################################## ### 42" Steel 48" max coil to coil Bmax 11.3 KG 1700 Amps 56 Turns ############### #### <4" Gap....3.5" is what is listed in Excel ##################################### ## Excel file shows an M1 and M2 magnet 1180 Amps and .934 T #################### ## Excel file shows M3 with same length and number of coils at 2274 Amps and 1.8 T #### #### Good field 4 inch wide or 101.6 mm I think this could be better???? ##### #### Gap is 3.5 inches or 88.9 mm ########### #### Steel length is 42 inches 1066.8 mm ########### #### Field length is iron length plus 1/2 gap on each end????? ######## ##### Iron height is ???? 30 inches (762 mm) is a guess need to measure ######## ##### Iron Width is a guess 36 inches (914.4) is a guess need to measure ######## ## Lariat magnets are at z= 4092.2 and z=4776.9 thus ours are initial (4092.2+4776.9)/2 # ######################################################################################### #param zPosition=((4092.2+4776.9)/2) #param halfSteel=(1066.8/2) param scale=10.0/16.0 genericbend sweeper fieldWidth=101.6 fieldHeight=88.9 fieldLength=1066.8+88.9/2+88.9/2 kill=$killshld \ ironColor=0,0,1 ironWidth=533.4 ironHeight=474.675 ironLength=1066.8 place sweeper By=-1.8*$scale z=$zPosition x=65 rotation=Y8 cornerarc z=$zPosition-$halfSteel angle=16*$scale centerRadius=3600.0/$Bscale ######################################################################################### ######################################################################################### ############################################################################## ######### FTBF Drift Chambers 3 and 4 ################################## ###### FTBF MWPC has 128 wires 1 mm apart for X and Y planes ############# ############################################################################## virtualdetector detector3 width=128.0 height=128.0 color=0,0,1 length=25.0 virtualdetector detector4 width=128.0 height=128.0 color=0,0,1 length=25.0 place detector3 z=$zPosition+$halfSteel+100 material=POLYSTYRENE format=ascii place detector4 z=$zPosition+$halfSteel+100+1000 material=POLYSTYRENE format=ascii # place DetBeam3 z=$zPosition+$halfSteel+100 ##I put these into to be more visible # place DetBeam3 z=$zPosition+$halfSteel+100+1000 ############################################################################## ############################################################################## ############################################################################################ ####### Nova Test Beams Detector ################################################### ############################################################################################ #put the detector 4 meters downstream of the final ftbf detector box nova height=2600 length=2600 width=2000 material=POLYSTYRENE color=.9,.6,1 kill=$killshld place nova y=0 z=$zPosition+$halfSteel+100+1000+4000+2000/2 place nova y=0 z=$zPosition+$halfSteel+100+1000+4000+2000/2+4000 ############################################################### ## wire chamber is 128 mm draw a triangle from ########## ## center of first chamber to edge of second ############# ## Trace this solid angle 4 meters to the detector ######## ## 128*4=512 mm so the good beam will only by in a 512 mm ### # portion of the detector ######## ############################################################### virtualdetector novaDetectorRegion width=512.0 height=512 color=0,0,1 length=25.0 place novaDetectorRegion z=$zPosition+$halfSteel+100+1000+4000+25/2-500 material=POLYSTYRENE format=ascii ############################################################################################ ############################################################################################ survey ########################################################################## ########################### Food For Thought ####################### ########################################################################## ### 53.1 MHz Buckets ### 4 second spill ### 53*10**6 cycles/second * 4 second = 212*10**6 cycles per spill ### 1,000,000 counts per spill on MC7SC1 is about the maximum (uneless I write why more is OK) ### 1,000,000 particles/spill divided by 212*10**6 buckets per spill is .005 particles per bucket ########################################################################## ##########################################################################