ILD Detector Optimization and Benchmarking at Tsinghua University 12-January-2009
ILD Introduction ILD origins in the European and Asian based Large Detector study. ILC Reference Design (RDR) in 2007 GLD Detector Outline Document (DOD) arxiv:physics/0607154 LDC DOC http://www.ilcldc.org/ Common feature: Tracker(Pixel & Silicon & Gas) + PFA calorimeter + LDC GLD At LCWS2007, we agreed to work together for a joint LOI GLD (B=3T, R ECAL =2.1m) + LDC(B=4T, R ECAL =1.6m ) ILD 2
International Large Detector ( ILD ) LDC and GLD had a common future; Pixel vertex detector placed very close to the beam pipe. Gaseous tracker, TPC, for highly efficient and precise track measurements, supplemented by silicon trackers. EM and HD calorimeters are placed inside a solenoid field and read our by very small sensors to achieve a good energy measurement by Particle Flow Analysis (PFA). But in detail: B=4 Tesla(LDC) vs 3 Tesla(GLD) ECAL radius: 1.6m(LDC) vs 2.1m(GLD) Sub Detector technologies. Simulation studies of physics performances are used to reach agreement of detector parameters. 3
How we optimize Optimization tools GLD Jupiter/Sattelites, LDC Mokka/MarlinReco intermediate detector models were introduced for comparison GLDPrim by Jupiter, and LDCPrim by Mokka, both having B=3.5T and R ECAL =1.85m. Performances have been studies as a function of major parameters. Reached a consensus on the ILD reference detector for LOI benchmark studies at Cambridge (Sep. 2008). Physics performance studies have been performed based on ILD model. 4
Jupiter/Satellites for Full Simulation Studies : GLD Tools for simulation Tools JUPITER JLC Unified Particle Interaction and Tracking EmulatoR Geant4 based Simulator MC truth generator Satellites IO LEDA Input/Output module set Monte-Calro Exact hits To Intermediate Simulated output Library Extention for Data Analysis METIS URANUS JSF/ROOT based Framework Event Reconstruction For real data Unified Reconstruction and ANalysis Utility Set JSF: the analysis flow controller based on ROOT The release includes event generators, Quick Simulator, and simple event display 5
Mokka Mokka is a full simulation using Geant4 and a realistic description of a detector for the future linear collider. LDC Home page: http://polzope.in2p3.fr:8081/mokka Mokka is now a part of the ilcsoft, http://ilcsoft.desy.de/portal/software_packages/ Detector Geometry: managed by MySQL data base and CGA (Common Geometry Access) API. LDC and other variants are prepared and used for ILD optimization. Implementation of detailed detector model based on engineering studies is in progress. ex. ECAL structure 6
PandoraPFA LCFIVertex
GLD + LDC Combined Framework Whizard Physsim LDC MOKKA StdHep Jupiter GLD StdHep: Same generator data LCIO: Common IO format GLDPrim/LDCPrim: Similar detector model Marlin LCIO Sattelites LCIO helps to collaborative works for detector optimization LCIO DST and Analysis After the LOI, two frameworks will be merged to a single framework. 8
Detector Parameters for Opt. studies GLD/GLDPrim/J4LDC prepared for Jupiter LDC/LDCPrim/LDCGLD prepared for Mokka Physics performance was compared between different geometries Jupiter Mokka GLD GLDPrim J4LDC LDCGLD LDCPrim LDC B(T) 3.0 3.5 4.0 3.0 3.5 4.0 VTX Rmin (cm) 1.75 1.6 1.5 1.65 1.50 1.4 # VTX layers 3 x double super layers 5 layers # IT layers 4 layers 2 layers TPC Rmin(cm) 43.7 43.5 34.0 37.1 ECAL Rmin(cm) 210 185 160 202 182.5 161 HCAL Thick. (Int.L) 6.79 6.29 5.67 5.86 Geometryiesin Mokka and Jupiter are similar, but there are many small differences in geometry and assumed detector technologies 9
Pt resolution Single muon, produced at cosq=0. by Jupiter+Satellites: TPC+IT+VTX fitting LDC : ~5% worse at high Pt Shorter Lever arm GLD/GLD : ~10%worse at low Pt Lower B 10
GLDPrim - LDCPrim 4mm LDCPrim(Mokka+Pandora) is better than GLDPrim(Jupiter+Sattelites) by 15~30%. Possible source: s rf (IT) 4mm(LDCPrim) 10mm(GLDPrim) Silicon External Tracker in Mokka 3x10-5 Sub-detector technology is more important than geometry 11
GLDPrim vs LDCPrim (s rf (IP)) GLDPrim is better than LDCPrim ; 3 double layers vs 5 layers? s rf =s Z =2.8mm Fast sim. study by M.Berggren 12
kaon_0l Energy Resolution Hadron Model: LCPhysics HCAL response is not smooth around 13 GeV - LE/HE behaviour ECAL resolution: same 13
Jet measurement: Particle Flow Analysis PFA: Charged particles by Tracker Neutral particles by Calorimeter, remove charged particle energies Performance studies depend on shower simulation; longitudinal, lateral, and tof distribution, neutron response, etc. 14
Jupiter data analyzed by PandoraPFA Pandora PFA: Sophisticated algorithm tuned to Geant4 shower shape has achieved the performance goal of ILC, DE/E ~ 30%/ E Z pole uds-pair events:gldprim Ejet(GeV) 15
Jet Energy Resolution Same trend is seen by analysis of Jupiter models, though performance is slightly worse than Mokka model 16
ECAL Seg. and HCAL thickness by M. THomson ILD: ECAL+HCAL= 6.8 Int. L.(48layers) 6.8 Int. L look OK, but worse resolution is seen for 90 o jets. Performance is strong function of ECAL seg. size. 2x2cm 2 too large, 1x1cm 2 would be ok for jets with E < 100 GeV 17
Physics Benchmark Studies ILC goal precise studies of Tera scale physics. According to the request by ILC Research Director (RD) and International Detector Advisory Group(IDAG), simulation studies for LOI should based on a realistic Monte Calro program based on a realistic reconstruction program include backgrounds by physics processes and those caused by accelerator. Signal processes: the minimum set. Recoil mass measurement by e + e - ZH e + e - /m + m - + H H c cbar decay in e + e - ZH process e + e - t tbar 6 jets and t (tbar) charge ID for A FB meas. e + e - t + t- and t pol. measurement. Separate WW and ZZ in Chargino/Neutrino pair production process 18
Higgs recoil mass meas. e + e - ZH e + e - X / m + m - X, Ecm=250 GeV, 250 fb-1 Analysis. by Itoh Kazutoshi Select e + e - / m + m - consistent with Z and study recoil mass Precise track meas. is a key for Compare 3 geometries GLDprim case, with backgound m + m - X e + e - Channel m + m - channel Differences are small. 19
Benchmark study: Example Using several detector models, performance to separate W/Z in jet mode have been studied using SUSY processes e e W W W W + - + - 0 0 + - 1 1 e e + - by Taikan Suehara ZZ 0 0 0 0 2 2 1 1 20
e + e - t + t -, t rn t is polarized probe New Phyaics by Taikan Suehara SM+NP SM 21
ILD reference detector model At 2 nd ILD WS at Cambridge, we agreed to created the new model, ILD reference design model for LOI, in Mokka: Model parameters, B=3.5 Tesla Rin ECAL=185cm, TPC: halfz=230cm VTX three double layers. Silicon trackers: ( SIT, FTD, SET, SOT) Calorimeters (ECAL 22X0, 0.5x0.5cm, HCAL ). for the sake of simulation, some detector technologies are assumed in Mokka. But as ILD, many detector technologies are open and not selected at the time of LOI. By the time of LOI, we have no time to merge Jupiter/Sattelites and Mokka/Marlin framework. A work to merge two framework for ILD Software will come after LOI. 22
Mokka model CAD Model
ILD_00 MC/DST production ILD performance are expected to be similar to GLDPrim/LDCPrim. But for consistent and complete study, new MC&DST production has been lunched with an improved software. Started since Dec. last year, using GRID Goal: 250 fb-1 @ 250 GeV, 500 fb-1@500 GeV, Signal + SM background StdHep (@SLAC) Sim(Mokka), reconstruction and DST maker. DST contents: lcio format contains : Tracks, PFOs, [23456]-Jets, LCFIVertex, MCParticls,.. Production profile: Typical CPU time: ~0.5 min.(mm) to 4min.(6f) Typical event size ( for uds-pair @ 500 GeV ) Sim. ~950kb, Rec.~1800kb, DST ~ 23kb 25
ILD_00 MC/DST production Ecm Signal Events NEvents L [1/fb] 250 Sig1: ZH, Z ee/mm 105k 4624 Sig2: ZH, Z nn, H qq 194k 1000 Sig3: ZH, Z qq, H qq 567k 1000 500 Sig4: ee tt 2385k 517 Sig5: ee tt bbqqqq 1012k 3737 Sig6: ee + 1 + 1, 0 2 0 2 137k 678 @ last week Rough summary: ~ 23M events O(50) TB Sim/Rec. files ~ 0.5 TB DSTs so far Production continues SM (250GeV) NEvents L [1/fb] 2f (w/o ee) 1314k 3.8 4f 10631k 772 6f 200k 3753k ee 95k 0.0014 gg X 0 0 eg eg 0 0 nn+ng 0 0 gg+ng 0 0 SM (500GeV) NEvents L [1/fb] 2f (w/o ee, tt) 760k 14.3 4f 2610k 42.4 6f (w/o bbqqqq) 624k 1221 ee 0 0.0 gg X 544k 3560 eg eg 0 0 nn+ng 0 0 gg+ng 0 0 26
GRID for MC production GRID provides Huge CPU and storage resources A way to communicate world wide VO ILC is hosted by DESY, based on LHC Computing GRID MC production and production are running on GRID Simulated, Reconstructed, and DST are placed on GRID. DST: 20~50 MB x O(10k) files or more. In Japan, Replications of DST to KEK/Tohoku/Kobe U. sites are in progress in parallel to the production. Resources in KEK will be increased in near feature. IP2P3 Grid more DESY KEK UK 27
Summary ILD has been optimized ILD MC and DST production is in progress, for performance studies of LOI 3 rd Workshop will be held at Seoul in Feb 16-18, LOI is due March 31. Presented at TILC09 ( 17-21, April ) Detector TDR phase will follow. Many physics channels are yet to be analyzed. Your participations are welcomed. 28
Backup Slides