Dear Colleagues,

You might be interested in this snapshot of the JLC simulation tools.

                               Michael Peskin


---------- Forwarded message ----------
Date: Fri, 22 Jan 1999 10:10:03 +0900 (JST)
From: [log in to unmask]
To: [log in to unmask]
Subject: A proposal of data format for ACFA JLC Physics study


Dear Colleague


>From several ACFA members, I have requested to consider simple data  
format for easy usage of  detector simulator.  Please send me comments about
the format.  Postscript version of this document is at
http://acfahep.kek.jp/subg/sim/lib/simdst-1-00.ps



regards
Akiya Miyamoto
KEK
E-mail:[log in to unmask]

*********************************************************************************


A proposal of data format for ACFA JLC physics study
-- version 1.00 --

We have Quick Simulator and JIM Full Simulator for the JLC physics
study.  However, output data format of them are slightly different
and they seem to be complicated for simple minded user.
In addition, combined track information output by Quick Simulator
is not sufficient for analysis using the vertex detector information.

In this report, new data format is proposed, which is intended to satisfy 
following conditions
1) Both Quick Simulator and JIM Full Simulator can output same format.
2) Data is written as a Fortran unformated file.
3) Converter to NTUPLE and/or ROOT is also provided.

The information saved in a file are,
1) CDC track parameter and error matrix
   hit points of vertex detector associated to the CDC track
   Track position and direction at the surface of EMC/HDC and 
   associated calorimeter information.
2) Calorimeter information not associated to the track.
3) Generator track information 


Fortran record format:

  CHARACTER*4 STRID, PRODUC
  INTEGER*4   IVERS
  INTEGER*4   NWTOT, NGEN, NTRK, NEMC, NHDC
  REAL*4      HEAD(2)
  REAL*4      GEN(11,n), TRKF(38,j),EMC(5,n), HDC(5,n)
  REAL*8      TRKD(15,j)
  INTEGER*4   NVTX(j)
  REAL*8      VTXD(3,k,j)

  READ(IU) STRID, PRODUC, IVERS, NWTOT, NGEN, NTRK, NEMC, NHDC, 
 >         (HEAD(I),I=1,2),
 >         ((GEN(K,I),K=1,11),I=1,NGEN),
 >         ((TRKF(K,I),K=1,38),I=1,NTRK),
 >         ((TRKD(K,I),K=1,15),I=1,NTRK),
 >         ((EMC(K,I),K=1,5),I=1,NEMC),
 >         ((HDC(K,I),K=1,5),I=1,NHDC),
 >         (NVTX(I),I=1,NTRK),
 >         (((VTXD(K1,K2,I),K1=1,3),K2=1,NVTX(I)),I=1,NTRK),

  STRID = 'LCSM'           ! Helpful to know byte order
  PRODUC = 'QIK' or 'JIM'  ! Data generator
  IVERS =  100               ! Data format version.

  NWTOT = 7  + 11*NWGEN+1 + 38*NTRK+1 + 2*15*NTRK+1 
          + 5*NEMC+1 + 5*NHDC+1 + NTRK + 3*Sum_of_NVTX

  HEAD(1) = Event Number
  HEAD(2) = Sum of initial e- and e+ energy (GeV)

  GEN( ) ! Generator track information (Those used as simulator input)
  GEN(1,i) = Px (GeV/c)
  GEN(2,i) = Py (GeV/c)
  GEN(3,i) = Pz (GeV/c)
  GEN(4,i) = E  (GeV/c)
  GEN(5,i) = vertex X  (cm)  
  GEN(6,i) = vertex Y  (cm)
  GEN(7,i) = vertex Z  (cm)
  GEN(8,i) = Track generated time(nsec)
  GEN(9,i) = Particle ID ( a la PDG)
  GEN(10,i) = Particle mass (GeV/c^2)
  GEN(11,i) = Charge

  TRKF() ! CDC Track parameter, togather with matched VTX hit and Cal information
  TRKF(1:3,j)  = Px, Py, Pz at closest approach
  TRKF(4,j)    = Energy assuming mass less particle
  TRKF(5:7,j)  = x,y,z of closest approach to the beam axis
  TRKF(8,j)    = apparent charge (+/-)
  TRKR(9,j)    = Generator track number i corresponding to the track.
  TRKF(10:14,j)  = Track Parameter, delta-r, phi, kappa, delta-z, tan-lambda
  TRKF(15:17,j) = X0, Y0, and Z0 of Pivot.
  TRKF(18,j)    = Number of measured CDC layer
  TRKF(19:23,j) = (r, theta, phi, delta-theta, delta-phi ) of 
                 the track at EMC entrance ( For end cap, r=+-(10000+ z location))
  TRKF(24,j)    = Number of counters containing 90% of EMC energy 
	( -1 when no EMC association)
  TRKF(25:27,j) = (r, theta, phi ) of EMC cluster
  TRKF(28,j)    = Energy of Associated EMC cluster
  TRKF(29:33,j) = (r, theta, phi, delta-theta, delta-phi ) of 
                 the track at HDC entrance ( For end cap, r=+-(10000+ z location))
  TRKF(34,j)    = Number of counters containing 90% of HDC energy 
	( -1 when no HDC association)
  TRKF(35:37,j) = (r, theta, phi ) of HDC cluster
  TRKF(38,j)    = Energy of Associated HDC cluster

  TRKD(1:15,j) = Upper half of the error matrix of the track parameter

  VTX() !  VTX hit information matched with CDC track
  NVTX(j)      = Number of VTX Hit associated to the track.
  VTXD(1,k,j)  = r of k-th VTX hit (cm)
  VTXD(2,k,j)  = azimuthal angle (radian) of k-th VTX hit
  VTXD(3,k,j)  = Z of k-th VTX hit


  EMC(), HDC() : EMC and HDC array is those without CDC track association.

  EMC(1,j)    = Number of counters containing 90% of EMC energy 
	( -1 when no EMC association)
  EMC(2:4,j) = (r, theta, phi ) of EMC cluster
  EMC(5,j)    = Energy of Associated EMC cluster

  HDC(1,j)    = Number of counters containing 90% of HDC energy 
	( -1 when no HDC association)
  HDC(2:4,j) = (r, theta, phi ) of HDC cluster
  HDC(5,j)    = Energy of Associated HDC cluster


About Helix Parameterization please see,
Postscript version of this document, http://acfahep.kek.jp/subg/sim/lib/simdst-1-00.ps
or "Introduction to Helical Track Manupilation" 
at http://www-jlc.kek.jp/subg/offl/lib/docs/helix_manip.ps.gz 
or http://www-jlc.kek.jp/subg/offl/lib/docs/helix_manip/main.html