Code indexing in gitaly is broken and leads to code not being visible to the user. We work on the issue with highest priority.

Skip to content
Snippets Groups Projects
Commit 4303afeb authored by kraus's avatar kraus
Browse files

have to set P0 to correct energy

parent 8871e3d4
No related branches found
No related tags found
No related merge requests found
......@@ -91,13 +91,18 @@ REAL global_shift = 1.0;
//================================================
// drift before bend (m).
REAL drift_before_bend = 0.1;
REAL drift_before_bend = 0.2;
// Bend angle (radians).
REAL bend_angle = 0.5235987755982988;
// Bend design energy (MeV).
REAL bend_energy = 7.0;
REAL bend_energy = 6.50762633;
REAL gamma = (bend_energy * 1e-3 + EMASS)/EMASS;
REAL beta = sqrt(1-(1/gamma^2));
REAL gambet = gamma*beta;
REAL P0 = gamma*beta*EMASS;
// Bend length (m).
REAL bend_length = 0.129409522551;
......@@ -186,7 +191,7 @@ if (CSR_ON > 0) {
L = bend_length,
GAP = 0.02,
WAKEF = FS_CSR_WAKE;
drift2:DRIFT, L=0.4, ELEMEDGE=bend_length + drift_before_bend, WAKEF = FS_CSR_WAKE;
drift2:DRIFT, L=0.4, ELEMEDGE=0.4 + bend_length + drift_before_bend, WAKEF = FS_CSR_WAKE;
} else {
bend1: SBend, ANGLE = bend_angle,
E1 = 0.08726646259971647,
......@@ -196,7 +201,7 @@ if (CSR_ON > 0) {
DESIGNENERGY = bend_energy,
L = bend_length,
GAP = 0.02;
drift2:DRIFT, L=0.4, ELEMEDGE=bend_length + drift_before_bend;
drift2:DRIFT, L=0.4, ELEMEDGE=0.4 + bend_length + drift_before_bend;
}
//******************************************************************************************************************
......
  • kraus @kraus

    mentioned in issue #61 (closed)

    ·

    mentioned in issue #61 (closed)

    Toggle commit list
  • Author Developer

    @adelmann I've shifted the drift in the input file. Please try to run the simulation after undoing this.

  • Author Developer

    @adelmann Please have a look at the changes in the lines 194 and 204. I've added 0.4 m. I've used the same input file with both 1.6 (only the necessary adaptions) and 1.9. Neither showed any effect in the drift, see posted results.

  • adelmann :reminder_ribbon: @adelmann ·
    Owner

    YES need to set ELEMEDGE to 0.225 then the reults agree :thumbsup:

    or remove 0.4 i.e. ELEMEDGE=bend_length + drift_before_bend;

    Edited by adelmann
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment