GAUSS - GEANT4 based simulat ion f or LHCb GEANT4 Workshop 2 Oct - PowerPoint PPT Presentation

GAUSS - GEANT4 based simulat ion f or LHCb GEANT4 Workshop 2 Oct ober 2002 W. Pokor ski / CERN

Cont ent s I nt r oduct ion t o LHCb sof t war e Overview of Gauss proj ect GiGa – Gaudi int erf ace t o GEANT4 Summary 02/10/2002 Gauss Application 2

GAUDI – LHCb sof t ware f ramework all of t he LHCb event processing sof t ware is built wit hin a f ramework – Gaudi f ramework separat ion bet ween dat a and algorit hms dat a st ore-cent red archit ect ural st yle separat ion bet ween t ransient and persist ent dat a t ransient as a bridge bet ween ot her represent at ions isolat ion of user codes/ algorit hms f rom underlying persist ency t echnologies component s int eract t hrough t heir abst ract int erf aces algorit hms manipulat e dat a conversion services change dat a represent at ions 02/10/2002 Gauss Application 3

GAUDI - Obj ect Diagram Convert er Applicat ion I nt er act ive Convert er Event Convert er Ser vice Manager Select or Message Per sist ency Event Dat a Dat a Ser vice Ser vice Ser vice Files Transient J obOpt ions Event St or e Algorit hm Ser vice Algorit hm Algorit hm Transient Par t icle Pr op. Per sist ency Det ec. Dat a Dat a Det ect or Ser vice Ser vice Ser vice Files St ore Ot her Transient Services Per sist ency Hist ogr am Dat a Hist ogram Ser vice Ser vice Files St ore 02/10/2002 Gauss Application 4

GAUDI - Applicat ions Gaudi Applicat ions Gaudi Applicat ions P luggable Component s Component s P luggable Gauss Simulat ion Gauss Simulat ion Event Gener at ion Event Gener at ion Brunel Brunel Reconst r uct ion Reconst r uct ion Gaudi GiGa Det ect or Simulat ion GiGa Det ect or Simulat ion Gaudi Gauss Gauss DaVinci Analysis DaVinci Analysis Digit izat ion Digit izat ion Panor amix Panor amixVisualizat ion Visualizat ion Embedded Geant 4 Geant 4 Embedded Det ect or Simulat ion GiGa GiGa Det ect or Simulat ion Geant 4 cont rol and Geant 4 cont rol and conf igurat ion, Act ions conf igurat ion, Act ions GiGa GiGa GiGa GiGa Conversion of t ransient t ransient Conversion of obj ect s t o/ f rom Geant 4 Geant 4 obj ect s t o/ f rom GiGa Conversion GiGa Conversion represent at ion represent at ion 02/10/2002 Gauss Application 5

Gauss applicat ion J obOpt s J obOpt s J obOpt s I nt .f ace GiGa Digi MCPar t icle Digit Geant 4 P yt hia Cnv Cnv HepMC Geant 4 MCVer t ex Alg MCDigit (GiGa) et c MCHit Cnv Geomet ry Generat or Det ect or Simulat ion Digit izat ion 02/10/2002 Gauss Application 6

GiGa overview GEANT4 I nt er f ace f or Gaudi Applicat ions or Gaudi I nt er f ace t o GEANT4 Applicat ions GEANT4 callable and cont rollable f rom wit hin GAUDI envir onment common det ect or geomet r y sour ce used by ot her applicat ions (r econst r uct ion, visualisat ion) communicat ion via Tr ansient St or es (Event , Det ect or Dat a) as any ot her ser vice or algor it hm in Gaudi use of common services (Par t iclePr oper t ySvc, RandomNumber Svc, Magnet icFieldSvc, et c.) 02/10/2002 Gauss Application 7

GiGa st ruct ure Per sist ency Dat a Applicat ion Ser vice Files Manager GiGaKine G4 Transient Conversion Event Kine Event Ser vice Ser vice St ore GiGaHit s Conversion G4 Hit s Convert er Cnv GiGa Ser vice Algorit hm Cnv Algorit hm Algorit hm Ser vice Geant 4 Transient GiGaGeom G4 Act ion Det ec. Act ion Det ect or Conversion Geom Ser vice St ore Ser vice Ot her Services Per sist ency Dat a Ser vice Files 02/10/2002 Gauss Application 8

