OPAL-Cyc - single particle tracking: not store right statistics in .h5 file
Simulation:
- single particle tracking in OPAL-Cyc (1.5.002)
- PSDUMPLOCALFRAME = FALSE
- Input file: /afs/psi.ch/project/SCGantry/CCT_Concept/OPAL_Test
OPAL outputs
-
trackOrbit: IDO in global frame. The particle coordinates are right
The six-dimensional phase space data in the global Cartesian coordinates
Part. ID x [mm] beta_xgamma y [mm] beta_ygamma z [mm] beta_z*gamma
ID0 -1.30000000e+03 0.00000000e+00 1.00000000e-01 6.56544376e-01 1.50000000e+02 0.00000000e+00 ID0 -1.30000000e+03 0.00000000e+00 3.39069250e+00 6.56544376e-01 1.50000000e+02 0.00000000e+00 ID0 -1.30000000e+03 0.00000000e+00 6.68138500e+00 6.56544376e-01 1.50000000e+02 0.00000000e+00 ID0 -1.30000000e+03 0.00000000e+00 9.97207750e+00 6.56544376e-01 1.50000000e+02 0.00000000e+00 ID0 -1.30000000e+03 0.00000000e+00 1.32627700e+01 6.56544376e-01 1.50000000e+02 0.00000000e+00 ID0 -1.30000000e+03 0.00000000e+00 1.65534625e+01 6.56544376e-01 1.50000000e+02 0.00000000e+00
* **h5**: the particle coordinates seem to be not stored. They are always zero. Here an example for different steps (0, 1, 10, 100, 120)
root [3] TNtupleD *opalC = dataC.GetNtuple(0); root [4] opalC->Scan()
- Row * x * y * z * px * py * pz * phiz * enez *
-
0 * 1.63e-322 * 6.90e-310 * 6.90e-310 * 0 * 0.6565443 * 0 * 0 * 0 *
root [5] TNtupleD *opalC = dataC.GetNtuple(1); root [6] opalC->Scan()
- Row * x * y * z * px * py * pz * phiz * enez *
-
0 * 1.63e-322 * 6.90e-310 * 6.90e-310 * 0 * 0.6565443 * 0 * 0 * 0 *
root [7] TNtupleD *opalC = dataC.GetNtuple(10); root [8] opalC->Scan()
- Row * x * y * z * px * py * pz * phiz * enez *
-
0 * 1.63e-322 * 6.90e-310 * 6.90e-310 * 0 * 0.6565443 * 0 * 0 * 0 *
root [9] TNtupleD *opalC = dataC.GetNtuple(100); root [10] opalC->Scan()
- Row * x * y * z * px * py * pz * phiz * enez *
-
0 * 1.63e-322 * 6.90e-310 * 6.90e-310 * 0.6408104 * -0.142872 * 0 * 0 * 0 *