Marc Labiche ASW CSNSM, 25 th -26 th Nov. 2013

Transcription

1 Marc Labiche ASW CSNSM, 25 th -26 th Nov. 2013

2 À la mémoire de Enrico Enrico Farnea ( )

3 Some generalities AGATA Code is based fully on GEANT4 transport code. Compile on different GEANT4 version From Geant4.7 (although not tested recently) up to Geant4.9.6 Since Geant4.9.5, compilation possible with cmake or make GEANT4.10 is coming soon ~(next month) make may no longer be available for compilation cmake only

4 The AGATA Simulation team Core team: Cesar Domingo-Pardo, Joa Ljungvall, Johan Nyberg, and myself (convener) Extra support for the different campaigns/ancillary detectors: - At GANIL: Marcin Palazc (N-Wall-NEDA), Grzegorz Jaworski (EXOGAM), Bart Bruyneel (FATIMA), - At GSI: Dan Bloor (LYCCA+AGATA), M. Reese (Ext evt generator for plunger), - At Legnaro: Dino Bazzacco,

5 The team roles/activities: Develop: - Code general improvement - Integration of Ancillary detectors - Event generators - Background - Preliminary analysis codes Assist: - Installation - Running - Bug fixing Maintain: - Compatibility with GEANT4 distributions Perform simulations: - Of general interests (Source or Commissioning runs) - Of experiment proposals, Occasionally

6 ASC Distribution Can be checked out here: Subversion is used for the control version, so to check out type: svn co (for the full version) svn co (for the original version) svn co (for GANIL version) svn co (for GSI version) To received updates, the command is: cd../trunk/ svn update To commit changes, you ll need login & password (just ask me)

7 AGATA repository structure Agata/ trunk/ branches/ GANIL/ trunk/ GSI/ a_event_generator/ b_event_builder/ trunk/ c_event_reconstruction/ Other branches can be added (LYCCA, ) trunk/ Currently each of the trunk/ directories (from top to bottom) are copies of Enrico s In preparation

8 Also distributed: The repository contains an INSTALL file For both GEANT4 and ASC. For GEANT4, the easiest way is to download the virtual machine from VMWare player or VirtualBox a new directory with some external basic generators. Produce ASCII input for ASC with gamma-emitter momentum and gamma-decay cascade Basic reaction available: Coulex, knockout (fragmentation), fusion-evaporation, transfer, elastic/inelastic Also valid for multi-layer target (plunger experiment). See: README, HOWTO

