Page MenuHomeHEPForge
Files F10881123

No OneTemporary
Actions

View Options

diff --git a/DIPSY/Bin.h b/DIPSY/Bin.h
new file mode 100644
--- /dev/null
+++ b/DIPSY/Bin.h
@@ -0,0 +1,170 @@
+// -*- C++ -*-
+#ifndef ThePEG_Bin_H
+#define ThePEG_Bin_H
+//
+// This is the declaration of the Bin class.
+//
+
+#include "ThePEG/Config/ThePEG.h"
+
+namespace ThePEG {
+
+/**
+ * The Bin class is used to valculate the average of several measurements.
+ */
+template <typename T, typename T2 = T>
+class Bin {
+
+public:
+
+ /** @name Standard constructors and destructors. */
+ //@{
+ /**
+ * The default constructor.
+ */
+ inline Bin(): sum(T()), sum2(T2()), sumw(0.0), n(0) {}
+
+ /**
+ * The constructor fo an initial value.
+ */
+ inline Bin(const T & t): sum(t), sum2(t*t), sumw(0.0), n(0) {}
+
+ /**
+ * The copy constructor fo an initial value.
+ */
+ inline Bin(const Bin<T,T2> & b): sum(b.sum), sum2(b.sum2),
+ sumw(b.sumw), n(b.n) {}
+
+ /**
+ * The assignment operator.
+ */
+ Bin<T,T2> & operator=(const Bin<T,T2> & b) {
+ sum = b.sum;
+ sum2 = b.sum2;
+ sumw = b.sumw;
+ n = b.n;
+ return *this;
+ }
+ //@}
+
+ /**
+ * Add a value
+ */
+ Bin<T,T2> & operator+=(const T & t) {
+ return fill(t, 1.0);
+ }
+
+ Bin<T,T2> & fill(const T & t, double w) {
+ if ( n >= 0 ) {
+ sum += t*w;
+ sum2 += t*t*w;
+ sumw += w;
+ ++n;
+ }
+ return *this;
+ }
+
+ T var() const {
+ return n > 0? sum2/sumw - average()*average(): sum2;
+ }
+
+ T2 err2() const {
+ return n > 0? var()/double(abs(n)): sum2;
+ }
+
+ T err() const {
+ return sqrt(err2());
+ }
+
+ T average() const {
+ return n > 0? sum/sumw: sum;
+ }
+
+ T operator()() const {
+ return average();
+ }
+
+ Bin<T,T2> & operator+=(const Bin<T,T2> & b) {
+ T av = average() + b.average();
+ T2 e2 = err2() + b.err2();
+ sum = av;
+ sum2 = e2;
+ sumw = 1.0;
+ n = -1;
+ return *this;
+ }
+
+ Bin<T,T2> & operator-=(const Bin<T,T2> & b) {
+ T av = average() - b.average();
+ T2 e2 = err2() + b.err2();
+ sum = av;
+ sum2 = e2;
+ sumw = 1.0;
+ n = -1;
+ return *this;
+ }
+
+ Bin<T,T2> & operator*=(const Bin<T,T2> & b) {
+ T av = average()*b.average();
+ T2 e2 = av*av*(err2()/(average()*average) +
+ b.err2()/(b.average()*b.average()));
+ sum = av;
+ sum2 = e2;
+ sumw = 1.0;
+ n = -1;
+ return *this;
+ }
+
+ Bin<T,T2> sqr() const {
+ Bin<T,T2> ret = *this;
+ ret.sum = average()*average();
+ ret.sum2 = 4.0*err2()*ret.sum;
+ ret.sumw = 1.0;
+ ret.n = -1;
+ return ret;
+ }
+
+
+ Bin<T,T2> & operator/=(const Bin<T,T2> & b) {
+ T av = average()/b.average();
+ T2 e2 = av*av*(err2()/(average()*average) +
+ b.err2()/(b.average()*b.average()));
+ sum = av;
+ sum2 = e2;
+ sumw = 1.0;
+ n = -1;
+ return *this;
+ }
+
+ Bin<T,T2> operator+(const Bin<T,T2> & b) const {
+ Bin<T,T2> ret = *this;
+ return ret += b;
+ }
+
+ Bin<T,T2> operator-(const Bin<T,T2> & b) const {
+ Bin<T,T2> ret = *this;
+ return ret -= b;
+ }
+
+ Bin<T,T2> operator*(const Bin<T,T2> & b) const {
+ Bin<T,T2> ret = *this;
+ return ret *= b;
+ }
+
+ Bin<T,T2> operator/(const Bin<T,T2> & b) const {
+ Bin<T,T2> ret = *this;
+ return ret /= b;
+ }
+
+private:
+
+ T sum;
+ T2 sum2;
+ double sumw;
+ long n;
+
+};
+
+}
+
+#endif /* ThePEG_Bin_H */
diff --git a/DIPSY/PAXSec.in b/DIPSY/PAXSec.in
--- a/DIPSY/PAXSec.in
+++ b/DIPSY/PAXSec.in
@@ -1,684 +1,701 @@
cd /DIPSY
## First we setup some tuned parameters
#read CurrentTune.in
read Tune31.in
## Now we set up an event generator. We start with running pp rahter than pA.
cp EventHandler PAEventHandler
set PAEventHandler:WFL stdProton
set PAEventHandler:WFR stdProton
set PAEventHandler:ConsistencyLevel 0
set PAEventHandler:XSecFn:CheckOffShell false
set PAEventHandler:CascadeHandler NULL
set PAEventHandler:HadronizationHandler NULL
set PAEventHandler:DecayHandler NULL
create ThePEG::LuminosityFunction PALumi
set PAEventHandler:LuminosityFunction PALumi
cp Generator PAGenerator
set PAGenerator:EventHandler PAEventHandler
set PAGenerator:NumberOfEvents 0
set PAGenerator:EventHandler:EventFiller:PTCut 0.6
set PAGenerator:Path PAXSec
erase PAGenerator:AnalysisHandlers[0]
set PAEventHandler:BGen:Width 10
set PAEventHandler:EffectivePartonMode Colours
## These are the analysess we will run
## First the Glauber analyses
erase PAEventHandler:AnalysisHandlers[0]
create DIPSY::GlauberAnalysis Glauber GlauberAnalysis.so
set Glauber:tMax 4
set Glauber:Nt 200
insert PAEventHandler:AnalysisHandlers[0] Glauber
## Some semi-inclusive cross section for DIPSY which need at least
## four combinations of left- and right-moving cascades.
create DIPSY::SemiInclusiveXSecAnalysis SemiIncl SemiInclusiveXSecAnalysis.so
insert PAEventHandler:AnalysisHandlers[0] SemiIncl
set PAEventHandler:PreSampleL 2
set PAEventHandler:PreSampleR 2
## This is just to keep track of the progress of a run
create DIPSY::AnalysisProgress AnaLog AnalysisProgress.so
set AnaLog:Interval 600
insert PAEventHandler:AnalysisHandlers[0] AnaLog
## Set the interaction frame
set PAEventHandler:YFrametest 0.5
# set PAEventHandler:YFrametest 0.9
## The sample rates need to be adjusted so that we get a reasonable
## statistics in a reasonable time. It is typically efficient to
## sample a number of impact parameter values for each pair of DIPSY
## cascades.
set PAEventHandler:PreSampleB 10
## But we need a good sample of cascades
set PAEventHandler:PreSamples 100000
## Set the pp energy we want to run with
set PALumi:BeamEMaxA 100
set PALumi:BeamEMaxB 100
## We run pp to get the nucleon--nucleon cross sections
saverun PA01pp0 PAGenerator
set PAEventHandler:BGen:Width 5
saverun PA01pp0t PAGenerator
set PAEventHandler:BGen:Width 10
saverun PA01pp0o PAGenerator
saverun PA01pp0c PAGenerator
# set PAEventHandler:EffectivePartonMode Relatives
set PALumi:BeamEMaxA 5000
set PALumi:BeamEMaxB 5000
saverun PA50pp1 PAGenerator
set PALumi:BeamEMaxA 200
set PALumi:BeamEMaxB 200
saverun PA02pp0 PAGenerator
set PALumi:BeamEMaxA 273
set PALumi:BeamEMaxB 273
saverun PAU2pp0 PAGenerator
set PAEventHandler:BGen:Width 1
saverun PAU2pp01 PAGenerator
set PAEventHandler:BGen:Width 10
# playground ...
set PAEventHandler:PreSampleB 4
set PAEventHandler:PreSamples 10000
saverun PAU2pp1 PAGenerator
saverun PAU2pp12 PAGenerator
set PAEventHandler:BGen:Width 4
saverun PAU2pp4 PAGenerator
set PAEventHandler:BGen:Width 10
set stdProton:R0 2.8
set PAEventHandler:RMax 3.1
set PAEventHandler:CoherenceRange 3.0
set stdEmitter:PMinusOrdering 1.0
saverun PAU2pp2 PAGenerator
set stdProton:R 0.001
+set PAEventHandler:PreSamples 100000
saverun PAU2pp3 PAGenerator
+set stdProton:R 0.5
+saverun PAU2pp5 PAGenerator
set stdProton:R 0.1
set stdEmitter:PMinusOrdering 0.8
set PAEventHandler:CoherenceRange 2.9
set PAEventHandler:RMax 2.9
set PAEventHandler:PreSampleB 10
set PAEventHandler:PreSamples 100000
# ... end
set stdProton:R0 2.9
set PALumi:BeamEMaxA 315
set PALumi:BeamEMaxB 315
saverun PAU3pp0 PAGenerator
set PALumi:BeamEMaxA 500
set PALumi:BeamEMaxB 500
saverun PA05pp0 PAGenerator
set PALumi:BeamEMaxA 900
set PALumi:BeamEMaxB 900
saverun PA09pp0 PAGenerator
set PALumi:BeamEMaxA 1000
set PALumi:BeamEMaxB 1000
saverun PA10pp0 PAGenerator
set PALumi:BeamEMaxA 2500
set PALumi:BeamEMaxB 2500
saverun PA25pp0 PAGenerator
set PALumi:BeamEMaxA 3500
set PALumi:BeamEMaxB 3500
saverun PA35pp0 PAGenerator
set PALumi:BeamEMaxA 5000
set PALumi:BeamEMaxB 5000
saverun PA50pp0 PAGenerator
set PAEventHandler:PreSampleB 1
set PAEventHandler:PreSamples 10000
create ThePEG::LWHFactory HFac LWHFactory.so
set HFac:StoreType flat
set HFac:Suffix dat
set PAGenerator:HistogramFactory HFac
create DIPSY::DipoleDensityAnalysis Density DipoleDensityAnalysis.so
insert PAEventHandler:AnalysisHandlers[0] Density
saverun Density70 PAGenerator
set PAEventHandler:Swinger NULL
saverun Density70ns PAGenerator
set PAEventHandler:Swinger stdSwinger
set PALumi:BeamEMaxA 1750
set PALumi:BeamEMaxB 1750
saverun Density35 PAGenerator
set PALumi:BeamEMaxA 875
set PALumi:BeamEMaxB 875
saverun Density18 PAGenerator
erase PAEventHandler:AnalysisHandlers[0]
set PAEventHandler:PreSampleB 10
set PAEventHandler:PreSamples 100000
set PALumi:BeamEMaxA 3500
set PALumi:BeamEMaxB 3500
## Now we take these cross sections and feed them into the Glauber analysis
## My numbers for 100 GeV (in millibarns)
set Glauber:TotalnnXSec 49.71
set Glauber:ElasticnnXSec 9.10
set Glauber:InElasticnnXSec 35.80
set PALumi:BeamEMaxA 100
set PALumi:BeamEMaxB 100
## Run again to make sure that the cross sections are reproduced, at least by "grey3 disc"
saverun PA01pp1 PAGenerator
set PAEventHandler:PreSamples 1000000
saverun PA01pp2 PAGenerator
saverun PA01pp3 PAGenerator
## Interlude for gamma-p Q^2 = 5 GeV^2
set virtualPhoton:Q2 5.0
set PAEventHandler:WFL virtualPhoton
set PAEventHandler:PreSamples 100000
saverun PA01gp1 PAGenerator
set ElXSec:NBins 200
set ElXSec:DeltaB 0.25
set ElXSec:NqBins 400
set ElXSec:Deltaq 0.02
set PAEventHandler:PreSamples 10000
saverun PA01gp2 PAGenerator
set Glauber:TotalnnXSec 0.0268
set Glauber:ElasticnnXSec 0.00000531
set Glauber:InElasticnnXSec 0.0227
set PAEventHandler:PreSamples 1000000
+set virtualPhoton:Q2 45.0
saverun PA01gp3 PAGenerator
+set virtualPhoton:Q2 15.0
+saverun PA01gp4 PAGenerator
+set virtualPhoton:Q2 3.5
+saverun PA01gp5 PAGenerator
+set virtualPhoton:Q2 90.0
+saverun PA01gp6 PAGenerator
+set PALumi:BeamEMaxA 50
+set PALumi:BeamEMaxB 50
+set virtualPhoton:Q2 45.0
+saverun PA01gp7 PAGenerator
+set PALumi:BeamEMaxA 100
+set PALumi:BeamEMaxB 100
+
set PAEventHandler:WFL stdProton
## My numbers for 200 GeV
set Glauber:TotalnnXSec 58.52
set Glauber:ElasticnnXSec 11.64
set Glauber:InElasticnnXSec 41.05
set PALumi:BeamEMaxA 200
set PALumi:BeamEMaxB 200
set PAEventHandler:PreSamples 100000
saverun PA02pp1 PAGenerator
set PAEventHandler:PreSamples 1000000
saverun PA02pp2 PAGenerator
saverun PA02pp3 PAGenerator
## My numbers for 500 GeV
set Glauber:TotalnnXSec 70.53
set Glauber:ElasticnnXSec 15.36
set Glauber:InElasticnnXSec 47.88
set PALumi:BeamEMaxA 500
set PALumi:BeamEMaxB 500
set PAEventHandler:PreSamples 100000
saverun PA05pp1 PAGenerator
set PAEventHandler:PreSamples 1000000
saverun PA05pp2 PAGenerator
saverun PA05pp3 PAGenerator
## My numbers for 1000 GeV
set Glauber:TotalnnXSec 79.63
set Glauber:ElasticnnXSec 18.35
set Glauber:InElasticnnXSec 52.85
set PALumi:BeamEMaxA 1000
set PALumi:BeamEMaxB 1000
set PAEventHandler:PreSamples 100000
saverun PA10pp1 PAGenerator
set PAEventHandler:PreSamples 1000000
saverun PA10pp2 PAGenerator
saverun PA10pp3 PAGenerator
## My numbers for 2500 GeV
set Glauber:TotalnnXSec 91.43
set Glauber:ElasticnnXSec 22.35
set Glauber:InElasticnnXSec 59.12
set PALumi:BeamEMaxA 2500
set PALumi:BeamEMaxB 2500
set PAEventHandler:PreSamples 100000
saverun PA25pp1 PAGenerator
saverun PA25pp2 PAGenerator
set PAEventHandler:PreSamples 1000000
saverun PA25pp3 PAGenerator
set Glauber:TotalnnXSec 95.86
set Glauber:ElasticnnXSec 23.88
set Glauber:InElasticnnXSec 61.41
set PALumi:BeamEMaxA 3500
set PALumi:BeamEMaxB 3500
set PAEventHandler:PreSamples 100000
saverun PA35pp1 PAGenerator
saverun PA35pp2 PAGenerator
set PALumi:BeamEMaxA 5000
set PALumi:BeamEMaxB 5000
set PAEventHandler:PreSamples 100000
saverun PA50pp2 PAGenerator
## Finally it's time to run with a heavy ion.
create DIPSY::GaussianImpactGenerator BGauss GaussianImpactGenerator.so
set PAEventHandler:BGen BGauss
set BGauss:Width 15
set PAEventHandler:WFR Oxygen
set PALumi:BeamEMaxA 100
set PALumi:BeamEMaxB 1600
set Glauber:TotalnnXSec 49.71
set Glauber:ElasticnnXSec 9.10
set Glauber:InElasticnnXSec 35.80
set Copper:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 1000
saverun PA01pO1 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA01pO2 PAGenerator
set Oxygen:Rn 0
saverun PA01pO2a PAGenerator
set Oxygen:Rn -1.3
saverun PA01pO2b PAGenerator
set Oxygen:Rn 1.3
saverun PA01pO2c PAGenerator
set Oxygen:w 0.051
saverun PA01pO2d PAGenerator
set Oxygen:Rn 0
saverun PA01pO2e PAGenerator
set PAEventHandler:WFR Copper
# set PAEventHandler:WFR Lead
set PALumi:BeamEMaxA 100
set PALumi:BeamEMaxB 6300
set Glauber:TotalnnXSec 49.71
set Glauber:ElasticnnXSec 9.10
set Glauber:InElasticnnXSec 35.80
set Copper:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 1000
saverun PA01pCu1 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA01pCu2 PAGenerator
set PALumi:BeamEMaxA 200
set PALumi:BeamEMaxB 12600
set Glauber:TotalnnXSec 58.52
set Glauber:ElasticnnXSec 11.64
set Glauber:InElasticnnXSec 41.05
set Copper:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 1000
saverun PA02pCu1 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA02pCu2 PAGenerator
set PALumi:BeamEMaxA 500
set PALumi:BeamEMaxB 31500
set Glauber:TotalnnXSec 70.53
set Glauber:ElasticnnXSec 15.36
set Glauber:InElasticnnXSec 47.88
set Copper:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 1000
saverun PA05pCu1 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA05pCu2 PAGenerator
set Copper:InterNucleonSwing On
do Copper:SetNucleus Cu GLISSANDO
set Copper:Recenter true
saverun PA05pCu3 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA05pCu5 PAGenerator
set Copper:InterNucleonSwing On
set Copper:Recenter false
saverun PA05pCu4 PAGenerator
do Copper:SetNucleus Cu
set PALumi:BeamEMaxA 1000
set PALumi:BeamEMaxB 63000
set Glauber:TotalnnXSec 79.63
set Glauber:ElasticnnXSec 18.35
set Glauber:InElasticnnXSec 52.85
set Copper:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 500
saverun PA10pCu1 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA10pCu2 PAGenerator
saverun PA10pCu3 PAGenerator
set PALumi:BeamEMaxA 2500
set PALumi:BeamEMaxB 157500
set Glauber:TotalnnXSec 91.43
set Glauber:ElasticnnXSec 22.35
set Glauber:InElasticnnXSec 59.12
set Copper:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 500
saverun PA25pCu1 PAGenerator
set Copper:InterNucleonSwing Off
saverun PA25pCu2 PAGenerator
# Let's do lead as well
set PAEventHandler:WFR Lead
# Just testing Glauber first
set Lead:Rn -0.9
set Lead:R 6.407
set Lead:a 0.459
set PALumi:BeamEMaxA 10
set PALumi:BeamEMaxB 2080
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 500
set Glauber:TotalnnXSec 49.71
set Glauber:ElasticnnXSec 9.10
set Glauber:InElasticnnXSec 35.80
saverun PA01pPbG PAGenerator
set Glauber:TotalnnXSec 58.52
set Glauber:ElasticnnXSec 11.64
set Glauber:InElasticnnXSec 41.05
saverun PA02pPbG PAGenerator
set Glauber:TotalnnXSec 70.53
set Glauber:ElasticnnXSec 15.36
set Glauber:InElasticnnXSec 47.88
saverun PA05pPbG PAGenerator
set Glauber:TotalnnXSec 79.63
set Glauber:ElasticnnXSec 18.35
set Glauber:InElasticnnXSec 52.85
saverun PA10pPbG PAGenerator
set Glauber:TotalnnXSec 91.43
set Glauber:ElasticnnXSec 22.35
set Glauber:InElasticnnXSec 59.12
set PAEventHandler:WFR stdProton
set PAEventHandler:PreSamples 5000
saverun PA10ppG PAGenerator
set PAEventHandler:PreSamples 500
set PAEventHandler:WFR Lead
saverun PA25pPbG PAGenerator
set Glauber:TotalnnXSec 95.6
set Glauber:ElasticnnXSec 23.82
set Glauber:InElasticnnXSec 61.22
saverun PA35pPbG PAGenerator
set Glauber:TotalnnXSec 95.86
set Glauber:ElasticnnXSec 23.88
set Glauber:InElasticnnXSec 61.41
saverun PA35pPbG2 PAGenerator
set Glauber:TotalnnXSec 95.79
set Glauber:ElasticnnXSec 23.86
set Glauber:InElasticnnXSec 61.36
saverun PA35pPbG3 PAGenerator
set Glauber:TotalnnXSec 99.8
set Glauber:ElasticnnXSec 25.1
set Glauber:InElasticnnXSec 63.4
saverun PA50pPbG PAGenerator
set PAEventHandler:PreSampleB 1
set PAEventHandler:PreSampleL 1
set PAEventHandler:PreSamples 20000
set Lead:R 6.407
set Lead:a 0.459
set Lead:Rn -0.9
set Lead:w 0
set Lead:NTry 0
saverun PA35pPb9GL0 PAGenerator
set Lead:NTry 1
saverun PA35pPb9GL1 PAGenerator
set Lead:NTry 100
saverun PA35pPb9GL100 PAGenerator
set Lead:R 6.491
set Lead:a 0.54
set Lead:Rn 0.0
set Lead:w 0
set Lead:NTry 0
saverun PA35pPb0GL0 PAGenerator
set Lead:NTry 1
saverun PA35pPb0GL1 PAGenerator
set Lead:NTry 100
saverun PA35pPb0GL100 PAGenerator
set Lead:R 6.531
set Lead:a 0.46
set Lead:Rn 0.0
set Lead:w 0
set Lead:NTry 0
saverun PA35pPb1GL0 PAGenerator
set Lead:NTry 1
saverun PA35pPb1GL1 PAGenerator
set Lead:NTry 100
saverun PA35pPb1GL100 PAGenerator
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 500
do Lead:SetNucleus Pb
set Lead:Rn 1.3
# Now include also DIPSY
set PALumi:BeamEMaxA 100
set PALumi:BeamEMaxB 20800
set Glauber:TotalnnXSec 49.71
set Glauber:ElasticnnXSec 9.10
set Glauber:InElasticnnXSec 35.80
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 500
saverun PA01pPb1 PAGenerator
saverun PA01pPb1o PAGenerator
set Lead:Rn -1.3
saverun PA01pPb1c PAGenerator
set Lead:Rn 1.3
saverun PA01pPb1C PAGenerator
set Lead:Rn 0.0
saverun PA01pPb1h PAGenerator
set Lead:Rn -0.9
set Lead:R 6.407
set Lead:a 0.459
saverun PA01pPb1g PAGenerator
set PAEventHandler:PreSamples 50
saverun PA01pPb1g2 PAGenerator
set PAEventHandler:PreSamples 500
do Lead:SetNucleus Pb
set Lead:Rn 1.3
set Lead:InterNucleonSwing Off
saverun PA01pPb2 PAGenerator
set PALumi:BeamEMaxA 200
set PALumi:BeamEMaxB 41600
set Glauber:TotalnnXSec 58.52
set Glauber:ElasticnnXSec 11.64
set Glauber:InElasticnnXSec 41.05
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 250
saverun PA02pPb1 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA02pPb2 PAGenerator
set Lead:InterNucleonSwing On
set Lead:Rn -0.9
set Lead:R 6.407
set Lead:a 0.459
saverun PA02pPb3 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA02pPb4 PAGenerator
set Lead:InterNucleonSwing On
do Lead:SetNucleus Pb
set Lead:Rn 1.3
set PALumi:BeamEMaxA 500
set PALumi:BeamEMaxB 10400
set Glauber:TotalnnXSec 70.53
set Glauber:ElasticnnXSec 15.36
set Glauber:InElasticnnXSec 47.88
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 250
saverun PA05pPb1 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA05pPb2 PAGenerator
do Lead:SetNucleus Pb GLISSANDO
set Lead:Recenter true
saverun PA05pPb3 PAGenerator
set Lead:Recenter false
saverun PA05pPb4 PAGenerator
do Lead:SetNucleus Pb
set PALumi:BeamEMaxA 1000
set PALumi:BeamEMaxB 208000
set Glauber:TotalnnXSec 79.63
set Glauber:ElasticnnXSec 18.35
set Glauber:InElasticnnXSec 52.85
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 250
saverun PA10pPb1 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA10pPb2 PAGenerator
set PALumi:BeamEMaxA 2500
set PALumi:BeamEMaxB 416000
set Glauber:TotalnnXSec 91.43
set Glauber:ElasticnnXSec 22.35
set Glauber:InElasticnnXSec 59.12
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 100
saverun PA25pPb1 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA25pPb2 PAGenerator
set Lead:InterNucleonSwing On
set Lead:Rn -0.9
set Lead:R 6.407
set Lead:a 0.459
saverun PA25pPb3 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA25pPb4 PAGenerator
set Lead:InterNucleonSwing On
do Lead:SetNucleus Pb GLISSANDO
saverun PA25pPb5 PAGenerator
set Lead:Recenter true
saverun PA25pPb7 PAGenerator
set Lead:Recenter false
set Lead:InterNucleonSwing Off
saverun PA25pPb6 PAGenerator
set Lead:Recenter true
saverun PA25pPb8 PAGenerator
set Lead:Recenter false
set PALumi:BeamEMaxA 10
set PALumi:BeamEMaxB 2080
set PAEventHandler:PreSamples 5000
saverun PA25pPbGG PAGenerator
set Lead:Recenter true
saverun PA25pPbGR PAGenerator
set Lead:Recenter false
set PAEventHandler:PreSamples 100
set Lead:InterNucleonSwing On
do Lead:SetNucleus Pb
set Lead:Rn 1.3
set PALumi:BeamEMaxA 3500
set PALumi:BeamEMaxB 728000
set Glauber:TotalnnXSec 95.79
set Glauber:ElasticnnXSec 23.86
set Glauber:InElasticnnXSec 61.36
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 40
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 100
saverun PA35pPb1 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA35pPb2 PAGenerator
set Lead:InterNucleonSwing On
set Lead:Rn -0.9
set Lead:R 6.407
set Lead:a 0.459
saverun PA35pPb3 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA35pPb4 PAGenerator
set Lead:InterNucleonSwing On
do Lead:SetNucleus Pb
set Lead:Rn 1.3
set PALumi:BeamEMaxA 5000
set PALumi:BeamEMaxB 832000
set Glauber:TotalnnXSec 99.54
set Glauber:ElasticnnXSec 25.12
set Glauber:InElasticnnXSec 63.2
set Lead:InterNucleonSwing On
set PAEventHandler:PreSampleB 20
set PAEventHandler:PreSampleL 10
set PAEventHandler:PreSamples 100
saverun PA50pPb1 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA50pPb2 PAGenerator
set Lead:InterNucleonSwing On
set Lead:Rn -0.9
set Lead:R 6.407
set Lead:a 0.459
saverun PA50pPb3 PAGenerator
set Lead:InterNucleonSwing Off
saverun PA50pPb4 PAGenerator
set Lead:InterNucleonSwing On
do Lead:SetNucleus Pb GLISSANDO
saverun PA50pPb5 PAGenerator
set Lead:Recenter true
saverun PA50pPb7 PAGenerator
set Lead:Recenter false
set Lead:InterNucleonSwing Off
saverun PA50pPb6 PAGenerator
set Lead:Recenter true
saverun PA50pPb8 PAGenerator
set Lead:Recenter false
set PALumi:BeamEMaxA 10
set PALumi:BeamEMaxB 2080
set PAEventHandler:PreSamples 5000
saverun PA50pPbGG PAGenerator
set Lead:Recenter true
saverun PA50pPbGR PAGenerator
set Lead:Recenter false
# set PALumi:BeamEMaxB 12600
# set PALumi:BeamEMaxB 31500
# set PALumi:BeamEMaxB 63000
# set PALumi:BeamEMaxB 157500
## We may need to adjust the sampling. Especially since the nucleus
## cascade is time consuming, it may be worth running several protons
## per nucleus.
# set PAEventHandler:PreSampleB 10
set PAEventHandler:PreSampleL 10
set PAEventHandler:PreSamples 2000
## Let's go
# saverun PA01pCu2 PAGenerator
## And let's run one more time with the swing between nucleons turned
## off
set Copper:InterNucleonSwing Off
# saverun PA01pCu3 PAGenerator
set PAEventHandler:BGen:Width 10
# saverun PA01pCu4 PAGenerator
set PAEventHandler:BGen:Width 15
# saverun PA01pCu5 PAGenerator
set PAEventHandler:PreSamples 200
create DIPSY::GaussianImpactGenerator BGauss GaussianImpactGenerator.so
# saverun PA01pCu6 PAGenerator
set BGauss:Width 21
set PAEventHandler:BGen BGauss
# saverun PA01pCu7 PAGenerator
set BGauss:Width 15
set Copper:InterNucleonSwing On
set PAEventHandler:PreSamples 1000
# saverun PA01pCu80 PAGenerator
set PAEventHandler:PreSampleL 20
set PAEventHandler:PreSamples 500
# saverun PA01pCu81 PAGenerator
set PAEventHandler:PreSampleR 4
set PAEventHandler:PreSamples 250
# saverun PA01pCu82 PAGenerator

File Metadata

Mime Type
text/x-diff
Expires
Sat, May 3, 5:54 AM (17 h, 17 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4982863
Default Alt Text
(23 KB)

Event Timeline