Detector simulation withMokka/Geant4

Gabriel Musat L.L.R. – Ecole polytechniqueECFA-DESY Linear Collider Workshop

Amsterdam, April 1-4, 2003

What is Mokka? Developped at L.L.R. since the end of 1999.

Became « de facto » the reference for the calorimetry studies shown in the T.D.R.

Great news : became the Geant4 detailed simulation tool for all detector modules for the F.L.C. AND calorimeter prototypes.

Gabriel Musat L.L.R. – Ecole polytechnique

Anna Lupi, Andrea Dotti – INFN Pisa

●G4 better simulates pions●G3 and G4 do not agree●Impact position dependencymuch less evident: hadronicshowers much wider thanTileCal sampling period

Pions 20 GeV h=0.25 (run 200599)

Resolution

Data = (17.7±0.1) %G4 = (17.1±0.2)%G3 = (14.0±0.1)%

2.83 ±0.042.42±0.022.79 ±0.01s

15.80 ±0.0317.30±0.0215.71±0.01Mean

G4G3DATA

Bug fix in Geant4-05

• PI Reconstructed Energy less than 0.8 * Real Energy:

–Before : 10%

–After: 14.7%

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka Mokka kernel - kernel - Common GeometryCommon GeometryAccess API (F77, C++,C, Java)Access API (F77, C++,C, Java)

Simulation(Mokka)

Reconstruction(user code)

Analysis(user code)

Geometry drivers

Run Time

CGA APIGeometryDatabase

Implements some reconstruction utilities.

Gabriel Musat L.L.R. – Ecole polytechnique

The CGA API and the reconstruction

Example, a 50 GeV pi in ECAL : the total number of X0 from the entrance point in ECAL is required

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka kernel - Common Geometry Access API

• New java/CGA sub-directory– CGARunManager class implemented using JNI– CGATest - example

• New examples/CGA sub-directory– F77 and C examples

Gabriel Musat L.L.R. – Ecole polytechnique

MokkaMokka kernel - extending the kernel - extending thepersistency modelpersistency model

Mokka

ASCIIfiles

Persistency Abstraction subsystem

LCIO

SIOfiles

Beta version inmokka-01-05

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka detectors - new mask

Karsten Buesser, January 2003

• New sub-detector mask03 and new detector

model D09M1

• The mask doesn’t collect hits

• LAT module implemented with a material

density equal to half tungsten density

• LCAL module not implemented

Mokka detectors - new mask

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka detectors - New model “ProtoEcalHcalRPC”

Ecal prototype:old prototype already available in Mokka

Hcal prototype :a Hcal module in a cube format (~1 m3), built with RPCs

Y

XBeam

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka detectors - in development (not yet integrated)

• Detector model DT01 reflecting the BRAHMS detector

(Ties Behnke)

• NLC (XML) detector model (Ties Behnke)

• CMOS model of VxD (Damien Grandjean)

• SET implementation (LPNHE Paris)

• Hcal GEM (Venkat, University of Texas at Arlington)

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka kernel- future developments I

• Modify the registration scheme for new detector drivers

– implement the Control class as a singleton

– add a static member to the driver class

• Put reference material definitions in the DB

• Put specific materials in the sub-detector driver – add a virtual method

– put material definitions in the sub-detector DB

• Improve the Event Display

Gabriel Musat L.L.R. – Ecole polytechnique

• Put MySQL scripts in the DB

– a) identify sub-detector master parameters

– b) compute size and placements of each detector piece as function of a)

– c) translate b) into G4 objects

• Eventually put XML files in the DB

Mokka kernel- future developments II

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka kernel- future developments III

• Add “cellIndex’’ method to CGA

– (P S M I J K) => X Y Z

– P, S, M = detector piece, stave, module numbers

– I,J = the cell coordinates in the cells matrix

– K = Sensitive (Si or scintillator or gas) layer number

– X,Y,Z = the cell centre in world coordinates

Gabriel Musat L.L.R. – Ecole polytechnique

Mokka - collaboration support

• Detector development integration– DB access in Write mode via accounts

– Install Mokka CVS repository on pollin1.in2p3.fr• global access

– sub-directories of Mokka

• Data server for simulated events on pollinclb.in2p3.fr

• Software documentation

Gabriel Musat L.L.R. – Ecole polytechnique