Layout with 9 frames 24 January 2018 Saba
Layout with 9 frames 24 January 2018 Saba Parsa
side Layout driven by physics considerations d: scintillator detector module a: ARMCO steel module d 2 d 1 [a 1 -3] d 3 d 4 [a 4] Case 1: Gap 20 cm Case 2: Gap 30 cm d 5 [a 5] [a 6] d 8 d 7 d 6 [a 7 -8] Gap 20 cm Gap 10 cm [a 9 -10] 10 cm d 9 d 11 d 10 [a 11 -12] [a 13 -14] 10 cm d 12 [a 15 -17] d 13 10 cm d 17 d 16 [a 26 -29] [a 22 -25] [a 18 -21] 10 cm d 15 d 14 10 cm d 18 [a 30 -33] 10 cm 2
GDML files: Case 1: MIND_jp_23 Jan_2 gaps. gdml Case 2: MIND_jp_23 Jan_1 gap. gdml • "Fe. Width" value="3500”, "Fe. Height" value="2100”, " Fe. Depth" value="30” 5 5 • "Al. Width" value="5" Simulated Armco • "Bar. Depth" value="7. 5"/>, "Sci. Depth" value="4*Bar. Depth"/> 30 5 5 50 7. 5 • "Gap 0" value="5" between Fe and Al • "Gap. Fe. Sci" value="15" between Fe and Sci • "Gap. Fe" value="20” between Fe and Fe Simulated Scintillator (In reality is it 40 mm) • "Gap 1" value="200"/> Large Gap • "Gap 2" value="100"/> 2 Frames • "Gap 3" value="50"/> Last. Frame 7. 5 30 • "Sci. Width" value="2900”, "Sci. Height" value="2000”, • "Sci. X 1 Width" value="2720”, " Sci. X 2 Width" value="2841. 0”, "Sci. Y 1 Width" value="2900", "Sci. Y 2 Width" value="2900”, • "Sci. X 1 Height" value="1972. 4, "Sci. Y 1 Height" value="1972. 4”, "Sci. Y 2 Height" value="1972. 4"/> • "Bar. Width. X 1" value="210. 6"/> , "Bar. Width. X 2" value="350. 7"/> , "Bar. Width. Y 1" value="31. 3"/> , "Bar. Width. Y 2" value="41. 3"/> , 3
Simulations with Saroman in vitual box!! Kinetic 0. 5 energy Ge. V 0. 6 0. 8 1 1. 5 2 5 Layout 2 gap Rec fails Done Layout 1 gap Rec fails Done The reconstruction fails for low momentum. The Reconstruction is written in a way that can deal only with a certain geometry. 4
Layout Case 1 1. 5 Ge. V 2 Ge. V 5
Layout Case 2 1 Ge. V 2 Ge. V 1. 5 Ge. V 6
Particle Generator (pythonlib/print_config. py) • Particle generator is by default GAUSSIAN • (I don’t know vertex location and size of the gaussian) • When changed to FIXED (vertex: 0, 0, -7000 , costh_min=1 costh_max=1) • Events map look exactly similar to gaussian FIXED target Z = -7000 Z = -5000 Z = -4000 Z = -3000 Z= -2900 Z = -2500 Costh min/max = 1 Done Rec fails Fixed target Costhmin/m ax =1 Costhmin/m ax =0. 99 Costhmin/ max =0. 98 Costhmin = 1/ Costhmax =0. 95 Costhmin = 1/ Costhmax =0. 9 Costhmin = 0. 95 Costhmax = 0. 9 Z = -3000 Done Rec fail Root empty 7
Particle Generator (pythonlib/print_config. py) FIXED target Z = -7000 Z = -5000 Z = -4000 Z = -3000 Z= -2900 Z = -2500 Costh min/max = 1 Done Rec fails Z = -7000 , Costh = 1 Done Rec fails Z = -3000 , Costh = 1 8
Particle Generator (pythonlib/print_config. py) Fixed target Costhmin/m ax =1 Costhmin/m ax =0. 99 Costhmin/ max =0. 98 Costhmin = 1/ Costhmax =0. 95 Costhmin = 1/ Costhmax =0. 9 Costhmin = 0. 95 Costhmax = 0. 9 Z = -3000 Done Rec fail Root empty Z = -3000 , Costh = 1 Z = -3000 , Costh min= 1 / Costh max =0. 9 9
Particle Generator (pythonlib/print_config. py) • It seems the settings for angle of generated particles does not function. • If we do not provide costh_min =1 the Reconstruction fails. • For any value of costh_max (0. 9, 0. 95, 0. 98) the distribution is identical -> this parameter is not used • It seems that with the Gaussain vertex, we get the same distribution as with FIXED vertex Is there any Questions? 10
Back up slide: Material in gdml • /> isotope< • > isotope N="54" Z="26" name="Fe 540 x 1931 f 340<" • > element name="C 0 x 19316 fa 0<" • > • <materials> • > • /> isotope< • > • /> • atom unit="g/mole" value="53. 9396</" isotope N="14" Z="7" name="N 140 x 19317540<" atom unit="g/mole" value="14. 0031</" • > atom unit="g/mole" value="37. 9627</" • /> isotope< isotope N="40" Z="18" name="Ar 400 x 19317 d 00<" • > fraction n="0. 9893" ref="C 120 x 19317200</" • /> isotope< • > fraction n="0. 0107" ref="C 130 x 193172 b 0</" • > isotope N="15" Z="7" name="N 150 x 193175 e 0<" • > isotope N="56" Z="26" name="Fe 560 x 1931 f 3 b 0<" • /> element< • > • /> isotope< • > isotope N="1" Z="1" name="H 10 x 19323 de 0<" • /> isotope< • > element name="Ar 0 x 19317990<" • > atom unit="g/mole" value="1. 00782503081372</" • > element name="N 0 x 19317310<" • > fraction n="0. 003365" ref="Ar 360 x 19317800</" isotope N="57" Z="26" name="Fe 570 x 1931 f 420<" • /> isotope< • > fraction n="0. 99632" ref="N 140 x 19317540</" • > fraction n="0. 000632" ref="Ar 380 x 19317 c 10</" • > isotope N="2" Z="1" name="H 20 x 193245 e 0<" • > fraction n="0. 00368" ref="N 150 x 193175 e 0</" • > fraction n="0. 996003" ref="Ar 400 x 19317 d 00</" • > atom unit="g/mole" value="2. 01410199966617</" • /> element< • /> isotope N="58" Z="26" name="Fe 580 x 193228 d 0<" • /> isotope< • > isotope N="16" Z="8" name="O 160 x 19317500<" • • > element name="H 0 x 19324360<" • > isotope< • > fraction n="0. 999885" ref="H 10 x 19323 de 0</" • /> element name="Fe 0 x 19322640<" • > fraction n="0. 000115" ref="H 20 x 193245 e 0</" • > fraction n="0. 05845" ref="Fe 540 x 1931 f 340</" • /> • > > fraction n="0. 91754" ref="Fe 560 x 1931 f 3 b 0</" • • /> isotope< • > fraction n="0. 02119" ref="Fe 570 x 1931 f 420</" <material name="G 4_POLYSTYRENE 0 x 19324240" state="solid"> isotope N="18" Z="8" name="O 180 x 19317950<" fraction n="0. 00282" ref="Fe 580 x 193228 d 0</" MEE unit="e. V" value="68. 7</" > • D unit="g/cm 3" value="1. 06</" > • /> • • > fraction n="0. 922577329294378" ref="C 0 x 19316 fa 0</" • /> isotope< • > element name="O 0 x 19317620<" • atom unit="g/mole" value="55. 9349</" atom unit="g/mole" value="56. 9354</" atom unit="g/mole" value="57. 9333</" element< <material name="G 4_Fe 0 x 1931 f 240" state="solid"> element< atom unit="g/mole" value="15. 0001</" MEE unit="e. V" value="85. 7</" isotope< • > D unit="g/cm 3" value="0. 00120479</" isotope N="17" Z="8" name="O 170 x 193178 c 0<" • > fraction n="0. 000124000124" ref="C 0 x 19316 fa 0</" • > fraction n="0. 755267755" ref="N 0 x 19317310</" • > fraction n="0. 231781232" ref="O 0 x 19317620</" • > fraction n="0. 0128270128" ref="Ar 0 x 19317990</" • /> atom unit="g/mole" value="16. 9991</" atom unit="g/mole" value="17. 9992</" • > fraction n="0. 99757" ref="O 160 x 19317500</" D unit="g/cm 3" value="7. 874</" • /> • > fraction n="0. 00038" ref="O 170 x 193178 c 0</" fraction n="1" ref="Fe 0 x 19322640</" • > material Z="13" name="G 4_Al 0 x 19323860" state="solid<" • > fraction n="0. 00205" ref="O 180 x 19317950</" • /> element< • > isotope N="36" Z="18" name="Ar 360 x 19317800<" • > • /> isotope< • > isotope N="38" Z="18" name="Ar 380 x 19317 c 10<" MEE unit="e. V" value="286</" • > • /> material< • > isotope N="12" Z="6" name="C 120 x 19317200<" • > • /> isotope< • > isotope N="13" Z="6" name="C 130 x 193172 b 0<" • > atom unit="g/mole" value="12</" atom unit="g/mole" value="13. 0034</" material< • > MEE unit="e. V" value="166</" • > D unit="g/cm 3" value="2. 699</" • > atom unit="g/mole" value="26. 9815</" • /> material< <material name="G 4_AIR 0 x 19316 e 30" state="gas"> > > fraction n="0. 0774226707056217" ref="H 0 x 19324360</" > element< • atom unit="g/mole" value="15. 9949</" • • atom unit="g/mole" value="39. 9624</" atom unit="g/mole" value="35. 9675</" • material< <material name="G 4_Galactic 0 x 19324780" state="gas"> • > T unit="K" value="2. 73</" • > P unit="pascal" value="3 e-18</" • > MEE unit="e. V" value="21. 8</" • > D unit="g/cm 3" value="1 e-25</" • > fraction n="1" ref="H 0 x 19324360</" • /> • material< </materials> 11
- Slides: 11