170 likes | 258 Vues
Certifying Geant4-based calorimeter simulations for the LCD. Dhiman Chakraborty, Guilherme Lima, Jeremy McCormick, Vishnu Zutshi NICADD, NIU ALCWG-Cal Meeting November 24, 2003. Why comparing Mokka and LCDG4?. Previous LCD studies based on Gismo
E N D
Certifying Geant4-basedcalorimeter simulationsfor the LCD Dhiman Chakraborty, Guilherme Lima, Jeremy McCormick, Vishnu Zutshi NICADD, NIU ALCWG-Cal Meeting November 24, 2003
Why comparing Mokka and LCDG4? • Previous LCD studies based on Gismo • Geant4-based LCD simulations are rather new, they need to be certified (then Gismo should be dropped) • (LCD)Mokka and LCDG4 were developed independently. Both are based on Geant4, so they should provide compatible results • XML-based geometry representation, like in Gismo
Fair comparison • Geant4 version 5.2 • SDJan03 geometry (cylindrical layers with virtual cells) • Physics list from Mokka • Range cut of 0.1mm • Text output implemented into projective LCDG4 • Same events are processed in both detector simulators (single particles: 50 GeV e, μ, π) input from binary stdhep file: θ = 90°, flat in φ • Same materials in sub-detectors (look at X0, λI )
Distributions used for comparisons • Energy depositions per layer • Energy depositions per cell • Number of hits per layer • Dependence of Nhits with threshold
Ecal: Energies per layer R A N G E O U T S MIP peaks
Hcal: energies per layer Significant Disagreements?
Ecal: energies per cell Discrepancies! ECal threshold at 0.04 MeV
Hcal: energies per cell Slightly different slopes… HCal threshold at 0.7 MeV
Ecal: energies in absorbers LCDG4 only, Mokka does not provide energies in absorbers MIP peaks
HCal: energies in absorber LCDG4 only, Mokka does not provide energies in absorbers MIP peaks
Cross checking: Ecal+Hcal, cell+abs LCDG4 only
ECal Nhits x threshold – 50 GeV muons ECal threshold
Ecal Nhits x threshold – 50 GeV pions ECal threshold
Ecal Nhits x threshold – 50 GeV positrons ECal threshold
Nhits on Hcal – Dependence on thresholds HCal threshold HCal threshold
Conclusions • Energy deposition: very good agreement on both layer and cell distributions • Number of hits: Good agreement on shapes, small disagreement on normalization • LCDG4 and Mokka give compatible results for the calorimeter simulations of SD detectors • MC particles contributing to hits: expect a bug fix quite soon! • Requests for event processing are welcome, pleasesend a message to jeremy@nicadd.niu.edu.