260 likes | 392 Views
First simulations for ACTAR . H. Alvarez, J. Benlliure, M. Caamaño, D. Cortina, I. Durán Universidad de Santiago de Compostela. EURONS ACTAR JRA 1 Ganil, November 2005. Before to start. Kinematics of key reactions elastic scattering transfer reactions
E N D
First simulations for ACTAR H. Alvarez, J. Benlliure, M. Caamaño, D. Cortina, I. Durán Universidad de Santiago de Compostela EURONS ACTAR JRA 1 Ganil, November 2005
Before to start • Kinematics of key reactions elastic scattering transfer reactions • Comparison GEANT4 versus SRIM Range Longitudinal straggling Documentation in Santiago webpage .... GANIL 2005
Key reactions 6he p GANIL 2005
Key reactions 132Sn p GANIL 2005
Key reactions 8n 8n 25Al 22Mg GANIL 2005
Interaction of HI with gas 8He(p,p) @ 50 MeV Range Calculation Isobutane r at STP Range Range GEANT GEANT SRIM SRIM proton 8He Elab Elab Bordeaux 2005
Interaction of HI with gas Range Calculation SRIM GEANT GANIL 2005
Interaction of HI with gas Range Calculation SRIM GEANT GANIL 2005
Interaction of HI with gas Range Calculation SRIM GEANT GANIL 2005
Interaction of HI with gas Range Calculation SRIM GEANT GANIL 2005
Interaction of HI with gas Range Calculation SRIM GEANT GANIL 2005
Interaction of HI with gas Range Calculation SRIM GEANT GANIL 2005
Interaction of HI with gas Straggling SRIM GEANT GANIL 2005
Before to start • Kinematics of key reactions elastic scattering transfer reactions • Comparison GEANT4 versus SRIM Range Longitudinal straggling Other gases? Higher preassure?? GANIL 2005
Basic Ingredients • Simulation package GEANT 4 geometry description cubic volume (1m3) filled with gas any gas, pressure and geometry is in principle possible “pseudo-event generator” kinematics of different reactions output of CINE.FOR description of light charged particles and heavy ions in gas includes the effect of a magnetic field perpendicular to the symmetry plane rough pad pixelation histogram binning--> effect on final resolution output analysis possibilities--> ascii files histogrames Root Trees ( data base) Cosmetics Documentation Internal More general Bordeaux 2005
What is next ? • Simulation package GEANT 4 geometry description cubic volume (1m3) filled with gas any gas, pressure and geometry is in principle possible “pseudo-event generator” kinematics of different reactions output of CINE.FOR description of light charged particles and heavy ions in gas includes the effect of a magnetic field perpendicular to the simetry plane rough pad pixelation histogram binning--> effect on final resolution output analysis possibilities--> ascii files histogrames Root Trees ( data base) Implementation of other geometries Together with Liverpool: Incorporate realistic event generators Work on track reconstruction GANIL 2005
Inputs 8He(p,t) @ 50 MeV f=0o isobutane1 atm 3 atm By= 2T 8He(p,t) @ 100 MeV f=0o 1 atm 3 atm Bordeaux 2005
Event Generator Kinematics **** CINEMATIQUE **** 8He(p,p)8He @ 50 MeV THLAB3 ELAB3 THLAB4 ELAB4 0.000 50.000 0.00 0.000 0.000 30.207 0.00 19.793 1.000 49.877 85.46 0.123 1.000 30.281 3.49 19.719 2.000 49.504 80.84 0.496 2.000 30.509 7.06 19.491 3.000 48.860 76.04 1.140 3.000 30.911 10.81 19.089 4.000 47.908 70.93 2.092 4.000 31.524 14.88 18.476 5.000 46.565 65.27 3.435 5.000 32.432 19.50 17.568 6.000 44.631 58.48 5.369 6.000 33.837 25.24 16.163 7.000 41.106 47.76 8.894 :::::::::::::::::::::::::::::::::::::::::::::::: • Kinematics fixed by an external program • Fixed reaction point • No energy loss of primary beam in the gas target • selection of the azimuthal angle ( or random) Bordeaux 2005
Magnetic Field 8He(p,p) @ 50 MeV 8He Elab=47.908 MeV qlab=4.0oflab= -125o p E=2.092 MeV qlab=70.93o flab= 55o Isobutane STP By=2T By=4T Range 8He ~ 319 mm Range p ~ 26 mm Bordeaux 2005
Magnetic Field 8He(p,p) @ 50 MeV 8He Elab= 33.8 MeV qlab=6oflab= -29.3o p E=16.2 MeV qlab=25.2o flab=150.6o Isobutane STP By=4T By=2T Range 8He ~171.0 mm Range p ~ 964.6 mm Bordeaux 2005
Magnetic Field Isobutane STP By=2T 8He(p,t) @ 100 MeV 8He proton Isobutane 3 atm By=2T 8He proton Bordeaux 2005
Track reconstruction 8He(p,t) @ 100 MeV Isobutane r at STP Y By= 2T Z X 6He Elab= 72.897 MeV qlab=16oflab= -29o t E=33.53 MeV qlab=25oflab=150o Bordeaux 2005
Track reconstruction 8He(p,t) @ 100 MeV Isobutane r at STP Y By= 4T X Z 6He Elab= 72.897 MeV qlab=16oflab= -29o t E=33.53 MeV qlab=25oflab=150o Bordeaux 2005
Track reconstruction 8He(p,t) @ 100 MeV By= 3T Isobutane r at STP Pad resolution 4mm Pad resolution 8mm 6He Elab= 72.897 MeV qlab=16oflab= -29o t E=33.53 MeV qlab=25oflab=150o Pad resolution 12mm Bordeaux 2005
Track reconstruction 8He(p,t) @ 100 MeV By= 3T Isobutane r at STP Pad resolution 4mm Pad resolution 8mm 6He Elab=42.09 MeV qlab=17oflab= 72o t E=64.3 MeV qlab=19.5oflab=-107o Pad resolution 12mm Bordeaux 2005
Simulation Output Bordeaux 2005