ILD Detector Optimization and Benchmarking Akiya Miyamoto KEK

ILD Detector Optimization and Benchmarking Akiya Miyamoto, KEK at Tsinghua University 12 -January-2009

ILD Introduction n ILD origins in the European and Asian based Large Detector study. n ILC Reference Design (RDR) in 2007 u GLD Detector Outline Document (DOD) ar. Xiv: physics/0607154 u LDC DOC http: //www. ilcldc. org/ Common feature: Tracker(Pixel & Silicon & Gas) + PFA calorimeter + … LDC GLD n At LCWS 2007, we agreed to work together for a joint LOI u GLD (B=3 T, RECAL=2. 1 m) + LDC(B=4 T, RECAL=1. 6 m ) ILD Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 2

International Large Detector ( ILD ) n LDC and GLD had a common future; u Pixel vertex detector placed very close to the beam pipe. u Gaseous tracker, TPC, for highly efficient and precise track measurements, supplemented by silicon trackers. u 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). n But in detail: u u B=4 Tesla(LDC) vs 3 Tesla(GLD) ECAL radius: 1. 6 m(LDC) vs 2. 1 m(GLD) Sub Detector technologies …. n Simulation studies of physics performances are used to reach agreement of detector parameters. Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 3

How we optimize n Optimization tools u GLD Jupiter/Sattelites, LDC Mokka/Marlin. Reco u intermediate detector models were introduced for comparison l GLDPrim by Jupiter, and LDCPrim by Mokka, both having B=3. 5 T and RECAL=1. 85 m. n Performances have been studies as a function of major parameters. n Reached a consensus on the ILD reference detector for LOI benchmark studies at Cambridge (Sep. 2008). n Physics performance studies have been performed based on ILD model. Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 4

Jupiter/Satellites for Full Simulation Studies : GLD For real data Tools for simulation Tools JUPITER JLC Unified Particle Interaction and Tracking Emulato. R Satellites IO Input/ Output module set METIS Monte-Calro Exact hits To Intermediate Simulated output URANUS Unified Reconstruction and ANalysis Utility Set LEDA Geant 4 based Simulator MC truth generator Library Extention for Data Analysis JSF/ROOT based Framework Event Reconstruction JSF: the analysis flow controller based on ROOT The release includes event generators, Quick Simulator, and simple event display Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 5

Mokka LDC − Mokka is a full simulation using Geant 4 and a realistic description of a detector for the future linear collider. − Home page: http: //polzope. in 2 p 3. fr: 8081/MOKKA Mokka is now a part of the ilcsoft, http: //ilcsoft. desy. de/portal/software_packages/ − Detector Geometry: − managed by My. SQL 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 Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 6

Pandora. PFA LCFIVertex

GLD + LDC Combined Framework Whizard Physsim Std. Hep LDC MOKKA GLD Jupiter LCIO Marlin Sattelites LCIO DST and Analysis Akiya Miyamoto, KEK üStd. Hep: Same generator data üLCIO: Common IO format ü GLDPrim/LDCPrim: Similar detector model LCIO helps to collaborative works for detector optimization After the LOI, two frameworks will be merged to a single framework. Seminar at Tsuingha Univ. , 12 -Jan-2009 8

Detector Parameters for Opt. studies l GLD/GLDPrim/J 4 LDC prepared for Jupiter l LDC/LDCPrim/LDCGLD prepared for Mokka Physics performance was compared between different geometries Jupiter Mokka GLDPrim J 4 LDC LDCGLD LDCPrim LDC B(T) 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 4 layers 2 layers # IT layers TPC Rmin(cm) 43. 7 43. 5 34. 0 ECAL Rmin(cm) 210 185 160 HCAL Thick. (Int. L) 6. 79 6. 29 5. 67 37. 1 202 182. 5 161 5. 86 Geometryiesin Mokka and Jupiter are similar, but there are many small differences in geometry and assumed detector technologies Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 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 Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 10

GLDPrim - LDCPrim 4 mm LDCPrim(Mokka+Pandora) is better than GLDPrim(Jupiter+Sattelites) by 15~30%. Possible source: srf(IT) 4 mm(LDCPrim) 10 mm(GLDPrim) Silicon External Tracker in Mokka 3 x 10 -5 Sub-detector technology is more important than geometry Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 11

GLDPrim vs LDCPrim (srf(IP)) n. GLDPrim is better than LDCPrim ; ü 3 double layers vs 5 layers ? srf=s. Z=2. 8 mm Fast sim. study by M. Berggren Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 12

kaon_0 L Energy Resolution Hadron Model: LCPhysics HCAL response is not smooth around 13 Ge. V - LE/HE behaviour ECAL resolution: same Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 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. Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 14

Jupiter data analyzed by Pandora. PFA Pandora PFA: Sophisticated algorithm tuned to Geant 4 shower shape has achieved the performance goal of ILC, DE/E ~ 30%/√E Z pole uds-pair events: GLDPrim Ejet(Ge. V) Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 15

Jet Energy Resolution Same trend is seen by analysis of Jupiter models, though performance Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-200916 is slightly worse than Mokka model

