Instead of zeroing fcyc outside the cavity, zero the impulse amplitud…

…e. So the phase advancement is preserved, but without driving the cavity ring-up.
project8 · Jul 25, 2023 · 938eea2 · 938eea2
1 parent 9c795d6
commit 938eea2
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/Source/Kassiopeia/LMCFieldCalculator.cc b/Source/Kassiopeia/LMCFieldCalculator.cc
@@ -397,11 +397,11 @@ namespace locust
     	double zLocation = tKassParticleXP[2];
     	double dimZ = fInterface->fField->GetDimL();
 		double orbitPhase = tKassParticleXP[6];  // radians
-
-		double cycFrequency = 0.;
+		double cycFrequency = tKassParticleXP[7];
+		double amplitude = 0.;
 		if (fabs(zLocation) < dimZ/2.) // Check whether the electron is inside the cavity length.
 		{
-			cycFrequency = tKassParticleXP[7];  // rad/s
+			amplitude = 1.;
 		}
 
     	if ( !BypassTF )
@@ -421,7 +421,7 @@ namespace locust
 
     			if (*it != 0.)
     			{
-    				fFIRBuffer.push_back(cos(orbitPhase));
+    				fFIRBuffer.push_back(amplitude * cos(orbitPhase));
     			}
     			else
     			{