macros/REST_Axion_FieldIntegrationTests.C Adding Z-cutoff plot genera…

…tion

macros/REST_Axion_FieldIntegrationTests.C

@@ -18,6 +18,8 @@
 //***
 //*** Author: Javier Galan
 //*******************************************************************************************************
+Bool_t performCutOffTest = false; // To be enabled in case we need to reproduce the Z-cutoff curves published in the RayTracing paper
+
 int REST_Axion_FieldIntegrationTests(Double_t sX = 10, Double_t sY = 10, Double_t sZ = 50, Double_t dm = 0.01,
                                      Double_t tolerance = 0.1, Double_t Ea = 4.2) {
     /// Setting up magnetic field and track to evaluate
@@ -46,5 +48,26 @@ int REST_Axion_FieldIntegrationTests(Double_t sX = 10, Double_t sY = 10, Double_
     std::pair<double, double> prob =
         ax->GammaTransmissionFieldMapProbability(Ea, ma - dm, tolerance, 100, 25);
 
+	if( performCutOffTest )
+	{
+		ax->SetDebug(false);
+		for( Double_t zStart = -10000; zStart < -4000; zStart += 100 )
+		{
+			magneticField.SetUserTrack(TVector3(0, 0, zStart), TVector3(0, 0, -zStart));
+
+			std::pair<double, double> prob = ax->GammaTransmissionFieldMapProbability(Ea, ma - dm, tolerance, 100, 25);
+
+			std::cout << "zStart: " << zStart << " Prob: " << prob.first << std::endl;
+		}
+
+		for( Double_t zStart = -10000; zStart < -4000; zStart += 100 )
+		{
+
+			Double_t field = magneticField.GetTransversalComponent(TVector3(0,0,zStart), TVector3 (0,0,1));
+
+			std::cout << "zStart: " << zStart << " Field: " << field << std::endl;
+		}
+	}
+
     return 0;
 }