Skip to Content.
Sympa Menu

opal - [Opal] Opal-T Energy Issue

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

[Opal] Opal-T Energy Issue


Chronological Thread 
  • From: Ryan Roussel <roussel AT g.ucla.edu>
  • To: opal AT lists.psi.ch
  • Subject: [Opal] Opal-T Energy Issue
  • Date: Thu, 15 Mar 2018 00:09:59 +0000

Hello all,
I am trying to understand how OPAL handles particle energy and momentum. I've seen two different methods for specifying a beam distribution with a non-zero initial energy. One can specify the central momentum/energy/gamma in the BEAM command or I have seen examples where a non-zero offset can be specified in the DISTRIBUTION command. 

I ran the same input file using both of these methods and I have been unable to interpret the results. In the case where I specify the energy only in the BEAM command I see the energy correctly specified in the BEAM portion of the OPAL-T output. But when I look at the final bunch dump BUNCH region, the kinetic energy is not correct. Furthermore it takes orders of magnitude more time to travel 0.1 m in a simple drift than it should.

Conversely, when the energy is specified in both the DISTRIBUTION command and the BEAM command the kinetic energy in the BUNCH dumps is correct and the amount of time to travel 0.1 m seems correct but the "spos" and reference trajectory are wrong.

Part of my confusion could come from units issues, where the energy and momentum specifications are in eV,Mev,or GeV in both the program and documentation. As a new user, this is difficult to work with and would be a request of mine to improve the code unless there is a major reason for this.

Regardless, this is something that I could use help with. I have attached the input file and the outputs from each run that I discussed before. The beamline is a simple drift with a screen.

Thanks for you time,
Ryan Roussel
Ippl> CommMPI: Parent process waiting for children ...
Ippl> CommMPI: Initialization complete.
> ____ _____ ___
> / __ \| __ \ /\ | |
> | | | | |__) / \ | |
> | | | | ___/ /\ \ | |
> | |__| | | / ____ \| |____
> \____/|_| /_/ \_\______|
OPAL>
OPAL> This is OPAL (Object Oriented Parallel Accelerator Library) Version
1.9.0
OPAL>
OPAL> (c) PSI, http://amas.web.psi.ch
OPAL> The optimiser (former opt-Pilot) is integrated
OPAL>
OPAL> Please send cookies, goodies or other motivations (wine and beer ... )
OPAL> to the OPAL developers opal AT lists.psi.ch
OPAL>
OPAL> Time: 15:22:44 date: 14/03/2018
OPAL>
OPAL> Couldn't find startup file "/home/roussel/init.opal".
OPAL> Note: this is not mandatory for an OPAL simulation!
OPAL>
OPAL> * Reading input stream "quickpic_test.tmpl".
OPAL>
OPAL> value: {EDES,P0}={60,60.0005}
OPAL>
OPAL> *
**********************************************************************************

OPAL> * Selected Tracking Method == PARALLEL-T, NEW TRACK
OPAL> *
**********************************************************************************

