1 / 41

Marc Labiche (marc.labiche@stfc.ac.uk) AGATA –GRETINA Workshop, ANL 5 th -8 th Dec 2016

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

veronicathomas
Download Presentation

Marc Labiche (marc.labiche@stfc.ac.uk) AGATA –GRETINA Workshop, ANL 5 th -8 th Dec 2016

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. AGATA Simulation Code(ASC) Overview Marc Labiche (marc.labiche@stfc.ac.uk) AGATA –GRETINA Workshop, ANL 5th -8th Dec 2016

  2. 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

  3. ASC Generalities • ASC is a Native GEANT4 simulation program • Almost as old as GEANT4 • Initially developed & maintained by E. Farnea (University of Padova) • Current ASC version compatible with Geant4.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 & SToGS (GANIL) • Slides: http://npg.dl.ac.uk/ICC-simulation/ICCWorkshopSim2015

  4. 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 G4 features (GPS, GDML) , + ancillaries (FATIMA, AIDA …) .

  5. A glimpse at the SVN Repository Agata/ trunk/ Most up-to-date branch branches/ GANIL/ trunk/ GSI/ a_event_generator/ trunk/ b_event_builder/ c_event_reconstruction/ ext_generators/ Svn co http://npg.dl.ac.uk/svn/agata or Svn co http://npg.dl.ac.uk/svn/agata/branches/GANIL/

  6. Agata + Ancillary Detectors Simulation of standalone ancillary detectors or combination of them can can also be carried out. LYCCA PARIS + Fast Timing Array Fast Timing Array Command: ./Agata –a Nbanc Idanc

  7. Current Ancillary Detectors

  8. Other detectors available: GRETINA 8 Miniball Defined as different Ge crystal shapes ORGAM (Eurogam) ./Agata –g 4 GALILEO ./Agata –g 5 Defined as standalone detector but not as AGATA ancillaries

  9. 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

  10. Ex: Effect of ancillaries on AGATA tracked spectrum ? Courtesy of Philipp R. John & Alain. Goasduff 190W 138Ba +2p channel in 136Xe+192Os at 900MeV Doppler corrected using 138Ba recoil, 138Ba g rays: 1435.8 and 462.8 keV, 190W (partner) g rays: 207, 357, 485, 591, 695 keV, + 20 g rays for low energy background

  11. Ex: Decay lifetime effect on FaTimA efficiency v/c=10% Note: Chamber present but not shown

  12. Ex: Simulated decay curves Acceptance effects on the decay curve when source is moving

  13. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Full energy peak efficiency at 1172 KeV when gating on 1332 KeV in Euroball Gated on 1332 keV: So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency

  14. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Full energy peak efficiency at 1172 KeV when gating on 1332 KeV in Euroball Gated on 1332 keV: Actually, 23 crystals were considered in the simulations So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency

  15. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Full energy peak efficiency at 1172 KeV when gating on 1332 KeV in Euroball Gated on 1332 keV: Now for 21 crystals ! So: ~7% discrepancy on the Core Common Efficiency ~12% discrepancy on the Calorimeter Efficiency

  16. Recent ASC developments GEANT4 General Particle Source (GPS) • GPS allows the user to simulate realistic sources (ex: 60Co, 152Eu, …) • Now also available in the AGATA code: ./Agata –gps • Ex. Macros also available in the GANIL branch svnrepostory: agata/branch/GANIL/trunk/macros/gps 152Eu However: Issues with some intensities !

  17. Recent ASC developments GEANT4 GDML capability tested • Implement CAD files of Mechanical Structure • Easy to include as “ancillary”: ./Agata –a N id • HoneyComb is invoked as “ancillary” with the id number: 26 • GANIL vacuum chamber + target framewith the id number: 27 • OUPS plunger • Import/export detector geometry from/to other frameworks

  18. Recent ASC developments From a step file converted into a gdml file CAD to GDML format with FastRad software (free or commercial versions) http://www.fastrad.net/ Also: FreeCaD at http://cad-gdml.in2p3.fr/ GDML files can be quite big ! All GDML files are available on this git repository: https://github.com/malabi/gdml-files All 60 ATC + honeycomb structure Note: 28MB gdml file

  19. Simulated Efficiency & P/T curves (after tracking with mgt) • Setup: - Nominal • - 10ATC+1ADC • at backward angles • Isotropic source (Mg=1) • + Vacuum Chamber

  20. Simulated Efficiency curves Isotropic source with and without Chamber 32 crystals - Nominal - Similar effect of absorption/scattering in the Chamber in both operation mode (Calorimeter or after Tracking) - Large attenuation in chamber at low energy ( ~ -20 % at 121 keV) - Expected but now better estimate with realistic chamber

  21. 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

  22. Recent ASC developments Simulation with timestamp information Led by JoaLjungvall (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 Geant4 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 t0 (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 AgataAlternativeGenerator class)

  23. Time Propagation: • For a realistic source: ./build/Agata –Gen –b macros_adf/Co60Source/sourcerun.mac • This macro uses the commands /Agata/file/info/enableTime /Agata/generator/emitter/SetGammaRaySourceaaEl Activity(kBq) Ex: /Agata/generator/emitter/SetGammaRaySource60Co1000 = This opens a user define file called “60Co.g4srcdata” ” and define the source activity to 1000kBq

  24. “60Co.g4srcdata” format explained:

  25. g-g matrices for 60Co sources Courtesy of J. Ljungvall Gated on 1332 keV

  26. Time propagation • For a beam, use the commands: • /Agata/file/info/enableTime • Intensity: pps, (if set ≤ 0 no time propagation) /Agata/generator/emitter/SetParticlePerSecondspps • Beam bunch frequency: HF, (if set ≤ 0 Continuous beam) /Agata/generator/emitter/SetAcceleratorHF HF units • Bunch length: Bt, should be >0 /Agata/generator/emitter/SetWidthBeamPulsBtunits • See examples in: macros_adf/Simulated experiments/

  27. Also available with the –Gen option: • Addition of g-ray backgrounds • Discrete g rays: /Agata/generator/emitter/AddDiscreteGamma E N E: energy in MeV N: number per event • Exponential energy distribution: /Agata/generator/emitter/setSlopeGammaBackground S /Agata/generator/emitter/setMaxEGammaBackground E /Agata/generator/emitter/setNumberofGammaBackgroundN

  28. Ex: Effect of ancillaries on AGATA tracked spectrum ? Courtesy of Philipp R. John & Alain. Goasduff 190W 138Ba +2p channel in 136Xe+192Os at 900MeV Doppler corrected using 138Ba recoil, 138Ba g rays: 1435.8 and 462.8 keV, 190W (partner) g rays: 207, 357, 485, 591, 695 keV, + 20 g rays for low energy background

  29. Future developments • Apply GEANT4 Multi-threading feature. • Develop ASC existing features to facilitate their transfer to other ( & more recent) frameworks • SToGs , NPTool (SPIRAL 2) – AGATA GDML geometry yet implemented. • ENSARRoot (FAIRRoot) All of these other frameworks are well connected to ROOT

  30. Summary • After the Concept and design of the AGATA array, the Geant4 based ASC is used more and more to prepare and and helps in the analysis of the experiment • The Geant4 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 geant4 GDML and GPS features • Significant impact on Efficiency & P/T ratio at low energy are observed when CAD GDML converted geometry files are used.

  31. Fin

  32. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency

  33. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Now with 21 crystals

  34. Recent developments GEANT4 General Particle Source • GPS allows the user to simulate realistic sources (ex: 60Co, 152Eu, …) • Now also available in the AGATA code: ./Agata –gps • Ex. Macros also available in the GANIL branch svnrepostory: agata/branch/GANIL/trunk/macros/gps GEANT4 GDML capability added • Implement Mechanical structure • Easy to include as ancillary: ./Agata –a N id • HoneyComb 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

  35. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Full energy peak efficiency So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency

  36. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Full energy peak efficiency Actually, 23 crystals were considered in the simulations So: ~20% discrepancy on the Core Common Efficiency ~33% discrepancy on the Calorimeter Efficiency

  37. Simulation validation • GSI source run with 21 crystals N. Lalovic et al. NIM A 806 (2016) 258-260 Full energy peak efficiency Now for 21 crystals ! So: ~7% discrepancy on the Core Common Efficiency ~12% discrepancy on the Calorimeter Efficiency

  38. Time propagation with beam time structure • Intensity: pps, (if set ≤ 0 no time propagation) /Agata/generator/emitter/SetParticlePerSecondspps • Beam bunch frequency: HF, (if set ≤ 0 Continuous beam) /Agata/generator/emitter/SetAcceleratorHF HF units • Bunch length: Bt, should be >0 /Agata/generator/emitter/SetWidthBeamPulsBt units

  39. Time propagation with beam time structure Ex: beam structure for 1010pps, beam repetition rate of 100MHz and beam bunches of 1ns(red) and 2ns(green) Courtesy of J. Ljungvall

  40. Recent ASC developments Ex: g-g matrices for different source activities 10 MBq60Co 100 KBq60Co Eg[keV] Eg[keV] Eg[keV] Eg[keV]

  41. Event generators included: • Fusion-evaporation /Agata/generator/emitter/BeamOut/SetPfe P1 (Pi = probability) • Coulomb excitation: /Agata/generator/emitter/BeamOut/SetPclx P2 • Fusion-Fission excitation: /Agata/generator/emitter/BeamOut/SetPff P3 • Multi-nucleon transfer: /Agata/generator/emitter/BeamOut/SetPtr P4

More Related