diff --git a/DIPSY/PAXSec.in b/DIPSY/PAXSec.in
--- a/DIPSY/PAXSec.in
+++ b/DIPSY/PAXSec.in
@@ -1,862 +1,871 @@
 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]
+# 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 1000000
 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
+create DIPSY::OldStyleEmitter oldEmitter OldStyleEmitter.so
+set PAEventHandler:Emitter oldEmitter
+set PAEventHandler:PreSamples 10000
+saverun PA35pp0OLD PAGenerator
+set PAEventHandler:Emitter stdEmitter
+saverun PA35pp0NEW PAGenerator
+set PAEventHandler:PreSamples 100000
+
+
 set PALumi:BeamEMaxA 5000
 set PALumi:BeamEMaxB 5000
 saverun PA50pp0 PAGenerator
 
 set PAEventHandler:PreSampleB 1
 set PAEventHandler:PreSamples 10000
 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 Glauber:SDnnXSec 3.75
 set PALumi:BeamEMaxA 200
 set PALumi:BeamEMaxB 200
 set PAEventHandler:PreSamples 100000
 saverun PA02pp1 PAGenerator
 saverun PA02pp4 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 Glauber:SDnnXSec 7.13
 set PALumi:BeamEMaxA 2500
 set PALumi:BeamEMaxB 2500
 set PAEventHandler:PreSamples 100000
 saverun PA25pp1 PAGenerator
 saverun PA25pp2 PAGenerator
 saverun PA25pp5 PAGenerator
 saverun PA25pp6 PAGenerator
 set PAEventHandler:PreSampleL 20
 set PAEventHandler:PreSamples 100
 saverun PA25pp7 PAGenerator
 set PAEventHandler:PreSampleL 200
 saverun PA25pp8 PAGenerator
 set PAEventHandler:PreSampleL 100
 set PAEventHandler:PreSampleB 100
 saverun PA25pp9 PAGenerator
 set PAEventHandler:PreSampleR 10
 set PAEventHandler:PreSamples 100
 set PAEventHandler:PreSampleL 20
 set PAEventHandler:PreSampleB 20
 saverun PA25pp22 PAGenerator
 set PAEventHandler:PreSampleL 20
 set PAEventHandler:PreSampleB 80
 saverun PA25pp28 PAGenerator
 set PAEventHandler:PreSampleL 20
 set PAEventHandler:PreSampleB 40
 saverun PA25pp24 PAGenerator
 set PAEventHandler:PreSampleL 40
 set PAEventHandler:PreSampleB 20
 saverun PA25pp42 PAGenerator
 set PAEventHandler:PreSampleL 80
 set PAEventHandler:PreSampleB 20
 saverun PA25pp82 PAGenerator
 set PAEventHandler:PreSampleL 40
 set PAEventHandler:PreSampleB 40
 saverun PA25pp44 PAGenerator
 set PAEventHandler:PreSampleR 2
 set PAEventHandler:PreSamples 10
 saverun PA25pp442 PAGenerator
 
 set PAEventHandler:PreSampleR 10
 set PAEventHandler:PreSamples 100
 set PAEventHandler:PreSampleL 80
 set PAEventHandler:PreSampleB 80
 saverun PA25pp88 PAGenerator
 
 set PAEventHandler:PreSampleR 2
 set PAEventHandler:PreSamples 100
 set PAEventHandler:PreSampleL 200
 set PAEventHandler:PreSampleB 200
 saverun PA25pp20 PAGenerator
 
 set PAEventHandler:PreSampleR 2
 set PAEventHandler:PreSamples 25
 set PAEventHandler:PreSampleL 400
 set PAEventHandler:PreSampleB 400
 saverun PA25pp40 PAGenerator
 
 set PAEventHandler:PreSampleR 1
 set PAEventHandler:PreSamples 20
 set PAEventHandler:PreSampleL 800
 set PAEventHandler:PreSampleB 800
 saverun PA25pp80 PAGenerator
 
 set PAEventHandler:PreSampleR 2
 set PAEventHandler:PreSampleL 2
 set PAEventHandler:PreSampleB 10
 set PAEventHandler:PreSamples 1000000
 saverun PA25pp3 PAGenerator
 saverun PA25pp4 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 Glauber:SDnnXSec 7.13
 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 Glauber:SDnnXSec 7.13
 set PAEventHandler:WFR stdProton
 set PAEventHandler:PreSamples 5000
 saverun PA10ppG PAGenerator
 set PAEventHandler:PreSamples 100000
 set PAEventHandler:PreSampleL 2
 saverun PA10ppG1 PAGenerator
 set PAEventHandler:PreSamples 100
 saverun PA10ppG10 PAGenerator
 set PAEventHandler:PreSamples 1000
 saverun PA10ppG11 PAGenerator
 set PAEventHandler:PreSamples 10000
 saverun PA10ppG12 PAGenerator
 set PAEventHandler:PreSamples 100000
 set PAEventHandler:PreSampleL 20
 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 520000
 set Glauber:TotalnnXSec 91.43
 set Glauber:ElasticnnXSec 22.35
 set Glauber:InElasticnnXSec 59.12
 set Glauber:SDnnXSec 7.13
 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 Glauber:TotalnnXSec 90.9
 set Glauber:ElasticnnXSec 22.15
 set Glauber:InElasticnnXSec 58.75
 set Glauber:SDnnXSec 7.16
 saverun PA25pPb20 PAGenerator
 set PAEventHandler:PreSampleB 20
 set PAEventHandler:PreSampleL 10
 set Lead:InterNucleonSwing On
 saverun PA25pPb1500 PAGenerator
 set Glauber:ParentalMode 1
 saverun PA25pPb1501 PAGenerator
 set Glauber:ParentalMode 2
 saverun PA25pPb1502 PAGenerator
 set Glauber:ParentalMode 0
 set Lead:InterNucleonSwing Off
 set PAEventHandler:YFrametest 0.99
 saverun PA25pPb299 PAGenerator
 set PAEventHandler:PreSamples 200
 saverun PA25pPb2990 PAGenerator
 set Glauber:ParentalMode 1
 saverun PA25pPb2991 PAGenerator
 set Glauber:ParentalMode 2
 saverun PA25pPb2992 PAGenerator
 set Glauber:ParentalMode 0
 set Glauber:TotalnnXSec 91.43
 set Glauber:ElasticnnXSec 22.35
 set Glauber:InElasticnnXSec 59.12
 set Glauber:SDnnXSec 7.13
 set PAEventHandler:PreSamples 100
 set PAEventHandler:YFrametest 0.5
 set PAEventHandler:PreSampleB 40
 set PAEventHandler:PreSampleL 20
 
 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
 saverun PA25pPbGR2 PAGenerator
 set PAEventHandler:PreSamples 500
 saverun PA25pPbGR3 PAGenerator
 set Glauber:TotalnnXSec 89.78
 set Glauber:ElasticnnXSec 21.89
 set Glauber:InElasticnnXSec 57.98
 set Glauber:SDnnXSec 7.07
 saverun PA25pPbGR4 PAGenerator
 set PAEventHandler:PreSamples 5000
 saverun PA25pPbGR5 PAGenerator
 set PAEventHandler:PreSamples 500
 set Glauber:TotalnnXSec 86.01
 set Glauber:ElasticnnXSec 18.13
 set Glauber:InElasticnnXSec 48.02
 set Glauber:SDnnXSec 12.13
 saverun PA25pPbGP4 PAGenerator
 saverun PA25pPbGP40 PAGenerator
 set Lead:NTry 1
 saverun PA25pPbGP41 PAGenerator
 set PAEventHandler:PreSamples 5000
 saverun PA25pPbGP51 PAGenerator
 set Lead:NTry 100
 saverun PA25pPbGP5 PAGenerator
 # sigtot taken from ?
 set Glauber:TotalnnXSec 93.76
 # sigel taken from ?
 set Glauber:ElasticnnXSec 23.21
 # should be (sigtot - sigel)*(1 - 20% (SD) - 12% (DD))
 set Glauber:InElasticnnXSec 47.97
 # should be (sigtot - sigel)*20%
 set Glauber:SDnnXSec 14.11
 saverun PA25pPbGP6 PAGenerator
 saverun PA25pPbGP6S PAGenerator
 set PAEventHandler:PreSamples 500
 saverun PA25pPbGP7S PAGenerator
 set PAEventHandler:PreSamples 500
 set Lead:NTry 1
 saverun PA25pPbGP61S PAGenerator
 saverun PA25pPbGP71S PAGenerator
 saverun PA25pPbGP71s PAGenerator
 set PAEventHandler:WFR stdProton
 set PAEventHandler:PreSamples 100
 saverun PA25ppGP61S PAGenerator
 saverun PA25ppGP71S PAGenerator
 saverun PA25ppGP71s PAGenerator
 set PAEventHandler:WFR Lead
 set PAEventHandler:PreSamples 500
 
 
 set Lead:Recenter false
 set Lead:NTry 100
 
 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 1040000
 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
 # 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