frib-high-level-controls/FLAME

Input file for FRIB (LS1+FS1)

Closed this issue · 1 comments

kryv commented

I'm trying to debug main.cpp with FRIB lattice file.
The input lattice file already exist at FLAME/python/uscsi/test/latticeout_IMP_withPV_consolidate.lat .
And executable compiled file 'main' also requires detailed data from TLM_Lite source codes.
So I have executed "main" with 2 arguments.

$ ./main /${path}/latticeout_IMP_withPV_consolidate.lat /${path}/tlm_lite
This returns
Parsing succeeds
Main exception: NCharge

Accordingly, I add
NCharge = [10111.0, 10531.0];
in latticeout_IMP_withPV_consolidate.lat .

And then, I was able to run the "main" file.

But, the calculated result of the beam parameters (eg. BeamRMS.out etc.) is differ greatly from TLM Java code (tlm_lite/example/runtlm).
Do you have any idea? or Is there any wrong command in above?
(Now I'm checking physics variables in the cpp source code also.)

kryv commented

Add information about the result.

After I add
NCharge = [10111.0, 10531.0]; ,

$ ./main /${path}/latticeout_IMP_withPV_consolidate.lat /${path}/tlm_lite
returns below

Parsing succeeds

Ion charge states:
  0.13866  0.14286

Ion charge amount:
 10111.0 10531.0

Barycenter: 
-----------
omission
-----------

Beam envelope:
-----------
omission
-----------

InitLong:
  IonZ = 0.13866  IonEs [Mev/u] = 931.49432, IonEk [Mev/u] = 0.50000, IonW [Mev/u] = 931.99432

InitLong: 0.01832 sec

InitLattice:
  IonZ = 0.13866  IonEs [Mev/u] = 931.49432, IonEk [Mev/u] = 0.50000, IonW [Mev/u] = 931.99432
  IonZ = 0.14286  IonEs [Mev/u] = 931.49432, IonEk [Mev/u] = 0.50000, IonW [Mev/u] = 931.99432
*** GetTransitFac: beta out of range 0.0327619

and
tail -2 BeamRMS.out

9.855276000000000e+00  7.023517629255475e+08  2.697435992480285e+05  7.325989069535193e+08  2.812918864926295e+05  3.013622197197078e-01  8.764218641711942e-04
9.855276000000000e+00

The main code stopped near s~10.0 [m] and the x RMS reached to 7.0e8 [mm].