OPAL> *
**********************************************************************************
OPAL> * Write initial distribution to file "data/quickpic_test_DIST.dat"
OPAL> *
**********************************************************************************
OPAL> * ************* D I S T R I B U T I O N
********************************************
OPAL> *
OPAL> * Number of particles: 1000
OPAL> *
OPAL> * Distribution type: GAUSS
OPAL> *
OPAL> * Sigma Time Rise = 7.000000e-12 [sec]
OPAL> * TPULSEFWHM = 0.000000e+00 [sec]
OPAL> * Sigma Time Fall = 7.000000e-12 [sec]
OPAL> * Longitudinal cutoff = 3.000000e+00 [units of Sigma Time]
OPAL> * Flat top modulation amplitude = 0.000000e+00 [Percent of distribution
amplitude]
OPAL> * Flat top modulation periods = 0.000000e+00
OPAL> * SIGMAX = 7.000000e-05 [m]
OPAL> * SIGMAY = 7.000000e-05 [m]
OPAL> * SIGMAPX = 1.000000e-06 [Beta Gamma]
OPAL> * SIGMAPY = 1.000000e-06 [Beta Gamma]
OPAL> * CORRX = 0.000000e+00
OPAL> * CORRY = 0.000000e+00
OPAL> * CUTOFFX = 3.000000e+00 [units of SIGMAX]
OPAL> * CUTOFFY = 3.000000e+00 [units of SIGMAY]
OPAL> * CUTOFFPX = 3.000000e+00 [units of SIGMAPX]
OPAL> * CUTOFFPY = 3.000000e+00 [units of SIGMAPY]
OPAL> *
OPAL> * Number of energy bins = 1
OPAL> * Distribution is emitted.
OPAL> * Emission time = 2.551626e-11 [sec]
OPAL> * Time per bin = 2.551626e-11 [sec]
OPAL> * Delta t during emission = 2.551626e-11 [sec]
OPAL> *
OPAL> * ------------- THERMAL EMITTANCE MODEL
--------------------------------------------
OPAL> * THERMAL EMITTANCE in NONE MODE
OPAL> * Kinetic energy added to longitudinal momentum = 1.000000e+00 [eV]
OPAL> *
----------------------------------------------------------------------------------
OPAL> *
OPAL> *
OPAL> *
**********************************************************************************
OPAL>
OPAL> * ************* B E A M
************************************************************
OPAL> * BEAM BEAM1
OPAL> * PARTICLE ELECTRON
OPAL> * CURRENT 1.300000e+06 A
OPAL> * FREQUENCY 1.300000e+09 MHz
OPAL> * CHARGE -e * 1.000000e+00
OPAL> * REST MASS 5.109989e-04 GeV
OPAL> * MOMENTUM 6.000051e+01
OPAL> * NPART 1.000000e+03
OPAL> *
**********************************************************************************

OPAL>
OPAL> * ************* F I E L D S O L V E R
**********************************************
OPAL> * FIELDSOLVER FS_SC
OPAL> * TYPE NONE
OPAL> * N-PROCESSORS 1
OPAL> * MX 1.600000e+01
OPAL> * MY 1.600000e+01
OPAL> * MT 3.200000e+01
OPAL> * BBOXINCR 1.000000e+00
OPAL> * XDIM serial
OPAL> * YDIM serial
OPAL> * Z(T)DIM parallel
OPAL> *
**********************************************************************************

OPAL>
OPAL>
OPAL> Size of map 3 sections
OPAL> Key: ( 0.000000- 0.095014) number of overlapping elements 0
OPAL> Key: ( 0.095014- 0.104978) number of overlapping elements 1
OPAL> YAG1
OPAL> Key: ( 0.104978- 0.300047) number of overlapping elements 0
OPAL>
OPAL>
OPAL> Track start at: 15:22:44, t= -12.758 [ps]; zstart at: 0.000 [um]
ParallelTTracker > Change dT from 1e-10 [s] to 2.55163e-11 [s] during
emission
OPAL> Executing ParallelTTracker, initial dt= 100.000 [ps];
OPAL> max integration steps 100000, next step= 0
ParallelTTracker > 15:22:45 Step 999 at 59.302 [mm], t= 99.987 [ns],
E=1.000 [eV]
ParallelTTracker > 15:22:45 Step 1999 at 118.612 [mm], t= 199.987 [ns],
E=1.000 [eV]
ParallelTTracker > 15:22:45 Step 2999 at 177.921 [mm], t= 299.987 [ns],
E=1.000 [eV]
ParallelTTracker > 15:22:45 Step 3999 at 237.231 [mm], t= 399.987 [ns],
E=1.000 [eV]
ParallelTTracker > 15:22:45 Step 4999 at 296.540 [mm], t= 499.987 [ns],
E=1.000 [eV]
OPAL >
OPAL > * ************** B U N C H
*********************************************************
OPAL > * NP = 1000
OPAL > * Qtot = 1.000 [nC] Qi = 1000.000 [fC]
OPAL > * Ekin = 1.000 [eV] dEkin = 0.000 [eV]
OPAL > * rmax = ( 512.58152 , 554.88105 , -26.88487 ) [um]
OPAL > * rmin = ( -476.49934 , -444.88622 , -50.86824 ) [um]
OPAL > * rms beam size = ( 156.96350 , 159.99970 , 4.20843 ) [um]
OPAL > * rms momenta = ( 9.45226e-07 , 9.63083e-07 , 0.00000e+00 )
[beta gamma]
OPAL > * mean position = ( 9.21950 , -1.37665 , -39.33849 ) [um]
OPAL > * mean momenta = ( 4.12256e-08 , 1.49518e-08 , 1.97836e-03 )
[beta gamma]
OPAL > * rms emittance = ( 3.19397e-08 , 3.26566e-08 , 0.00000e+00 )
(not normalized)
OPAL > * rms correlation = ( 9.04773e-01 , 9.07862e-01 , 0.00000e+00 )
OPAL > * hr = ( -1.00000 , -1.00000 , -1.00000 ) [m]
OPAL > * dh = 1.00000e-10 [%]
OPAL > * t = 505.787 [ns] dT = 100.000 [ps]
OPAL > * spos = 300.000 [mm]
OPAL > *
**********************************************************************************

