[Rivet-svn] r3834 - in schools/2010-Lauterbad/setups: . day1 day1/herwig day1/pythia day1/shared day1/sherpa day2 day2/herwig day2/herwig/ME day2/herwig/QED day2/herwig/shower day2/pythia day2/pythia/ME day2/pythia/QED day2/pythia/shower day2/shared day2/sherpa day2/sherpa/ME day2/sherpa/QED day2/sherpa/merging day3 day3/herwig day3/pythia day3/shared day3/shared/publications day3/shared/rivet day3/sherpa

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Tue Jul 17 14:44:01 BST 2012


Author: hoeth
Date: Tue Jul 17 14:44:00 2012
New Revision: 3834

Log:
Lauterbad setups

Added:
   schools/2010-Lauterbad/setups/
   schools/2010-Lauterbad/setups/day1/
   schools/2010-Lauterbad/setups/day1/herwig/
   schools/2010-Lauterbad/setups/day1/pythia/
   schools/2010-Lauterbad/setups/day1/shared/
   schools/2010-Lauterbad/setups/day1/sherpa/
   schools/2010-Lauterbad/setups/day1/sherpa/Run.15-3.dat
   schools/2010-Lauterbad/setups/day1/sherpa/Run.20-3.dat
   schools/2010-Lauterbad/setups/day1/sherpa/Run.dat
   schools/2010-Lauterbad/setups/day2/
   schools/2010-Lauterbad/setups/day2/herwig/
   schools/2010-Lauterbad/setups/day2/herwig/ME/
   schools/2010-Lauterbad/setups/day2/herwig/ME/Base.in
   schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-0jet.in
   schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-1jet.in
   schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-Powheg.in
   schools/2010-Lauterbad/setups/day2/herwig/QED/
   schools/2010-Lauterbad/setups/day2/herwig/QED/Base.in
   schools/2010-Lauterbad/setups/day2/herwig/QED/TVT-0jet-QED-off.in
   schools/2010-Lauterbad/setups/day2/herwig/QED/TVT-0jet-QED-on.in
   schools/2010-Lauterbad/setups/day2/herwig/shower/
   schools/2010-Lauterbad/setups/day2/herwig/shower/Base.in
   schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-0jet-MEC.in
   schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-0jet-soft.in
   schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-1jet.in
   schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-Powheg.in
   schools/2010-Lauterbad/setups/day2/pythia/
   schools/2010-Lauterbad/setups/day2/pythia/ME/
   schools/2010-Lauterbad/setups/day2/pythia/ME/0jet.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/ME/1jet.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/ME/pythia_rivet   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/pythia/QED/
   schools/2010-Lauterbad/setups/day2/pythia/QED/QED.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/QED/noQED.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/QED/pythia_rivet   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/pythia/Zpythia.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/shower/
   schools/2010-Lauterbad/setups/day2/pythia/shower/0jet.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/shower/1jet.cmnd
   schools/2010-Lauterbad/setups/day2/pythia/shower/pythia_rivet   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/shared/
   schools/2010-Lauterbad/setups/day2/shared/CDF_2008_S7540469.pdf   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/shared/D0_2007_S7075677.pdf   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/shared/D0_2008_S7554427.pdf   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/shared/D0_2008_S7863608.pdf   (contents, props changed)
   schools/2010-Lauterbad/setups/day2/sherpa/
   schools/2010-Lauterbad/setups/day2/sherpa/ME/
   schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.0jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.1jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.2jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.3jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/QED/
   schools/2010-Lauterbad/setups/day2/sherpa/QED/Run.ME.dat
   schools/2010-Lauterbad/setups/day2/sherpa/QED/Run.Off.dat
   schools/2010-Lauterbad/setups/day2/sherpa/merging/
   schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.0jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.1jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.2jet.dat
   schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.3jet.dat
   schools/2010-Lauterbad/setups/day3/
   schools/2010-Lauterbad/setups/day3/herwig/
   schools/2010-Lauterbad/setups/day3/herwig/Minbias-1960.in
   schools/2010-Lauterbad/setups/day3/herwig/Minbias-900.in
   schools/2010-Lauterbad/setups/day3/pythia/
   schools/2010-Lauterbad/setups/day3/pythia/minbias1960.cmnd
   schools/2010-Lauterbad/setups/day3/pythia/minbias900.cmnd
   schools/2010-Lauterbad/setups/day3/shared/
   schools/2010-Lauterbad/setups/day3/shared/publications/
   schools/2010-Lauterbad/setups/day3/shared/publications/ATLAS_2010_S8591806.pdf   (contents, props changed)
   schools/2010-Lauterbad/setups/day3/shared/publications/CDF_2009_S8233977.pdf   (contents, props changed)
   schools/2010-Lauterbad/setups/day3/shared/rivet/
   schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.aida
   schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.cc
   schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.info
   schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.plot
   schools/2010-Lauterbad/setups/day3/shared/rivet/CDF_2009_S8233977.cc
   schools/2010-Lauterbad/setups/day3/shared/rivet/Makefile
   schools/2010-Lauterbad/setups/day3/shared/setup.sh
   schools/2010-Lauterbad/setups/day3/sherpa/
   schools/2010-Lauterbad/setups/day3/sherpa/Run.1960.dat
   schools/2010-Lauterbad/setups/day3/sherpa/Run.900.dat

Added: schools/2010-Lauterbad/setups/day1/sherpa/Run.15-3.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day1/sherpa/Run.15-3.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,129 @@
+(run){
+  !-------------------------------------------------------
+  !--- Run parameters ------------------------------------
+  !-------------------------------------------------------
+  EVENTS          = 10000  	! Number of Events
+  OUTPUT          = 2      	! Output level:
+  !                        	! 0=errors,1=events,2=run info,
+  !                 			 	! 3=events+information
+  !
+  ANALYSIS        = Rivet
+  ANALYSIS_OUTPUT = Analysis.15-3
+}(run)
+
+(beam){
+  !-------------------------------------------------------
+  !-- Beam parameters ------------------------------------
+  !-------------------------------------------------------
+  BEAM_1          = 2212                 ! possible beam particles:  P+, P-, e+, e-
+  BEAM_ENERGY_1   = 980.                 ! in GeV
+  K_PERP_MEAN_1   = 0.33
+  K_PERP_SIGMA_1  = 0.8
+  !
+  BEAM_2          = -2212                ! possible beam particles:  P+, P-, e+, e- 
+  BEAM_ENERGY_2   = 980.                 ! in GeV
+  K_PERP_MEAN_2   = 0.33
+  K_PERP_SIGMA_2  = 0.8
+}(beam)
+
+(processes){
+  !-------------------------------------------------------
+  !-- Processes to calculate -----------------------------
+  !-------------------------------------------------------
+  !
+  ! jet jet -> electron neutrino + up to 3 jets
+  !
+  Process : 93 93 -> 11 -12 93{3}
+  Integration_Error 0.0025 {3}    ! Accuracy to which the different final
+  Integration_Error 0.005 {4}     ! state multiplicities are integrated
+  Integration_Error 0.01 {5}
+  Order electroweak : 2
+  CKKW sqr(15/E_CMS)
+  End process
+  !
+  ! jet jet -> positron neutrino + up to 3 jets
+  !
+  Process : 93 93 -> -11 12 93{3}
+  Integration_Error 0.0025 {3}    ! Accuracy to which the different final
+  Integration_Error 0.005 {4}     ! state multiplicities are integrated
+  Integration_Error 0.01 {5}
+  Order electroweak : 2
+  CKKW sqr(15/E_CMS)
+  End process
+}(processes)
+
+(selector){
+  !-------------------------------------------------------
+  !-- Parton level selectors -----------------------------
+  !-------------------------------------------------------
+  Mass 11 -12 1 14000
+  Mass -11 12 1 14000
+}(selector)
+
+
+(me){
+  !-------------------------------------------------------
+  !-- ME generators --------------------------------------
+  !-------------------------------------------------------
+  RESULT_DIRECTORY = Results.15-3
+  ME_SIGNAL_GENERATOR = Internal Comix
+  EVENT_GENERATION_MODE = Unweighted
+}(me)
+
+(integration){
+  !-------------------------------------------------------
+  !-- Phase space setup ----------------------------------
+  !-------------------------------------------------------
+  ERROR =  1.e-3   ! Error by calculating matrix-elements 
+}(integration)
+
+(shower){
+  !-------------------------------------------------------
+  !-- Parton showers -------------------------------------
+  !-------------------------------------------------------
+  SHOWER_GENERATOR = CSS
+  CSS_PT2MIN = 1.
+}(shower)
+
+(fragmentation){
+  !-------------------------------------------------------
+  !-- Fragmentation parameters----------------------------
+  !-------------------------------------------------------
+  FRAGMENTATION = Ahadic       ! Off, Ahadic, Lund
+  DECAYMODEL    = Hadrons      ! Lund or Hadrons
+  YFS_MODE      = 2            ! Soft QED on hadron decays
+}(fragmentation)
+
+(model){
+  !-------------------------------------------------------
+  !-- Model parameters -----------------------------------
+  !-------------------------------------------------------
+  MODEL         = SM           ! Model 
+  STABLE[15] = 0            ! Set tau lepton (PDG-Id 15) unstable
+}(model)
+
+(isr){
+  PDF_LIBRARY = CTEQ6Sherpa
+  PDF_SET = cteq6.6m
+}(isr)
+
+(mi){
+  !-------------------------------------------------------
+  !-- Underlying Event Setup file ------------------------
+  !-------------------------------------------------------
+  !
+  ! general parameters
+  !
+  MI_HANDLER    = None                  ! Amisic / None
+  !
+  ! hard underlying event parameters
+  !
+  SCALE_MIN         = 2.225             ! minimum scale
+  RESCALE_EXPONENT    = 0.16            ! rescaling exponent
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_WJETS
+  } END_RIVET 
+}(analysis)

