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 + +