Transcript G4Beamline for KEK test K. Yonehara Illinois Institute of Tech 3/31/2005 KEK, Japan p2 beam line 3/31/2005 KEK, Japan.

G4Beamline for KEK test
K. Yonehara
Illinois Institute of Tech
3/31/2005
KEK, Japan
1
p2 beam line
3/31/2005
KEK, Japan
2
3/31/2005
KEK, Japan
3
3/31/2005
KEK, Japan
4
Initial beam
*
ppro.in 3/10/05 KY
*
* Pion production input file
# QGSP is the "default" physics use-case for High Energy Physics
physics QGSP doStochastics=1
trackcuts keep=proton,pi+,mu+,e+ kineticEnergyCut=10 maxTime=10000
particlecolor proton=0,0.8,0 pi+=1,0,0 mu+=0,0,1 plus=1,0,1 minus=1,1,0 \
neutral=0,1,1
# the beam is nominally headed in the +Z direction
beam gaussian particle=proton nEvents=1000000 beamZ=0.0 \
sigmaX=0.10 sigmaY=0.10 sigmaXp=0.00 sigmaYp=0.00 \
meanMomentum=12000.0 sigmaP=10.0 meanT=0.0 sigmaT=0.0
# inner target
tubs InnerT innerRadius=0.0 outerRadius=0.5 length=20.0 material=W \
color=0.0,0.0,0.5
# slit
tubs Col innerRadius=200.0 outerRadius=1600.0 length=1.0 material=Pb \
color=0.5,0.0,0.0 kill=1
# define the detector
asciifile Det0 file=piontrack.dat radius=1000.0 color=0,1,0 length=0.001 material=Vacuum
# place four detectors, putting their number into their namesp
#place Det0 z=-1000.0 rename=Det#
place BeamVis z=0
place InnerT z=20.0
#place Det0 z=1020.0 rename=Det#
#place Det0 z=2020.0 rename=Det#
place Col z=3000.0 x=528.7
place Det0 z=3020.0 rename=Det#
3/31/2005
KEK, Japan
5
Produced Pion & Muon Spectra
in the dynamic aperture of p2 line
Inner target
3/31/2005
KEK, Japan
6
Beam transport (1)
#
kekb.in KY 03/10/05
# fringe=0 for no fringe fields, fringe='' to enable the default computation
param fringe=''
param B1=-1.770 B2=-1.770 B3=-2.001 B4=-2.001
param pionP=4300.0 muP=4300.0 protonP=12000.0
param Q1=-0.8356 Q2=0.7153 Q3=0.6363 Q4=-0.5666
physics LHEP_BIC doStochastics=1
trackcuts keep=pi+,mu+,e+ kineticEnergyCut=10 maxTime=50000
particlecolor proton=0,0.8,0 pi+=1,0,0 mu+=0,0,1 plus=1,0,1 minus=1,1,0 \
neutral=0,1,1
beam ascii filename=rotpion.dat beamZ=2800.0 format=BLTrackFile renumber=1 nEvents=14180
###
### VirtualDetector used throughout the beamline
###
virtualdetector BeamDet radius=200 noSingles=1 color=1,1,1
###
### define the beamline magnets
###
genericquad QuadTypeIV fieldLength=1074.0 apertureRadius=100.0 ironRadius=381 \
ironLength=1100 ironColor=0,.6,0 kill=1 fringe=$fringe maxStep=20.0
# Type I bend shimmed to 200 mm gap
genericbend BendTypeI fieldWidth=1060 fieldHeight=200 fieldLength=1069 \
ironColor=1,0,0 ironWidth=2028 ironHeight=1320 ironLength=1100 \
fringe=$fringe maxStep=20.0
tubs Col innerRadius=300.0 outerRadius=2000.0 length=20.0 material=W \
color=0.5,0.0,0.0 kill=1
asciifile Det0 file=pionatbl.dat radius=1000.0 color=0,1,0 \
length=0.001 material=Vacuum
3/31/2005
KEK, Japan
7
Beam transport (2)
###
### Lay out the beamline
###
place BeamDet z=1
place InnerT z=100
#cornerarc z=50 angle=10.0 centerRadius=200.0
place Det0 z=2800.1 rename=Det#
place QuadTypeIV rename=Q1 gradient=$Q1 z=3467
place QuadTypeIV rename=Q2 gradient=$Q3 z=5367
place BendTypeI rename=B1 By=$B1 z=7174 x=0 rotation=Y3.75
cornerarc z=6663 angle=7.5 centerRadius=1038
place BendTypeI rename=B2 By=$B2 z=8973 x=0 rotation=Y3.75
cornerarc z=8463 angle=7.5 centerRadius=2005.3
place BendTypeI rename=B3 By=$B3 z=17003 x=0 rotation=Y4.24
cornerarc z=16493 angle=8.48 centerRadius=2005.3
place QuadTypeIV rename=Q3 gradient=$Q3 x=-100.0 z=18967
place QuadTypeIV rename=Q4 gradient=$Q4 x=-100.0 z=20867
place BendTypeI rename=B4 By=$B4 z=23305 x=0 rotation=Y4.24
cornerarc z=22795 angle=8.48 centerRadius=2005.3
virtualdetector TOF height=400 width=400 length=1 color=1,1,1
tubs TrackerPipe innerRadius=150 outerRadius=250 length=1260 color=1,1,1 kill=1
virtualdetector TrackerDet radius=150 length=1.0 material=Vacuum color=1,1,1
place TrackerPipe z=38531.6
#place TrackerDet z=38531.6 x=-100.0 rename=Tracker1
place Col z=29000.0
place Det0 z=30000.0 rename=Det#
3/31/2005
KEK, Japan
8
3/31/2005
KEK, Japan
9
Pion spectrum at tracker position
# of proton = 106
3/31/2005
KEK, Japan
10
Experiment & Simulation
3/31/2005
KEK, Japan
11
Conclusion
► Beam
simulation study for KEK beam line
has been done in g4beamline.
► The pion and muon spectra in the dynamic
aperture of beam line were extracted.
► Track pions and muons in the beam line.
► Extract spectra of pion and muon at the end
of beam line.
3/31/2005
KEK, Japan
12
Next to do
► Optimize
the field parameters to extract the
low momentum muon.
 Is there any limitation to adjust parameters?
► Compare
the simulation results with the
experimental one.
 Slit size?
 Field parameters?
3/31/2005
KEK, Japan
13