Added: schools/2010-Lauterbad/setups/day1/sherpa/Run.20-3.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day1/sherpa/Run.20-3.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,129 @@
+(run){
+  !-------------------------------------------------------
+  !--- Run parameters ------------------------------------
+  !-------------------------------------------------------
+  EVENTS          = 10000  	! Number of Events
+  OUTPUT          = 2      	! Output level:
+  !                        	! 0=errors,1=events,2=run info,
+  !                 			 	! 3=events+information
+  !
+  ANALYSIS        = Rivet
+  ANALYSIS_OUTPUT = Analysis.20-3
+}(run)
+
+(beam){
+  !-------------------------------------------------------
+  !-- Beam parameters ------------------------------------
+  !-------------------------------------------------------
+  BEAM_1          = 2212                 ! possible beam particles:  P+, P-, e+, e-
+  BEAM_ENERGY_1   = 980.                 ! in GeV
+  K_PERP_MEAN_1   = 0.33
+  K_PERP_SIGMA_1  = 0.8
+  !
+  BEAM_2          = -2212                ! possible beam particles:  P+, P-, e+, e- 
+  BEAM_ENERGY_2   = 980.                 ! in GeV
+  K_PERP_MEAN_2   = 0.33
+  K_PERP_SIGMA_2  = 0.8
+}(beam)
+
+(processes){
+  !-------------------------------------------------------
+  !-- Processes to calculate -----------------------------
+  !-------------------------------------------------------
+  !
+  ! jet jet -> electron neutrino + up to 3 jets
+  !
+  Process : 93 93 -> 11 -12 93{3}
+  Integration_Error 0.0025 {3}    ! Accuracy to which the different final
+  Integration_Error 0.005 {4}     ! state multiplicities are integrated
+  Integration_Error 0.01 {5}
+  Order electroweak : 2
+  CKKW sqr(20/E_CMS)
+  End process
+  !
+  ! jet jet -> positron neutrino + up to 3 jets
+  !
+  Process : 93 93 -> -11 12 93{3}
+  Integration_Error 0.0025 {3}    ! Accuracy to which the different final
+  Integration_Error 0.005 {4}     ! state multiplicities are integrated
+  Integration_Error 0.01 {5}
+  Order electroweak : 2
+  CKKW sqr(20/E_CMS)
+  End process
+}(processes)
+
+(selector){
+  !-------------------------------------------------------
+  !-- Parton level selectors -----------------------------
+  !-------------------------------------------------------
+  Mass 11 -12 1 14000
+  Mass -11 12 1 14000
+}(selector)
+
+
+(me){
+  !-------------------------------------------------------
+  !-- ME generators --------------------------------------
+  !-------------------------------------------------------
+  RESULT_DIRECTORY = Results.20-3
+  ME_SIGNAL_GENERATOR = Internal Comix
+  EVENT_GENERATION_MODE = Unweighted
+}(me)
+
+(integration){
+  !-------------------------------------------------------
+  !-- Phase space setup ----------------------------------
+  !-------------------------------------------------------
+  ERROR =  1.e-3   ! Error by calculating matrix-elements 
+}(integration)
+
+(shower){
+  !-------------------------------------------------------
+  !-- Parton showers -------------------------------------
+  !-------------------------------------------------------
+  SHOWER_GENERATOR = CSS
+  CSS_PT2MIN = 1.
+}(shower)
+
+(fragmentation){
+  !-------------------------------------------------------
+  !-- Fragmentation parameters----------------------------
+  !-------------------------------------------------------
+  FRAGMENTATION = Ahadic       ! Off, Ahadic, Lund
+  DECAYMODEL    = Hadrons      ! Lund or Hadrons
+  YFS_MODE      = 2            ! Soft QED on hadron decays
+}(fragmentation)
+
+(model){
+  !-------------------------------------------------------
+  !-- Model parameters -----------------------------------
+  !-------------------------------------------------------
+  MODEL         = SM           ! Model 
+  STABLE[15] = 0            ! Set tau lepton (PDG-Id 15) unstable
+}(model)
+
+(isr){
+  PDF_LIBRARY = CTEQ6Sherpa
+  PDF_SET = cteq6.6m
+}(isr)
+
+(mi){
+  !-------------------------------------------------------
+  !-- Underlying Event Setup file ------------------------
+  !-------------------------------------------------------
+  !
+  ! general parameters
+  !
+  MI_HANDLER    = None                  ! Amisic / None
+  !
+  ! hard underlying event parameters
+  !
+  SCALE_MIN         = 2.225             ! minimum scale
+  RESCALE_EXPONENT    = 0.16            ! rescaling exponent
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_WJETS
+  } END_RIVET 
+}(analysis)

Added: schools/2010-Lauterbad/setups/day1/sherpa/Run.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day1/sherpa/Run.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,127 @@
+(run){
+  !-------------------------------------------------------
+  !--- Run parameters ------------------------------------
+  !-------------------------------------------------------
+  EVENTS          = 10000       ! Number of Events
+  OUTPUT          = 2           ! Output level:
+  !                             ! 0=errors,1=events,2=run info,
+  !                             ! 3=events+information
+  !
+  ANALYSIS        = Rivet
+  ANALYSIS_OUTPUT = Analysis
+}(run)
+
+(beam){
+  !-------------------------------------------------------
+  !-- Beam parameters ------------------------------------
+  !-------------------------------------------------------
+  BEAM_1          = 2212                 ! possible beam particles:  P+, P-, e+, e-
+  BEAM_ENERGY_1   = 980.                 ! in GeV
+  K_PERP_MEAN_1   = 0.33
+  K_PERP_SIGMA_1  = 0.8
+  !
+  BEAM_2          = -2212                ! possible beam particles:  P+, P-, e+, e- 
+  BEAM_ENERGY_2   = 980.                 ! in GeV
+  K_PERP_MEAN_2   = 0.33
+  K_PERP_SIGMA_2  = 0.8
+}(beam)
+
+(processes){
+  !-------------------------------------------------------
+  !-- Processes to calculate -----------------------------
+  !-------------------------------------------------------
+  !
+  ! jet jet -> electron neutrino + up to 0 jets
+  !
+  Process : 93 93 -> 11 -12 93{0}
+  Integration_Error 0.0025 {3}    ! Accuracy to which the different final
+  Integration_Error 0.005 {4}     ! state multiplicities are integrated
+  Order electroweak : 2
+  CKKW sqr(20/E_CMS)
+  End process
+  !
+  ! jet jet -> positron neutrino + up to 0 jets
+  !
+  Process : 93 93 -> -11 12 93{0}
+  Integration_Error 0.0025 {3}    ! Accuracy to which the different final
+  Integration_Error 0.005 {4}     ! state multiplicities are integrated
+  Order electroweak : 2
+  CKKW sqr(20/E_CMS)
+  End process
+}(processes)
+
+(selector){
+  !-------------------------------------------------------
+  !-- Parton level selectors -----------------------------
+  !-------------------------------------------------------
+  Mass 11 -12 1 14000
+  Mass -11 12 1 14000
+}(selector)
+
+
+(me){
+  !-------------------------------------------------------
+  !-- ME generators --------------------------------------
+  !-------------------------------------------------------
+  RESULT_DIRECTORY = Results
+  ME_SIGNAL_GENERATOR = Internal Comix
+  EVENT_GENERATION_MODE = Unweighted
+}(me)
+
+(integration){
+  !-------------------------------------------------------
+  !-- Phase space setup ----------------------------------
+  !-------------------------------------------------------
+  ERROR =  1.e-3   ! Error by calculating matrix-elements 
+}(integration)
+
+(shower){
+  !-------------------------------------------------------
+  !-- Parton showers -------------------------------------
+  !-------------------------------------------------------
+  SHOWER_GENERATOR = CSS
+  CSS_PT2MIN = 1.
+}(shower)
+
+(fragmentation){
+  !-------------------------------------------------------
+  !-- Fragmentation parameters----------------------------
+  !-------------------------------------------------------
+  FRAGMENTATION = Ahadic       ! Off, Ahadic, Lund
+  DECAYMODEL    = Hadrons      ! Lund or Hadrons
+  YFS_MODE      = 2            ! Soft QED on hadron decays
+}(fragmentation)
+
+(model){
+  !-------------------------------------------------------
+  !-- Model parameters -----------------------------------
+  !-------------------------------------------------------
+  MODEL      = SM           ! Model 
+  STABLE[15] = 0            ! Set tau lepton (PDG-Id 15) unstable
+}(model)
+
+(isr){
+  PDF_LIBRARY = LHAPDFSherpa
+  PDF_SET = cteq61.LHpdf
+}(isr)
+
+(mi){
+  !-------------------------------------------------------
+  !-- Underlying Event Setup file ------------------------
+  !-------------------------------------------------------
+  !
+  ! general parameters
+  !
+  MI_HANDLER    = Amisic                  ! Amisic / None
+  !
+  ! hard underlying event parameters
+  !
+  SCALE_MIN         = 2.225             ! minimum scale
+  RESCALE_EXPONENT    = 0.16            ! rescaling exponent
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_WJETS
+  } END_RIVET 
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/herwig/ME/Base.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/ME/Base.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,100 @@
+##################################################
+# Global generator settings
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 1000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 1000
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+
+cp MEqq2gZ2ff MEZe
+set MEZe:GammaZ All
+cp MEZe MEZmu
+set MEZe:Process Electron
+set MEZmu:Process Muon
+
+
+cp MEZJet MEZJetE
+set MEZJetE:GammaZ All
+cp MEZJetE MEZJetMu
+set MEZJetE:ZDecay Electron
+set MEZJetMu:ZDecay Muon
+
+
+cp PowhegMEqq2gZ2ff PowhegMEZe
+set PowhegMEZe:GammaZ All
+cp PowhegMEZe PowhegMEZmu
+set PowhegMEZe:Process Electron
+set PowhegMEZmu:Process Muon
+
+##################################################
+#  Setup the POWHEG shower
+##################################################
+cd /Herwig/Shower
+create Herwig::PowhegEvolver PowhegEvolver HwPowhegShower.so
+set PowhegEvolver:ShowerModel ShowerModel
+set PowhegEvolver:SplittingGenerator SplittingGenerator
+set PowhegEvolver:MECorrMode No
+set PowhegEvolver:IntrinsicPtGaussian 1.9*GeV
+
+# create and use the Drell-yan hard emission generator
+create Herwig::DrellYanHardGenerator DrellYanHardGenerator
+set DrellYanHardGenerator:ShowerAlpha AlphaQCD
+insert PowhegEvolver:HardGenerator 0 DrellYanHardGenerator
+
+
+##################################################
+# Tevatron physics parameters (override defaults) 
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+set LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
+
+##################################################
+# Control flags
+##################################################
+
+cd /Herwig/Shower
+set SplittingGenerator:ISR No
+set SplittingGenerator:FSR No
+
+cd /Herwig/EventHandlers
+set LHCHandler:CascadeHandler:MPIHandler NULL
+set LHCHandler:HadronizationHandler NULL
+set LHCHandler:DecayHandler NULL
+set /Herwig/Analysis/Basics:CheckQuark No
+
+# disable QED radiation
+erase LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+#  Rivet analyses
+##################################################
+cd /Herwig/Generators
+create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+insert Rivet:Analyses 0 MC_ZJETS
+insert Rivet:Analyses 0 CDF_2008_S7540469
+insert Rivet:Analyses 0 D0_2008_S7863608
+insert Rivet:Analyses 0 D0_2007_S7075677
+insert Rivet:Analyses 0 D0_2008_S7554427
+#set Rivet:Debug Yes
+
+erase LHCGenerator:AnalysisHandlers 0
+insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+set /Herwig/Cuts/MassCut:MinM 66*GeV
+set /Herwig/Cuts/MassCut:MaxM 116*GeV
+
+#insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+#set /Herwig/Analysis/HepMCFile:PrintEvent 10000
+#set /Herwig/Analysis/HepMCFile:Format GenEvent
+#set /Herwig/Analysis/HepMCFile:Units GeV_mm
+
+# set /Herwig/Model:EW/Sin2ThetaW 0.2224