OPAL >
OPAL> done executing ParallelTTracker at 15:22:45
OPAL> * End of input stream "quickpic_test.tmpl".
Timings> -----------------------------------------------------------------
Timings> Timing results for 1 nodes:
Timings> -----------------------------------------------------------------
Timings> mainTimer........... Wall tot = 0.93762, CPU tot = 0.83
Timings>
Timings> Binaryrepart........ Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Boundingbox......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Boundingbox-bounds.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Boundingbox-update.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Compute Statistics.. Wall max = 0.02055, CPU max = 0
Timings> Wall avg = 0.02055, CPU avg = 0
Timings> Wall min = 0.02055, CPU min = 0
Timings>
Timings> Create Distr........ Wall max = 0.002984, CPU max = 0
Timings> Wall avg = 0.002984, CPU avg = 0
Timings> Wall min = 0.002984, CPU min = 0
Timings>
Timings> External field eval. Wall max = 0.066062, CPU max = 0.07
Timings> Wall avg = 0.066062, CPU avg = 0.07
Timings> Wall min = 0.066062, CPU min = 0.07
Timings>
Timings> Fast inside test.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Histogram........... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Initialize geometry. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Inside test......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Load Distr.......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Particle Inside..... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Ray tracing......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Secondary emission.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SelfField total..... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> TIntegration1....... Wall max = 0.081502, CPU max = 0.11
Timings> Wall avg = 0.081502, CPU avg = 0.11
Timings> Wall min = 0.081502, CPU min = 0.11
Timings>
Timings> TIntegration2....... Wall max = 0.250465, CPU max = 0.18
Timings> Wall avg = 0.250465, CPU avg = 0.18
Timings> Wall min = 0.250465, CPU min = 0.18
Timings>
Timings> WakeField........... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Write H5-File....... Wall max = 0.031713, CPU max = 0.01
Timings> Wall avg = 0.031713, CPU avg = 0.01
Timings> Wall min = 0.031713, CPU min = 0.01
Timings>
Timings> Write Stat.......... Wall max = 0.077891, CPU max = 0.06
Timings> Wall avg = 0.077891, CPU avg = 0.06
Timings> Wall min = 0.077891, CPU min = 0.06
Timings>
Timings> -----------------------------------------------------------------
Ippl> CommMPI: Parent process waiting for children ...
Ippl> CommMPI: Initialization complete.
> ____ _____ ___
> / __ \| __ \ /\ | |
> | | | | |__) / \ | |
> | | | | ___/ /\ \ | |
> | |__| | | / ____ \| |____
> \____/|_| /_/ \_\______|
OPAL>
OPAL> This is OPAL (Object Oriented Parallel Accelerator Library) Version
1.9.0
OPAL>
OPAL> (c) PSI, http://amas.web.psi.ch
OPAL> The optimiser (former opt-Pilot) is integrated
OPAL>
OPAL> Please send cookies, goodies or other motivations (wine and beer ... )
OPAL> to the OPAL developers opal AT lists.psi.ch
OPAL>
OPAL> Time: 15:24:22 date: 14/03/2018
OPAL>
OPAL> Couldn't find startup file "/home/roussel/init.opal".
OPAL> Note: this is not mandatory for an OPAL simulation!
OPAL>
OPAL> * Reading input stream "quickpic_test.tmpl".
OPAL>
OPAL> value: {EDES,P0}={60,60.0005}
OPAL>
OPAL> *
**********************************************************************************

