A Novel Fast Track Engine for the ATLAS Experiment on behalf of the ATLAS Collaboration Physikalisches Institut University of Bonn February 26th 2010 ACAT 2010 Jaipur, India
1 The ATLAS detector Fast detector simulations Why? 2 Track simulation in Digitisation Combined simulation: Tracking and calorimetry Comparison to other simulation 3 Comparison with full simulation Comparison with collision data Special use-cases and applications 4 22
The ATLAS detector Muon Spectrometer Hadronic Calorimeter Electromagnetic Calorimeter Solenoid 3 Tracking System TRT Pixel/SCT
The ATLAS Inner Detector TRT drift tubes 300 000 straw tubes 130 µm resolution (Rφ) Xe/CO 2 /O 2 about 30 measurements per track SCT silicon strips 6 million Si strips resolution: 17 µm (Rφ)/580 µm (Z) 4 (double) measurements / track silicon Pixel 80 million Si pixels resolution: 10 µm (Rφ)/115 µm (Z) 3 measurements / track 2 T solenoidal field 44
Why do we need fast detector simulations? Monte Carlo simulation of detector response is needed to compare theoretical predictions (by Monte Carlo event generators) to data Detailed simulation of particles penetrating the detector material is CPU-time consuming of a single t t event in full Geant4 simulation takes about 30 ksi2kminutes Fast simulation techniques can increase the amount of simulated events 55
Track simulation 4-momentum, particle type (PDG ID) si standard high energy physics libraries (e.g. PYTHIA) Geant4 Detector, Full physics list Digitisation Reconstruction c Track Material effects Particle decay Photon conversions Digitisation ri ATLFAST Track represen h smearing Track detector simulation simulates passage of particles through the detector, charge deposition in sensitive detectors, showers, material interactions digitisation emulates readout reconstruction pattern recognition, track finding, track fitting 66
The simulation scheme of uses the component model of the ATLAS software take extrapolation engine of the ATLAS track reconstruction reconstruction modules, such as the estimation of energy loss, are replaced by Monte Carlo implementations event data objects identical to full simulation or real data Where feasible, Geant 4 modules are used, such as particle decays Nearly all effects are estimated from first principles (Bethe-Bloch, etc.) no parametrisations used, despite hadronic interactions 77
ATLAS Inner Detector Geometry Geant 4 r [mm] 0.0 0.25 0.5 0.75 1.0 1.5 2.0 Tracking Geometry r [mm] 0.0 0.25 0.5 0.75 1.0 1.5 2.0 1000 1000 800 600 400 2.5 3.0 800 600 400 2.5 3.0 200 Geant4 0-4000 -3000-2000 -1000 0 1000 2000 3000 4000 z [mm] 200 TrackingGeometry 0-4000 -3000-2000 -1000 0 1000 2000 3000 4000 z [mm] Track reconstruction in ATLAS uses greatly simplified detector description Sensitive elements identical to full Geant 4 description! uses the same! 88
Track simulation in MC Event Generator 99
Track simulation in MC Event Generator µ 99
Track simulation in MC Event Generator K 0 S µ 99
Track simulation in MC Event Generator e K 0 S µ 99
Track simulation in MC Event Generator π + e K 0 S µ 99
Track simulation in MC Event Generator γ π + e K 0 S µ 99
Track simulation in MC Event Generator e γ π + e K 0 S µ 99
Track simulation in MC Event Generator γ π + e KS 0 µ e Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ π + e KS 0 µ Extrapolation to next surface Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ γ secondaries, brem photons, conversion electrons Material Effects Extrapolation to next surface π + e K 0 S µ Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ γ secondaries, brem photons, conversion electrons Material Effects Extrapolation to next surface Decay point reached? π + e K 0 S µ Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ γ secondaries, brem photons, conversion electrons Material Effects Extrapolation to next surface no Decay point reached? π + e K 0 S µ Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ γ secondaries, brem photons, conversion electrons Material Effects Extrapolation to next surface no Decay point reached? π + e K 0 S µ Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ γ secondaries, brem photons, conversion electrons Material Effects Extrapolation to next surface no Decay point reached? π + e K 0 S µ Particle unstable? Simulate decay length 99
Track simulation in MC Event Generator γ γ π + secondaries, brem photons, conversion electrons Material Effects Extrapolation to next surface no Decay point reached? yes Create decay products e K 0 S µ Particle unstable? Simulate decay length 99
Track simulation in Postprocessing 10 10
Track simulation in Postprocessing Extract measurements from simulated tracks 10 10
Track simulation in Postprocessing Extract measurements from simulated tracks 10 10
Track simulation in Postprocessing Extract measurements from simulated tracks 10 10 Add noise
Track simulation in Postprocessing Extract measurements from simulated tracks 10 10 Add noise and merge clusters
Digitisation in Silicon clusterisation with geometric approach d track intersection with surface exit of sensor material ( ) pixel position (vetoed) cluster position A B C p y p x A B C Without Lorentz angle: Q 1 Q 2 Q 3 ixeldetector With Lorentz angle: h of the track Q1 Trackel is o the charge Q2 Track 11 11 θl u θl Lorentz Angle ring tion is random landau mbers shift to θl
Combined and FastCaloSim 0 Y (cm) 100 Calo cells and FastCaloSim interfaced: Conversion point Electron tracks LAr Calorimeter Presampler Brem point Solenoid coil Elect simulates the Inner Detector including secondaries Calorimeter deposits simulated by FastCaloSim Muon System simulated by 0 X (cm) 12 12
The Fast and the Furious Atlfast I Atlfast II Geant4/full Atlfast IIF Atlfast II ID track per FatrasID digitisatio reconstruc digitisation reconstructio Calo cluster FastCaloS muons: fu digitisation reconstructio MS track perig FatrasM digitisatio reconstruc digitisation reconstructio rel. gain - 13 13
Material effects Comparison of and Geant4 2 10 Multiple scattering Energy loss 10 Geant4 2 GeV Single µ - Gaussian traversing 1mm (Highland) Si 1 (~1.06 Gaussian % X Mixture 0) (Frühwirth et al.) -0.004-0.002 0 0.002 Radiation of brem photons Entries/bin 5 10 4 10 Hadronic interactions Entries 307120 Geant4 Mean 4 111 Entries 195013 Mean 5 874 Entries/bin 3 10 2 10 14 14 10 1 0 20 40 60 80 p [GeV] 100 Energy 120of secondary particles
Digitisation in Tuning of the silicon clusterisation Pixels and SCT clusterisation tuned by adapting minimal required path length in the Pixels cell strength of Landau smearing Tuning can be done within 24h 15 15 Cluster width in Pixels Measurement residual in Pixels
Track parameter resolutions Single muon events with p T =1 GeV, 5 GeV, 100 GeV In general good agreement, but still some parameters to tune in the digitisation In particular tails better described than in ultra-fast sim Entries/bin offline/newt 3 ATLFAST 10 10 2 1 GeV 5 GeV 100 GeV Geant 4 1 GeV 5 GeV 100 GeV 10 1-1 -0.5 0 0.5rec 1 d true -d (mm) 0 0 rec true 16 16 d 0 - d 0 [mm] d 0 -d 0 (mm) Transverse impact parameter d 0 Inverse track momentum q/p T
Reconstructed tracks in minimum bias events at s = 900 GeV on minimum bias Monte Carlo (Pythia) at the center of mass energy s = 900 GeV Still some discrepancies, but not yet tuned to data and misalignments of detector modules in data 17 17
Reconstructed tracks in minimum bias events Average number of Pixels hits per track 18 18 Pixels hits / track vs η Average Number of Pi 4 3.5 3 2.5 ATLAS Preliminary Fatras 2-2.5-2 -1.5-1 -0.5 0 0.5 1 1.5 2 sinoidal structure due to inactive Pixels modules in b-layer folded with the z-position of the primary vertex Pixels hits / track vs φ Average Number of Pixe 3.5 3 2.5 4 Fatras ATLAS Preliminary 2-3 -2-1 0 Detector conditions like inactive modules automatically taken into account in Precise description of the detector geometry
Reconstructed tracks in minimum bias events at s = 900GeV Pixel cluster size mean cluster size (Δ η 1.8 1.6 1.4 1.2 1 0.8 0.6 0.4 0.2 0 Q 1 Fatras Without Lorentz angle: Q 2 Q 3 u ATLAS Preliminary θ L -2-1 0 1 2 Track incident angl Cluster size depends on incident angle, because of sensor thickness With Lorentz angle: ing Q 1 Q 2 Track ion is random landau mbers θ L shift to Track θ L Lorentz Angle Mean cluster size (in η direction) vs incident angle (η) on the Pixel module Very sensitive test of the clusterisation model 19 19
Reconstructed tracks in minimum bias events at s = 900GeV Transverse impact parameter wrt. primary vertex Arbitrary -1 10 Fatras -2 10-3 10 ATLAS Preliminary -2-1.5-1 -0.5 0 0.5 1 1.5 d 0 [m Longitudinal impact parameter Arbitrary 10 Fatras -2 10-3 10 ATLAS Preliminary -2-1.5-1 -0.5 0 0 Position and size of the beam spot in the simulation taken from detector conditions data base 20 20
for Super-LHC upgrade studies Occupancy 0.3 0.25 10 Pix Pix Pix Pix 0.2 was used for SLHC upgrade studies of the ATLAS tracker Allows easy testing of various geometries at a reasonable time scale 0.15 0.1 0.05 0 0 50 21 21 Detector occupancies can be derived reliably Reconstruction effects are included in momentum resolutions, etc.
Stress tests for track fitters simulation with high noise levels Detailed truth information by allows to evaluate performance of track fitters Quick simulation of arbitrary noise levels ) [MeV] T σ(p 85 80 75 70 65 60 55 Kalman DAF Silicon-only 0.0 0.1 0.2 0.3 0.4 0.5 0.6 TRT noise occupany 22 22 5 GeV single muon events Example: Study of adaptive track fitter (Deterministic Annealing Filter) High detector occupancy Solution of left-right ambiguities in the TRT
Conclusions is a new track simulation concept between full Geant4 simulation and conventional fast detector simulations The full reconstruction chain can be run on output Speed improvement mostly due to simplified Tracking Geometry and extrapolation (Nearly) no parametrisations needed All important physics effects included, like multiple scattering, brem, conversions, particle decays, hadronic interactions Allows studies to be performed that cannot easily be done either with full simulation or conventional fast simulations Currently in the tuning phase Validation with collision data has started 23 23
Acknowledgements Special thanks to Andreas Salzburger (CERN), Sharka Todorova (Tufts University) and Simone Zimmermann (Bonn) and PHYSICS AT THE TERA SCALE Helmholtz Alliance 24 24
25 25
The ATLAS Muon System 26 26 Monitored Drift Tubes (MDT) 354k straw tubes barrel and forward region 80 µm straw resolution (Z) 20 measurements / track Resistive Plate Chambers (RPC) barrel region chamber resolution: 10 mm (Z)/10 mm (φ) 6 measurements / track trigger (+ 2 nd coordinate) Thin Gap Chambers (TGC) end-cap chamber resolution: 2 6 mm (R)/3 7 mm (φ) 9 measurements / track trigger (+ 2 nd coordinate) barrel toroid: 1.5 5.5 Tm bending power, end-cap toroids: 1 7.5 Tm Cathode-Strip Chambers (CSC) forward region multi-wire prop. chambers plane resolution: 60µm (R)/5mm (φ) 4 measurements / track Muons with momenta of 4 GeV and 20 GeV in the bending plane of the barrel muon spectrometer.
The green simulation time/event, ksi2kseconds Minimum Bias t t Jets W ± e ± ν e Heavy Ion Full Sim 551 1990 2640 1150 56, 000 Fast G4 Sim 246 757 832 447 21, 700 ATLFAST-II 31.2 101 93.6 57.0 3050 ATLFAST-IIF 2.13 7.41 7.68 4.09 203 ATLFAST-I 0.029 0.097 0.084 0.050 6 Hard to estimate how much CO 2 is emitted for one CPU second by the grid Taking 0.01 grams/s: Simulating 100k t t events with ATLFAST-IIF ( + FastCaloSim) instead of Full saves about 2 tons of CO 2 27 27
Digitisation in of hits with transition radiation in the Transition Radiation Tracker transition radiation in the TRT produces hits with stronger signal and is used for particle ID probability of transition radiation depends on relativistic γ factor measured with test beams and cosmic ray data fit of turn-on curve has been fed into 28 28 γ
Combined simulation Comparison of reconstructed muon p T resolution T central muons ( η < 1.2) in Z µ + µ events Muon System stand-alone Resolution = σ Muon p 0.08 0.06 0.04 A F F Combined Muon Spectrometer / Inner Detector reconstruction 0.02 0 0 10 20 30 40 50 60 7 29 29
Reconstruction efficiencies Single electrons and muons with transverse momentum p T = 5 GeV Shape roughly reproduced, but electrons Needs some extra fudge factors too perfect for Efficiency 1. 0. 0 0. 0 0. 0 30 30