Added: schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-0jet.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-0jet.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,8 @@
+read Base.in
+
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZe
+insert SimpleQCD:MatrixElements[0] MEZmu
+
+cd /Herwig/Generators
+saverun TVT-0jet LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-1jet.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-1jet.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,8 @@
+read Base.in
+
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZJetE
+insert SimpleQCD:MatrixElements[0] MEZJetMu
+
+cd /Herwig/Generators
+saverun TVT-1jet LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-Powheg.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/ME/TVT-Powheg.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,34 @@
+read Base.in
+
+##################################################
+#  Need to use an NLO PDF
+##################################################
+set /Herwig/Particles/p+:PDF    /Herwig/Partons/MRST-NLO
+set /Herwig/Particles/pbar-:PDF /Herwig/Partons/MRST-NLO
+##################################################
+#  and NLO strong coupling
+##################################################
+cd /Herwig/Generators
+create Herwig::O2AlphaS O2AlphaS 
+set LHCGenerator:StandardModelParameters:QCD/RunningAlphaS O2AlphaS
+##################################################
+#  Setup the POWHEG shower
+##################################################
+cd /Herwig/Shower
+# use the general recon for now
+set KinematicsReconstructor:ReconstructionOption General
+# create the Powheg evolver and use it instead of the default one
+set ShowerHandler:Evolver PowhegEvolver
+
+
+
+
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] PowhegMEZe
+insert SimpleQCD:MatrixElements[0] PowhegMEZmu
+
+
+
+cd /Herwig/Generators
+
+saverun TVT-Powheg LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/QED/Base.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/QED/Base.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,83 @@
+##################################################
+# Global generator settings
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 1000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 1000
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+
+cp MEqq2gZ2ff MEZe
+set MEZe:GammaZ All
+cp MEZe MEZmu
+set MEZe:Process Electron
+set MEZmu:Process Muon
+
+
+cp MEZJet MEZJetE
+set MEZJetE:GammaZ All
+cp MEZJetE MEZJetMu
+set MEZJetE:ZDecay Electron
+set MEZJetMu:ZDecay Muon
+
+
+cp PowhegMEqq2gZ2ff PowhegMEZe
+set PowhegMEZe:GammaZ All
+cp PowhegMEZe PowhegMEZmu
+set PowhegMEZe:Process Electron
+set PowhegMEZmu:Process Muon
+
+##################################################
+#  Setup the POWHEG shower
+##################################################
+cd /Herwig/Shower
+create Herwig::PowhegEvolver PowhegEvolver HwPowhegShower.so
+set PowhegEvolver:ShowerModel ShowerModel
+set PowhegEvolver:SplittingGenerator SplittingGenerator
+set PowhegEvolver:MECorrMode No
+set PowhegEvolver:IntrinsicPtGaussian 1.9*GeV
+
+# create and use the Drell-yan hard emission generator
+create Herwig::DrellYanHardGenerator DrellYanHardGenerator
+set DrellYanHardGenerator:ShowerAlpha AlphaQCD
+insert PowhegEvolver:HardGenerator 0 DrellYanHardGenerator
+
+
+##################################################
+# Tevatron physics parameters (override defaults) 
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+set LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
+
+##################################################
+#  Rivet analyses
+##################################################
+cd /Herwig/Generators
+create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+insert Rivet:Analyses 0 MC_ZJETS
+insert Rivet:Analyses 0 MC_ZJETS_NOCLUS
+insert Rivet:Analyses 0 CDF_2008_S7540469
+insert Rivet:Analyses 0 D0_2008_S7863608
+insert Rivet:Analyses 0 D0_2007_S7075677
+insert Rivet:Analyses 0 D0_2008_S7554427
+#set Rivet:Debug Yes
+
+insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+set /Herwig/Cuts/MassCut:MinM 66*GeV
+set /Herwig/Cuts/MassCut:MaxM 116*GeV
+
+#insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+#set /Herwig/Analysis/HepMCFile:PrintEvent 10000
+#set /Herwig/Analysis/HepMCFile:Format GenEvent
+#set /Herwig/Analysis/HepMCFile:Units GeV_mm
+
+# set /Herwig/Model:EW/Sin2ThetaW 0.2224

Added: schools/2010-Lauterbad/setups/day2/herwig/QED/TVT-0jet-QED-off.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/QED/TVT-0jet-QED-off.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,27 @@
+##################################################
+# Global generator settings
+##################################################
+read Base.in
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZe
+insert SimpleQCD:MatrixElements[0] MEZmu
+
+##################################################
+# Control flags
+##################################################
+cd /Herwig/EventHandlers
+set LHCHandler:CascadeHandler:Evolver:MECorrMode Yes
+
+# disable QED radiation
+erase LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-0jet-QED-off LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/QED/TVT-0jet-QED-on.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/QED/TVT-0jet-QED-on.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,24 @@
+##################################################
+# Global generator settings
+##################################################
+read Base.in
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZe
+insert SimpleQCD:MatrixElements[0] MEZmu
+
+##################################################
+# Control flags
+##################################################
+cd /Herwig/EventHandlers
+set LHCHandler:CascadeHandler:Evolver:MECorrMode Yes
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-0jet-QED-on LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/shower/Base.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/shower/Base.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,89 @@
+##################################################
+# Global generator settings
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 1000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 1000
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+
+cp MEqq2gZ2ff MEZe
+set MEZe:GammaZ All
+cp MEZe MEZmu
+set MEZe:Process Electron
+set MEZmu:Process Muon
+
+
+cp MEZJet MEZJetE
+set MEZJetE:GammaZ All
+cp MEZJetE MEZJetMu
+set MEZJetE:ZDecay Electron
+set MEZJetMu:ZDecay Muon
+
+
+cp PowhegMEqq2gZ2ff PowhegMEZe
+set PowhegMEZe:GammaZ All
+cp PowhegMEZe PowhegMEZmu
+set PowhegMEZe:Process Electron
+set PowhegMEZmu:Process Muon
+
+##################################################
+#  Setup the POWHEG shower
+##################################################
+cd /Herwig/Shower
+create Herwig::PowhegEvolver PowhegEvolver HwPowhegShower.so
+set PowhegEvolver:ShowerModel ShowerModel
+set PowhegEvolver:SplittingGenerator SplittingGenerator
+set PowhegEvolver:MECorrMode No
+set PowhegEvolver:IntrinsicPtGaussian 1.9*GeV
+
+# create and use the Drell-yan hard emission generator
+create Herwig::DrellYanHardGenerator DrellYanHardGenerator
+set DrellYanHardGenerator:ShowerAlpha AlphaQCD
+insert PowhegEvolver:HardGenerator 0 DrellYanHardGenerator
+
+
+##################################################
+# Tevatron physics parameters (override defaults) 
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+set LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
+
+##################################################
+# Control flags
+##################################################
+cd /Herwig/EventHandlers
+erase LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+#  Rivet analyses
+##################################################
+cd /Herwig/Generators
+create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+insert Rivet:Analyses 0 MC_ZJETS
+insert Rivet:Analyses 0 CDF_2008_S7540469
+insert Rivet:Analyses 0 D0_2008_S7863608
+insert Rivet:Analyses 0 D0_2007_S7075677
+insert Rivet:Analyses 0 D0_2008_S7554427
+#set Rivet:Debug Yes
+
+insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+set /Herwig/Cuts/MassCut:MinM 66*GeV
+set /Herwig/Cuts/MassCut:MaxM 116*GeV
+
+
+#insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+#set /Herwig/Analysis/HepMCFile:PrintEvent 10000
+#set /Herwig/Analysis/HepMCFile:Format GenEvent
+#set /Herwig/Analysis/HepMCFile:Units GeV_mm
+
+# set /Herwig/Model:EW/Sin2ThetaW 0.2224

Added: schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-0jet-MEC.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-0jet-MEC.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,24 @@
+##################################################
+# Global generator settings
+##################################################
+read Base.in
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZe
+insert SimpleQCD:MatrixElements[0] MEZmu
+
+##################################################
+# Control flags
+##################################################
+cd /Herwig/EventHandlers
+set LHCHandler:CascadeHandler:Evolver:MECorrMode Yes
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-0jet-MEC LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-0jet-soft.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-0jet-soft.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,24 @@
+##################################################
+# Global generator settings
+##################################################
+read Base.in
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZe
+insert SimpleQCD:MatrixElements[0] MEZmu
+
+##################################################
+# Control flags
+##################################################
+cd /Herwig/EventHandlers
+set LHCHandler:CascadeHandler:Evolver:MECorrMode Soft
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-0jet-soft LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-1jet.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-1jet.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,18 @@
+##################################################
+# Global generator settings
+##################################################
+read Base.in
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEZJetE
+insert SimpleQCD:MatrixElements[0] MEZJetMu
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-1jet LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-Powheg.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/herwig/shower/TVT-Powheg.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,40 @@
+##################################################
+# Global generator settings
+##################################################
+read Base.in
+
+##################################################
+#  Need to use an NLO PDF
+##################################################
+set /Herwig/Particles/p+:PDF    /Herwig/Partons/MRST-NLO
+set /Herwig/Particles/pbar-:PDF /Herwig/Partons/MRST-NLO
+##################################################
+#  and NLO strong coupling
+##################################################
+cd /Herwig/Generators
+create Herwig::O2AlphaS O2AlphaS 
+set LHCGenerator:StandardModelParameters:QCD/RunningAlphaS O2AlphaS
+
+##################################################
+#  Setup the POWHEG shower
+##################################################
+cd /Herwig/Shower
+# use the general recon for now
+set KinematicsReconstructor:ReconstructionOption General
+# create the Powheg evolver and use it instead of the default one
+set ShowerHandler:Evolver PowhegEvolver
+
+##################################################
+# Matrix Elements for hadron-hadron collisions 
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] PowhegMEZe
+insert SimpleQCD:MatrixElements[0] PowhegMEZmu
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+
+saverun TVT-Powheg LHCGenerator