OPAL> * Selected Tracking Method == PARALLEL-T, NEW TRACK
OPAL> *
**********************************************************************************

OPAL> *
**********************************************************************************
OPAL> * Write initial distribution to file "data/quickpic_test_DIST.dat"
OPAL> *
**********************************************************************************
OPAL> * ************* D I S T R I B U T I O N
********************************************
OPAL> *
OPAL> * Number of particles: 1000
OPAL> *
OPAL> * Distribution type: GAUSS
OPAL> *
OPAL> * Sigma Time Rise = 7.000000e-12 [sec]
OPAL> * TPULSEFWHM = 0.000000e+00 [sec]
OPAL> * Sigma Time Fall = 7.000000e-12 [sec]
OPAL> * Longitudinal cutoff = 3.000000e+00 [units of Sigma Time]
OPAL> * Flat top modulation amplitude = 0.000000e+00 [Percent of distribution
amplitude]
OPAL> * Flat top modulation periods = 0.000000e+00
OPAL> * SIGMAX = 7.000000e-05 [m]
OPAL> * SIGMAY = 7.000000e-05 [m]
OPAL> * SIGMAPX = 1.000000e-06 [Beta Gamma]
OPAL> * SIGMAPY = 1.000000e-06 [Beta Gamma]
OPAL> * CORRX = 0.000000e+00
OPAL> * CORRY = 0.000000e+00
OPAL> * CUTOFFX = 3.000000e+00 [units of SIGMAX]
OPAL> * CUTOFFY = 3.000000e+00 [units of SIGMAY]
OPAL> * CUTOFFPX = 3.000000e+00 [units of SIGMAPX]
OPAL> * CUTOFFPY = 3.000000e+00 [units of SIGMAPY]
OPAL> *
OPAL> * Number of energy bins = 1
OPAL> * Distribution is emitted.
OPAL> * Emission time = 2.551626e-11 [sec]
OPAL> * Time per bin = 2.551626e-11 [sec]
OPAL> * Delta t during emission = 2.551626e-11 [sec]
OPAL> *
OPAL> * ------------- THERMAL EMITTANCE MODEL
--------------------------------------------
OPAL> * THERMAL EMITTANCE in NONE MODE
OPAL> * Kinetic energy added to longitudinal momentum = 1.000000e+00 [eV]
OPAL> *
----------------------------------------------------------------------------------
OPAL> *
OPAL> *
OPAL> *
**********************************************************************************
OPAL>
OPAL> * ************* B E A M
************************************************************
OPAL> * BEAM BEAM1
OPAL> * PARTICLE ELECTRON
OPAL> * CURRENT 1.300000e+06 A
OPAL> * FREQUENCY 1.300000e+09 MHz
OPAL> * CHARGE -e * 1.000000e+00
OPAL> * REST MASS 5.109989e-04 GeV
OPAL> * MOMENTUM 6.000051e+01
OPAL> * NPART 1.000000e+03
OPAL> *
**********************************************************************************

OPAL>
OPAL> * ************* F I E L D S O L V E R
**********************************************
OPAL> * FIELDSOLVER FS_SC
OPAL> * TYPE NONE
OPAL> * N-PROCESSORS 1
OPAL> * MX 1.600000e+01
OPAL> * MY 1.600000e+01
OPAL> * MT 3.200000e+01
OPAL> * BBOXINCR 1.000000e+00
OPAL> * XDIM serial
OPAL> * YDIM serial
OPAL> * Z(T)DIM parallel
OPAL> *
**********************************************************************************