GiGa f eat ures it ’s a Gaudi ser vice pr ovides access t o int er nal G4 event loop via GiGaRunManager all int eract ions wit h Geant 4 only t hrough abst r act int erf aces of GiGa Service minimizes t he couplings t o Geant 4 allows loading ext er nal physics list s inst ant iat es (using Abst r act Fact or y pat t er n) dif f er ent “act ions” (makes t hem t o be plugable component s) 02/10/2002 Gauss Application 9

GiGa - Geomet ry Conversion Xml descr ipt ion Mat er ials Mat er ials Volumes Volumes Xml Cnv Geo Conver sion Ser vice Convert er Convert er Convert er Convert er Gaudi t ransient st ore Geant 4 Geant 4 Geant 4 Geant 4 Geant 4 Geant 4 Mat erials Mat erials Volumes Volumes Mat erials Volumes 02/10/2002 Gauss Application 10

Geomet ry conversion ex. (RI CH 1) Xml � GiGa � G4 � OpenGL 02/10/2002 Gauss Application 11

GiGa - Kinemat ics Conversion HepMC Event LHCb LHCb LHCb LHCb Ver t ices Ver t ices Part icles Part icles Kine Conver sion Ser vice Convert er Convert er Convert er Convert er Convert er Convert er Geant 4 Geant 4 Geant 4 Geant 4 Tr aj ect or ies Tr aj ect or ies Primary Primary Ver t ices Ver t ices 02/10/2002 Gauss Application 12

GiGa – Hit s Conversion LHCb LHCb Ecal Hit s Ecal Hit s LHCb LHCb Muon Hit s Muon Hit s LHCb LHCb Velo Velo Hit s Hit s Convert er Convert er Hit s Conver sion Hit s Conver sion Convert er Convert er Convert er Convert er Ser vice Ser vice Geant 4 Geant 4 Geant 4 Geant 4 Geant 4 Hit s Geant 4 Geant 4 Geant 4 Hit s Hit s Geant 4 Hit s Hit s Hit s Hit s Hit s Hit s 02/10/2002 Gauss Application 13

Det ect or Simulat ion – “physics list s” physics list s: crucial part of t he whole simulat ion program –will cert ainly require several t uning it erat ions most of t he st uf f already implement ed in Geant 4 some specif ic processes needed implement at ion f or RI CH: phot oelect ric process (creat ion of phot oelect rons in HP Ds), energy loss: in t he silicon of HP Ds new f eat ure recent ly added t o GiGa : modular physics list s allows dynamic loading (via j obOpt ions) of part icular physics “sublist s” expect ed t o increase f lexibilit y and t o make validat ion easier 02/10/2002 Gauss Application 14

Sensit ive Det ect ors & Hit s GiGaSensDetTracker G4TrackerHits creat es GiGaTrackerHitsCnv ProcessHit() (Geant4) Geant 4 world conver t s t o: invoked when part icle passed t hrough t he sensit ive volume MCHits lvVolume Gaudi world (/Event/MC/OT/Hits) (XmlDDDB) <logvol name="lvU_ActiveLayer" … sensdet="GiGaSensDetTracker/myDet"> 02/10/2002 Gauss Application 15

RI CH1 wit h SinglePart icleGun RICH1 Event Pion with 7 GeV/c. Cherenkov Photons In Aerogel and C 4 F 10 . Rayleigh scattering Switched off for Illustration. S. Easo 02/10/2002 Gauss Application 16

RI CH1 Hit s 02/10/2002 Gauss Application 17

Panoramix view of MCHit s 02/10/2002 Gauss Application 18

Panoramix view of MCHit s (2) 02/10/2002 Gauss Application 19

Summary Current St at us: we can run GEANT4 simulat ion wit h all t he subdet ect ors included we are able t o produce MCPar t icles, MCHit s (f rom LHCb event model) and save t hem in a ROOT f ile we ar e t est ing/ validat ing subdet ect or specif ic code and physics pr ocesses we have st ar t ed adding digit izat ion algor it hms t o Gauss 02/10/2002 Gauss Application 20

Summary (2) Planning: by end of t he year: t o have a complet e GEANT4 based simulat ion applicat ion by t he summer : t o per f or m some “t est pr oduct ions” by … : t o def init ely move f r om t he old GEANT3 based simulat ion pr ogr am t o t he new one, GEANT4 based 02/10/2002 Gauss Application 21