9 Also distributed: The two tracking algorithms: BUT no manuals run with default parameters. OFT (contact: MGT (contact: Dino.Bazzacco:pd.infn.it) Both algorithms are in the directory: trunk/analysis Many examples to run simulation are available in: trunk/macros/ tunk/macros_ganil/ trunk/macros_pespec/

10 The ASC manual In each trunk/docs/ - the ASC manual is available: agatamanual.pdf (Enrico) It needs to be updated for more recent developments but still completely valid.

11 The agata/ repository

12 The branches

13 The Agata/branches/GANIL/trunk/ directory Geant4-based AGATA core program (Leave it to the developpers) Compilation files for make (GNUMakefile) and cmake (CMakeLists.txt, build/)

14 The AGATA Cmakefile.list Produces the executable file Agata and saves a copy of the configurations files & macros in a build directory (e.g: build1). build can be put any where on your system. you can then work in this directory: create new macros or config. Files run your simulations in that build directory. Suggestion: when a new version of the AGATA code comes along (ie: after an, svn update) You can create a new build directory (e.g build2) and compile the new version there. If satisfied, you can delete build1!! The current Cmakelist.txt works but still under construction Not all CPPFLAGS are set yet. The config. & macros files not yet copied to the build directory

15 The Agata/branches/GANIL/trunk/ directory Geometry config files for AGATA

16 The files A180/A180eulerxxx.list Probably the only file you may want/need to modify To add new Triple clusters (ATC) or double clusters (ADC) in the geometry And to position them as you which Note: For full array use: A180euler.list Ex: A180/A180euler5D5T.list # The Euler angles (degree) and shifts (mm) of the 15 clusters # cl cl# psi(rz) theta(ry) phi(rz) dx dy dz ATC ADC

17 The Agata/branches/GANIL/trunk/ directory Geometry config files for AGATA Geometry files for some ancillary det.

18 The Agata/branches/GANIL/trunk/ directory Geometry config files for AGATA Geometry files for some ancillary det. EXOGAM input configuration file

19 The EXOGAM input file In the file Exogam_input :. # CONFIGURATION X CONFIGURATION B # # Id Large/Small In/Out T->Ge can Theta Phi # Clover (cm) (deg.) (deg.) # #Flange 2: 0 L #Flange 4: 1 L #Flange 3: 2 L #Flange 5: 3 L #Flange 10: 4 L #Flange 11: 5 L #Flange 12: 6 L #Flange 13: 7 L #Flange 6: 8 L #Flange 7: 9 L #Flange 8: 10 L #Flange 9: 11 L #Flange 15: 12 L #Flange 1: 13 L

20 The Agata/branches/GANIL/trunk/ directory Geometry config files for AGATA EXOGAM input configuration file Geometry files for some ancillary det. Examples of macros

21 The Agata/branches/GANIL/trunk/ directory mgt and OFT tracking code

22 The Agata/branches/GANIL/trunk/ directory mgt and OFT tracking code Event input file examples

23 The Agata/branches/GANIL/trunk/ directory mgt and OFT tracking code Event input file examples Agata code manual

24 The Agata/branches/GANIL/trunk/ directory mgt and OFT tracking code Event input file examples Agata code manual

25 On the Geant4 VM only Compilations with GNUMakefile In a terminal window or in your.cshrc file, run the commands: source /usr/local/geant4.9.6.p02/share/geant /geant4make/geant4make.csh /usr/local/geant4.9.6.p02/share/geant /geant4make/ (Note:there is a space between../geant4make.csh and /usr/local/.) unsetenv G4UI_USE_QT unsetenv G4VIS_USE_OPENGLQT unsetenv G4UI_USE_XM unsetenv G4VIS_USE_OPENGLXM In the trunk/ directory, compile agata with the command: make clean make The executable Agata should then be created in the directory $G4WORKDIR/bin/Linuxg++/ To execute the program run: $G4WORKDIR/bin/Linux-g++/Agata

26 Compilations with CMakefile on VM In trunk directory, create a build directory and enter it: mkdir Build_Agata cd Build_Agata Build with cmake: cmake DGeant4_DIR=$G4LIB/Geant / Compile with make: make The executable Agata should then be created in this Build_Agata/ directory. To run the code, from the trunk directory:./build_agata/agata

27 Time for a break

28 Running the Agata code The options available The AGATA commands With built-in generators What generators? Simulation Output file description The other generator options Simulation Input file description Using Ancillary detectors

29 The options >$G4WORKDIR/bin/Linux-g++/Agata -h (or./build/agata -h) -g index select geometry (0=AGATA, 1=SHELL, 2=Simple, 3=Geom for test of Physics 4=Orgam (eurogam)) -a N_anc n_1... n_n select ancillary devices -n use also neutrons, protons,... -nole use the standard treatment of the interactions for gamma/e-/e+ -nolowh use the standard treatment of the interactions for hadrons -p consider the linear polarization of the photons -Vol enables the use of geantinos for volume and solid angle calculation -Ext reads the event sequence from file -Gen uses user-supplied event generation -Test Event generator using geant4 gun commands -b file.mac runs in batch mode executing file.mac -Path <dir> sets base directory to <dir> -seed uses current time as seed for random number engine -run r sets first run number to r -SN Use G4ScreenedNuclear instead of multiplescattering for Generic Ions -decay zzz.aaaa make radioactive decayfiles -AddOn library_name Add some kind of det. via shared lib -h help, print this list

30 The AGATA commands >$G4WORKDIR/bin/Linux-g++/Agata (or./build/agata) Idle> ls Idle> ls /Agata/ Idle> ls /Agata/generator/ Idle> ls /Agata/generator/gamma/ List of all interactive commands (IC) List of Agata/ IC

31 Generators Built-in >./path_to/agata External From event input files >./path_to/agata Ext From user-defined generators (on-the-fly) >./path_to/agata -Gen

32 Built-in generators (1) Many options (including v/c recoil & polarisation) For the energy of the emitted g, in CM: 7 types (0 to 6), set with the command: /Agata/generator/gamma/gunType type Defined in class src/agataemitted.cc type =0 : monochromatic (=default) type =1 : rotational band type =2 : discreet energies read from file type =3 : flat energy distribution type =4 : sampled from a continuous spectrum type =5 : same as 2, with intensity normalisation to 1 (100%) type =6 : monochromatic + E2 angular distribution

33 Built-in generators (2) Recoil velocity: v/c value (in %) set with the command: /Agata/generator/recoil/beta value Other parameters for /Agata/generator/recoil/ are: sigma * Define sigma of the gaussian distribution of recoil v/c. openingangle * Define opening angle for the recoil cone. gaussangle * Set opening angle for the recoil cone to a gaussian distribution. uniformangle * Set opening angle for the recoil cone to a uniform distribution. direction * Define a fixed direction of the recoil. thetaphi * Define a fixed direction of the recoil. sourcetype * Define behaviour of the source. positionofsource * Define initial position of the source.

34 Output file description (1) The header file: AGATA OUTPUT_MASK # ndet Energy AbsolutePosition(x y z) RelativePosition(x' y' z') # RelativePosition for mgs (x'' y'' z'') nseg time interaction # 0--> disabled; 1--> enabled DISTFACTOR 1 ENERFACTOR 1 G4TRACKING 0 DATE Sat Nov 23 21:00: GEOMETRY AGATA SUMMARY TRANSFORMATION PASSIVE 1 CAPSULES 1 ENDGEOMETRY GENERATOR 0 ENDGENERATOR GAMMA RECOIL SOURCE $ = Information written for each hits = Some scaling factors = Tracking method = Date & time BEGINNIN OF GEOMETRY = Detector name = Rmax, Rmin, Ncrist., Ncl, Nslice(Ncl), Nseg(Ncl),.. = Translation & Rotation parameters = Passive material implemented = Capsule included END OF GEOMETRY = Generator type = emitted particle multiplicity = kinetic energy = beta; beta sigma; direction(x,y,z); Op. Angle = type; motion; position(x,y,z)

35 Output file description (2) The output events in GammaEvents.xxxx gamma Crystal # Incident energy and direction (dx,dy,dz) Evt # Segment # Deposited energy and absolute Interaction Position

36 Using external event input file Can be produced by hand or via gammaware, Ext_generators/, ASCII format: Example: Events00 file in Agata/events/: FORMAT 0 0 = Emitter info, Emitted info # (Level of detailed information: 0 to 4) # REACTION = Z beam, A beam, Z target, A Target, E beam # (Default= 1, 1, 6, 12, 0) # EMITTED = N total, Type i=0,, Type i=n-1 # $ = Beginning of event = Emitter: z, A, E, x Dir, y Dir,z Dir, x S, y S,z S = 1 st emitted: E, x Dir, y Dir,z Dir, x S, y S,z S, [t, P] = 2 nd emitted: E, x Dir, y Dir,z Dir, x S, y S,z S, [t, P] index Type 1 gamma 2 neutron 3 proton 4 deuterium 5 Triton 6 3He 7 alpha 8 Generic ion 97 Electron 98 Positron 99 Geantino

37 AGATA Ancillary detectors RFD EUCLIDES DIAMANT AIDA EXOGAM FATIMA

38 The list of Ancillaries so far Ancillaries Index Sensitive Det. instances Offset Koeln Si det Shell (default) DANTE EUCLIDES ,5000 Brick (PRISMA dipole ) N-Wall DIAMANT EXOGAM HELENA RFD TRACE? ,13000 CUP? GASPARD? Cassandra AIDA FATIMA

39 Multi-detector example By running the command: Agata a Nb of ancillaries EXOGAM id DIAMANT id N-Wall id

40 Exogam event output description Incoming gamma 9000: ExoGam offset 400: BGO shield offset 36: Flange_ID *4 + Cryst_ID Energy Dep. position(x,y,z) Seg_ID (00 to 03)

41 The EXOGAM input file Flange ID of the first clover activated. # CONFIGURATION X CONFIGURATION B # # Id Large/Small In/Out T->Ge can Theta Phi # Clover (cm) (deg.) (deg.) # #Flange 2: 0 L #Flange 4: 1 L #Flange 3: 2 L #Flange 5: 3 L #Flange 10: 4 L #Flange 11: 5 L #Flange 12: 6 L #Flange 13: 7 L #Flange 6: 8 L #Flange 7: 9 L #Flange 8: 10 L #Flange 9: 11 L #Flange 15: 12 L #Flange 1: 13 L

42 MINIBALL Also, some special ones Not defined as other ancillaries Close to AGATA design Just different cluster size & segmentation Thus = modified AGATA clusters Use mgt tracking code to analyse MINIBALL data PRISMA Merging of PRISMA code output and ASC output LYCCA grazing code the input file (emitters + g s) Same input file used for both simulation codes Coincidences with acceptance effects + Doppler correction Merging of LYCCA Sim. package + ASC code Contact: Dan Bloor db654@york.ac.uk