Status of Software Tools Akiya Miyamoto KEK 8

Status of Software Tools Akiya Miyamoto KEK 8 th ACFA, EXCO, Daegu 12 July 2005 Short summary of work done by K. Fujii(KEK), T. Fujikawa(Tohoku), K. Kawagoe(Kobe), T. Takeshita(Shinshu), H. Matsunaga(Tsukuba), A. Yamaguchi(Tsukuba), Y. Yamaguchi(Tsukuba), T. Nagamine(Tohoku), S. Yamamoto(GUAS), H. Ono(Niigata), T. Yoshioka(Tokyo), A. L. C. Sanchez(MISIIT), … 1

Overview of our tools n lcbase : configuration files n Leda : Analysis tools (Kalman fitter, 4 vector and jet findinder utilities ) n jsf : Root-based framework n lclib : Quick. Sim and other fortran based utilities n physsim : Helas-based generator n Jupiter : Full simulation based on Geant 4 n Uranus : Data analysis packages n Satellites : Data analysis packages for MC data Ø We use only C++, except old fortran tools. A. Miyamoto , 8 th ACFA, 12 July 2005 2

Quick. Sim overview n Purpose: simulate detector effects for physics study n Components: w VTX, IT, TPC, CAL n Model for tracker w circular trajectory parabolic trajectory w With multiple scattering, without energy loss w Equally spacing sampling n Model for Calorimeter w EM signal by e/g, HD signal by hadron, muon no signal w Segmented calorimeter. Lateral spreads are generated by an analytic form. A. Miyamoto , 8 th ACFA, 12 July 2005 3

Quick. Sim update n Change in source code w Gas density in tracker, Cal. Shower shapre can be changed by input data card. n GLD-V 2 : new parameter file for GLD configuration w Beam pipe: r=1. 5 cm, Be, 25 mm w VTX : (r, thickness)=(1. 9 cm, 20 mm), (4. 1, 200), (6. 3, 200) ß not consistent with Sugimoto san’s parameter yet. ß Beam pipe and VTX 1 st layer parameter are not realistic to reproduce performances. u u TPC: r=40~205 cm, 200 sampling, u srf=150 mm, sz=1 mm CAL : EM segmentation 2 cm equivalent HD segmentation 6 cm equivalent Parameters for lateral spread have not been changed. A. Miyamoto , 8 th ACFA, 12 July 2005 4

Quick. Sim Tracker performance Momentum resolution similar to the Jupiter/Satellites result IP resolution at 1 Ge. V is about factor 2 worse than the Jupiter/Satellites result A. Miyamoto , 8 th ACFA, 12 July 2005 5

Higgs signals by Quick. Sim Mh=118. 4 Ge. V/c 2 s. M=4. 02 Ge. V/c 2 Mh=117. 6 Ge. V/c 2 s. M=4. 808 Ge. V/c 2 Criteria for event selections are same Mh=116. 2 Ge. V/c 2 s. M=3. 81 Ge. V/c 2 A. Miyamoto , 8 th ACFA, 12 July 2005 Possible improvments: - parameters of track-cluster matching - detector configuration dependant event selection - …… 6

Jupiter/Satellites Concepts For real data Tools for simulation Tools Satellites JUPITER IO JLC Unified Particle Interaction and Tracking Emulato. R 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 A. Miyamoto , 8 th ACFA, 12 July 2005 7

Jupiter feature - 1 n Based on Geant 4 7. 0 p 1 n Modular organization of source codes for easy installation of sub-detectors n Geometry w parameters ( size, material, etc ) can be modified by input ASCII file. w Only simple geometries are implemented ( OK for the detector optimization ) n Input: w Std. Hep file(ASCII, Binary), Hep. Evt, CAIN, or any generators implemented in JSF. A. Miyamoto , 8 th ACFA, 12 July 2005 8

Jupiter feature - 2 n Run mode: w A standalone Geant 4 application w JSF application to output a ROOT file. n Output: w Exact Hits of each detectors (Smearing in Satellites) w Pre- and Post- Hits at before/after Calorimeter Used to record true track information which enter CAL. w Format: ROOT (standard), ASCII(debug), LCIO n Physics List w LCPhysics. List w QGSP, LHEP, … Calorimeter resolution : Depends on physics list, and not as good as we expected for single particle A. Miyamoto , 8 th ACFA, 12 July 2005 9

Geometry in Jupiter Muon/Iron Solenoid TPC Hadron Calorimeter Elemag. Calorimeter A. Miyamoto , 8 th ACFA, 12 July 2005 QC 1 IT Forward Cal. VTX 10

Calorimeter Geometry Full One Tower 　EM + HD 6. 1λ 210 cm 40 cm 27 X 0 270 cm n. Default sensor size: ØEM: 4 cmx 1 mm, 38 layers ØHD: 12 cmx 2 mm, 130 layers n. Replica ØPhi direction : Tower and mini-tower n. Sandwich structure of X/Y scinti structure can be defined. A. Miyamoto , 8 th ACFA, 12 July 2005 11

