ACTARSim: a simulation code for ACtive TARgets H. Alvarez-Pol - - PowerPoint PPT Presentation

actarsim a simulation code for active targets
SMART_READER_LITE
LIVE PREVIEW

ACTARSim: a simulation code for ACtive TARgets H. Alvarez-Pol - - PowerPoint PPT Presentation

ACTARSim: a simulation code for ACtive TARgets H. Alvarez-Pol GENP, Univ. de Santiago de Compostela on behalf of the ACTAR-TPC collaboration ICC/SPIRAL2 Meeting on Common Simulation & Data Analysis Framework, Orsay, 26th November 2014


slide-1
SLIDE 1

ACTARSim: a simulation code for ACtive TARgets

  • H. Alvarez-Pol

GENP, Univ. de Santiago de Compostela

  • n behalf of the ACTAR-TPC collaboration

ICC/SPIRAL2 Meeting on Common Simulation & Data Analysis Framework, Orsay, 26th November 2014

slide-2
SLIDE 2

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

2

Context ACTARSim: a simulation code for ACTAR TPC

Topics covered:

  • Active Targets.
  • The ACTAR TPC program.
  • ActarSim code in detail.
  • Tests and status.
  • Ongoing development.

Work supported by Xunta de Galicia “Proxectos Plan Galego IDT”, project 2013-PG015: “Física de núcleos exóticos con Detectores Activos”

