|
[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/sherpablackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue 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 |