ilc @ kek ipns and grid akiya miyamoto kek at 28-feburary-2007 fjppl meeting
TRANSCRIPT
ILC @ KEK IPNSand
GRIDAkiya Miyamoto
KEKAt 28-Feburary-2007
FJPPL meeting
Physics Scenario at ILC
Detector Concepts SiD
– Silicon tracker, 5T field("small" radius)– SiW ECAL– http://www-sid.slac.stanford.edu/
LDC – TPC, 4T field– SiW ECAL (“medium” radius)– http://www.ilcldc.org/
GLD – TPC (+Silicon IT), 3T field– W/Scintillator ECAL (“large” radius)– http://ilcphys.kek.jp/gld/
4th Concept – Dual Readout(Scintillator+Cristal) Calorimeter– Dual Air Core Coil – http://www.4thconcept.org/
GLD Design Concepts
Separate neutral and charged particle PFA == Key for the good energy resolution
Segmentation of calorimeter : do to the limit of Moliere Radius Large radius calorimeter : spatially separate particle
( also good for tracker ) High B field : separate charged and neutral
Figure of Merit
2
2 2M
BR
R
: CALgranularity
:EffectiveMoliere radiusMR
Neutral energy inside certain distance from cahrgedscale ~ 1/R2
GLD Configuration
GLD Side view Moderate B Field : 3T R(ECAL) ~ 2.1m
ECAL: 33 layers of 3mmt W/2mmt Scint./1mmt Gap HCAL: 46 layers of20mmt Fe/5mmt Scint./1mmt Gap
Photon sensor: MPPC ~O(10M) ch. Configuration of sensor is one of the R&D item
s
GLD Configuration - 2
TPC:R: 0.452.0m, ~200 radial sampleHalf Z: 2.3mMPGD readout: r<150m
SIT: Silicon Strip Barrel/Endcap
VTX:Fine Pixel CCD: ~5x5mm2
2 layers x 3 Super Layers
cos 0.9
GLD organizationMember :16 countries, 77 Univ./Inst. 224 members
Contact Persons H.Yamamoto, H.B.Park (Asia), G.Wilson(NA) R.Settles, M.Thomson(EU)
UK 5Germany 3Italy 2Netherlands 1Rusia 1
Japan 28 Philipine 2Korea 8 Australia 2China 5 India 4Singapole 1 Vietnum 1
USA 11Canada 1
# inst.
Executive boardS.Yamashita - BenchmarkA.Miyamoto - SoftwareY.Sugimoto - Vertex DetectorH.J.Kim - Intermediate TrackerA.Sugiyama/R.Settles – TPCT.Takeshita - Calorimeter/MuonT.Tauchi - Interaction RegionH.Yamaoka - Coil & StructureP.Ledu - DAQM.Tomson - Space
GLD Concepts has been developed through E-mails and TV meetings discussion http://ilcphys.kek.jp/gld [email protected]
GLD DOD: physics/0607154
Task forces (since March 2006) IR (T.Tauchi ) PFA (T.Yoshioka) Tracking ( to be decided)
ILC crossing angle, Detector hall, push/pull options, etc are hot topics in recent meetings
GLD is a team for detector concept studies
Not a collaboration
ILC @ KEK IPNS
Hardware studies Vertex Detector Fine Pixel CCD (~5x5m2 pixel size )for ILC
Vertex Detector
TPC with MPGD readoutAs a member of LCTPC Collaborationparticipate beam tests at DESY
Software study Detector simulator based on Geant4 Event reconstruction tools Physics studies
Simulation
ROOT objects : Event Tree & Configuration
Our software tools
BeamtestAnalysis
EventReconstruction
Digitizer Finder Fitter
DetectorSimulator
QuickSim FullSim
EventGenerator
Pythia CAIN StdHep
PhysicsAnalysis
Jet finder
Link to various tools at http://acfahep.kek.jp/subg/sim/soft GLD Software at http://ilcphys.kek.jp/soft All packages are kept in the CVS. Accessible from http://jlccvs.kek.jp/
JSF
Framework: JSF = Root based application All functions based on C++, compiled or through CINT Provides common framework for event generations, detector simulations,
analysis, and beam test data analysis Unified framework for interactive and batch job: GUI, event display Data are stored as root objects; root trees, ntuples, etc
Release includes other tools QuickSim, Event generators, beamstrahlung spectrum generator, etc.
Example of QuickSim Study; ,e e X e
Incl. beamstrahlung350GeV, nominal(Mh)~109MeV
Incl. beamstrahlung350GeV, high-lum(Mh)~164MeV
Incl. beamstrahlung250GeV, nominal(Mh)~27MeV
E/E(beam)~0.1%Differential Luminosity(500GeV)
min/ no als s
Jupiter/Satellites for Full Simulation Studies
JUPITERJLC Unified
Particle Interactionand
Tracking EmulatoR
IOInput/Outputmodule set
URANUS
LEDA
Monte-Calro Exact hits ToIntermediate Simulated output
Unified Reconstructionand
ANalysis Utility Set
Library Extention for
Data Analysis
METISSatellites
Geant4 basedSimulator
JSF/ROOT basedFramework
JSF: the analysis flow controller based on ROOT The release includes event generators, Quick Simulator, and simple event display
MC truth generator Event Reconstruction
Tools for simulation Tools For real data
Jupiter feature - 1
Modular structure easy installation of sub-
detectors
Currently using Geant4 8.0p1
Geometry Simple geometries are
implemented ( enough for the detector optimization )
parameters ( size, material, etc ) can be modified by input ASCII file.
Parameters are saved as a ROOT object for use in Satellites later
Jupiter feature - 2
Run mode: A standalone Geant4 application JSF application to output a ROOT file.
Output: Exact Hits of each detectors (Smearing in Satellites) Pre- and Post- Hits at before/after Calorimeter
Used to record true track information which enter CAL/FCAL/BCAL.
Break points in tracking volume Interface to LCIO format is prepared
in the JSF framework Compatibility is yet to be tested.
Break point
Post-hits
Input: StdHep file(ASCII), HepEvt, CAIN, or any generators implemented in JSF. Binary StdHep file interface was implemented, but yet to be tested.
GLD Geometry in Jupiter
FCAL
BCAL
IT
VTX
CH2mask1 module
Include 10cm air gap as a readout space
DID Field and Backgrounds
High energy
Low energy
Exit hole
cm
Detector Installed Dipole-magentic Field : reduce background hits
Jupiter Simulation
FCALBCAL
CH2Mask
e+e- hit distribution at BCAL
4m入射ビーム
Typical Event Display
- ZH → h : Two jets from Higgs can be seen.
Jet Measurements in ILC Det. Particle reconstruction
Charged particles in tracking DetectorPhotons in the ECALNeutral hadrons in the HCAL (and possibly ECAL)b/c ID: Vertex Detector
Large detector – spatially separate particles High B-field – separate charged/neutrals High granularity ECAL/HCAL – resolve particles
For good jet erngy resolution Separate energy deposits from different particles
e+
e-
Realistic PFA Critical part to complete detector design
Large R & medium granularity vs small R & fine granularity Large R & medium B vs small R & high B Importance of HD Cal resolution vs granuality …
Algorithm developed in GLD: Consists of several steps Small-clustering Gamma Finding Cluster-track matching Neutral hadron clustering
Red : pionYellow :gammaBlue : neutron
- Performance in the EndCap region is remarkably improved recently.- Almost no angular dependence : 31%/√E for |cos|<0.9.
All angle
- Z → uds @ 91.2GeV, tile calorimeter, 2cm x 2cm tile size
Jet Energy Resolution (Z-pole)
T.Yoshioka (Tokyo)
Next step !
Detector configuration optimization.
Other energy points & physics channels
Benchmark Processes
Benchmark processes recommended by the Benchmark Panel.
GRID for ILC studies
ILC GRID in Japan has just begun
ILC Computing requirements in coming years Full simulation based detector studies
Detector optimization Background studies & IR designs, etc These studies will be based on many bench mark processes. Cross checking of analysis codes
Sharing and access to beam test data
GRID will be an infrastructure for these studies. Data sharing Utilize CPU resources Among domestic/regional colleagues – easier access to codes &
data
We are beginner. as a first step, First: minimum resources Get familiar tools and data sharing
CALICEVO & ILCVO
GRID Configuration: JPY2006
KEKCC LCG environment
LCG resource outside KEK
GRID-LAN F/W
KEK-LAN F/W
NFSSLC4
Existing KEK-ILC group resource
ILC NFS
SE
LCG
LCG
LCG/Grid Protocol
SE: Storage Element
LCG/Grid Protocol
University
CPUServers
LCG UI
SLC3
SE2
LCG UI
SLC3 TohokuKobe
GRID : Future Direction
KEKCC LCG environment
LCG resource outside KEK
GRID-LAN F/W
KEK-LAN F/W
NFSSE
LCG
LCG
LCG/Grid Protocol
SE: Storage Element
LCG/Grid Protocol
Universities
SE2
No GRID univ’s
TohokuKobe
ILC Resource
KEK User PC’s
GIRD sites
ILC VO, CALICE VO, JHEP VO, …
Backup slides
Description Detector Language IO-Format Region
Simdet Fast Monte Carlo TeslaTDR Fortran Stdhep/LCIO EU
SGV Fast Monte Carlo flexible C++ None(LCIO) EU
Lelaps Fast Monte Carlo SiD, flexible C++ SIO, LCIO US
QuickSim Fast Monte Carlo GLD Fortran ROOT Asia
Brahms-Sim Full sim. - Geant3 TeslaTDR C++ ASCII, LCIO EU
Mokka Full sim. – Geant4 TeslaTDR, LDC C++ LCIO EU
SLIC Full sim. – Geant4 SiD C++ LCIO US
ILC-ROOT Full sim. – Geant4 4th C++ ROOT US+EU
Jupiter Full sim. – Geant4 GLD C++ ROOT, LCIO Asia
Brahms-Reco Reconstruction framework TeslaTDR Fortran LCIO EU
Marlin Reconstruction Analysis framework
Flexible,LDC C++ LCIO EU
Org-lcsim Reconstruction packages SiD(flexible) Java LCIO US
Satellites Reconstruction packages GLD C++ ROOT Asia
LCCD Conditiions data toolkit LDC, SiD, .. C++ MySQL, LCIO EU
GEAR Geometry Description Flexible C++ XML EU
LCIO Persistency/Datamodel All C++,Java,Fortran
- EU,US,Asia
JAS3/WIRED Analysis tool/Event display LDC, SiD … Java XML,LCIO,stdhep, heprep, US, EU
JSF Analysis framework All C++ ROOT/LCIO Asia
Software tools in the world
FJPPL
• France Japan Particle Physics Laboratory• France-(CNRS-CEA: LAPP, LLR, LPNHE,
LAL, DAPNIA); Japan-KEK • http://acpp.in2p3.fr/cgi-bin/twiki/bin/view/FJH
EPL/WebHome– You have to register yourself to access internal
info.
• ILC detector R&D is one of the main projects of FJPPL.
K.Kawagoe, Sep. 2006
“A common R&D on the new generation detector for the ILC”
• Members– J.C. Brient, H. Videau, J.C. Vanel, C. de la Thaille, R. Po
eschl, D. Boutigni– K.Kawagoe, T. Takeshita, S. Yamashita, T. Yoshioka, A.
Miyamoto, S. Kawabata, T.Sasaki, G. Iwai• Goals(1) Design and development of reliable Particle Flow
Algorithm (PFA) which allows studying the design and the geometry of the future detector for the ILC,
(2) Design and development of a DAQ system compatible with the new generation of calorimeter currently in design and prototype for the ILC, and
(3) Design and development of the optimized detector for the final ILC project, leading the participation of the LOI and TDR document for the ECAL point of view.
K.Kawagoe, Sep. 2006
CPU time/Data size
Data size
Xenon 3 GHz(32bit)
for 500 1/fb
Process MB/ev CPU sec/ev. #ev. GB cpu day
qq 91 GeV ~1.5 ~150
qq 350 GeV ~3.0 ~270
eeZHH 350GeV ~2.0 ~300 14k 28 48.6
eeH 350GeV ~1.9 ~170 15k 29 29.5
eeeeH 350GeV ~2.3 ~380 1.5k
3.5 15.4
eeZZqq 350GeV ~1.7 ~200 110k 187 255.6
eeeWeqq 350GeV ~1.6 ~240 1000k 16000 2777.8
ZHqqH 350GeV ~3.8 ~300 49k 186 170.1
ZZqqqq 350GeV ~3.3 ~340 434k 1432 1707.8( cpu time is about half at KEKCC, AMD 2.5GHz 64bit )
Cross sections
5k events/4y
SM processes+ New physics