Slic A Geant4-based detector simulation package Norman Graf, Jeremy McCormick SLAC October 15, 2009.

Presentation on theme: "Slic A Geant4-based detector simulation package Norman Graf, Jeremy McCormick SLAC October 15, 2009."— Presentation transcript:

1 slic A Geant4-based detector simulation package Norman Graf, Jeremy McCormick SLAC October 15, 2009

2 2 Mission Statement Provide full detector response simulation capabilities for Linear Collider physics program. Need flexibility for new detector geometries and readout technologies. The system should be flexible, powerful, yet simple to install and maintain. Limited resources demand efficient solutions, focused effort.

3 3 Full Detector Response Simulation Use Geant4 toolkit to describe interaction of particles with matter. Thin layer of non-G4 C++ provides access to: –Event Generator input –Detector Geometry description input –Detector Hits output Geometries fully described at run-time! –In principle, as fully detailed as desired. –Uses lcdd, an extension of GDML. Solution is applicable beyond LC problem domain.

4 Geometry Definition Goal was to free the end user from having to write C++ code to define the detector. All of the detector properties should be definable at runtime with an easy to use format. Selected xml, and extended the existing GDML format for pure geometry description. 4

5 5 Why XML? Simplicity: Rigid set of rules Extensibility: easily add custom features, data types Interoperability: OS, languages, applications Self-describing data, validate against schema Hierarchical structure  OOP, detector/subdetector Open W3 standard, lingua franca for B2B Many tools for validating, parsing, translating Automatic code-generation for data-binding Plain text: easily edited, cvs versioning

6 LCDD and GDML GDML LCDD expressions (CLHEP) materials solids volume definitions geometry hierarchy detector info identifiers sensitive detectors regions physics limits & cuts visualization magnetic fields Adopted GDML as base geometry definition, then extended it to incorporate missing detector elements.

7 LCDD Structure LCDD Root Element Information about the Detector Identifier Specifications Detector Readouts Physics Limits Regions (sets of volumes) Visualization Attributes GDML Root Element Constants, Positions, Rotations Material Definitions Solid Definitions Volume Hierarchy Magnetic Field

8 lcdd Features Regions: production cuts Physics limits: track length, step length, etc. Visualization: color, level of detail, wireframe/solid Sensitive detectors –calorimeter, optical calorimeter, tracker –segmentation IDs –volume identifiers (physvolid) Magnetic fields –dipole, solenoid, field map utilities –information on Geant4 stores –GDML load/dump 8

9 “Compact” Description The lcdd file is very descriptive, but therefore also very verbose. Can be written by hand, but prone to human error. –Also, just specific to the simulation and not easily accessible to reconstruction and visualization. Developed a “compact” detector description which encapsulates the basic properties of a detector and which is further processed by code to produce the input specific to different clients. 9

10 Compact Detector Description A number of generally useful detector types (at least for HEP collider detectors) have been developed, such as: – sampling calorimeters –TPCs –Silicon trackers (microstrip as well as pixel) –Generic geometrical support structures Can also incorporate GDML snippets –Allows inclusion of more complicated volumes derived for instance from engineering (CAD) drawings. 10

11 11 GeomConverter Compact Description GeomConverter LCDD GODL org.lcsim Analysis & Reconstruction HEPREP Small Java program for converting from compact description to a variety of other formats slic lelaps wired lcio

12 Compact Description - Example <detector id=“3” name=“HADBarrel” type=“CylindricalBarrelCalorimeter” readout=“HcalBarrHits” vis=“HADVis”> global unique name global unique identifier detector type readout collection layering sensitive layer absorber visualization settings

13 13 xml: Defining a Tracker Module <module_component width="7.6" length="125.0" thickness="0.26" material="CarbonFiber" sensitive="false"> <module_component width="7.6" length="125.0" thickness="0.05" material="Epoxy" sensitive="false"> <module_component width="9.6" length="125.0" thickness="0.1" material="Silicon" sensitive="true">

14 xml: Placing the modules

15 15 A Barrel Vertex Detector

16 Example Vertex Detector CAD Drawing GEANT Volumes LCIO Hits

17 17 Barrel Outer Tracker Composed of overlapping silicon wafers.

18 Complete Silicon Tracker 18

19 19 Generic Hits Problem Statement We wish to define a generic output hit format for full simulations of the response of detector elements to physics events. Want to preserve the “true” Monte Carlo track information for later comparisons. Want to defer digitization as much as possible to allow various resolutions, readout technologies, etc. to be efficiently studied.

20 20 Types of Hits “Tracker” Hits –Position sensitive. –Particle unperturbed by measurement. –Save “ideal” hit information. “Calorimeter” Hits –Energy sensitive. –Enormous number of particles in shower precludes saving of each “ideal” hit. –Quantization necessary at simulation level.

21 21 Tracker Hit MC Track handle Encoded detector ID (detector dependent ) Hit position in sensitive volume Track momentum at hit position. Energy deposited in sensitive volume. Time of track's crossing. Path length in sensitive volume.

22 22 Calorimeter Hit Encoded detector ID (detector dependent) MC IDs for tracks contributing to this cell. Energy deposited. Time of energy deposition. Repeated for each energy contribution. Support recently added for optical calorimeters –Can store Cerenkov and scintillation light. –See other talks at this workshop.

23 23 LCIO Persistency framework for LC simulations. Currently uses SIO: Simple Input Output –on the fly data compression –some OO capabilities, e.g. pointers –C++, Java (and FORTRAN!) implementations available Changes in IO engine designed for (e.g. root). Extensible event data model –Generic Tracker and Calorimeter Hits. –Monte Carlo particle heirarchy.

24 24 slic: The Executable Build static executable on Linux, Windows, Mac. –SimDist build kit uses GNU autoconf Commandline or G4 macro control. Only dependence is local detector description file. –Trivial Grid usage (no database call-backs, etc.) –Grid ready, Condor and lsf scripts available. Event input via stdhep, particle gun, … Detector input via GDML, lcdd Response output via LCIO using generic hits.

25 25 Detector Full Simulation MC Event (stdhep) Geometry (lcdd) Raw Event (slcio) GEANT4 slic Compact Geometry Description (xml) Reconstruction, Visualization, …

26 26 Detector Variants Runtime XML format allows variations in detector geometries to be easily set up and studied: –Sampling calorimeters: absorber materials, dimensions Readout technologies, e.g. RPC, scintillator Layering (radii, number, composition) Readout segmentation (size, projective vs. nonprojective) –Total absorption crystal calorimeters Optical properties –Tracking detector technologies & topologies TPC, silicon microstrip, silicon pixels

27 ILC Full Detector Concepts SiD GLD LDC

28 28 slic & lcdd: Summary Provides a complete and flexible detector simulation package capable of simulating arbitrarily complex detectors with runtime detector description. Being used by ILC detector community for simultaneous and iterative evolution of different detector concepts and their variations. Being used for ATLAS upgrade tracking studies Has been applied to CPT simulations. Could be used by other communities (astro, medical) for rapid prototyping or simulation.

29 Additional Information ILC Detector Simulation http://www.lcsim.orghttp://www.lcsim.org ILC Forum http://forum.linearcollider.orghttp://forum.linearcollider.org SLIC http://www.lcsim.org/software/slichttp://www.lcsim.org/software/slic LCDD http://www.lcsim.org/software/lcddhttp://www.lcsim.org/software/lcdd GDML http://gdml.web.cern.ch/GDML/http://gdml.web.cern.ch/GDML/ LCIO http://lcio.desy.dehttp://lcio.desy.de