A Monte Carlo code for AGATA based on Geant4 E.Farnea INFN Sezione di Padova, Italy
� ✁ � � ✁ Why Geant4? The code is well mantained and in widespread use Object-oriented, suitable for big projects C++ based Leaves the user responsibility on geometry definition, detector response and read out, event generation Possibility to import complex geometries from CAD systems?
✂ ✂ ✂ ✂ ✂ ✂ Philosophy of the program Command-line UI based on tcsh Built-in commands to change simulation parameters without recompiling Additional possibilities through switches at the start of the program Sequences of commands automatized through macro files Graphics enabled only when needed Concentrate on the production of list-mode output files rather than making on-line analysis
Class structure of the program Agata *Agata *Agata Agata Agata Agata SteppingAction RunAction EventAction PhysicsList VisManager *Agata *Agata Agata GeneratorEmitter Analysis Agata SteppingOmega *Agata GeneratorOmega GeneratorAction CSpec1D CSpec2D *Agata *Agata *Agata GeneratorGamma GeneratorNeutron Detector *Agata Shell *Agata Detector *Agata DetectorAncillary Construction *Agata DetectorArray CConvex Detector * Possibility to Polyhedron Simple Agata *Agata change parameters DetectorReadOut SensitiveDetector via a messenger class Agata Agata Messenger classes are not shown! DummySD HitDetector
✄ ✄ ✄ AgataDetectorConstruction Generate here only material definition, experimental hall and other passive objects (target, reaction chamber) Handles actual detector arrangement using auxiliary classes Implemented geometries selectable via switch at the start of the main program: single germanium detector ( AgataDetectorSimple ), germanium 4 π shell ( AgataDetectorShell ) and the AGATA array ( AgataDetectorArray )
☎ ✖ ✔ ✒ ✆ ☛ ✝ ✕✖ ✗ ✖ ✘ ✘ ✘✙ ☎ ✖ ✘ ✚ ✖ ✘ ✗✖ ✙ ✖ ✖ ✙ ✛ ☛ ✏ ✆ ✓ ☎✝ ✞ ✡ ✆ ☞ ☎✌ ☎ ✍ ✝ ✞ ☛ ✎ ✏ ☎ ✑ ✒ ✆ ☛ ✝ ✒ ☛ ✖ ✞✠☛ ✞✠✟
✜ ✪ ✾ ✪ ✰ ✰ ✭✧ ✭ ✽ ✮ ✻ ✻ ✧❯ ✮ ✮ ✪ ✲ ✰ ✧ ✳ ✲ ✸ ✲ ✮ ✻ ✪ ❅ ✰ ✧ ✫ ❇ ✬ ✳ ✭ ✬ ❀ ❅ ❆ ❅ ✯ ✴ ✴ ❂ ✧ ❁ ✯ ✵ ✸ ✧ ✦ ✣❙✥ ❅ ❅ ✷ P ✰ ✧ ✳ ❅ ✩ ✪ ✸ ✾ ❅ ❅ ✷ ✷ ❃ ✜ ✧ ✪ ✸ ❍ ✰ ✳ ✷ ❚ ✲ ✧ ✪✵ ✦ ✲ ✲ ✫ ✮ ✪ ✳ ✲ ✱ ✰ ✪ ✭ ✬ ✭ ✻ ✫ ● ✧ ▼ ❂ ✷ ❍ ▲ ❑ ❀ ❈ ❉ ✴ ❂ ❃ ▲ ❑ ❀ ● ✴ ✶ ✪ ◆❖ ✯ ❃ ✷ ❍ ▲ ❑ ❀ ● P ✶ ❂ ▲ ✬ ❑ ❀ ● P ❈ ❏ ✮ ✪ ✳ ✱❏ ✭ ❅ ✪ ✰ ✧ ❋ ✬ ✯ ✧ ✮ ✫ ✭ ✾ ✫ ✸ ✮ ✬ ❊ ✫ ❉ ✷ ❚ ❍ ✪ ✧ ✧ ✧ ✸ ✯ ✸ ❇ ❇ ■ ✲ ✳ ✭ ❅ ❄ ✫ ✧ ✲ ✷ P ✶ ❚ ● ❱ ✳ ✦ ✭ ✬ ✰ ✧ ✫ ❇ ✬ ✳ ✭ ✻ ❅ ✪ ❆ ❅ ✯ ❄ ❃ ❂ ✧ ❁ ❀ ❅ ✭ ✾ ✸ ❋ ✧ ✯ ✧ ✮ ✫ ✭ ✪ ✾ ✮ ✸ ✬ ❊ ✫ ❉ ✷ ❈ ✴ ❅ ✻ ✲ ✪ ✰ ✪ ✲ ✧ ✪✵ ✦ ✲ ✴ ✲ ✮ ✳ ✷ ✲ ✱ ✪ ✭ ✪✫✬ ✩ ✣✤✥ ✜ ✢ ✶ ✳ ✰ ✲ ✰ ✧ ✭ ✭ ✽ ✲ ✧ ✪✵ ✦ ✴ ✰ ✲ ✰ ✧ ✳ ✲ ✯ ✮ ✵ ✸ ✧ ✬ ✪ ✮ ❃ ▲ ❑ ❀ ● P ❈ ❏ ✪ ✶ ✳ ✬ ◆❖ ❉ ✪ ✫ ▼ ❂ ❂ ❈ ❍ ✲ ❘ ❏ ✧ ✲ ✲ ✮ ✪ ✳ ✱ ● ✯✰ ✧ ❁ ❂ ✷ ❍ ▲ ❑ ❀ ✷ ✧ ▲ ❅ ✯ ✸ ❇ ❇ ■ ✲ ✳ ✭ ✧ ✫ ✭ ✧ ✲ ✷ ✴ ❍ ✶ ● ✴ ❄ ✧ ✧ ✸ ✯ ❀ ❑ ❀ ● ✷ ✴ ❂ ❃ ❑ ▲ ● ❍ ❄ ✱❏ ❂ γ γ γ γ ✮★✯ ✮✼✻ ✧✺✹ ✸✿◗ γ γ γ γ ✮★✯✰ ✮✼✻ ✸✿✻ ✧✺✹ ✦★✧
❳ ❳ ❳ ❲ ❳ AgataDetectorArray Irregular polyhedra generated with the CConvexPolyhedron class (D.Bazzacco) Actual detector shape can consider the intersection of such polyhedra with a closed-end cylinder Vertexes of the polyhedra calculated with an external program ( MarsView by D.Bazzacco) Available data files for the geometries with 180 crystals and 120 crystals (grouped in triple or quadruple clusters, or with the shape used in GRETA) Possibility to add extra passive materials to emulate an ancillary device ( AgataDetectorAncillary )
❨ ❩ ❩ ❩ ❨ ❨ ❩ ❨ ❩ ❩ ❨ ❨ Start with a On its faces, draw a regular platonic solid pattern of triangles grouped Project the faces on e.g. an icosahedron as hexagons and pentagons. the enclosing sphere; E.g. with 110 hexagons and flatten the hexagons. (always) 12 pentagons
❬ ❬ ❬ ❭ ❭ ❬ ❭ ❬ ❬ ❭ ❭ ❭ Al capsules 0.7 mm spacing 0.8 mm thick Al canning 2 mm spacing 2 mm thick A radial projection of the spherical tiling generates Space for encapsulation and the shapes of the detectors. canning obtained cutting the Ball with 180 hexagons . Add encapsulation and crystals. In the example 3 part of the cryostats for crystals form a triple cluster realistic MC simulations
Configuration A=180
Configuration A=180 – solids Solid 1 Solid 2 Solid 3
Configuration A=120
Configuration A=120 – solids Solid 1 Solid 2
❪ ② ❡ ❡ ❢ ➇ t ♦ s ✈ ♣ q ❸ s ✉ ✉ s ➈ ④ ❸ ❸ ♦ t ✉ ⑨ ⑥ ❤❢ ❡ t ❤ t ➂ ➃ ⑧ ➄ ➅ q t ② ➆ ❷ ❸ ③ ♦ ④ ➁ ➂ ➃ ❤ ⑩ ❢ ♠ q ♣ ✇ ✉ ✈ ⑥ ⑨ ⑤ ❥ ♦ ✉ ① ④ s ♦ ➁ s ♠ ➃ ✐ ❣ ♣ ⑨ ⑩ ❺ ❻ ❣ ❼❽ ④ ✈ ❤ ➁ ➂ ➃ ♦ r ⑥ q ❤ ⑨ ⑥ ❤ ❢ ⑥ ⑨ ➉ ② ❸ ③ ✉ t ♦♣ ❧ ✉ ✈ ♦ ♣ ✉ ❥ ❤ ➊ ➋ ❸ ❸ ❢ ➁ ✈ ❧ ① ♦ ✉ ✐ ✐ ✐ ❢ ② ♣ ♣ ♦③ ♣ ✈ ⑧ ⑥ ❤❢ ⑥ ❤ ❢ ⑥ ⑧ ⑥ ⑤ ❡ ❥ ✇ ❣ ❣ ❣ ❛ ❜ ❝ ❛ ❫ ❞ ❪ ❜ ❫ ❛ ❴ ❡❢ ❤ ✉ ✐ ❢ ❣ ❥ ♠ q t ❧ ✉ ✈ ♦ ♣ ❢ ❤ ✐⑨ ② s q ⑤⑩ ♦ ♣ ④ ③ ♦ q r ♣ ❸ ④ ③ ✉ ❥ ⑨ ❣ ✐ ⑩ ❺ ❻ ❣ ❼❽ ④ ♣ ❶ t ❤ ✇ ❧ ♣ ❶ ❢ ♦ ❥ ③ ❧ ♦ ❥ t ⑧ ④ ⑩ ♣ ⑨ ❣ ③ ♣ Efficiency reduces also but all nice symmetries remain; smaller crystals simplify PSA . To reduce cost of germanium, A-180 could be squeezed to similar size as A-120. t★④ ⑤⑦⑥ ❫❵❴ ❾➀❿ ❾➀❿ r✼s r✼s r✼s ♥★♦♣ ♥★♦♣ ❦✿❧ ❦✿❧ ❷❹❸ ε ε
Recommend
More recommend