for the study of 20Mg(p,p) and 20Mg(p,p') resonant reactions for the spectroscopy of 21Al, exotic nuclei in the rp-process by transfer (57Cu, 61Ga, 73Rb, 77Y), and fission after transfer reactions.

slide-3
SLIDE 3

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

3

The concept behind ACtive TARgets In short: the detection gas is the target itself.

  • High geometric efficiency.
  • High thickness. High luminosity.
  • Control on the reaction energy. Very low threshold.
  • Full 3D tracking of participant particles.
  • In many cases, close to zero background. Ideal for low statistics

reactions and short-lived species. Physics cases particularly suited:

  • Resonant reactions.
  • One- and multi-nucleon transfer on light nuclei reactions or fission.
  • Inelastic scattering to GRs (GMR, GQR, ...).
  • Nuclear astrophysics.
  • Search for exotic decays, delayed modes, ...

ACtive TARget introduction

slide-4
SLIDE 4

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

4

ACtive TARget precursors

CENBG TPC (2005):

  • For two-proton radioactivity.
  • Energies and angles between protons.
  • GEM's + 200 μm pitch μ-grooves.

TACTIC (2007):

  • Astrophysical interest reactions.
  • Ready for high intensities (107 ions/s).
  • Target gas isolated and interaction

vertex ionization not reachable.

ACtive TARget precursors in nuclear physics

K.A.Chipps et al. Proceedings of Science (NIC XI) 261 (2010) B.Blank et al. Nucl. Instrum. Meth. Phys. Res. A 613, 65 (2010)

slide-5
SLIDE 5

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

5

MAYA Active Target

MAYA main characteristics:

  • Box-like geometry with ancillary detectors (Si and CsI[Tl]) in the front wall.
  • 1024 hexagonal pads with Frisch grid and proportional wire amplification.
  • Versatile multipurpose design.

ACtive TARget precursors in nuclear physics

C.E.Demonchy et al. Nucl. Instrum. Meth. Phys. Res. A 583, 341 (2007)

slide-6
SLIDE 6

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

6

MAYA Active Target:

Physics highlights:

  • Structure for unbound isotopes: analysis of 7H (most exotic nuclei ever).

ACtive TARget precursors in nuclear physics

M.Caamano et al. Phys. Rev. Lett. 99, 062502 (2007)

slide-7
SLIDE 7

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

7

MAYA Active Target:

Physics highlights:

  • Structure for unbound isotopes: analysis of 7H (most exotic nuclei ever).
  • First study of 11Li 2-neutron halo via transfer reactions at low energy.

ACtive TARget precursors in nuclear physics

I.Tanihata et al. Phys. Rev. Lett. 100, 192502 (2008) T.Roger et al. Phys. Rev. C. 79, 031603 (2009)

slide-8
SLIDE 8

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

8

ACtive TARget precursors in nuclear physics

C.Monrozeau et al. Phys. Rev. Lett. 100, 042501 (2008)

  • M. Vandebrouck et al., Accepted in Phys. Rev. Lett (2014)

MAYA Active Target:

Physics highlights:

  • Structure for unbound isotopes: analysis of 7H (most exotic nuclei ever).
  • First study of 11Li 2-neutron halo via transfer reactions at low energy.
  • Giant Resonances analysis for Ni isotopes

(56Ni(α,α’), 68Ni(α, α’), 68Ni(d,d’)).

  • ...
  • Future measurements: transfer induced fission (starting next week).
slide-9
SLIDE 9

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

9

The ACTAR TPC project The ACTAR-TPC detector

Leading a new generation of Active target devices:

  • Overcomes many of the limitations with present devices.
  • Funded via an ERC Starting grant (2014-2019).

Physics cases:

  • One and two nucleon transfer reactions.
  • Rare and exotic nuclear decay (2p, β-2p, …).
  • Transfer-induced fission.
  • Inelastic scattering and giant resonances.
  • Resonant scattering and astrophysics.

Detector Design:

  • Amplification = MICROMEGAS (+ GEMs).
  • Pad sizes of 2x2 mm2: 16384 channels.
  • ANR General Electronics for TPC’s (GET).
  • Improved data throughput + internal trigger.
slide-10
SLIDE 10

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

10

ACTARSim: simulation code for ACTAR- TPC ACTARSim (since 2005)

ActarSim is a Geant4 + ROOT application for ACtive TARgets simulation.

  • Developed for the ACTAR TPC design and MAYA analysis comparison.
  • Initial development at USC, maintained and extended in GANIL since 2008.
  • New developments ongoing during ACTAR-TPC construction period,

responsabilities back to USC (2013). GEANT4 is used for the production and tracking of primaries and secondaries above Ecut = 1 keV. ROOT is used for the calculation of the drift and diffusion

  • f the electronic clouds, the induction in the pad plane, the visualization, ...
slide-11
SLIDE 11

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

11

ACTARsim main characteristics:

Detectors geometry:

  • Cube (MAYA or ACTAR-TPC).
  • Cylinder (drift on endcaps or sides).
  • Easy size modification for design.

ACTARSim in detail

slide-12
SLIDE 12

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

12

ACTARsim main characteristics:

Detectors geometry:

  • Cube (MAYA or ACTAR-TPC).
  • Cylinder (drift on endcaps or sides).
  • Easy size modification for design.

Gas and pressure:

  • Several gas predefined.
  • Presure, density, … given at runtime.

Required for the Geant4 initialization.

ACTARSim in detail

slide-13
SLIDE 13

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

13

ACTARSim in detail ACTARsim main characteristics:

Detectors geometry:

  • Cube (MAYA or ACTAR-TPC).
  • Cylinder (drift on endcaps or sides).
  • Easy size modification for design.

Gas and pressure:

  • Several gas predefined.
  • Presure, density, … given at runtime.

Required for the Geant4 initialization. Ancillary detectors:

  • Configurable silicon DSSD layers.
  • Configurable CsI[Tl] detectors.
slide-14
SLIDE 14

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

14

ACTARSim in detail ACTARsim Geant4 part:

Kinematics calculator integrated:

  • Event generator with proper vertex

calculation based on kinematics.

slide-15
SLIDE 15

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

15

ACTARSim in detail ACTARsim Geant4 part:

Kinematics calculator integrated:

  • Event generator with proper vertex

calculation based on kinematics. Physics interaction from Geant4:

  • Realistic description of discrete and

continuum processes.

  • User selectable physics packages.
  • Realistic beam-gas interaction with

gaussian beam profile.

  • Beam shielding tube possible.
slide-16
SLIDE 16

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

16

ACTARSim in detail ACTARsim Geant4 part:

Kinematics calculator integrated:

  • Event generator with proper vertex

calculation based on kinematics. Physics interaction from Geant4:

  • Realistic description of discrete and

continuum processes.

  • User selectable physics packages.
  • Realistic beam-gas interaction with

gaussian beam profile.

  • Beam shielding tube possible.

That's all for the Geant4 part!

  • The result is a TTree (TClonesArray) of strides (groups of steps) per particle.
  • Energy loss used for the calculation of the ionization in the gas.
  • The TTree contains the Hits in the ancillary detectors.
slide-17
SLIDE 17

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

17

ACTARSim in detail ACTARsim ROOT part:

Drift in electric and magnetic fields:

  • Constant electric field included.
  • Drift parameters selected by the user.
  • Drift and difussion is calculated using

the electric maps.

slide-18
SLIDE 18

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

18

ACTARSim in detail ACTARsim ROOT part:

Drift in electric and magnetic fields:

  • Constant electric field included.
  • Drift parameters selected by the user.
  • Drift and difussion is calculated using

the electric maps. Amplification and pad induction:

  • GEM, micromegas and wires

amplification schemes are included, as well as the induction in the pad plane.

  • Variable size, hexagonal- or square-

shaped pad planes included.

slide-19
SLIDE 19

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

19

ACTARSim in detail ACTARsim ROOT part:

Drift in electric and magnetic fields:

  • Constant electric field included.
  • Drift parameters selected by the user.
  • Drift and difussion is calculated using

the electric maps. Amplification and pad induction:

  • GEM, micromegas and wires

amplification schemes are included, as well as the induction in the pad plane.

  • Variable size, hexagonal- or square-

shaped pad planes included. Induction calculation and output:

  • The induced charge in each pad is calculated including timing.
  • Results are stored in a TTree (TClonesArray) of pads with signal.
slide-20
SLIDE 20

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

20

ACTARSim visualization ACTARsim visualization and analysis

77Ni (d,p) 78Ni, few events

slide-21
SLIDE 21

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

21

ACTARSim visualization ACTARsim visualization and analysis

A high resolution case with 1mm square pad in a box-like active target

slide-22
SLIDE 22

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

22

ACTARSim current work Sistematic energy loss evaluation:

Evaluation of the energy loss properties of the physics libraries in Geant4:

  • Energy loss obtained from the step energy loss sum.
  • Energy straggling as RMS of the energy loss distribution.
  • Angular straggling given by the step angle at the end of the absorber.
  • Lateral spread given by the sigma of the position distribution (x-position in

the XZ plane) at the end of the absorber.

  • D2, H2, iC4H10 as target gases.
  • p, 4He, 12C, 24Mg, 56Fe as projectiles.
  • Geant4 libraries: emstandard_opt3 and ionGasModels.

Energy loss and energy straggling match pretty good in all cases, angular straggling and lateral spread match better at higher energy for heavier ions.

Evaluation performed by Piotr Konczykowski (USC)

slide-23
SLIDE 23

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

23

ACTARSim current work Sistematic energy loss evaluation:

12C on deuterium gas STP:

Evaluation performed by Piotr Konczykowski (USC)

slide-24
SLIDE 24

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

24

ACTARSim current work Sistematic energy loss evaluation:

Proton on deuterium gas STP:

Evaluation performed by Piotr Konczykowski (USC)

slide-25
SLIDE 25

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

25

ACTARSim current work Sistematic energy loss evaluation:

Proton on iC4H10 gas STP:

Evaluation performed by Piotr Konczykowski (USC)

slide-26
SLIDE 26

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

26

ACTARSim current work Sistematic energy loss evaluation:

56Fe on iC4H10 gas STP:

Evaluation performed by Piotr Konczykowski (USC)

slide-27
SLIDE 27

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

27

ACTARSim current work ACTARSim ongoing studies:

  • ACTAR-TPC demonstrator.
  • Range analysis and comparison for ions at low energies in deuterium and
  • ther gases of interest for the Physical Cases.
  • Complete tests of the Geant4 physics libraries and comparison of the

longitudinal range, straggling and radial range.

  • Merge code for MAIKoSim (Osaka Active Target) (MAIKoSim was a branch

from ActarSim previous version made by Yassid Ayyad, RCNP Osaka Univ.).

  • Integration of cosmic rays generator for the analysis of ACTAR-TPC

demonstrator laboratory data.

  • Documentation.
slide-28
SLIDE 28

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

28

ACTARSim current work ACTAR-TPC demonstrator in ACTARSim:

slide-29
SLIDE 29

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

29

ACTARSim GitHub repository

https://github.com/hapol/ActarSim/

slide-30
SLIDE 30

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

30

People involved ACTAR-TPC collaborators:

slide-31
SLIDE 31

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

31

ADDITIONAL SLIDES

slide-32
SLIDE 32

ACTARSim: a simulation for ACTAR TPC - H. Alvarez Pol Orsay, 26th November 2014

32

Additional information ACTARSim and Physics Simulations

Reactions (as described in ERC Starting Grant):

  • Commissioning: single-nucleon transfer reactions with 19Ne (start of rp

process producing heavy neutron-deficient nuclei) at GANIL:

  • 19Ne(d,3He)18F: test reconstruction capabilities for 18F excited states

energies, spins and parities though the angular distributions.

  • 19Ne(d,t)18Ne: test of triton energy deposition (~1 keV/mm). Aditional

interest of 18Ne as 2-proton emitter.

  • 19Ne(d,3He)18F: software (GET) trigger test, due to kinematics, allowing

external forward detection to test triggers efficiency.

  • Two-proton decay of 54Zn at LISE.
  • Single nucleon transfer in neutron-rich:80Zn(d,p)81Zn as a preparation for

78Ni(d,p)79Ni at HIE-ISOLDE.