OPAL>
OPAL>
OPAL> Size of map 3 sections
OPAL> Key: ( 0.000000- 0.095014) number of overlapping elements 0
OPAL> Key: ( 0.095014- 0.104978) number of overlapping elements 1
OPAL> YAG1
OPAL> Key: ( 0.104978- 0.300047) number of overlapping elements 0
OPAL>
OPAL>
OPAL> Track start at: 15:24:22, t= -12.758 [ps]; zstart at: 0.000 [um]
ParallelTTracker > Change dT from 1e-10 [s] to 2.55163e-11 [s] during
emission
OPAL> Executing ParallelTTracker, initial dt= 100.000 [ps];
OPAL> max integration steps 100000, next step= 0
OPAL >
OPAL > * ************** B U N C H
*********************************************************
OPAL > * NP = 1000
OPAL > * Qtot = 1.000 [nC] Qi = 1000.000 [fC]
OPAL > * Ekin = 30.154 [MeV] dEkin = 3.835 [eV]
OPAL > * rmax = ( 0.20371 , 0.18810 , 341.58498 ) [mm]
OPAL > * rmin = ( -0.20523 , -0.19590 , 329.46378 ) [mm]
OPAL > * rms beam size = ( 0.06686 , 0.06709 , 2.12694 ) [mm]
OPAL > * rms momenta = ( 9.45226e-07 , 9.63083e-07 , 7.50633e-06 )
[beta gamma]
OPAL > * mean position = ( 0.00297 , -0.00364 , 335.29091 ) [mm]
OPAL > * mean momenta = ( 4.12256e-08 , 1.49518e-08 , 6.00025e+01 )
[beta gamma]
OPAL > * rms emittance = ( 1.05309e-12 , 1.07673e-12 , 2.66081e-10 )
(not normalized)
OPAL > * rms correlation = ( -1.93590e-02 , -1.13376e-02 , -2.73438e-06 )
OPAL > * hr = ( -1.00000 , -1.00000 , -1.00000 ) [m]
OPAL > * dh = 1.00000e-10 [%]
OPAL > * t = 1.187 [ns] dT = 100.000 [ps]
OPAL > * spos = 704.149 [um]
OPAL > *
**********************************************************************************

OPAL >
OPAL> done executing ParallelTTracker at 15:24:22
OPAL> * End of input stream "quickpic_test.tmpl".
Timings> -----------------------------------------------------------------
Timings> Timing results for 1 nodes:
Timings> -----------------------------------------------------------------
Timings> mainTimer........... Wall tot = 0.183193, CPU tot = 0.13
Timings>
Timings> Binaryrepart........ Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Boundingbox......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Boundingbox-bounds.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Boundingbox-update.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Compute Statistics.. Wall max = 0.00014, CPU max = 0
Timings> Wall avg = 0.00014, CPU avg = 0
Timings> Wall min = 0.00014, CPU min = 0
Timings>
Timings> Create Distr........ Wall max = 0.00324, CPU max = 0
Timings> Wall avg = 0.00324, CPU avg = 0
Timings> Wall min = 0.00324, CPU min = 0
Timings>
Timings> External field eval. Wall max = 0.000196, CPU max = 0
Timings> Wall avg = 0.000196, CPU avg = 0
Timings> Wall min = 0.000196, CPU min = 0
Timings>
Timings> Fast inside test.... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Histogram........... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Initialize geometry. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Inside test......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Load Distr.......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Particle Inside..... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Ray tracing......... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Secondary emission.. Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> SelfField total..... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> TIntegration1....... Wall max = 0.000185, CPU max = 0
Timings> Wall avg = 0.000185, CPU avg = 0
Timings> Wall min = 0.000185, CPU min = 0
Timings>
Timings> TIntegration2....... Wall max = 0.000692, CPU max = 0
Timings> Wall avg = 0.000692, CPU avg = 0
Timings> Wall min = 0.000692, CPU min = 0
Timings>
Timings> WakeField........... Wall max = 0, CPU max = 0
Timings> Wall avg = 0, CPU avg = 0
Timings> Wall min = 0, CPU min = 0
Timings>
Timings> Write H5-File....... Wall max = 0.031428, CPU max = 0
Timings> Wall avg = 0.031428, CPU avg = 0
Timings> Wall min = 0.031428, CPU min = 0
Timings>
Timings> Write Stat.......... Wall max = 0.003441, CPU max = 0
Timings> Wall avg = 0.003441, CPU avg = 0
Timings> Wall min = 0.003441, CPU min = 0
Timings>
Timings> -----------------------------------------------------------------

Attachment: quickpic_test.tmpl
Description: Binary data




Archive powered by MHonArc 2.6.19.

Top of Page