Added: schools/2010-Lauterbad/setups/day2/pythia/ME/0jet.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/ME/0jet.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,47 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 50000       ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 4) Z0 production at lowest order (+ shower with ME correction).
+WeakSingleBoson:ffbar2gmZ = on    ! switch on f fbar -> gamma/Z0
+
+! 6) A few ways to change shower activity for lowest-order Z0 production.
+SpaceShower:MEcorrections = off   ! switch off matrix corrections
+SpaceShower:pTmaxMatch = 1        ! no radiation above factorization scale
+PartonLevel:ISR = off             ! switch off all initial-state radiation
+#BeamRemnants:primordialKT = off   ! switch off primordial kT
+
+! 7) Some other switches, less important/interesting (but save time).
+PartonLevel:FSR = off             ! no final-state radiation
+PartonLevel:MI = off              ! no multiple parton interactions
+HadronLevel:all = off             ! no hadronization/decay
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/ME/1jet.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/ME/1jet.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,49 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 50000       ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 5) Alternative production process Z0 + jet (+ shower).
+WeakBosonAndParton:qqbar2gmZg = on ! switch on q qbar -> gamma/Z0 + g
+WeakBosonAndParton:qg2gmZq = on    ! switch on q g    -> gamma/Z0 + q
+PhaseSpace:pTHatMin = 20.          ! set minimum pT for 2 -> 2 processes
+
+! 6) A few ways to change shower activity for lowest-order Z0 production.
+SpaceShower:MEcorrections = off   ! switch off matrix corrections
+SpaceShower:pTmaxMatch = 1        ! no radiation above factorization scale
+PartonLevel:ISR = off             ! switch off all initial-state radiation
+#BeamRemnants:primordialKT = off   ! switch off primordial kT
+
+! 7) Some other switches, less important/interesting (but save time).
+PartonLevel:FSR = off             ! no final-state radiation
+PartonLevel:MI = off              ! no multiple parton interactions
+HadronLevel:all = off             ! no hadronization/decay
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/ME/pythia_rivet
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/ME/pythia_rivet	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,14 @@
+#! /bin/bash
+
+if [ $# != 1 ]; then 
+    echo "$0: Need a config file name as argument." >&2
+    exit 1
+fi
+
+NAME=$(basename $1 .cmnd)
+PIPENAME=$NAME.pipe
+
+mkfifo $PIPENAME
+main32.exe "$1" $PIPENAME &
+rivet -a MC_ZJETS -a CDF_2008_S7540469 -a D0_2008_S7863608 -a D0_2007_S7075677 -a D0_2008_S7554427 -H$NAME.aida  $PIPENAME
+rm -f $PIPENAME

Added: schools/2010-Lauterbad/setups/day2/pythia/QED/QED.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/QED/QED.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,36 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 10000       ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 4) Z0 production at lowest order (+ shower with ME correction).
+WeakSingleBoson:ffbar2gmZ = on    ! switch on f fbar -> gamma/Z0
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+#TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/QED/noQED.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/QED/noQED.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,36 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 10000       ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 4) Z0 production at lowest order (+ shower with ME correction).
+WeakSingleBoson:ffbar2gmZ = on    ! switch on f fbar -> gamma/Z0
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/QED/pythia_rivet
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/QED/pythia_rivet	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,14 @@
+#! /bin/bash
+
+if [ $# != 1 ]; then 
+    echo "$0: Need a config file name as argument." >&2
+    exit 1
+fi
+
+NAME=$(basename $1 .cmnd)
+PIPENAME=$NAME.pipe
+
+mkfifo $PIPENAME
+main32.exe "$1" $PIPENAME &
+rivet -a MC_ZJETS -a MC_ZJETS_NOCLUS -a CDF_2008_S7540469 -a D0_2008_S7863608 -a D0_2007_S7075677 -a D0_2008_S7554427 -H$NAME.aida  $PIPENAME
+rm -f $PIPENAME

Added: schools/2010-Lauterbad/setups/day2/pythia/Zpythia.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/Zpythia.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,52 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 1000        ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 4) Z0 production at lowest order (+ shower with ME correction).
+WeakSingleBoson:ffbar2gmZ = on    ! switch on f fbar -> gamma/Z0
+
+! 5) Alternative production process Z0 + jet (+ shower).
+#WeakBosonAndParton:qqbar2gmZg = on ! switch on q qbar -> gamma/Z0 + g
+#WeakBosonAndParton:qg2gmZq = on    ! switch on q g    -> gamma/Z0 + q
+#PhaseSpace:pTHatMin = 20.          ! set minimum pT for 2 -> 2 processes
+
+! 6) A few ways to change shower activity for lowest-order Z0 production.
+#SpaceShower:MEcorrections = off   ! switch off matrix corrections
+#SpaceShower:pTmaxMatch = 1        ! no radiation above factorization scale
+#PartonLevel:ISR = off             ! switch off all initial-state radiation
+#BeamRemnants:primordialKT = off   ! switch off primordial kT
+
+! 7) Some other switches, less important/interesting (but save time).
+#PartonLevel:FSR = off             ! no final-state radiation
+#PartonLevel:MI = off              ! no multiple parton interactions
+#HadronLevel:all = off             ! no hadronization/decay
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+#TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/shower/0jet.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/shower/0jet.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,42 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 10000       ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 4) Z0 production at lowest order (+ shower with ME correction).
+WeakSingleBoson:ffbar2gmZ = on    ! switch on f fbar -> gamma/Z0
+
+! 6) A few ways to change shower activity for lowest-order Z0 production.
+#SpaceShower:MEcorrections = off   ! switch off matrix corrections
+#SpaceShower:pTmaxMatch = 1        ! no radiation above factorization scale
+#PartonLevel:ISR = off             ! switch off all initial-state radiation
+#BeamRemnants:primordialKT = off   ! switch off primordial kT
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/shower/1jet.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/shower/1jet.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,38 @@
+! File: Zpythia.cmnd
+! This file contains commands to be read in for a Pythia8 run. 
+! Lines not beginning with a letter or digit are comments.
+! Names are case-insensitive  -  but spellings-sensitive!
+! The changes here are illustrative of different Z production scenarios.
+
+! 1) Settings that will be used in the main program.
+Main:numberOfEvents = 10000       ! number of events to generate
+Main:timesToShow = 10             ! show how far along run is this many times
+Main:timesAllowErrors = 3         ! abort run after this many flawed events
+Main:showChangedSettings = on     ! print changed flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+
+! 2) Beam parameter settings for LHC at 7 TeV.
+Beams:idA = 2212                  ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                 ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960.                 ! CM energy of collision
+
+!) 3) Restrict Z0 mass range, to avoid low-mass gamma* region.
+23:mMin = 66.                     ! lowest allowed mass (default 10 GeV)
+23:mMax = 116.                    ! highest allowed mass (default unbounded)
+
+! 5) Alternative production process Z0 + jet (+ shower).
+WeakBosonAndParton:qqbar2gmZg = on ! switch on q qbar -> gamma/Z0 + g
+WeakBosonAndParton:qg2gmZq = on    ! switch on q g    -> gamma/Z0 + q
+PhaseSpace:pTHatMin = 20.          ! set minimum pT for 2 -> 2 processes
+
+! 8) Let Z0 decays to leptons. Optionally switch off QED radiation.
+23:onMade = off                   ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 13 -13      ! ...and back on Z0 -> e- e+
+TimeShower:QEDshowerByL = off     ! switch off photon radiation off leptons
+
+! 9) Optional choice of random number sequences
+#Random:setSeed = on               ! allow to override default seed
+#Random:seed = 123456789           ! any number between 1 and 900000000
+
+! 10) Optionally choose between a few different (draft) tunes
+#Tune:pp = 3                        ! options 2, 3 and 4 are most interesting

