Transcript Using Jupiter and Satellites Akiya Miyamoto KEK Feb 2007 Target of Jupiter/Satellites/Uranus  Jupiter and Satellites are tools for detector optimization based on Geant4 Full detector.

Using
Jupiter and Satellites
Akiya Miyamoto
KEK
Feb 2007
1
Target of Jupiter/Satellites/Uranus
 Jupiter and Satellites are tools for detector
optimization based on Geant4 Full detector simulation.
 Implement an ideal “GLD” geometry in Jupiter
( other geometry is also possible )
 Study PFA performances by an ultimate condition
 Implement “tower” calorimeter
 Develop analysis tools
 Study physics performance vs. detector choice.
 Not for simulation studies for beam test studies
 A study by detector simulation is also crucial for
the design of the IR design of ILC.
2
Jupiter/Satellites Concepts
For real data
Tools for simulation Tools
JUPITER
JLC Unified
Particle Interaction
and
Tracking EmulatoR
Satellites
IO
Input/Output
module set
METIS
Monte-Calro Exact hits To
Intermediate Simulated output
URANUS
Unified Reconstruction
and
ANalysis Utility
Set
LEDA
Geant4 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
Jupiter Directory Structure
Jupiter
bin
config
include
lib
macros
Readme
Linux-g++
Jupiter
doc
sol
lcexp
kern
cal
cd
c
Jupiter.cc
Detector components to use
are switched by CPP flags in
Jupiter.cc or parameters in JSFJ4
sources
tmp
Sample G4
scripts
common.gmk
main
scripts
src
ct
tpc
vtx
it
bd
ir
include
Sub-Detector files
4
Standalone Jupiter
$ Jupiter [options] [g4macro file]
options:
-v N : N is a message level.
> 10 ; Write message every time new geometry is created.
5-10 ; Write message of geometry whose rank is less than N-4.
< 5 ; No message is written at geometry creation.
-f FILE : Read parameters from the file FILE
More than one -f option can be specified. In this case
a value in the first file is used.
-w
: Write default parameters defined in the program
to the file, jupiter.defaults
-help : Print help message
Useful to display events and debug
5
G4macro file
 Commands to Geant4 and Jupiter are provided as a
g4macro file. A list of commands can be found at
http://www-jlc.kek.jp/~miyamoto/Jupiter/html/_.html
 When Jupiter is executed as a standalone application
interactively, help command can be used to get help
information.
6
Jupiter parameter list file
 Parameters for detector configuration can be
modified by a parameter file given by –f option.
 $JUPITERROOT/doc/HowToRun for more details
 Standard configuration file as of Jan 2006 is at
http://ilcphys.kek.jp/soft/samples_dod/dec05/index.html
7
Satellites/Uranus
8
JSFJ4
 Interface to run Jupiter in JSF environment.
 Allows to write
 Jupiter output by ROOT or LCIO
 Use generators in JSF or StdHep format file.
 Satellites package reads ROOT file created by JSFJ4
9
Metis package
 Metis is a collection of reconstruction tools for Jupiter
data.
 Current aim is to prepare a minimum set of Metis
modules for studies of Particle Flow Algorithm.
 Novice users will be able to do physics analysis using
information of PFO classes.
 As a first step, a cheated track finder and a cluster
maker, etc are in preparation in order to know ultimate
performance.
 Each module is independent, thus shall be easy to
implement different reconstruction algorithm according
to interests
 Development of a real clustering code begins recently
10
Metis Analysis Flow
Jupiter result
make smeared TPC hits from exact hit
make tracks from TPC
make hybrid tracks ( TPC+IT+VTX)
make smeared/merged CAL hits from exact hit
make cluster from CAL hits
make Particle Flow Objects
jet clustering
Physics study
11
Satellites Directory Structure
Satellites
io
jsfj4
bin
kern
examples
include
mctruth
lib
src
j42lcio
test
cal
vtx
cd
c ….
Run Jupiter in JSF
to create a ROOT
file
S4xxxExactHit class
= J4xxxHit class
Output LCIO data
examples
macro
Leda
metis
12
Metis Directory Structure
metis
(JSF’s) Modules for MC data analysis
hitmaker
make smeared TPC hits from exact hit
trackmaker
make tracks from TPC
hybt
hybtmaker
make hybrid tracks ( TPC+IT+VTX)
cal
hitmaker
make smeared/merged CAL hits from exact hit
clustermaker
make cluster from CAL hits
pfo
pfomaker
make Particle Flow Objects
jet
jetmaker
make jet Objects
tpc
13
Uranus Directory Structure
Uranus
data
include
Packages for real data analysis
kern
lib
VTX
detconfig
src
IT
hitmaker
TPC
trackmaker
hybt
14
Relation among TrackHits
TObjArray
U4VTrack
U4VTrackHit
double fChi2
U4VTrackMaker* fU4V
JSFHelicalTrack* fHT
double fResidual
double fT0
U4VHit *fHitPointer
U4XXXTrack
U4XXXTrackHit
U4XXXHit
S4XXXTrack
S4XXXTrackHit
S4XXXHit
double fError
S4XXXExactHit *Hit
15
Relation among CAL Hits
TObject
TAttLockable
S4VHit
+GetExactHitsPtr()
+GetAddress()
+GetHitPosition()
S4CALHit
+ GetCALType()
+GetClusterType()
S4VExactHit
GetAddressPtr()
GetSmearedHitPtr()
S4CALExactHit
+ GetPreHitPtr()
16
Typical results
17
Momentum resolution
spt/pt2 (GeV -1)
Exact hit points created by single m were fitted by Kalman filter package
By A.Yamaguchi(Tsukuba)
18
EM Cal Performance
By A.L.C.Sanchez (Niigata U.)
Linearity is good, but to get energy resolution similar to beam test results,
Randge cut of O(1) mm is required. It is very small and we don’t know
why we need O(1) mm.
19
Cheated Particle Flow Analysis
A sample detector signal
 Key point of PFA is use Cal. signals only
for neutral particles
 Remove CAL signal if connected to a
track
 Use Simulation information to connect
track and cal. signals
 Cheated PFA
 What affects jet energy resolution
 Signal sampling fluctuation in Cal.
 Tracker resolution
 Treatment of V0, decays,and
interactions before Cal.
 … more
 Understand factors which affect
resolution
20
Cheated PFO analysis
ZH event at Ecm=500 GeV
- Exact hit points of TPC and CAL
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)
21
ROOT’s X3d view of the
same event
X3D
X3D-Jet
Same event, after a forced
4-jet clustering on PFObjects
e+e-  Z  q qbar at sqrt(E)=91.1GeV
Jet energy resolution
By Sumie Yamamoto
24
DE for Perfect CAL./Tracker
By Sumie Yamamoto
Jet energy resolution when resolution of CAL
and tracker are perfect.
Contribution of each detector
GeV/c2
Total
HD Cal
EM Cal
Tracker
Others
2.7 GeV
2.0
1.2
0.4
1.4 (under investigation)
25
Satellites examples









Exam01
Exam02
Exam03
Exam04
Exam05
Exam06
Exam07
Exam08
Exam09
 Analysis code examples
26