Detector Parameters n GLD_V 1: as an initial set of detector parameters è Beam Pipe : 1. 8 cmf, Be 500 mmt è VTX : 4 layers from R=2. 4 cm to 6. 0 cm, Si 330 mmt, |cosq|<0. 9 srj=sz=4 mm, Cylindrical shape, not FPCCD geometry yet è IT : 4 layers from R=9 cm to 30 cm, Si 600 mmt, |cosq|<0. 9 srj=sz=10 mm, Cylindrical shape è TPC : R=40 cm to 206 cm, |Z|<235 cm, # of sampling 200 srj=150 mm, sz=500 mm, Gas P 10 è Solenoid : 3 T è Muon : Barrel : 4 mud layers, total 270 cmt Fe Endcap : 5 mud layers, total 260 cm t Fe A. Miyamoto , 8 th ACFA, 12 July 2005 Signals are smeared in Satellites All parameters needs to be optimized 12

New IR/VTX geometry VTX Geometry for 2 mrad crossing angle. CH 2 mask QC 1 BCAL(W) W-Mask FCAL(W) IT Beam Pipe: r=1. 5 cm, 250 mmt Be VTX: r=2. 0, 2. 2, 3. 4, 4. 8, 5. 0 cm L=13, 40, 40, 40 (cm) Thickness=50 mm srf=sz=2 mm Work in progress: - 20 mrad crossing - Forward tracker(FT) - Sensitive detector for FCAL/BCAL/FT A. Miyamoto , 8 th ACFA, 12 July 2005 13

Jupiter for background studies n LCBDS w Developed at U. of Tokyo; Tanabe san, Abe san, and their colleagues. w Geant 4 based simulator for the beam delivery and the dump line of ILC ( from the exit of main linac to the beam dump) w Inputs to LCBDS : SAD data file are used to define geometry easy to follow changes of ILC parameters. Accelerator components are approximated by cubic and/or cylindical shape w Output from LCBDS: Synchrotron g’s, muons, off-momentum electrons, etc are created Particles information near the interaction point are saved in Std. Hep ASCII format Inputs to Jupiter w Status First version to write Std. Hep data were prepared. Waiting for the confirmation by Jupiter. A. Miyamoto , 8 th ACFA, 12 July 2005 14

Jupiter for background studies - 2 n CAIN w Simulator for beamstrahlung effects w Generates low energy e+/e-, g, and off-momentum electrons w Interface to CAIN data exists in Jupiter n IR geometry w Still old configuration for GLC. Urgent updates are required. A. Miyamoto , 8 th ACFA, 12 July 2005 15

Typical Background events 100 e+/e- CAIN events 1 bunch contains ~ 66 k events ILC 500 Ge. V nominal parameter TPC BCAL CH 4 mask FCAL Ready for background studies ! A. Miyamoto , 8 th ACFA, 12 July 2005 16

Metis package n Metis is a collection of reconstruction tools for Jupiter data. n Initial aim was to prepare a minimum set of Metis modules for studies of Particle Flow Algorithm. n Novice users will be able to do physics analysis using information of PFO classes. n As a first step, a cheated track finder and a cluster maker, etc are in preparation in order to know ultimate performance. n Each module is independent, thus it is easy to implement different reconstruction algorithms Real clustering codes are prepared using this framework A. Miyamoto , 8 th ACFA, 12 July 2005 17

Metis Analysis Jupiter result • Root file make smeared TPC hits from exact hit make tracks from TPC make hybrid tracks ( TPC+IT+VTX) • Khalman fitter • No track finder yet make smeared/merged CAL hits from exact hit make cluster from CAL hits • Strip/Tile • Cheated/Real Clustering/PFA make Particle Flow Objects jet clustering A. Miyamoto , 8 th ACFA, 12 July 2005 Physics study JSF Framework 18

Cheated PFO analysis ZH event at Ecm=500 Ge. V - Exact hit points of detectors are displayed. -Hits belong to the same PFO are shown with the same color -A framework of event display in JSF is used. By K. Fujii(KEK), S. Yamamoto(GUAS), A. Yamaguchi(Tsukuba) A. Miyamoto , 8 th ACFA, 12 July 2005 19

Web information n Our software tools are maintained in CVS server, jlccvs. kek. jp. n At http: //jlccvs. kek. jp/, w Description about how to download latest version. w Web interface to the CVS repository, http: //jlccvs. kek. jp/cgi-bin/cvsweb. cgi/ w Snap shot of source codes. http: //jlccvs. kek. jp/snapshots/ n Sim. Tools: binary codes of our tools w Web page: http: //acfahep. kek. jp/subg/simtools/index. html w Examples and documents are also prepared. w The latest version: Sim. Tools-1 -01 (May, 2005). New version will be prepared soon. n GLD Software page: w http: //ilcphys. kek. jp/soft/index. html A. Miyamoto , 8 th ACFA, 12 July 2005 20

Summary n Studies of detector performances based on Jupiter/Satellites have been performed by many people, which are covered by other parallel session talks. n Implementation of new geometry in Jupiter has started for studies of background and IR design. n Quick. Sim w Performance of tracker and Higgs signal were studied using GLD-V 2 parameter set. w Detector parameters can be changed easily, but consistency with full simulator is not automatic. Please use it with caution. Jupiter users/developers meet every Tuesday 13: 30 -15: 00 (JST) at TV meeting (KEK MCU, ID 31300). Please join. A. Miyamoto , 8 th ACFA, 12 July 2005 21