diff --git a/options/Gauss/PbAr-Beam2680GeV-0GeV-md100-expected-2023-0.37.py b/options/Gauss/PbAr-Beam2680GeV-0GeV-md100-expected-2023-0.37.py
new file mode 100644
index 0000000000000000000000000000000000000000..58b91d45f01c5b79bae2dd42b23dd888d050160d
--- /dev/null
+++ b/options/Gauss/PbAr-Beam2680GeV-0GeV-md100-expected-2023-0.37.py
@@ -0,0 +1,32 @@
+# Beam file for 2023 PbAr, specifying the FixTargetParticle and nu for PbAr case,
+# other beam parameters are the same as PbPb-Beam2680GeV-md100-expected-2023-fix1.py.
+# Thus, this beam file should be used together with PbPb-Beam2680GeV-md100-expected-2023-fix1.py,
+# And requires Gauss v56r0 or higher.
+
+# Syntax is:
+# gaudirun.py  $APPCONFIGOPTS/Gauss/PbPb-Beam2680GeV-md100-expected-2023-fix1.py
+#              $APPCONFIGOPTS/Gauss/PbAr-Beam2680GeV-0GeV-md100-expected-2023-0.37.py
+#              $GAUSSOPTS/BeforeVeloGeometry.py
+#              $APPCONFIGOPTS/Gauss/Smog2TriangularProfile.py
+#              $GAUSSOPTS/BeamGas.py #(only PbAr; if PbAr+PbPb, change to: $GAUSSOPTS/BeamGasWithBeamBeam.py)
+#              $LBCRMCROOT/options/EPOSBeamGas.py  #(EPOS for BeamGas)
+#              $DECFILESROOT/options/30000000.py #(i.e. event type, if embeding, add $GAUSSOPTS/Embedding.py)
+#              $APPCONFIGOPTS/Gauss/Run3-detector.py;
+#              $APPCONFIGOPTS/Gauss/DataType-2023.py;
+#              $APPCONFIGOPTS/Gauss/G4PL_FTFP_BERT_EmOpt2.py
+#              $GAUSSOPTS/GenStandAlone.py
+
+
+from Gauss.Configuration import *
+
+# Specify beam particle
+Gauss().B2Momentum = -2680*SystemOfUnits.GeV
+Gauss().B1Particle = 'Pb'
+Gauss().B2Particle = 'Pb'
+Gauss().FixedTargetParticle = 'Ar'
+
+#  assuming L(beam-gas)=1.8 x 10^31 cm-2 s-1 with 2600 beam1 bunches
+#  corresponds to nu(beam-gas)=0.37
+Gauss().FixedTargetLuminosity = 7*(10**27)/(SystemOfUnits.cm2*SystemOfUnits.s)
+Gauss().FixedTargetXSection = 600.*SystemOfUnits.millibarn
+
diff --git a/options/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py b/options/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py
new file mode 100644
index 0000000000000000000000000000000000000000..fb00c2349a5440e3b1da43017cc2d2fe2374486c
--- /dev/null
+++ b/options/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py
@@ -0,0 +1,34 @@
+# Beam file for 2022 pAr, specifying the FixTargetParticle and nu for pAr case, 
+# other beam parameters are the same as Beam6800GeV-md100-2022-nu2.1.py.
+# Thus, this beam file should be used together with Beam6800GeV-md100-2022-nu2.1.py,
+# And requires Gauss v56r0 or higher.
+#
+#  Syntax is: 
+#  gaudirun.py  $APPCONFIGOPTS/Gauss/Beam6800GeV-md100-2022-nu2.1.py
+#               $APPCONFIGOPTS/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py
+#               $GAUSSOPTS/BeforeVeloGeometry.py
+#       	$APPCONFIGOPTS/Gauss/Smog2TriangularProfile.py
+#		$GAUSSOPTS/BeamGas.py #(only pAr; if pAr+pp, change to: $GAUSSOPTS/BeamGasWithBeamBeam.py)
+#		$LBCRMCROOT/options/EPOSBeamGas.py  #(EPOS for BeamGas)
+#               $DECFILESROOT/options/30000000.py #(i.e. event type, if embeding, add $GAUSSOPTS/Embedding.py)
+#		$APPCONFIGOPTS/Gauss/Run3-detector-NoUT.py;
+#		$APPCONFIGOPTS/Gauss/DataType-2022.py;
+#		$APPCONFIGOPTS/Gauss/G4PL_FTFP_BERT_EmOpt2.py
+#		$GAUSSOPTS/GenStandAlone.py
+
+
+from Gauss.Configuration import *
+
+# Specify beam particle
+Gauss().B2Momentum = -6800*SystemOfUnits.GeV
+Gauss().B1Particle = 'p'
+Gauss().B2Particle = 'p'
+Gauss().FixedTargetParticle = 'Ar'
+
+
+#  assuming L(beam-gas)=1.8 x 10^31 cm-2 s-1 with 2600 beam1 bunches
+#  corresponds to nu(beam-gas)=0.37
+Gauss().FixedTargetLuminosity = 7*(10**27)/(SystemOfUnits.cm2*SystemOfUnits.s)
+Gauss().FixedTargetXSection = 600.*SystemOfUnits.millibarn
+
+
diff --git a/options/Gauss/pAr-Beam6800GeV-0GeV-md100-expected-2023-0.37.py b/options/Gauss/pAr-Beam6800GeV-0GeV-md100-expected-2023-0.37.py
new file mode 100644
index 0000000000000000000000000000000000000000..cf9ba4c9b6b850a8f16d5d43cd3c7fa3fcf5bb23
--- /dev/null
+++ b/options/Gauss/pAr-Beam6800GeV-0GeV-md100-expected-2023-0.37.py
@@ -0,0 +1,34 @@
+# Beam file for expected 2023 pAr, specifying the FixTargetParticle and nu for pAr case, taken from 2022 pAr beamfile.
+# other beam parameters are the same as Beam6800GeV-md100-expected-2023-nu1.4.py.
+# Thus, this beam file should be used together with Beam6800GeV-md100-expected-2023-nu1.4.py,
+# And requires Gauss v56r0 or higher.
+#
+#  Syntax is: 
+#  gaudirun.py  $APPCONFIGOPTS/Gauss/Beam6800GeV-md100-expected-2023-nu1.4.py
+#               $APPCONFIGOPTS/Gauss/pAr-Beam6800GeV-0GeV-md100-expected-2023-0.37.py
+#               $GAUSSOPTS/BeforeVeloGeometry.py
+#       	$APPCONFIGOPTS/Gauss/Smog2TriangularProfile.py
+#		$GAUSSOPTS/BeamGas.py #(only pAr; if pAr+pp, change to: $GAUSSOPTS/BeamGasWithBeamBeam.py)
+#		$LBCRMCROOT/options/EPOSBeamGas.py  #(EPOS for BeamGas)
+#               $DECFILESROOT/options/30000000.py #(i.e. event type, if embeding, add $GAUSSOPTS/Embedding.py)
+#		$APPCONFIGOPTS/Gauss/Run3-detector.py;
+#		$APPCONFIGOPTS/Gauss/DataType-2023.py;
+#		$APPCONFIGOPTS/Gauss/G4PL_FTFP_BERT_EmOpt2.py
+#		$GAUSSOPTS/GenStandAlone.py
+
+
+from Gauss.Configuration import *
+
+# Specify beam particle
+Gauss().B2Momentum = -6800*SystemOfUnits.GeV
+Gauss().B1Particle = 'p'
+Gauss().B2Particle = 'p'
+Gauss().FixedTargetParticle = 'Ar'
+
+
+#  assuming L(beam-gas)=1.8 x 10^31 cm-2 s-1 with 2600 beam1 bunches
+#  corresponds to nu(beam-gas)=0.37
+Gauss().FixedTargetLuminosity = 7*(10**27)/(SystemOfUnits.cm2*SystemOfUnits.s)
+Gauss().FixedTargetXSection = 600.*SystemOfUnits.millibarn
+
+