From 9f46f9dc80f339af7fc5f5dbb8bd5561d08dd08b Mon Sep 17 00:00:00 2001
From: qiuchan <qiuchan.lu@cern.ch>
Date: Wed, 2 Aug 2023 11:34:10 +0200
Subject: [PATCH 1/2] option files for Ar injected in Run3

---
 ...am2680GeV-0GeV-md100-expected-2023-0.37.py | 32 +++++++++++++++++
 ...> pAr-Beam6800GeV-0GeV-md100-2022-0.37.py} |  8 ++---
 ...am6800GeV-0GeV-md100-expected-2023-0.37.py | 34 +++++++++++++++++++
 3 files changed, 70 insertions(+), 4 deletions(-)
 create mode 100644 options/Gauss/PbAr-Beam2680GeV-0GeV-md100-expected-2023-0.37.py
 rename options/Gauss/{pArBeam6800GeV-0GeV-md100-2022-0.37.py => pAr-Beam6800GeV-0GeV-md100-2022-0.37.py} (75%)
 create mode 100644 options/Gauss/pAr-Beam6800GeV-0GeV-md100-expected-2023-0.37.py

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 00000000..58b91d45
--- /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/pArBeam6800GeV-0GeV-md100-2022-0.37.py b/options/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py
similarity index 75%
rename from options/Gauss/pArBeam6800GeV-0GeV-md100-2022-0.37.py
rename to options/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py
index 13fd3544..fb00c234 100644
--- a/options/Gauss/pArBeam6800GeV-0GeV-md100-2022-0.37.py
+++ b/options/Gauss/pAr-Beam6800GeV-0GeV-md100-2022-0.37.py
@@ -5,12 +5,12 @@
 #
 #  Syntax is: 
 #  gaudirun.py  $APPCONFIGOPTS/Gauss/Beam6800GeV-md100-2022-nu2.1.py
-#               $APPCONFIGOPTS/Gauss/pArBeam6800GeV-0GeV-md100-2022-0.37.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)
+#		$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
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 00000000..cf9ba4c9
--- /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
+
+
-- 
GitLab


From 445586d8fec9e7f89680111dda782288d104ae70 Mon Sep 17 00:00:00 2001
From: qiuchan <qiuchan.lu@cern.ch>
Date: Fri, 4 Aug 2023 14:18:50 +0200
Subject: [PATCH 2/2] keep original beamfile

---
 .../pArBeam6800GeV-0GeV-md100-2022-0.37.py    | 34 +++++++++++++++++++
 1 file changed, 34 insertions(+)
 create mode 100644 options/Gauss/pArBeam6800GeV-0GeV-md100-2022-0.37.py

diff --git a/options/Gauss/pArBeam6800GeV-0GeV-md100-2022-0.37.py b/options/Gauss/pArBeam6800GeV-0GeV-md100-2022-0.37.py
new file mode 100644
index 00000000..13fd3544
--- /dev/null
+++ b/options/Gauss/pArBeam6800GeV-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/pArBeam6800GeV-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)
+#		$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
+
+
-- 
GitLab