Added: schools/2010-Lauterbad/setups/day2/pythia/shower/pythia_rivet
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/pythia/shower/pythia_rivet	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,14 @@
+#! /bin/bash
+
+if [ $# != 1 ]; then 
+    echo "$0: Need a config file name as argument." >&2
+    exit 1
+fi
+
+NAME=$(basename $1 .cmnd)
+PIPENAME=$NAME.pipe
+
+mkfifo $PIPENAME
+main32.exe "$1" $PIPENAME &
+rivet -a MC_ZJETS -a CDF_2008_S7540469 -a D0_2008_S7863608 -a D0_2007_S7075677 -a D0_2008_S7554427 -H$NAME.aida  $PIPENAME
+rm -f $PIPENAME

Added: schools/2010-Lauterbad/setups/day2/shared/CDF_2008_S7540469.pdf
==============================================================================
Binary file. No diff available.

Added: schools/2010-Lauterbad/setups/day2/shared/D0_2007_S7075677.pdf
==============================================================================
Binary file. No diff available.

Added: schools/2010-Lauterbad/setups/day2/shared/D0_2008_S7554427.pdf
==============================================================================
Binary file. No diff available.

Added: schools/2010-Lauterbad/setups/day2/shared/D0_2008_S7863608.pdf
==============================================================================
Binary file. No diff available.

Added: schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.0jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.0jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,41 @@
+(run){
+  EVENTS = 100000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.0jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{0}
+  Order_EW 2;
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.0jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+  CSS_MAXEM = 0
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.1jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.1jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,43 @@
+(run){
+  EVENTS = 10000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.1jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{1}
+  Order_EW 2;
+  Integration_Error 0.0025 {3};
+  End process;
+}(processes)
+
+(selector){
+  JetFinder sqr(20/E_CMS) 1
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.1jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+  CSS_MAXEM = 0
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.2jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.2jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,44 @@
+(run){
+  EVENTS = 10000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.2jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{2}
+  Order_EW 2;
+  Integration_Error 0.0025 {3};
+  Integration_Error 0.005 {4};
+  End process;
+}(processes)
+
+(selector){
+  JetFinder sqr(20/E_CMS) 1
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.2jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+  CSS_MAXEM = 0
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.3jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/ME/Run.3jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,45 @@
+(run){
+  EVENTS = 10000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.3jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{3}
+  Order_EW 2;
+  Integration_Error 0.0025 {3};
+  Integration_Error 0.005 {4};
+  Integration_Error 0.01 {5};
+  End process;
+}(processes)
+
+(selector){
+  JetFinder sqr(20/E_CMS) 1
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.3jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+  CSS_MAXEM = 0
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/QED/Run.ME.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/QED/Run.ME.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,41 @@
+(run){
+  EVENTS = 100000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.ME
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{0}
+  Order_EW 2;
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY      = Results.ME
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = On
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+  CSS_EW_SCHEME = 0
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS MC_ZJETS_NOCLUS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/QED/Run.Off.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/QED/Run.Off.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,41 @@
+(run){
+  EVENTS = 100000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.Off
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{0}
+  Order_EW 2;
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY      = Results.Off
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = Off
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+  CSS_EW_SCHEME = 0
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS MC_ZJETS_NOCLUS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.0jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.0jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,40 @@
+(run){
+  EVENTS = 100000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.0jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{0}
+  Order_EW 2;
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.0jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.1jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.1jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,42 @@
+(run){
+  EVENTS = 10000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.1jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{1}
+  Order_EW 2;
+  CKKW sqr(20/E_CMS)
+  Integration_Error 0.0025 {3};
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.1jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.2jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.2jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,43 @@
+(run){
+  EVENTS = 10000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.2jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{2}
+  Order_EW 2;
+  CKKW sqr(20/E_CMS)
+  Integration_Error 0.0025 {3};
+  Integration_Error 0.005 {4};
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.2jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.3jet.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day2/sherpa/merging/Run.3jet.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,44 @@
+(run){
+  EVENTS = 10000
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.3jet
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+}(beam)
+
+(processes){
+  Process 93 93 -> 11 -11 93{3}
+  Order_EW 2;
+  CKKW sqr(20/E_CMS)
+  Integration_Error 0.0025 {3};
+  Integration_Error 0.005 {4};
+  Integration_Error 0.01 {5};
+  End process;
+}(processes)
+
+(selector){
+  Mass 11 -11 66 116
+}(selector)
+
+(me){
+  RESULT_DIRECTORY = Results.3jet
+  ME_SIGNAL_GENERATOR = Internal Comix
+  ME_QED = 0
+}(me)
+
+(shower){
+  SHOWER_GENERATOR = CSS
+}(shower)
+
+(mi){
+  MI_HANDLER = None  # None or Amisic
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+  -a MC_ZJETS CDF_2008_S7540469 D0_2007_S7075677 D0_2008_S7554427
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day3/herwig/Minbias-1960.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/herwig/Minbias-1960.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,31 @@
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+set LHCGenerator:EventHandler:BeamB /Herwig/Particles/pbar-
+
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 1.25*GeV
+
+cd /Herwig/MatrixElements/
+insert SimpleQCD:MatrixElements[0] MEMinBias
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+# min ibas cuts
+set /Herwig/Cuts/JetKtCut:MinKT  0.0*GeV
+set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
+set /Herwig/Cuts/QCDCuts:X1Min 0.01
+set /Herwig/Cuts/QCDCuts:X2Min 0.01
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+
+cd /Herwig/Generators
+create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+insert Rivet:Analyses 0 CDF_2009_S8233977
+insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+
+
+saverun Minbias-1960 LHCGenerator

Added: schools/2010-Lauterbad/setups/day3/herwig/Minbias-900.in
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/herwig/Minbias-900.in	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,30 @@
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 1
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 900.0
+
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 1.25*GeV
+
+cd /Herwig/MatrixElements/
+insert SimpleQCD:MatrixElements[0] MEMinBias
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+# min ibas cuts
+set /Herwig/Cuts/JetKtCut:MinKT  0.0*GeV
+set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
+set /Herwig/Cuts/QCDCuts:X1Min 0.01
+set /Herwig/Cuts/QCDCuts:X2Min 0.01
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+
+cd /Herwig/Generators
+create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+insert Rivet:Analyses 0 ATLAS_900_MINBIAS
+insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+
+
+saverun Minbias-900 LHCGenerator

Added: schools/2010-Lauterbad/setups/day3/pythia/minbias1960.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/pythia/minbias1960.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,21 @@
+! 1) Settings that will be used in a main program.
+Main:numberOfEvents = 100000       ! number of events to generate
+Main:timesToShow = 100             ! show how far along run is this many times
+Main:timesAllowErrors = 30         ! abort run after this many flawed events
+Main:showChangedSettings = on      ! print changed flags/modes/parameters
+#Main:showAllSettings = on         ! print all flags/modes/parameters
+Main:showChangedParticleData = on  ! print changed particle and decay data
+#Main:showAllParticleData = on     ! print all particle and decay data
+
+! 2) Beam parameter settings. Values below agree with default ones.
+Beams:idA = 2212                   ! first beam, p = 2212, pbar = -2212
+Beams:idB = -2212                  ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 1960                   ! CM energy of collision
+
+! 3) Pick processes and kinematics cuts
+SoftQCD:minBias on
+
+! Particles with lifetime c*tau > 10 mm are considered stable
+ParticleDecays:limitTau0 on
+ParticleDecays:tau0Max 10.0
+

Added: schools/2010-Lauterbad/setups/day3/pythia/minbias900.cmnd
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/pythia/minbias900.cmnd	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,21 @@
+! 1) Settings that will be used in a main program.
+Main:numberOfEvents = 100000       ! number of events to generate
+Main:timesToShow = 100             ! show how far along run is this many times
+Main:timesAllowErrors = 30         ! abort run after this many flawed events
+Main:showChangedSettings = on      ! print changed flags/modes/parameters
+#Main:showAllSettings = on         ! print all flags/modes/parameters
+Main:showChangedParticleData = on  ! print changed particle and decay data
+#Main:showAllParticleData = on     ! print all particle and decay data
+
+! 2) Beam parameter settings. Values below agree with default ones.
+Beams:idA = 2212                   ! first beam, p = 2212, pbar = -2212
+Beams:idB = 2212                   ! second beam, p = 2212, pbar = -2212
+Beams:eCM = 900                    ! CM energy of collision
+
+! 3) Pick processes and kinematics cuts
+SoftQCD:minBias on
+
+! Particles with lifetime c*tau > 10 mm are considered stable
+ParticleDecays:limitTau0 on
+ParticleDecays:tau0Max 10.0
+

Added: schools/2010-Lauterbad/setups/day3/shared/publications/ATLAS_2010_S8591806.pdf
==============================================================================
Binary file. No diff available.

Added: schools/2010-Lauterbad/setups/day3/shared/publications/CDF_2009_S8233977.pdf
==============================================================================
Binary file. No diff available.

Added: schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.aida	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,631 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">
+<aida version="3.3">
+  <implementation version="1.0" package="HepData"/>
+  <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/ATLAS_900_MINBIAS" title="(1/N)*(D(N)/DETARAP" >
+    <dataPoint>
+      <measurement value="0.0" errorPlus="0.2" errorMinus="0.2"/>
+      <measurement value="1.333" errorPlus="0.04011234224026316" errorMinus="0.04011234224026316"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/ATLAS_900_MINBIAS" title="(1/N)/DN(P=3)/DETARAP" >
+    <dataPoint>
+      <measurement value="-2.45" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/>
+      <measurement value="1.141" errorPlus="0.08728115489611718" errorMinus="0.08728115489611718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-2.3499999999999996" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/>
+      <measurement value="1.179" errorPlus="0.08229823813423955" errorMinus="0.08229823813423955"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-2.25" errorPlus="0.04999999999999982" errorMinus="0.04999999999999982"/>
+      <measurement value="1.212" errorPlus="0.0643816744112795" errorMinus="0.0643816744112795"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-2.1500000000000004" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/>
+      <measurement value="1.21" errorPlus="0.06239390995922599" errorMinus="0.06239390995922599"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-2.05" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/>
+      <measurement value="1.233" errorPlus="0.06338769596696192" errorMinus="0.06338769596696192"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.95" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.236" errorPlus="0.0643816744112795" errorMinus="0.0643816744112795"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.85" errorPlus="0.050000000000000044" errorMinus="0.04999999999999982"/>
+      <measurement value="1.254" errorPlus="0.06537583651472462" errorMinus="0.06537583651472462"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.3" errorPlus="0.06736467917239718" errorMinus="0.06736467917239718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.65" errorPlus="0.04999999999999982" errorMinus="0.050000000000000044"/>
+      <measurement value="1.323" errorPlus="0.06736467917239718" errorMinus="0.06736467917239718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.55" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.323" errorPlus="0.04751841748206689" errorMinus="0.04751841748206689"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.45" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.321" errorPlus="0.04751841748206689" errorMinus="0.04751841748206689"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999982"/>
+      <measurement value="1.317" errorPlus="0.046529560496527365" errorMinus="0.046529560496527365"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.25" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.341" errorPlus="0.0455411901469428" errorMinus="0.0455411901469428"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.15" errorPlus="0.04999999999999982" errorMinus="0.050000000000000044"/>
+      <measurement value="1.342" errorPlus="0.0455411901469428" errorMinus="0.0455411901469428"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-1.05" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.356" errorPlus="0.044553338819890924" errorMinus="0.044553338819890924"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.95" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="1.346" errorPlus="0.04257933771208754" errorMinus="0.04257933771208754"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.8500000000000001" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="1.344" errorPlus="0.04257933771208754" errorMinus="0.04257933771208754"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.341" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.6499999999999999" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="1.354" errorPlus="0.04257933771208754" errorMinus="0.04257933771208754"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.55" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="1.338" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="1.339" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.35" errorPlus="0.04999999999999999" errorMinus="0.050000000000000044"/>
+      <measurement value="1.342" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="1.342" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.15000000000000002" errorPlus="0.05000000000000002" errorMinus="0.04999999999999999"/>
+      <measurement value="1.333" errorPlus="0.040607881008493905" errorMinus="0.040607881008493905"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="-0.05" errorPlus="0.05" errorMinus="0.05"/>
+      <measurement value="1.342" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+      <measurement value="1.33" errorPlus="0.040607881008493905" errorMinus="0.040607881008493905"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.15000000000000002" errorPlus="0.04999999999999999" errorMinus="0.05000000000000002"/>
+      <measurement value="1.33" errorPlus="0.040607881008493905" errorMinus="0.040607881008493905"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.25" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="1.322" errorPlus="0.04044749683231337" errorMinus="0.04044749683231337"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.35" errorPlus="0.050000000000000044" errorMinus="0.04999999999999999"/>
+      <measurement value="1.335" errorPlus="0.04044749683231337" errorMinus="0.04044749683231337"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.45" errorPlus="0.04999999999999999" errorMinus="0.04999999999999999"/>
+      <measurement value="1.341" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="1.327" errorPlus="0.040607881008493905" errorMinus="0.040607881008493905"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.6499999999999999" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="1.336" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.348" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.8500000000000001" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="1.335" errorPlus="0.041593268686170845" errorMinus="0.041593268686170845"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.95" errorPlus="0.050000000000000044" errorMinus="0.04999999999999993"/>
+      <measurement value="1.34" errorPlus="0.04257933771208754" errorMinus="0.04257933771208754"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.05" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.349" errorPlus="0.043566041821583926" errorMinus="0.043566041821583926"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.15" errorPlus="0.050000000000000044" errorMinus="0.04999999999999982"/>
+      <measurement value="1.359" errorPlus="0.046529560496527365" errorMinus="0.046529560496527365"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.25" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.35" errorPlus="0.046529560496527365" errorMinus="0.046529560496527365"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.35" errorPlus="0.04999999999999982" errorMinus="0.050000000000000044"/>
+      <measurement value="1.338" errorPlus="0.04751841748206689" errorMinus="0.04751841748206689"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.45" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.325" errorPlus="0.04751841748206689" errorMinus="0.04751841748206689"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.55" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.316" errorPlus="0.04751841748206689" errorMinus="0.04751841748206689"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.65" errorPlus="0.050000000000000044" errorMinus="0.04999999999999982"/>
+      <measurement value="1.319" errorPlus="0.06736467917239718" errorMinus="0.06736467917239718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.75" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.303" errorPlus="0.06736467917239718" errorMinus="0.06736467917239718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.85" errorPlus="0.04999999999999982" errorMinus="0.050000000000000044"/>
+      <measurement value="1.281" errorPlus="0.06736467917239718" errorMinus="0.06736467917239718"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.95" errorPlus="0.050000000000000044" errorMinus="0.050000000000000044"/>
+      <measurement value="1.238" errorPlus="0.06537583651472462" errorMinus="0.06537583651472462"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.05" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/>
+      <measurement value="1.234" errorPlus="0.0643816744112795" errorMinus="0.0643816744112795"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.1500000000000004" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/>
+      <measurement value="1.237" errorPlus="0.06338769596696192" errorMinus="0.06338769596696192"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.25" errorPlus="0.04999999999999982" errorMinus="0.04999999999999982"/>
+      <measurement value="1.225" errorPlus="0.0643816744112795" errorMinus="0.0643816744112795"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.3499999999999996" errorPlus="0.050000000000000266" errorMinus="0.04999999999999982"/>
+      <measurement value="1.179" errorPlus="0.08130190649671139" errorMinus="0.08130190649671139"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.45" errorPlus="0.04999999999999982" errorMinus="0.050000000000000266"/>
+      <measurement value="1.157" errorPlus="0.08728115489611718" errorMinus="0.08728115489611718"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/ATLAS_900_MINBIAS" title="(1/N)*(1/(2*PI*PT))*D2(N)/DETARAP/DPT IN GEV**-2" >
+    <dataPoint>
+      <measurement value="0.55" errorPlus="0.04999999999999993" errorMinus="0.050000000000000044"/>
+      <measurement value="0.961" errorPlus="0.08241329019035705" errorMinus="0.08241329019035705"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.648" errorPlus="0.051999999999999935" errorMinus="0.04800000000000004"/>
+      <measurement value="0.59" errorPlus="0.02522023195769619" errorMinus="0.02522023195769619"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.748" errorPlus="0.052000000000000046" errorMinus="0.04800000000000004"/>
+      <measurement value="0.374" errorPlus="0.01581746025757612" errorMinus="0.01581746025757612"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.848" errorPlus="0.052000000000000046" errorMinus="0.04799999999999993"/>
+      <measurement value="0.245" errorPlus="0.010315320838442206" errorMinus="0.010315320838442206"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="0.948" errorPlus="0.052000000000000046" errorMinus="0.04799999999999993"/>
+      <measurement value="0.162" errorPlus="0.006813708535005001" errorMinus="0.006813708535005001"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.048" errorPlus="0.052000000000000046" errorMinus="0.04800000000000004"/>
+      <measurement value="0.111" errorPlus="0.004632347568997819" errorMinus="0.004632347568997819"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.148" errorPlus="0.052000000000000046" errorMinus="0.04799999999999982"/>
+      <measurement value="0.0762" errorPlus="0.0031911892767430765" errorMinus="0.0031911892767430765"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.248" errorPlus="0.052000000000000046" errorMinus="0.04800000000000004"/>
+      <measurement value="0.0537" errorPlus="0.002250294425180847" errorMinus="0.002250294425180847"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.348" errorPlus="0.051999999999999824" errorMinus="0.04800000000000004"/>
+      <measurement value="0.0382" errorPlus="0.0015994924194881325" errorMinus="0.0015994924194881325"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.448" errorPlus="0.052000000000000046" errorMinus="0.04800000000000004"/>
+      <measurement value="0.0275" errorPlus="0.0011488098188995426" errorMinus="0.0011488098188995426"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.549" errorPlus="0.051000000000000156" errorMinus="0.04899999999999993"/>
+      <measurement value="0.0203" errorPlus="8.531992733236475E-4" errorMinus="8.531992733236475E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.648" errorPlus="0.052000000000000046" errorMinus="0.04799999999999982"/>
+      <measurement value="0.0148" errorPlus="6.26647237287455E-4" errorMinus="6.26647237287455E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.748" errorPlus="0.052000000000000046" errorMinus="0.04800000000000004"/>
+      <measurement value="0.0111" errorPlus="4.6918572015780703E-4" errorMinus="4.6918572015780703E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.849" errorPlus="0.050999999999999934" errorMinus="0.04899999999999993"/>
+      <measurement value="0.00838" errorPlus="3.5771776584340905E-4" errorMinus="3.5771776584340905E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="1.949" errorPlus="0.050999999999999934" errorMinus="0.049000000000000155"/>
+      <measurement value="0.00638" errorPlus="2.7533014727777267E-4" errorMinus="2.7533014727777267E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.049" errorPlus="0.051000000000000156" errorMinus="0.04899999999999993"/>
+      <measurement value="0.00484" errorPlus="2.109621767047354E-4" errorMinus="2.109621767047354E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.149" errorPlus="0.051000000000000156" errorMinus="0.04899999999999993"/>
+      <measurement value="0.00375" errorPlus="1.665918365346874E-4" errorMinus="1.665918365346874E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.249" errorPlus="0.05099999999999971" errorMinus="0.04899999999999993"/>
+      <measurement value="0.00291" errorPlus="1.3126400877620643E-4" errorMinus="1.3126400877620643E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.349" errorPlus="0.05099999999999971" errorMinus="0.04900000000000038"/>
+      <measurement value="0.00224" errorPlus="1.0249297536904663E-4" errorMinus="1.0249297536904663E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.449" errorPlus="0.051000000000000156" errorMinus="0.04899999999999993"/>
+      <measurement value="0.00182" errorPlus="8.519776992386597E-5" errorMinus="8.519776992386597E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.617" errorPlus="0.133" errorMinus="0.11699999999999999"/>
+      <measurement value="0.0012" errorPlus="5.806375805956759E-5" errorMinus="5.806375805956759E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.867" errorPlus="0.133" errorMinus="0.11699999999999999"/>
+      <measurement value="7.0E-4" errorPlus="3.630440744592866E-5" errorMinus="3.630440744592866E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.116" errorPlus="0.1339999999999999" errorMinus="0.1160000000000001"/>
+      <measurement value="4.39E-4" errorPlus="2.4796614688299692E-5" errorMinus="2.4796614688299692E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.366" errorPlus="0.1339999999999999" errorMinus="0.1160000000000001"/>
+      <measurement value="2.6E-4" errorPlus="1.6581447463958023E-5" errorMinus="1.6581447463958023E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.72" errorPlus="0.2799999999999998" errorMinus="0.2200000000000002"/>
+      <measurement value="1.33E-4" errorPlus="1.0078397690109276E-5" errorMinus="1.0078397690109276E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.212" errorPlus="0.28800000000000026" errorMinus="0.21199999999999974"/>
+      <measurement value="5.72E-5" errorPlus="5.735477312308017E-6" errorMinus="5.735477312308017E-6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.726" errorPlus="0.274" errorMinus="0.22599999999999998"/>
+      <measurement value="2.55E-5" errorPlus="3.4504637369489916E-6" errorMinus="3.4504637369489916E-6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="5.407" errorPlus="0.593" errorMinus="0.40700000000000003"/>
+      <measurement value="9.53E-6" errorPlus="1.8775782273982621E-6" errorMinus="1.8775782273982621E-6"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="6.408" errorPlus="0.5919999999999996" errorMinus="0.40800000000000036"/>
+      <measurement value="3.1E-6" errorPlus="9.576137008209523E-7" errorMinus="9.576137008209523E-7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="7.425" errorPlus="0.5750000000000002" errorMinus="0.4249999999999998"/>
+      <measurement value="1.04E-6" errorPlus="4.894374321606389E-7" errorMinus="4.894374321606389E-7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="8.73" errorPlus="1.2699999999999996" errorMinus="0.7300000000000004"/>
+      <measurement value="3.46E-7" errorPlus="2.2541597104020824E-7" errorMinus="2.2541597104020824E-7"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="11.704" errorPlus="3.2959999999999994" errorMinus="1.7040000000000006"/>
+      <measurement value="5.01E-8" errorPlus="4.524522074208502E-8" errorMinus="4.524522074208502E-8"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="16.25" errorPlus="3.75" errorMinus="1.25"/>
+      <measurement value="4.57E-9" errorPlus="4.5744617169673645E-9" errorMinus="4.5744617169673645E-9"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/ATLAS_900_MINBIAS" title="(1/N)*D(N)/DN(P=3)" >
+    <dataPoint>
+      <measurement value="1.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.131" errorPlus="0.01772348001945442" errorMinus="0.01772348001945442"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.129" errorPlus="0.009343833046453686" errorMinus="0.009343833046453686"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.118" errorPlus="0.005118078154932767" errorMinus="0.005118078154932767"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0999" errorPlus="0.0032858052285550947" errorMinus="0.0032858052285550947"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="5.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0822" errorPlus="0.00258744835697256" errorMinus="0.00258744835697256"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="6.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0679" errorPlus="0.002213437146159791" errorMinus="0.002213437146159791"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="7.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0571" errorPlus="0.0019963990082145404" errorMinus="0.0019963990082145404"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="8.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0482" errorPlus="0.0018296447742663056" errorMinus="0.0018296447742663056"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="9.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0409" errorPlus="0.0016925826420000885" errorMinus="0.0016925826420000885"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="10.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0349" errorPlus="0.0015670548809789656" errorMinus="0.0015670548809789656"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="11.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0298" errorPlus="0.0014613541665181647" errorMinus="0.0014613541665181647"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="12.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0256" errorPlus="0.0013763051987113904" errorMinus="0.0013763051987113904"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="13.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.022" errorPlus="0.001310430845180317" errorMinus="0.001310430845180317"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="14.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0187" errorPlus="0.0012354901861204724" errorMinus="0.0012354901861204724"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="15.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0157" errorPlus="0.0011516179922179055" errorMinus="0.0011516179922179055"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="16.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0132" errorPlus="0.0010487783369234893" errorMinus="0.0010487783369234893"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="17.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.011" errorPlus="9.437255957109566E-4" errorMinus="9.437255957109566E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="18.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00937" errorPlus="8.67332116319925E-4" errorMinus="8.67332116319925E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="19.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0080" errorPlus="8.01870313704155E-4" errorMinus="8.01870313704155E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="20.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00678" errorPlus="7.364190383198957E-4" errorMinus="7.364190383198957E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="21.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.0057" errorPlus="6.767547561709486E-4" errorMinus="6.767547561709486E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="22.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00478" errorPlus="6.168646529020771E-4" errorMinus="6.168646529020771E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="23.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00395" errorPlus="5.513565089848854E-4" errorMinus="5.513565089848854E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="24.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00324" errorPlus="4.8654393429576324E-4" errorMinus="4.8654393429576324E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="25.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00264" errorPlus="4.2419924563817885E-4" errorMinus="4.2419924563817885E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="26.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00215" errorPlus="3.675210062023666E-4" errorMinus="3.675210062023666E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="27.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00177" errorPlus="3.216031871732617E-4" errorMinus="3.216031871732617E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="28.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00145" errorPlus="2.789134095019456E-4" errorMinus="2.789134095019456E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="29.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.00122" errorPlus="2.481328071819605E-4" errorMinus="2.481328071819605E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="30.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="9.95E-4" errorPlus="2.1524404753674373E-4" errorMinus="2.1524404753674373E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="31.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="7.82E-4" errorPlus="1.791439644531738E-4" errorMinus="1.791439644531738E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="32.846" errorPlus="1.6540000000000035" errorMinus="1.3459999999999965"/>
+      <measurement value="5.18E-4" errorPlus="1.3020096005790434E-4" errorMinus="1.3020096005790434E-4"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="36.623" errorPlus="2.8770000000000024" errorMinus="2.1229999999999976"/>
+      <measurement value="2.23E-4" errorPlus="6.688086422886593E-5" errorMinus="6.688086422886593E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="42.492" errorPlus="7.008000000000003" errorMinus="2.9919999999999973"/>
+      <measurement value="4.61E-5" errorPlus="1.9833809518093087E-5" errorMinus="1.9833809518093087E-5"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="52.869" errorPlus="7.631" errorMinus="3.3689999999999998"/>
+      <measurement value="4.07E-6" errorPlus="4.234548382059177E-6" errorMinus="4.234548382059177E-6"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/ATLAS_900_MINBIAS" title="MEAN(NAME=PT) IN GEV" >
+    <dataPoint>
+      <measurement value="1.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.753" errorPlus="0.007687424796380125" errorMinus="0.007687424796380125"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="2.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.767" errorPlus="0.007984747960956565" errorMinus="0.007984747960956565"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.785" errorPlus="0.008402142583888945" errorMinus="0.008402142583888945"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.802" errorPlus="0.008809250819451107" errorMinus="0.008809250819451107"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="5.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.819" errorPlus="0.009161528256792095" errorMinus="0.009161528256792095"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="6.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.835" errorPlus="0.009436699634935933" errorMinus="0.009436699634935933"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="7.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.849" errorPlus="0.009652398665616749" errorMinus="0.009652398665616749"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="8.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.862" errorPlus="0.009781492728617652" errorMinus="0.009781492728617652"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="9.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.874" errorPlus="0.009861465408345758" errorMinus="0.009861465408345758"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="10.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.884" errorPlus="0.009882550278141772" errorMinus="0.009882550278141772"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="11.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.893" errorPlus="0.00991717701767998" errorMinus="0.00991717701767998"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="12.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.901" errorPlus="0.009970842491986322" errorMinus="0.009970842491986322"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="13.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.908" errorPlus="0.010038451075738726" errorMinus="0.010038451075738726"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="14.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.915" errorPlus="0.010148201811158466" errorMinus="0.010148201811158466"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="15.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.922" errorPlus="0.010290330412576654" errorMinus="0.010290330412576654"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="16.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.928" errorPlus="0.010320096898769894" errorMinus="0.010320096898769894"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="17.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.934" errorPlus="0.010451717562200005" errorMinus="0.010451717562200005"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="18.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.939" errorPlus="0.010589693102257497" errorMinus="0.010589693102257497"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="19.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.945" errorPlus="0.010734784580977859" errorMinus="0.010734784580977859"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="20.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.95" errorPlus="0.010890459127144274" errorMinus="0.010890459127144274"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="21.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.956" errorPlus="0.010968085521183724" errorMinus="0.010968085521183724"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="22.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.96" errorPlus="0.011061663527697812" errorMinus="0.011061663527697812"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="23.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.965" errorPlus="0.011169874663576131" errorMinus="0.011169874663576131"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="24.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.969" errorPlus="0.011301433537388078" errorMinus="0.011301433537388078"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="25.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.973" errorPlus="0.011459410979627183" errorMinus="0.011459410979627183"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="26.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.977" errorPlus="0.01173718876051672" errorMinus="0.01173718876051672"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="27.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.98" errorPlus="0.0120491327488745" errorMinus="0.0120491327488745"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="28.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.984" errorPlus="0.012500575986729571" errorMinus="0.012500575986729571"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="29.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.988" errorPlus="0.012975684953018858" errorMinus="0.012975684953018858"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="30.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.993" errorPlus="0.01366238998125877" errorMinus="0.01366238998125877"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="31.0" errorPlus="0.5" errorMinus="0.5"/>
+      <measurement value="0.997" errorPlus="0.014211702923998938" errorMinus="0.014211702923998938"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="32.789" errorPlus="1.7109999999999985" errorMinus="1.2890000000000015"/>
+      <measurement value="1.01" errorPlus="0.01588112401563567" errorMinus="0.01588112401563567"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="36.442" errorPlus="3.058" errorMinus="1.9420000000000002"/>
+      <measurement value="1.02" errorPlus="0.021721012867727877" errorMinus="0.021721012867727877"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="42.393" errorPlus="7.106999999999999" errorMinus="2.8930000000000007"/>
+      <measurement value="1.04" errorPlus="0.05261796271236658" errorMinus="0.05261796271236658"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="52.75" errorPlus="7.75" errorMinus="3.25"/>
+      <measurement value="1.04" errorPlus="0.12028054705562326" errorMinus="0.12028054705562326"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Added: schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.cc	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,51 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Tools/Logging.hh"
+
+namespace Rivet {
+
+
+  /// @brief ATLAS minimum bias analysis at 900 GeV
+  class ATLAS_900_MINBIAS : public Analysis {
+  public:
+
+    ATLAS_900_MINBIAS() : Analysis("ATLAS_900_MINBIAS") {
+      setNeedsCrossSection(false);
+      setBeams(PROTON, PROTON);
+    }
+
+
+    void init() {
+      // This method is called once before the first event. You can book
+      // your histograms here and initialise a final state projection.
+    }
+
+
+    void analyze(const Event& event) {
+      // This method is called for each individual event. Analyse the
+      // event here and fill the the proper histograms. Be careful that
+      // generators might use weighted events, so you will have to
+      // propagate the weights through your code.
+    }
+
+
+    void finalize() {
+      // This method is called only once after analysing all events.
+      // That's the place where you need to take care of proper
+      // histogram normalisation.
+    }
+
+
+  private:
+    // Put your private variable declarations here.
+
+  };
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<ATLAS_900_MINBIAS> plugin_ATLAS_900_MINBIAS;
+
+}
+

Added: schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.info	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,40 @@
+Name: ATLAS_900_MINBIAS
+Year: 2010
+Summary: Charged particles at 900 GeV in ATLAS
+Experiment: ATLAS
+Collider: LHC 900GeV
+SpiresID: 8591806
+Status: PLAYGROUND
+Authors:
+ - MCnet school participants
+References:
+ - arXiv:1003.3124
+RunInfo:
+  pp QCD interactions at 900 GeV including diffractive events.
+NumEvents: 1000000
+Beams: [p+, p+]
+Energies: [900]
+PtCuts: [0]
+Description:
+  The first measurements with the ATLAS detector at the LHC.
+  Data were collected using a minimum-bias trigger in December 2009 during
+  proton-proton collisions at a centre of mass energy of 900 GeV. The charged-
+  particle density, its dependence on transverse momentum and pseudorapid-
+  ity, and the relationship between transverse momentum and charged-particle
+  multiplicity are measured for events with at least one charged particle in
+  the kinematic range $|\eta| < 2.5$ and $p_\perp > 500$ MeV.
+  All data is corrected to the particle level.
+BibKey: Aad:2010rd
+BibTeX: '@Article{Aad:2010rd,
+     author    = "Aad, G. and others",
+ collaboration = "ALTAS",
+     title     = "{Charged-particle multiplicities in pp interactions at
+                  sqrt(s) = 900 GeV measured with the ATLAS detector at the
+                  LHC}",
+     year      = "2010",
+     eprint    = "1003.3124",
+     archivePrefix = "arXiv",
+     primaryClass  =  "hep-ex",
+     doi       = "10.1016/j.physletb.2010.03.064",
+     SLACcitation  = "%%CITATION = 1003.3124;%%"
+}'

Added: schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/rivet/ATLAS_900_MINBIAS.plot	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,32 @@
+# BEGIN PLOT /ATLAS_900_MINBIAS/d02-x01-y01
+Title=Charged particle multiplicity as function of $\eta$
+XLabel=$\eta$
+YLabel=$1/N_\mathrm{ev} \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta$
+LogY=0
+YMin=0.75
+YMax=1.55
+# END PLOT
+
+# BEGIN PLOT /ATLAS_900_MINBIAS/d03-x01-y01
+Title=Charged particle multiplicity as function of $p_\perp$
+XLabel=$p_\perp$ [GeV]
+YLabel=$1/N_\mathrm{ev} \; 1/(2\pi p_\perp) \; \mathrm{d}^2N_\mathrm{ch}/\mathrm{d}\eta\mathrm{d}p_\perp$ [GeV$^{-2}$]
+YMin=2e-11
+YMax=10
+LogX=1
+# END PLOT
+
+# BEGIN PLOT /ATLAS_900_MINBIAS/d04-x01-y01
+Title=Charged particle density
+XLabel=$N_\mathrm{ch}$
+YLabel=$1/N_\mathrm{ev} \mathrm{d}N_\mathrm{ev}/\mathrm{d}N_\mathrm{ch}$
+# END PLOT
+
+# BEGIN PLOT /ATLAS_900_MINBIAS/d05-x01-y01
+Title=Average transverse momentum as function of $N_\mathrm{ch}$
+XLabel=$N_\mathrm{ch}$
+YLabel=$\langle p_\perp \rangle$ [GeV]
+YMin=0.55
+YMax=1.3
+LogY=0
+# END PLOT

Added: schools/2010-Lauterbad/setups/day3/shared/rivet/CDF_2009_S8233977.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/rivet/CDF_2009_S8233977.cc	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,122 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "Rivet/Projections/TriggerCDFRun2.hh"
+
+namespace Rivet {
+
+
+  /* @brief CDF Run II min-bias cross-section
+   *
+   * Measurement of \f$ \langle p_T \rangle \f$ vs. \f$ n_\text{ch} \f$,
+   * the track \f$ p_T \f$ distribution, and the \f$ \sum E_T \f$ distribution.
+   * Particles are selected within |eta|<1 and with pT>0.4 GeV.
+   * There is no pT cut for the \f$ \sum E_T \f$ measurement.
+   *
+   * @par Run conditions
+   *
+   * @arg \f$ \sqrt{s} = \f$ 1960 GeV
+   * @arg Run with generic QCD events.
+   * @arg Set particles with c*tau > 10 mm stable
+   *
+   */
+  class CDF_2009_S8233977 : public Analysis {
+  public:
+
+    /// Constructor
+    CDF_2009_S8233977()
+      : Analysis("CDF_2009_S8233977"),
+        _sumWeightSelected(0.0)
+    {
+      setBeams(PROTON, ANTIPROTON);
+      setNeedsCrossSection(true);
+    }
+ 
+ 
+    /// Book histograms and projections
+    void init() {
+      addProjection(TriggerCDFRun2(), "Trigger");
+
+      // -1 < eta < 1, pT cuts at 0.0 and 0.4, respectively
+      addProjection(FinalState(-1.0, 1.0, 0.0*GeV), "EtFS");
+      addProjection(ChargedFinalState(-1.0, 1.0, 0.4*GeV), "CFS");
+
+      _hist_pt = bookHistogram1D(1, 1, 1);               // d01-x01-x01
+      _hist_pt_vs_multiplicity = bookProfile1D(2, 1, 1); // d02-x01-x01
+      _hist_sumEt = bookHistogram1D(3, 1, 1);            // d03-x01-x01
+    }
+
+
+ 
+    /// Do the analysis
+    void analyze(const Event& evt) {
+      // MinBias Trigger
+      const bool trigger = applyProjection<TriggerCDFRun2>(evt, "Trigger").minBiasDecision();
+      if (!trigger) vetoEvent;
+
+      // Get the event weight
+      const double weight = evt.weight();
+
+      const ChargedFinalState& trackfs = applyProjection<ChargedFinalState>(evt, "CFS");
+      const size_t numParticles = trackfs.size();
+      foreach (const Particle& p, trackfs.particles()) {
+        const double pT = p.momentum().pT() / GeV;
+        _hist_pt_vs_multiplicity->fill(numParticles, pT, weight);
+     
+        // The weight for entries in the pT distribution should be weight/(pT*dPhi*dy).
+        //
+        // - dPhi = 2*PI
+        //
+        // - dy depends on the pT: They calculate y assuming the particle has the
+        //   pion mass and assuming that eta=1:
+        //   dy = 2 * 1/2 * ln [(sqrt(m^2 + (a+1)*pT^2) + a*pT) / (sqrt(m^2 + (a+1)*pT^2) - a*pT)]
+        //   with a = sinh(1).
+        //
+        // sinh(1) = 1.1752012
+        // m(charged pion)^2 = (139.57 MeV)^2 = 0.019479785 GeV^2
+        const double sinh1 = 1.1752012;
+        const double apT  = sinh1 * pT;
+        const double mPi = 139.57*MeV;
+        const double root = sqrt(mPi*mPi + (1+sinh1)*pT*pT);
+        const double dy = std::log((root+apT)/(root-apT));
+        const double dphi = TWOPI;
+        _hist_pt->fill(pT, weight/(pT*dphi*dy));
+      }
+
+      // Calc sum(Et) from calo particles
+      const FinalState& etfs = applyProjection<FinalState>(evt, "EtFS");
+      double sumEt = 0.0;
+      foreach (const Particle& p, etfs.particles()) {
+        sumEt += p.momentum().Et();
+      }
+      _hist_sumEt->fill(sumEt, weight);
+      _sumWeightSelected += evt.weight();
+    }
+
+ 
+    /// Normalize histos
+    void finalize() {
+      scale(_hist_sumEt, crossSection()/millibarn/(4*M_PI*_sumWeightSelected));
+      scale(_hist_pt, crossSection()/millibarn/_sumWeightSelected);
+      getLog() << Log::DEBUG << "sumOfWeights()     = " << sumOfWeights() << std::endl;
+      getLog() << Log::DEBUG << "_sumWeightSelected = " << _sumWeightSelected << std::endl;
+    }
+ 
+
+  private:
+
+    double _sumWeightSelected;
+    AIDA::IProfile1D *_hist_pt_vs_multiplicity;
+    AIDA::IHistogram1D *_hist_pt;
+    AIDA::IHistogram1D *_hist_sumEt;
+
+  };
+
+
+  // This global object acts as a hook for the plugin system
+  AnalysisBuilder<CDF_2009_S8233977> plugin_CDF_2009_S8233977;
+
+}

Added: schools/2010-Lauterbad/setups/day3/shared/rivet/Makefile
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/rivet/Makefile	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,13 @@
+
+shared_flags = -shared -fPIC
+# on Mac OSX you want to use these flags instead:
+#shared_flags = -undefined dynamic_lookup -bundle
+all: RivetATLAS_900_MINBIAS.so
+
+RivetATLAS_900_MINBIAS.so: ATLAS_900_MINBIAS.cc
+	g++ $(shared_flags) `rivet-config --cppflags` ATLAS_900_MINBIAS.cc -o RivetATLAS_900_MINBIAS.so
+
+install: all
+
+clean:
+	rm -f RivetATLAS_900_MINBIAS.so

Added: schools/2010-Lauterbad/setups/day3/shared/setup.sh
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/shared/setup.sh	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,5 @@
+
+export RIVET_ANALYSIS_PATH=/opt/lib:/home/mcnet/day3/shared/rivet
+export RIVET_REF_PATH=/opt/share/Rivet:/home/mcnet/day3/shared/rivet
+export RIVET_INFO_PATH=/opt/share/Rivet:/home/mcnet/day3/shared/rivet
+

Added: schools/2010-Lauterbad/setups/day3/sherpa/Run.1960.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/sherpa/Run.1960.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,46 @@
+(run){
+  EVENTS = 10000
+
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.1960
+
+  ERROR = 0.1
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 980;
+  BEAM_2 = -2212; BEAM_ENERGY_2 = 980;
+
+  K_PERP_MEAN_1 = 1.23
+  K_PERP_MEAN_2 = 1.23
+  K_PERP_SIGMA_1 = 0.908
+  K_PERP_SIGMA_2 = 0.908
+}(beam)
+
+(processes){
+  Process 93 93 -> 93 93
+  Order_EW 0
+  End process;
+}(processes)
+
+(selector){
+  NJetFinder  2  2.034354851  0.0  1.0
+}(selector)
+
+(me){
+  ME_SIGNAL_GENERATOR = Comix
+  EVENT_GENERATION_MODE = Weighted
+  RESULT_DIRECTORY = Results.1960
+}(me)
+
+(mi){
+  MI_HANDLER = Amisic
+  SCALE_MIN = 2.0
+  RESCALE_EXPONENT = 0.2
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+    -a CDF_2009_S8233977
+  } END_RIVET
+}(analysis)

Added: schools/2010-Lauterbad/setups/day3/sherpa/Run.900.dat
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ schools/2010-Lauterbad/setups/day3/sherpa/Run.900.dat	Tue Jul 17 14:44:00 2012	(r3834)
@@ -0,0 +1,46 @@
+(run){
+  EVENTS = 10000
+
+  ANALYSIS = Rivet
+  ANALYSIS_OUTPUT = Analysis.900
+
+  ERROR = 0.1
+}(run)
+
+(beam){
+  BEAM_1 =  2212; BEAM_ENERGY_1 = 450;
+  BEAM_2 =  2212; BEAM_ENERGY_2 = 450;
+
+  K_PERP_MEAN_1 = 1.23
+  K_PERP_MEAN_2 = 1.23
+  K_PERP_SIGMA_1 = 0.908
+  K_PERP_SIGMA_2 = 0.908
+}(beam)
+
+(processes){
+  Process 93 93 -> 93 93
+  Order_EW 0
+  End process;
+}(processes)
+
+(selector){
+  NJetFinder  2  1.74110112659  0.0  1.0
+}(selector)
+
+(me){
+  ME_SIGNAL_GENERATOR = Comix
+  EVENT_GENERATION_MODE = Weighted
+  RESULT_DIRECTORY = Results.900
+}(me)
+
+(mi){
+  MI_HANDLER = Amisic
+  SCALE_MIN = 2.0
+  RESCALE_EXPONENT = 0.2
+}(mi)
+
+(analysis){
+  BEGIN_RIVET {
+    -a ATLAS_900_MINBIAS
+  } END_RIVET
+}(analysis)


More information about the Rivet-svn mailing list