ECAL Seg. and HCAL thickness by M. THomson ILD: ECAL+HCAL= 6. 8 Int. L. (48 layers) 6. 8 Int. L look OK, but worse resolution is seen for 90 o jets. Akiya Miyamoto, KEK −Performance is strong function of ECAL seg. size. − 2 x 2 cm 2 too large, 1 x 1 cm 2 would be ok for jets with E < 100 Ge. V Seminar at Tsuingha Univ. , 12 -Jan-2009 17

Physics Benchmark Studies n ILC goal precise studies of Tera scale physics. n According to the request by ILC Research Director (RD) and International Detector Advisory Group(IDAG), simulation studies for LOI should u based on a realistic Monte Calro program u based on a realistic reconstruction program u include backgrounds by physics processes and those caused by accelerator. n Signal processes: the minimum set. u Recoil mass measurement by e+e- ZH e+e-/m+m- + H u H c cbar decay in e+e- ZH process u e+e- t tbar 6 jets and t (tbar) charge ID for AFB meas. u e+e- t+ t- and t pol. measurement. u Separate WW and ZZ in Chargino/Neutrino pair production process Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 18

Higgs recoil mass meas. n e+e - ZH e+e-X / m+m-X , Ecm=250 Ge. V, 250 fb-1 n Analysis. by Itoh Kazutoshi u Select e+e- / m+m - consistent with Z and study recoil mass u Precise track meas. is a key for GLDprim case, with backgound e+e - Channel Akiya Miyamoto, KEK m+m- channel Seminar at Tsuingha Univ. , 12 -Jan-2009 Compare 3 geometries m+m- X Differences are small. 19

Benchmark study: Example n Using several detector models, performance to separate W/Z in jet mode have been studied using SUSY processes Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 by Taikan Suehara 20

e+e- t+t-, t rn nt is polarized probe New Phyaics by Taikan Suehara SM+NP SM Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 21

ILD reference detector model n At 2 nd ILD WS at Cambridge, we agreed to created the new model, ILD reference design model for LOI, in Mokka: n Model parameters, u B=3. 5 Tesla u Rin ECAL=185 cm, TPC: half. Z=230 cm u VTX three double layers. u Silicon trackers: ( SIT, FTD, SET, SOT) u Calorimeters (ECAL 22 X 0, 0. 5 x 0. 5 cm, HCAL ) u …. 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. Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 22

Mokka model CAD Model

ILD_00 MC/DST production n 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. n Started since Dec. last year, using GRID u Goal: 250 fb-1 @ 250 Ge. V, 500 fb-1@500 Ge. V, Signal + SM background n Std. Hep (@SLAC) Sim(Mokka), reconstruction and DST maker. n DST contents: u lcio format u contains : Tracks, PFOs, [23456]-Jets, LCFIVertex, MCParticls, . . n Production profile: n Typical CPU time: ~0. 5 min. (mm) to 4 min. (6 f) n Typical event size ( for uds-pair @ 500 Ge. V ) n Sim. ~950 kb, Rec. ~1800 kb, DST ~ 23 kb Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 25

ILD_00 MC/DST production Ecm Signal Events 250 500 NEvents L [1/fb] Sig 1: ZH, Z ee/mm 105 k 4624 Sig 2: ZH, Z nn, H qq 194 k 1000 Sig 3: ZH, Z qq, H qq 567 k 1000 2385 k 517 1012 k 3737 137 k 678 Sig 4: ee tt Sig 5: ee tt bbqqqq Sig 6: ee c+ SM (250 Ge. V) 1 c+ 1, c 0 2 NEvents L [1/fb] SM (500 Ge. V) Production continues NEvents L [1/fb] 14. 3 2610 k 42. 4 624 k 1221 0 0. 0 544 k 3560 eg eg 0 0 0 nn+ng 0 0 0 gg+ng 0 0 3. 8 4 f 10631 k 772 6 f 200 k 3753 k 6 f (w/o bbqqqq) ee 95 k 0. 0014 ee gg X 0 0 gg X eg eg 0 0 nn+ng 0 gg+ng 0 Akiya Miyamoto, KEK Rough summary: ~ 23 M events O(50) TB Sim/Rec. files ~ 0. 5 TB DSTs so far 760 k 1314 k 2 f (w/o ee) @ last week 2 f (w/o ee, tt) 4 f Seminar at Tsuingha Univ. , 12 -Jan-2009 26

GRID for MC production n GRID provides u Huge CPU and storage resources u A way to communicate world wide n VO ILC is hosted by DESY, based on LHC Computing GRID u MC production and production are running on GRID u Simulated, Reconstructed, and DST are placed on GRID. DST: 20~50 MB x O(10 k) files or more … IP 2 P 3 n In Japan, DESY u Replications of DST to KEK/Tohoku/Kobe U. sites are in progress in parallel to the production. u Resources in KEK will be increased in near feature. UK Grid KEK Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 27

Summary n ILD has been optimized n ILD MC and DST production is in progress, for performance studies of LOI u 3 rd Workshop will be held at Seoul in Feb 16 -18, u LOI is due March 31. Presented at TILC 09 ( 17 -21, April ) u Detector TDR phase will follow. n Many physics channels are yet to be analyzed. Your participations are welcomed. Akiya Miyamoto, KEK Seminar at Tsuingha Univ. , 12 -Jan-2009 28

Backup Slides