AGATA Simulation Code ASC Overview Marc Labiche marc
AGATA Simulation Code (ASC) Overview Marc Labiche (marc. labiche@stfc. ac. uk) AGATA –GRETINA Workshop, ANL 5 th -8 th Dec 2016
Outline Ø ASC Generalities § Version / Distribution/ Working group § Update on Ancillary detectors § Validation status Ø Recent Developments § GPS & GDML capability § Simulated Time Stamped Data § g-ray b ackground simulation Ø Further developments in progress
ASC Generalities Ø ASC is a Native GEANT 4 simulation program Ø Almost as old as GEANT 4 Ø Initially developed & maintained by E. Farnea (University of Padova) Ø Current ASC version compatible with Geant 4. 10. 02. Ø Distribution: Ø SVN repository: http: //npg. dl. ac. uk/svn/agata Ø Schools/workshop Ø Nov. 2013 – AGATA Simulation mini-school (Orsay) Ø Slides: http: //npg. dl. ac. uk/AGATA/ Ø Dec. 2015 – ICC workshop on Simulation with AGATA & STo. GS (GANIL) Ø Slides: http: //npg. dl. ac. uk/ICC-simulation/ICCWorkshop. Sim 2015
ASC Generalities Ø Most recent contributors of the ASC Working Group: J. Ljungvall (CSNSM) Event Gen. and Time-Stamped Simulated Data C. Domingo-Pardo (IFIC Valencia) AGATA Double Cluster – Event Gen. (GSI) M. Ciemala (Krakow) PARIS & VAMOS G. Jaworski (Warsaw) EXOGAM & NEDA array B. Bruneels (CEA-IRFU) FATIMA D. Bloor (Uni. Of York) LYCCA detector + MOCADI external evt generator M. Reese (TU Darmstadt) External Evt Generator (plunger experiments) A. Goasduff (INFN - LNL) GALILEO M. Labiche (STFC Daresbury) Maintenance/Support & new G 4 features (GPS, GDML) , + ancillaries (FATIMA, AIDA …).
A glimpse at the SVN Repository Agata/ trunk/ branches/ Most up-to-date branch GANIL/ trunk/ GSI/ a_event_generator/ b_event_builder/ c_event_reconstruction/ ext_generators/ or Svn co http: //npg. dl. ac. uk/svn/agata/branches/GANIL/ trunk/
Agata + Ancillary Detectors Simulation of standalone ancillary detectors or combination of them can also be carried out. LYCCA Command: . /Agata –a Nbanc Idanc Fast Timing Array PARIS + Fast Timing Array
Current Ancillary Detectors Ancillaries Index Scoring ID Offset Koeln Si det. 1 1000 Shell (default) 2 2000 DANTE 3 3000 EUCLIDES 4 4000, 5000 Brick (PRISMA “dipole”) 6 6000 N-Wall 7 7000 DIAMANT 8 8000 EXOGAM 9 9000 HELENA 10 10000 RFD 11 11000 Cassandra 16 16000 AIDA 17 17000 Fast Timing Array 18 18000 PARIS 19 19000 LYCCA 22 22000
Other detectors available: GRETINA 8 Miniball Defined as different Ge crystal shapes ORGAM (Eurogam). /Agata –g 4 Defined as standalone detector but not as AGATA ancillaries GALILEO. /Agata –g 5
Why a Simulation Code ? Investigate Concept & Design of the Array Choose optimal design/configuration E. Farnea et al. NIM A 621 (2010) 331– 343 Develop and test analysis codes (tracking) Cost effective way to start and develop a project Preparation of experiments with different array configurations: Feasibility checks. Given the expect number of crystal available Number of crystals continuously increases but not yet reached 180. Nominal and Compact configuration Presence of ancillaries All has an impact on efficiency, P/T, tracking … Data Analyse & Interpretation of results
Ex: Effect of ancillaries on AGATA tracked spectrum ? Courtesy of Philipp R. John & Alain. Goasduff 190 W 138 Ba +2 p channel in 136 Xe+192 Os at 900 Me. V Doppler corrected using 138 Ba recoil, g rays: 1435. 8 and 462. 8 ke. V, (partner) g rays: 207, 357, 485, 591, 695 ke. V, + 20 g rays for low energy background 138 Ba 190 W
Ex: Decay lifetime effect on Fa. Tim. A efficiency v/c=10% Note: Chamber present but not shown
Ex: Simulated decay curves Acceptance effects on the decay curve when source is moving
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 Full energy peak efficiency at 1172 Ke. V when gating on 1332 Ke. V in Euroball Gated on 1332 ke. V: @ 1172 ke. V e. Exp (%) e. Sim (%) P/Texp (%) P/Tsim (%) Core Common 2. 38(2) 2. 84 18. 3(2) 22. 5 Calorimeter 3. 30(2) 4. 21 32. 2(3) 42. 5 So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 Full energy peak efficiency at 1172 Ke. V when gating on 1332 Ke. V in Euroball Gated on 1332 ke. V: @ 1172 ke. V e. Exp (%) e. Sim (%) P/Texp (%) P/Tsim (%) Core Common 2. 38(2) 2. 84 18. 3(2) 22. 5 Calorimeter 3. 30(2) 4. 21 32. 2(3) 42. 5 Actually, 23 crystals were considered in the simulations So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 Full energy peak efficiency at 1172 Ke. V when gating on 1332 Ke. V in Euroball Gated on 1332 ke. V: @ 1172 ke. V e. Exp (%) e. Sim (%) P/Texp (%) P/Tsim (%) Core Common 2. 38(2) 2. 55 18. 3(2) 22. 5 Calorimeter 3. 30(2) 3. 71 32. 2(3) 42. 5 Now for 21 crystals ! So: ~7% discrepancy on the Core Common Efficiency ~12% discrepancy on the Calorimeter Efficiency
Recent ASC developments GEANT 4 General Particle Source (GPS) § GPS allows the user to simulate realistic sources (ex: 60 Co, 152 Eu, § Now also available in the AGATA code: . /Agata –gps § Ex. Macros also available in the GANIL branch svn repostory: agata/branch/GANIL/trunk/macros/gps 152 Eu However: Issues with some intensities ! …)
Recent ASC developments GEANT 4 GDML capability tested § Implement CAD files of Mechanical Structure § Easy to include as “ancillary”: . /Agata –a N id § Honey. Comb is invoked as “ancillary” with the id number: 26 § GANIL vacuum chamber + target frame with the id number: 27 § OUPS plunger § Import/export detector geometry from/to other frameworks
Recent ASC developments CAD to GDML format with Fast. Rad software (free or commercial versions) http: //www. fastrad. net/ From a step file converted into a gdml file All 60 ATC + honeycomb structure Also: Free. Ca. D at http: //cad-gdml. in 2 p 3. fr/ GDML files can be quite big ! All GDML files are available on this git repository: https: //github. com/malabi/gdml-files Note: 28 MB gdml file
Simulated Efficiency & P/T curves (after tracking with mgt) Setup: - Nominal - 10 ATC+1 ADC at backward angles Ø Isotropic source (Mg=1) Ø + Vacuum Chamber
Simulated Efficiency curves Isotropic source with and without Chamber 32 crystals - Nominal - Large attenuation in chamber at low energy ( ~ -20 % at 121 ke. V) - Expected but now better estimate with realistic chamber - Similar effect of absorption/scattering in the Chamber in both operation mode (Calorimeter or after Tracking)
Simulated P/T curves Isotropic source with and without Chamber 32 crystals - Nominal - P/T ratio clearly also affected by presence of the chamber - Similar effect in both operation modes
Recent ASC developments Simulation with timestamp information Led by Joa Ljungvall (CSNSM) Why ? Need to take into account Source activity or Beam intensity and time structure: To estimate pile-up (and associated efficiency loss) and background in time spectrum To analyse Simulated and real data with same tools. How? As Geant 4 resets the clock to 0 at the beginning of each event we have to keep track of the time between events: We do it by setting the time to t 0 (Day, Hour, Min, …, nsec) at the beginning of the run, and propagate the time according to beam structure or source activity Choose option “–Gen”: . /Agata –Gen (see Agata. Alternative. Generator class)
Time Propagation: For a realistic source: . /build/Agata –Gen –b macros_adf/Co 60 Source/sourcerun. mac This macro uses the commands /Agata/file/info/enable. Time /Agata/generator/emitter/Set. Gamma. Ray. Source aa. El Activity(k. Bq) Ex: /Agata/generator/emitter/Set. Gamma. Ray. Source 60 Co 1000 = This opens a user define file called “ 60 Co. g 4 srcdata” ” and define the source activity to 1000 k. Bq
“ 60 Co. g 4 srcdata” format explained:
g-g matrices for 60 Co sources Courtesy of J. Ljungvall Gated on 1332 ke. V
Time propagation For a beam, use the commands: /Agata/file/info/enable. Time Intensity: pps, (if set ≤ 0 no time propagation) /Agata/generator/emitter/Set. Particle. Per. Seconds pps Beam bunch frequency: HF, (if set ≤ 0 Continuous beam) /Agata/generator/emitter/Set. Accelerator. HF HF units Bunch length: Bt, should be >0 /Agata/generator/emitter/Set. Width. Beam. Puls Bt units See examples in: macros_adf/Simulated experiments/
Also available with the –Gen option: Addition of g-ray backgrounds Discrete g rays: /Agata/generator/emitter/Add. Discrete. Gamma E N E: energy in Me. V N: number per event Exponential energy distribution: /Agata/generator/emitter/set. Slope. Gamma. Background S /Agata/generator/emitter/set. Max. EGamma. Background E /Agata/generator/emitter/set. Numberof. Gamma. Background N
Ex: Effect of ancillaries on AGATA tracked spectrum ? Courtesy of Philipp R. John & Alain. Goasduff 190 W 138 Ba +2 p channel in 136 Xe+192 Os at 900 Me. V Doppler corrected using 138 Ba recoil, g rays: 1435. 8 and 462. 8 ke. V, (partner) g rays: 207, 357, 485, 591, 695 ke. V, + 20 g rays for low energy background 138 Ba 190 W
Future developments Apply GEANT 4 Multi-threading feature. Develop ASC existing features to facilitate their transfer to other ( & more recent) frameworks STo. Gs , NPTool (SPIRAL 2) – AGATA GDML geometry yet implemented. ENSARRoot (FAIRRoot) All of these other frameworks are well connected to ROOT
Summary After the Concept and design of the AGATA array, the Geant 4 based ASC is used more and more to prepare and helps in the analysis of the experiment The Geant 4 based ASC continue to be well maintained and developed Main novelty: Timestamp added to the simulated data Simulated data real data after PSA Enriched by the addition of new “Ancillaries” and the new geant 4 GDML and GPS features Significant impact on Efficiency & P/T ratio at low energy are observed when CAD GDML converted geometry files are used.
Fin
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 @ 1172 ke. V e. Exp (%) e. Sim (%) P/Texp (%) P/Tsim (%) Core Common 2. 38(2) 2. 84(9) 18. 3(2) 22. 5(6) Calorimeter 3. 30(2) 4. 21(8) 32. 2(3) 42. 5(10) So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 @ 1172 ke. V e. Exp (%) e. Sim (%) P/Texp (%) P/Tsim (%) Core Common 2. 38(2) 2. 55(14) 18. 3(2) 23 Calorimeter 3. 30(2) 3. 71(17) 32. 2(3) 42 Now with 21 crystals
Recent developments GEANT 4 General Particle Source § GPS allows the user to simulate realistic sources (ex: 60 Co, 152 Eu, § Now also available in the AGATA code: . /Agata –gps § Ex. Macros also available in the GANIL branch svn repostory: agata/branch/GANIL/trunk/macros/gps GEANT 4 GDML capability added § Implement Mechanical structure § § Easy to include as ancillary: . /Agata –a N id § Honey. Comb is invoked as ancillary with id number: 26 § Ganil Vacuum chamber is invoked as ancillary with id number: 27 § OUPS (unique id not yet attributed) Import/export detector from/to other framework …)
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 Full energy peak efficiency @ 1172 ke. V e. Exp (%) e. Sim (%) Core Common 2. 38(2) 2. 84(9) Calorimeter 3. 30(2) 4. 21(8) So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 Full energy peak efficiency @ 1172 ke. V e. Exp (%) e. Sim (%) Core Common 2. 38(2) 2. 84(9) Calorimeter 3. 30(2) 4. 21(8) Actually, 23 crystals were considered in the simulations So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency
Simulation validation GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258 -260 Full energy peak efficiency @ 1172 ke. V e. Exp (%) e. Sim (%) Core Common 2. 38(2) 2. 55(14) Calorimeter 3. 30(2) 3. 71(17) Now for 21 crystals ! So: ~7% discrepancy on the Core Common Efficiency ~12% discrepancy on the Calorimeter Efficiency
Time propagation with beam time structure Intensity: pps, (if set ≤ 0 no time propagation) /Agata/generator/emitter/Set. Particle. Per. Seconds pps Beam bunch frequency: HF, (if set ≤ 0 Continuous beam) /Agata/generator/emitter/Set. Accelerator. HF HF units Bunch length: Bt, should be >0 /Agata/generator/emitter/Set. Width. Beam. Puls Bt units
Time propagation with beam time structure Ex: beam structure for 1010 pps, beam repetition rate of 100 MHz and beam bunches of 1 ns(red) and 2 ns(green) Courtesy of J. Ljungvall
Recent ASC developments Ex: g-g matrices for different source activities
Event generators included: Fusion-evaporation /Agata/generator/emitter/Beam. Out/Set. Pfe P 1 Coulomb excitation: /Agata/generator/emitter/Beam. Out/Set. Pclx P 2 Fusion-Fission excitation: /Agata/generator/emitter/Beam. Out/Set. Pff P 3 Multi-nucleon transfer: /Agata/generator/emitter/Beam. Out/Set. Ptr P 4 (Pi = probability)
- Slides: 41