setup for impacting muons

Author: philippwindischhofer

examples/dipole_ice/dipole_ice.cxx

@@ -58,11 +58,11 @@ int main(int argc, char* argv[]) {
 
   // CylinderGeometry geom(20, -15, 15, eps);
   // InfEDipoleAntenna dipole(0.0, 10.0, 0.0, impulse_response);
-  // scalar_t t_end = 25;
+  // scalar_t t_end = 250;
 
   CylinderGeometry geom(150, -150, 150, eps);
   InfEDipoleAntenna dipole(0.0, 10.0, -30.0, impulse_response);
-  scalar_t t_end = 150;
+  scalar_t t_end = 200;
 
   CylindricalWeightingFieldCalculator wfc(geom, dipole, t_end);
   wfc.Calculate(wf_path, "/scratch/midway3/windischhofer/eisvogel/");

examples/point_charge_cherenkov/point_charge_cherenkov.cxx

@@ -18,22 +18,23 @@ int main(int argc, char* argv[]) {
   // test trajectory: a point charge moving parallel to the z-axis 
   // with a constant impact parameter of 'b' along the x-axis
   scalar_t b = 50;
-  scalar_t tstart = -100, tend = 100;
+  scalar_t tstart = -120, tend = -1;
   scalar_t charge = 1;
-  scalar_t beta = 1.2;
+  scalar_t beta = 0.99;
 
   std::cout << "Building trajectory ..." << std::endl;
   Current0D track({
-      CoordUtils::MakeCoordVectorTXYZ(tstart, b, 0, beta * tstart),
-      CoordUtils::MakeCoordVectorTXYZ(tend, b, 0, beta * tend)
+      CoordUtils::MakeCoordVectorTXYZ(tstart, b, 0, -beta * tstart),
+      CoordUtils::MakeCoordVectorTXYZ(tend, b, 0, -beta * tend)
     },
     {charge}
     );
 
   std::cout << "Computing signal ..." << std::endl;
 
   std::vector<scalar_t> signal_times, signal_values;
-  for(scalar_t cur_t = 10; cur_t < 45; cur_t += 0.1) {
+  //for(scalar_t cur_t = 10; cur_t < 45; cur_t += 0.1) {
+  for(scalar_t cur_t = 40; cur_t < 110; cur_t += 0.1) {
     scalar_t cur_signal = calc.ComputeSignal(track, cur_t);
     signal_times.push_back(cur_t);
     signal_values.push_back(cur_signal);