|
|
[Rivet-svn] r3830 - in schools/2011-Kyoto/setups: . day1 day1/herwig day1/pythia 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/rivet-plugin day3/shared day3/sherpa day3/sherpa/ttbar day3/sherpa/wbbjets
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Tue Jul 17 14:18:40 BST 2012
Author: hoeth
Date: Tue Jul 17 14:18:39 2012
New Revision: 3830
Log:
Kyoto setups
Added:
schools/2011-Kyoto/setups/
schools/2011-Kyoto/setups/day1/
schools/2011-Kyoto/setups/day1/herwig/
schools/2011-Kyoto/setups/day1/pythia/
schools/2011-Kyoto/setups/day1/sherpa/
schools/2011-Kyoto/setups/day1/sherpa/Run.15-3.dat
schools/2011-Kyoto/setups/day1/sherpa/Run.20-3.NLO.dat
schools/2011-Kyoto/setups/day1/sherpa/Run.20-3.dat
schools/2011-Kyoto/setups/day1/sherpa/Run.dat
schools/2011-Kyoto/setups/day2/
schools/2011-Kyoto/setups/day2/herwig/
schools/2011-Kyoto/setups/day2/herwig/ME/
schools/2011-Kyoto/setups/day2/herwig/ME/Base.in
schools/2011-Kyoto/setups/day2/herwig/ME/TVT-0jet.in
schools/2011-Kyoto/setups/day2/herwig/ME/TVT-1jet.in
schools/2011-Kyoto/setups/day2/herwig/ME/TVT-Powheg.in
schools/2011-Kyoto/setups/day2/herwig/QED/
schools/2011-Kyoto/setups/day2/herwig/QED/Base.in
schools/2011-Kyoto/setups/day2/herwig/QED/TVT-0jet-QED-off.in
schools/2011-Kyoto/setups/day2/herwig/QED/TVT-0jet-QED-on.in
schools/2011-Kyoto/setups/day2/herwig/QED/extract_noclus.sh (contents, props changed)
schools/2011-Kyoto/setups/day2/herwig/shower/
schools/2011-Kyoto/setups/day2/herwig/shower/Base.in
schools/2011-Kyoto/setups/day2/herwig/shower/TVT-0jet-MEC.in
schools/2011-Kyoto/setups/day2/herwig/shower/TVT-0jet-soft.in
schools/2011-Kyoto/setups/day2/herwig/shower/TVT-1jet.in
schools/2011-Kyoto/setups/day2/herwig/shower/TVT-Powheg.in
schools/2011-Kyoto/setups/day2/pythia/
schools/2011-Kyoto/setups/day2/pythia/ME/
schools/2011-Kyoto/setups/day2/pythia/ME/0jet.cmnd
schools/2011-Kyoto/setups/day2/pythia/ME/1jet.cmnd
schools/2011-Kyoto/setups/day2/pythia/ME/pythia_rivet (contents, props changed)
schools/2011-Kyoto/setups/day2/pythia/QED/
schools/2011-Kyoto/setups/day2/pythia/QED/QED.cmnd
schools/2011-Kyoto/setups/day2/pythia/QED/noQED.cmnd
schools/2011-Kyoto/setups/day2/pythia/QED/pythia_rivet (contents, props changed)
schools/2011-Kyoto/setups/day2/pythia/Zpythia.cmnd
schools/2011-Kyoto/setups/day2/pythia/shower/
schools/2011-Kyoto/setups/day2/pythia/shower/0jet.cmnd
schools/2011-Kyoto/setups/day2/pythia/shower/1jet.cmnd
schools/2011-Kyoto/setups/day2/pythia/shower/pythia_rivet (contents, props changed)
schools/2011-Kyoto/setups/day2/shared/
schools/2011-Kyoto/setups/day2/shared/CDF_2008_S7540469.pdf (contents, props changed)
schools/2011-Kyoto/setups/day2/shared/D0_2007_S7075677.pdf (contents, props changed)
schools/2011-Kyoto/setups/day2/shared/D0_2008_S7554427.pdf (contents, props changed)
schools/2011-Kyoto/setups/day2/shared/D0_2008_S7863608.pdf (contents, props changed)
schools/2011-Kyoto/setups/day2/sherpa/
schools/2011-Kyoto/setups/day2/sherpa/ME/
schools/2011-Kyoto/setups/day2/sherpa/ME/Run.0jet.dat
schools/2011-Kyoto/setups/day2/sherpa/ME/Run.1jet.dat
schools/2011-Kyoto/setups/day2/sherpa/ME/Run.2jet.dat
schools/2011-Kyoto/setups/day2/sherpa/ME/Run.3jet.dat
schools/2011-Kyoto/setups/day2/sherpa/QED/
schools/2011-Kyoto/setups/day2/sherpa/QED/Run.ME.dat
schools/2011-Kyoto/setups/day2/sherpa/QED/Run.Off.dat
schools/2011-Kyoto/setups/day2/sherpa/merging/
schools/2011-Kyoto/setups/day2/sherpa/merging/Run.0jet.dat
schools/2011-Kyoto/setups/day2/sherpa/merging/Run.1jet.dat
schools/2011-Kyoto/setups/day2/sherpa/merging/Run.2jet.dat
schools/2011-Kyoto/setups/day2/sherpa/merging/Run.3jet.dat
schools/2011-Kyoto/setups/day3/
schools/2011-Kyoto/setups/day3/herwig/
schools/2011-Kyoto/setups/day3/herwig/TopBackground.in
schools/2011-Kyoto/setups/day3/herwig/TopSignal.in
schools/2011-Kyoto/setups/day3/pythia/
schools/2011-Kyoto/setups/day3/pythia/pythia_rivet (contents, props changed)
schools/2011-Kyoto/setups/day3/pythia/topbackground.cmnd
schools/2011-Kyoto/setups/day3/pythia/topsignal.cmnd
schools/2011-Kyoto/setups/day3/rivet-plugin/
schools/2011-Kyoto/setups/day3/rivet-plugin/MC_TOP.cc
schools/2011-Kyoto/setups/day3/rivet-plugin/MC_TOP.plot
schools/2011-Kyoto/setups/day3/rivet-plugin/Makefile
schools/2011-Kyoto/setups/day3/rivet-plugin/Sherpa_References.tex
schools/2011-Kyoto/setups/day3/shared/
schools/2011-Kyoto/setups/day3/shared/Wposbbjj.lhe.bz2 (contents, props changed)
schools/2011-Kyoto/setups/day3/sherpa/
schools/2011-Kyoto/setups/day3/sherpa/ttbar/
schools/2011-Kyoto/setups/day3/sherpa/ttbar/Run.dat
schools/2011-Kyoto/setups/day3/sherpa/wbbjets/
schools/2011-Kyoto/setups/day3/sherpa/wbbjets/Run.dat
Added: schools/2011-Kyoto/setups/day1/sherpa/Run.15-3.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day1/sherpa/Run.15-3.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,124 @@
+(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
+ !
+ BEAM_2 = -2212 ! possible beam particles: P+, P-, e+, e-
+ BEAM_ENERGY_2 = 980. ! in GeV
+}(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
+}(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.63 ! minimum scale
+ RESCALE_EXPONENT = 0.192 ! rescaling exponent
+}(mi)
+
+(analysis){
+ BEGIN_RIVET {
+ -a MC_WJETS
+ } END_RIVET
+}(analysis)
Added: schools/2011-Kyoto/setups/day1/sherpa/Run.20-3.NLO.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day1/sherpa/Run.20-3.NLO.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,131 @@
+(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.NLO
+}(run)
+
+(beam){
+ !-------------------------------------------------------
+ !-- Beam parameters ------------------------------------
+ !-------------------------------------------------------
+ BEAM_1 = 2212 ! possible beam particles: P+, P-, e+, e-
+ BEAM_ENERGY_1 = 980. ! in GeV
+ !
+ BEAM_2 = -2212 ! possible beam particles: P+, P-, e+, e-
+ BEAM_ENERGY_2 = 980. ! in GeV
+}(beam)
+
+(processes){
+ !-------------------------------------------------------
+ !-- Processes to calculate -----------------------------
+ !-------------------------------------------------------
+ !
+ ! jet jet -> electron neutrino + up to 3 jets
+ !
+ Process : 93 93 -> 11 -12 93{3}
+ NLO_QCD_Part BVIRS {2};
+ Loop_Generator Internal;
+ ME_Generator Amegic {2};
+ 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}
+ NLO_QCD_Part BVIRS {2};
+ Loop_Generator Internal;
+ ME_Generator Amegic {2};
+ 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 Amegic
+ EVENT_GENERATION_MODE = Weighted
+ NLO_Mode = 2
+}(me)
+
+(integration){
+ !-------------------------------------------------------
+ !-- Phase space setup ----------------------------------
+ !-------------------------------------------------------
+ ERROR = 1.e-3 ! Error by calculating matrix-elements
+}(integration)
+
+(shower){
+ !-------------------------------------------------------
+ !-- Parton showers -------------------------------------
+ !-------------------------------------------------------
+ SHOWER_GENERATOR = CSS
+}(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.63 ! minimum scale
+ RESCALE_EXPONENT = 0.192 ! rescaling exponent
+}(mi)
+
+(analysis){
+ BEGIN_RIVET {
+ -a MC_WJETS
+ } END_RIVET
+}(analysis)
Added: schools/2011-Kyoto/setups/day1/sherpa/Run.20-3.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day1/sherpa/Run.20-3.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,124 @@
+(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
+ !
+ BEAM_2 = -2212 ! possible beam particles: P+, P-, e+, e-
+ BEAM_ENERGY_2 = 980. ! in GeV
+}(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
+}(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.63 ! minimum scale
+ RESCALE_EXPONENT = 0.192 ! rescaling exponent
+}(mi)
+
+(analysis){
+ BEGIN_RIVET {
+ -a MC_WJETS
+ } END_RIVET
+}(analysis)
Added: schools/2011-Kyoto/setups/day1/sherpa/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day1/sherpa/Run.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,122 @@
+(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
+ !
+ BEAM_2 = -2212 ! possible beam particles: P+, P-, e+, e-
+ BEAM_ENERGY_2 = 980. ! in GeV
+}(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
+}(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.63 ! minimum scale
+ RESCALE_EXPONENT = 0.192 ! rescaling exponent
+}(mi)
+
+(analysis){
+ BEGIN_RIVET {
+ -a MC_WJETS
+ } END_RIVET
+}(analysis)
Added: schools/2011-Kyoto/setups/day2/herwig/ME/Base.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/ME/Base.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,90 @@
+##################################################
+# 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
+##################################################
+set /Herwig/Shower/Evolver:HardEmissionMode POWHEG
+create Herwig::O2AlphaS /Herwig/AlphaQCD_O2
+set /Herwig/Generators/LHCGenerator:StandardModelParameters:QCD/RunningAlphaS /Herwig/AlphaQCD_O2
+
+##################################################
+# 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_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/2011-Kyoto/setups/day2/herwig/ME/TVT-0jet.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/ME/TVT-0jet.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/ME/TVT-1jet.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/ME/TVT-1jet.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/ME/TVT-Powheg.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/ME/TVT-Powheg.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/QED/Base.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/QED/Base.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,73 @@
+##################################################
+# 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
+##################################################
+set /Herwig/Shower/Evolver:HardEmissionMode POWHEG
+create Herwig::O2AlphaS /Herwig/AlphaQCD_O2
+set /Herwig/Generators/LHCGenerator:StandardModelParameters:QCD/RunningAlphaS /Herwig/AlphaQCD_O2
+
+##################################################
+# 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_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/2011-Kyoto/setups/day2/herwig/QED/TVT-0jet-QED-off.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/QED/TVT-0jet-QED-off.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/QED/TVT-0jet-QED-on.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/QED/TVT-0jet-QED-on.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/QED/extract_noclus.sh
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/QED/extract_noclus.sh Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,11 @@
+#! /bin/bash
+
+for i in $@ ; do
+ aida2flat -m _NOCLUS $i | sed "s/_NOCLUS//g" > ${i/.aida/}_NOCLUS.dat
+ flat2aida ${i/.aida/}_NOCLUS.dat
+ rm ${i/.aida/}_NOCLUS.dat
+
+ aida2flat -M _NOCLUS $i > ${i/.aida/}_CLUS.dat
+ flat2aida ${i/.aida/}_CLUS.dat
+ rm ${i/.aida/}_CLUS.dat
+done
Added: schools/2011-Kyoto/setups/day2/herwig/shower/Base.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/shower/Base.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,79 @@
+##################################################
+# 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
+##################################################
+set /Herwig/Shower/Evolver:HardEmissionMode POWHEG
+create Herwig::O2AlphaS /Herwig/AlphaQCD_O2
+set /Herwig/Generators/LHCGenerator:StandardModelParameters:QCD/RunningAlphaS /Herwig/AlphaQCD_O2
+
+##################################################
+# 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_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/2011-Kyoto/setups/day2/herwig/shower/TVT-0jet-MEC.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/shower/TVT-0jet-MEC.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/shower/TVT-0jet-soft.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/shower/TVT-0jet-soft.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/shower/TVT-1jet.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/shower/TVT-1jet.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/herwig/shower/TVT-Powheg.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/herwig/shower/TVT-Powheg.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/pythia/ME/0jet.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/ME/0jet.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,48 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+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/2011-Kyoto/setups/day2/pythia/ME/1jet.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/ME/1jet.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,50 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+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/2011-Kyoto/setups/day2/pythia/ME/pythia_rivet
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/ME/pythia_rivet Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/pythia/QED/QED.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/QED/QED.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,37 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+#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/2011-Kyoto/setups/day2/pythia/QED/noQED.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/QED/noQED.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,37 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+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/2011-Kyoto/setups/day2/pythia/QED/pythia_rivet
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/QED/pythia_rivet Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/pythia/Zpythia.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/Zpythia.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,53 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+#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/2011-Kyoto/setups/day2/pythia/shower/0jet.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/shower/0jet.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,43 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+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/2011-Kyoto/setups/day2/pythia/shower/1jet.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/shower/1jet.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,39 @@
+! 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:onMode = off ! switch off all Z0 decay modes...
+23:onIfMatch = 11 -11 ! ...and back on Z0 -> e- e+
+23:onIfMatch = 13 -13 ! ...and back on Z0 -> mu- mu+
+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/2011-Kyoto/setups/day2/pythia/shower/pythia_rivet
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/pythia/shower/pythia_rivet Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/shared/CDF_2008_S7540469.pdf
==============================================================================
Binary file. No diff available.
Added: schools/2011-Kyoto/setups/day2/shared/D0_2007_S7075677.pdf
==============================================================================
Binary file. No diff available.
Added: schools/2011-Kyoto/setups/day2/shared/D0_2008_S7554427.pdf
==============================================================================
Binary file. No diff available.
Added: schools/2011-Kyoto/setups/day2/shared/D0_2008_S7863608.pdf
==============================================================================
Binary file. No diff available.
Added: schools/2011-Kyoto/setups/day2/sherpa/ME/Run.0jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/ME/Run.0jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/ME/Run.1jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/ME/Run.1jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/ME/Run.2jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/ME/Run.2jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/ME/Run.3jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/ME/Run.3jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/QED/Run.ME.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/QED/Run.ME.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/QED/Run.Off.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/QED/Run.Off.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/merging/Run.0jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/merging/Run.0jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/merging/Run.1jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/merging/Run.1jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/merging/Run.2jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/merging/Run.2jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day2/sherpa/merging/Run.3jet.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day2/sherpa/merging/Run.3jet.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -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/2011-Kyoto/setups/day3/herwig/TopBackground.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/herwig/TopBackground.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,85 @@
+##################################################
+# Example generator based on LHC parameters
+# usage: Herwig++ read LHC.in
+##################################################
+cd /Herwig/EventHandlers
+# create the object ot read to madgraph sample
+create ThePEG::LesHouchesFileReader TopReader LesHouches.so MadGraphReader.so
+set TopReader:WeightWarnings 0
+set TopReader:PDFA /Herwig/Partons/MRST
+set TopReader:PDFB /Herwig/Partons/MRST
+set TopReader:EBeamA 3500.*GeV
+set TopReader:EBeamB 3500.*GeV
+set TopReader:BeamA /Herwig/Particles/p+
+set TopReader:BeamB /Herwig/Particles/p+
+set TopReader:Cuts /Herwig/Cuts/QCDCuts
+set TopReader:PartonExtractor /Herwig/Partons/QCDExtractor
+
+# prepared samples
+set TopReader:FileName /home/rivet/school/day3/shared/Wposbbjj.lhe
+
+do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+;
+
+##################################################
+# Connect the LesHouches reader to the generator
+##################################################
+# create the event handler
+create ThePEG::LesHouchesEventHandler TopEventHandler LesHouches.so
+insert TopEventHandler:LesHouchesReaders[0] TopReader
+set TopEventHandler:CascadeHandler /Herwig/Shower/ShowerHandler
+set TopEventHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler
+set TopEventHandler:DecayHandler /Herwig/Decays/DecayHandler
+cd /Herwig/Generators
+set LHCGenerator:EventHandler /Herwig/EventHandlers/TopEventHandler
+
+##################################################
+# Technical parameters for this run
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 0
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+# Intrinsic pT tune extrapolated to LHC energy
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
+
+##################################################
+# Enable HepMC format output
+##################################################
+cd /Herwig/Generators
+
+# Either use block A -or- block B here
+
+
+###
+# A - save to HepMC file. Uses lots of space, good for debugging analysis
+###
+
+insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+set /Herwig/Analysis/HepMCFile:PrintEvent 5000
+set /Herwig/Analysis/HepMCFile:Format GenEvent
+set /Herwig/Analysis/HepMCFile:Units GeV_mm
+set /Herwig/Analysis/HepMCFile:Filename top-background.hepmc
+
+###
+# B - good for high statistics, use to get better plots
+###
+
+#create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+#insert Rivet:Analyses 0 MC_TOP
+#insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+
+
+
+
+get LHCGenerator:AnalysisHandlers
+erase LHCGenerator:AnalysisHandlers 1
+get LHCGenerator:AnalysisHandlers
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+saverun TopBackground LHCGenerator
Added: schools/2011-Kyoto/setups/day3/herwig/TopSignal.in
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/herwig/TopSignal.in Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,84 @@
+##################################################
+# Example generator based on LHC parameters
+# usage: Herwig++ read LHC.in
+##################################################
+
+##################################################
+# Technical parameters for this run
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:NumberOfEvents 10000000
+set LHCGenerator:RandomNumberGenerator:Seed 31122001
+set LHCGenerator:DebugLevel 0
+set LHCGenerator:PrintEvent 10
+set LHCGenerator:MaxErrors 10000
+
+##################################################
+# Changes to top decays to get semi-leptonic top
+# and hadronic tbar
+##################################################
+# change so t and tbar properties not automatically
+# charge conjugates of each other
+set /Herwig/Particles/t:Synchronized Not_synchronized
+set /Herwig/Particles/tbar:Synchronized Not_synchronized
+
+# top decay modes, only semi-leptonic on
+do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;
+
+# tbar decay modes, only hadronic on
+do /Herwig/Particles/tbar:SelectDecayModes tbar->b,bbar,cbar; tbar->bbar,cbar,d; tbar->bbar,cbar,s; tbar->bbar,s,ubar; tbar->bbar,ubar,d;
+
+##################################################
+# LHC physics parameters (override defaults here)
+##################################################
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 7000.0
+
+# Intrinsic pT tune extrapolated to LHC energy
+set /Herwig/Shower/Evolver:IntrinsicPtGaussian 2.2*GeV
+
+##################################################
+# Matrix Elements for hadron-hadron collisions
+# (by default only gamma/Z switched on)
+##################################################
+cd /Herwig/MatrixElements/
+
+# top-antitop production
+insert SimpleQCD:MatrixElements[0] MEHeavyQuark
+
+##################################################
+# Enable HepMC format output
+##################################################
+cd /Herwig/Generators
+
+###
+# Either use block A -or- block B here
+###
+
+###
+# A - save HepMC to file. Uses lots of space, good for debugging analysis
+###
+
+insert LHCGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile
+set /Herwig/Analysis/HepMCFile:PrintEvent 5000
+set /Herwig/Analysis/HepMCFile:Format GenEvent
+set /Herwig/Analysis/HepMCFile:Units GeV_mm
+set /Herwig/Analysis/HepMCFile:Filename top-signal.hepmc
+
+
+###
+# B - good for high statistics, use to get better plots
+###
+
+#create ThePEG::RivetAnalysis Rivet RivetAnalysis.so
+#insert Rivet:Analyses 0 MC_TOP
+#insert LHCGenerator:AnalysisHandlers 0 Rivet
+
+
+
+
+
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+saverun TopSignal LHCGenerator
Added: schools/2011-Kyoto/setups/day3/pythia/pythia_rivet
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/pythia/pythia_rivet Tue Jul 17 14:18:39 2012 (r3830)
@@ -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_TOP -H$NAME.aida $PIPENAME
+rm -f $PIPENAME
Added: schools/2011-Kyoto/setups/day3/pythia/topbackground.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/pythia/topbackground.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,23 @@
+! File: topbackground.cmnd
+! The intent is to study the top mass in top pair production,
+! with one W decaying leptonically and the other hadronically.
+! This file generates the background of W + 4 jets from a MadGraph LHEF,
+! alternatively fewer jets, with rest to come from showers etc.
+! At parton-level the lepton and all jets have pT > 20 and Delta-R > 0.7.
+
+! 1) Settings that will be used in a main program.
+Main:numberOfEvents = 2000 ! number of events to generate
+Main:timesToShow = 0 ! show how far along run is this many times
+Main:timesAllowErrors = 10 ! abort run after this many flawed events
+Main:showChangedSettings = on ! print changed flags/modes/parameters
+Main:showAllSettings = off ! print all flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+Main:showAllParticleData = off ! print all particle and decay data
+
+! 2) Beam and process selection from a Les Houches Event File.
+Beams:frameType = 4 ! read info from a LHEF
+
+! W b bbar jet jet with W+
+! Here jet = d, dbar, u, ubar, s, sbar, c, cbar, g.
+Beams:LHEF = /home/rivet/school/day3/shared/Wposbbjj.lhe
+
Added: schools/2011-Kyoto/setups/day3/pythia/topsignal.cmnd
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/pythia/topsignal.cmnd Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,35 @@
+! File: topsignal.cmnd
+! The intent is to study the top mass in top pair production,
+! with one W decaying leptonically and the other hadronically.
+
+! 1) Settings that will be used in a main program.
+Main:numberOfEvents = 2000 ! number of events to generate
+Main:timesToShow = 0 ! show how far along run is this many times
+Main:timesAllowErrors = 10 ! abort run after this many flawed events
+Main:showChangedSettings = on ! print changed flags/modes/parameters
+Main:showAllSettings = off ! print all flags/modes/parameters
+Main:showChangedParticleData = on ! print changed particle and decay data
+Main:showAllParticleData = off ! 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 = 7000. ! CM energy of collision
+
+! 3) Pick processes.
+Top:gg2ttbar = on ! g g -> t tbar
+Top:qqbar2ttbar = on ! q qbar -> t tbar
+
+! 4) Pick sign of lepton by forcing W decay modes.
+24:onMode = off ! all W decays are switched off
+
+! 4.1) Alternative 1: W+ -> e+, W- -> jets.
+24:onPosIfAny = 11 ! W+ can decay to e+
+24:onNegIfAny = 1 2 3 4 5 ! W- can decay to quark (d u s c b)
+
+! 4.2) Alternative 2: W+ -> jets, W- -> e-.
+#24:onPosIfAny = 1 2 3 4 5 ! W+ can decay to quark (d u s c b)
+#24:onNegIfAny = 11 ! W- can decay to e-
+
+! 5) Other settings.
+#6:m0 = 174.3 ! change the top mass
Added: schools/2011-Kyoto/setups/day3/rivet-plugin/MC_TOP.cc
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/rivet-plugin/MC_TOP.cc Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,233 @@
+#include "Rivet/Analysis.hh"
+#include "Rivet/Projections/FinalState.hh"
+#include "Rivet/Projections/ChargedLeptons.hh"
+#include "Rivet/Projections/MissingMomentum.hh"
+#include "Rivet/Projections/FastJets.hh"
+#include "Rivet/AnalysisLoader.hh"
+#include "Rivet/RivetAIDA.hh"
+
+namespace Rivet {
+
+
+ class MC_TOP : public Analysis {
+ public:
+
+ /// Minimal constructor
+ MC_TOP() : Analysis("MC_TOP")
+ {nevents=0; }
+
+
+ /// @name Analysis methods
+ //@{
+
+ /// Set up projections and book histograms
+ void init() {
+
+ // A FinalState is used to select particles within |eta| < 4.2 and with pT
+ // > 30 GeV, out of which the ChargedLeptons projection picks only the
+ // electrons and muons, to be accessed later as "LFS".
+ addProjection(ChargedLeptons(FinalState(-4.2, 4.2, 30*GeV)), "LFS");
+ // A second FinalState is used to select all particles in |eta| < 4.2,
+ // with no pT cut. This is used to construct jets and measure missing
+ // transverse energy.
+ FinalState fs(-4.2, 4.2, 0*GeV);
+ addProjection(FastJets(fs, FastJets::ANTIKT, 0.6), "Jets");
+ // EXERCISE
+ // If you decide to add a missing ET cut, register a MissingMomentum projection here
+
+ // Booking of histograms
+ _h_jet_1_pT = bookHistogram1D("jet_1_pT", 50, 0, 500);
+ _h_jet_2_pT = bookHistogram1D("jet_2_pT", 50, 0, 400);
+ _h_jet_3_pT = bookHistogram1D("jet_3_pT", 50, 0, 300);
+ _h_jet_4_pT = bookHistogram1D("jet_4_pT", 50, 0, 200);
+ //
+ _h_bjet_1_pT = bookHistogram1D("jetb_1_pT", 50, 0, 250);
+ _h_bjet_2_pT = bookHistogram1D("jetb_2_pT", 50, 0, 250);
+ //
+ _h_ljet_1_pT = bookHistogram1D("jetl_1_pT", 50, 0, 250);
+ _h_ljet_2_pT = bookHistogram1D("jetl_2_pT", 50, 0, 250);
+ //
+ _h_W_mass = bookHistogram1D("W_mass", 75, 30, 180);
+ _h_t_mass = bookHistogram1D("t_mass", 30, 130, 280);
+ _h_t_mass_W_cut = bookHistogram1D("t_mass_W_cut", 40, 130, 330);
+ _h_W_comb_mass = bookHistogram1D("W_comb_mass", 50, 30, 180);
+ _h_t_comb_mass = bookHistogram1D("t_comb_mass", 40, 130, 330);
+ }
+
+
+ void analyze(const Event& event) {
+ const double weight = event.weight();
+
+ // Use the "LFS" projection to require at least one hard charged
+ // lepton. This is an experimental signature for the leptonically decaying
+ // W. This helps to reduce pure QCD backgrounds.
+ const ChargedLeptons& lfs = applyProjection<ChargedLeptons>(event, "LFS");
+ MSG_DEBUG("Charged lepton multiplicity = " << lfs.chargedLeptons().size());
+ foreach (Particle lepton, lfs.chargedLeptons()) {
+ MSG_DEBUG("Lepton pT = " << lepton.momentum().pT());
+ }
+ if (lfs.chargedLeptons().empty()) {
+ MSG_DEBUG("Event failed lepton multiplicity cut");
+ vetoEvent;
+ }
+
+ // EXERCISE
+ // Here you could add use of a missing ET cut to favour the hard neutrino
+ // from a leptonically decaying W. This helps to reduce pure QCD
+ // backgrounds.
+
+ // Use the "Jets" projection to check that there are at least 4 jets of
+ // any pT. Getting the jets sorted by pT ensures that the first jet is the
+ // hardest, and so on. We apply no pT cut here only because we want to
+ // plot all jet pTs to help optimise our jet pT cut.
+ const FastJets& jetpro = applyProjection<FastJets>(event, "Jets");
+ const Jets jets = jetpro.jetsByPt(0*GeV);
+ if (jets.size() < 4) {
+ MSG_DEBUG("Event failed jet multiplicity cut");
+ vetoEvent;
+ }
+ nevents++;
+
+ // Fill all-jets histograms
+ _h_jet_1_pT->fill(jets[0].momentum().pT(), weight);
+ _h_jet_2_pT->fill(jets[1].momentum().pT(), weight);
+ _h_jet_3_pT->fill(jets[2].momentum().pT(), weight);
+ _h_jet_4_pT->fill(jets[3].momentum().pT(), weight);
+
+ // Insist that the 4 jets pass pT hardness cuts. If we don't find at least
+ // 4 such jets, we abandon this event.
+ // EXERCISE
+ // Vary these cuts to improve the analysis efficiency / purity.
+ if (jets[0].momentum().pT() < 45*GeV || jets[1].momentum().pT() < 35*GeV || jets[3].momentum().pT() < 30*GeV) {
+ MSG_DEBUG("Event failed jet pT cut");
+ vetoEvent;
+ }
+
+ // EXERCISE
+ // Add cuts on H_T and/or Centrality. You could either apply these to all
+ // jets or just your four identified jets. Exactly which choice you make affects
+ // where you should place this cut code.
+
+ // Sort the jets into b-jets and light jets. We expect one hard b-jet from
+ // each top decay, so our 4 hardest jets should include two b-jets. The
+ // Jet::containsBottom() method is equivalent to perfect experimental
+ // b-tagging, in a generator-independent way.
+ Jets bjets, ljets;
+ foreach (const Jet& jet, jets) {
+ if (jet.containsBottom()) {
+ bjets.push_back(jet);
+ } else {
+ ljets.push_back(jet);
+ }
+ }
+ MSG_DEBUG("Number of b-jets = " << bjets.size());
+ if (bjets.size() != 2) {
+ MSG_DEBUG("Event failed b-tagging cut");
+ vetoEvent;
+ }
+
+ // EXERCISE
+ // Modify the code above to insert a lepton--jet isolation cut. The point
+ // of this is to avoid jets being constructed with dominant contributions
+ // from the hard leptons: the jets should be hadronic. This isolation cut
+ // should, for example, take the form of filling the ljets and bjets
+ // vectors by iterating through all the hard jets in decreasing pT order
+ // and vetoing any jets which are within a delta(R) of the lepton, until
+ // the vectors each contain two appropriate jets. The choice of delta(R)
+ // will affect the purity and efficiency of the cut.
+
+ // Plot the pTs of the identified jets.
+ _h_bjet_1_pT->fill(bjets[0].momentum().pT(), weight);
+ _h_bjet_2_pT->fill(bjets[1].momentum().pT(), weight);
+ _h_ljet_1_pT->fill(ljets[0].momentum().pT(), weight);
+ _h_ljet_2_pT->fill(ljets[1].momentum().pT(), weight);
+
+ // Construct the hadronically decaying W momentum 4-vector from the two
+ // hardest non-b-tagged jets. There are two b-jets with which this can be
+ // combined to make the hadronically decaying top, one of which is correct
+ // and the other is not... but we have no way to identify which is which,
+ // so we construct both possible top momenta and fill the histograms with
+ // both.
+ // EXERCISE
+ // Instead make the W candidate from all pairs of light jets and pick the combination
+ // which has an invariant mass closest to the W mass.
+ const FourMomentum W = ljets[0].momentum() + ljets[1].momentum();
+ const FourMomentum t1 = W + bjets[0].momentum();
+ const FourMomentum t2 = W + bjets[1].momentum();
+ _h_W_mass->fill(W.mass(), weight);
+ _h_t_mass->fill(t1.mass(), weight);
+ _h_t_mass->fill(t2.mass(), weight);
+
+ // Placing a cut on the W mass (which is well known) helps to reduce
+ // backgrounds.
+ if (inRange(W.mass()/GeV, 30, 130)) {
+ MSG_DEBUG("W found with mass " << W.mass()/GeV << " GeV");
+ _h_t_mass_W_cut->fill(t1.mass(), weight);
+ _h_t_mass_W_cut->fill(t2.mass(), weight);
+ }
+
+ // Fill histograms for all of the combinatoric 2-jet masses
+ _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[2].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[0].momentum() + jets[3].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[1].momentum() + jets[2].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[1].momentum() + jets[3].momentum()), weight);
+ _h_W_comb_mass->fill(mass(jets[2].momentum() + jets[3].momentum()), weight);
+
+ // Fill histograms for all of the combinatoric 3-jet masses
+ _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum() + jets[2].momentum()), weight);
+ _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[1].momentum() + jets[3].momentum()), weight);
+ _h_t_comb_mass->fill(mass(jets[0].momentum() + jets[2].momentum() + jets[3].momentum()), weight);
+ _h_t_comb_mass->fill(mass(jets[1].momentum() + jets[2].momentum() + jets[3].momentum()), weight);
+ }
+
+
+ void finalize() {
+ // After the run, we normalise all of the histograms to the number of
+ // events supplied by the generator. This is hence not physically
+ // invariant, but will be scaled by the selection efficiency to help you
+ // optimise the analysis cuts.
+ scale(_h_jet_1_pT, 1/sumOfWeights());
+ scale(_h_jet_2_pT, 1/sumOfWeights());
+ scale(_h_jet_3_pT, 1/sumOfWeights());
+ scale(_h_jet_4_pT, 1/sumOfWeights());
+ scale(_h_bjet_1_pT, 1/sumOfWeights());
+ scale(_h_bjet_2_pT, 1/sumOfWeights());
+ scale(_h_ljet_1_pT, 1/sumOfWeights());
+ scale(_h_ljet_2_pT, 1/sumOfWeights());
+ scale(_h_W_mass, 1/sumOfWeights());
+ scale(_h_t_mass, 1/sumOfWeights());
+ scale(_h_W_comb_mass, 1/sumOfWeights());
+ scale(_h_t_comb_mass, 1/sumOfWeights());
+ scale(_h_t_mass_W_cut, 1/sumOfWeights());
+ std::cout << "Nevents = " << nevents << std::endl;
+ }
+
+ //@}
+
+
+ private:
+
+ // @name Histogram data members
+ //@{
+
+ AIDA::IHistogram1D *_h_jet_1_pT, *_h_jet_2_pT, *_h_jet_3_pT, *_h_jet_4_pT;
+ AIDA::IHistogram1D *_h_bjet_1_pT, *_h_bjet_2_pT;
+ AIDA::IHistogram1D *_h_ljet_1_pT, *_h_ljet_2_pT;
+ AIDA::IHistogram1D *_h_W_mass;
+ AIDA::IHistogram1D *_h_t_mass;
+ AIDA::IHistogram1D *_h_W_comb_mass;
+ AIDA::IHistogram1D *_h_t_comb_mass;
+ AIDA::IHistogram1D *_h_t_mass_W_cut;
+
+ int nevents;
+ //@}
+
+ };
+
+
+
+ // The hook for the plugin system
+ AnalysisBuilder<MC_TOP> plugin_MC_TOP;
+
+}
Added: schools/2011-Kyoto/setups/day3/rivet-plugin/MC_TOP.plot
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/rivet-plugin/MC_TOP.plot Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,4 @@
+# BEGIN PLOT /MC_TOP
+LogY=0
+#NormalizeToIntegral=1
+# END PLOT
Added: schools/2011-Kyoto/setups/day3/rivet-plugin/Makefile
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/rivet-plugin/Makefile Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,5 @@
+Rivet_MC_TOP.so : MC_TOP.cc
+ rivet-buildplugin $@ $^
+
+clean:
+ rm -f Rivet_MC_TOP.so
Added: schools/2011-Kyoto/setups/day3/rivet-plugin/Sherpa_References.tex
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/rivet-plugin/Sherpa_References.tex Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,83 @@
+%% Citation summary file generated by Sherpa 1.3.0
+%% PID 8973 on Wed Sep 7 16:27:52 2011
+
+\documentclass{article}
+
+\begin{document}
+
+The complete Sherpa package is published under \cite{Gleisberg:2008ta}.
+The Catani-Seymour subtraction based shower is published under \cite{Schumann:2007mg}.
+Amegic is published under \cite{Krauss:2001iv}.
+Comix is published under \cite{Gleisberg:2008fv}.
+
+\end{document}
+
+
+%% You have used the following configuration:
+%% SVN branch trunk/SHERPA, revision 17407 {
+%% 38ec49fe9a854afe050e1f7c68bf25b8 AHADIC++/Decays
+%% f60dc4791463c0f2005fa59b008b251c AHADIC++/Formation
+%% e102b2dec4ee1b4045de32b42e15500b AHADIC++/Main
+%% d71182438874e54e6ca2e3560b3953ef AHADIC++/Tools
+%% 7d6e949b5df087f33a2924ad958c9547 AMEGIC++/Amplitude
+%% c97533fb683025800ca91c8953ddff12 AMEGIC++/Amplitude/Zfunctions
+%% c9afff32100a50ab40864776956d2b57 AMEGIC++/Cluster
+%% a9a5a146d4092fb0edad95a57f1634e3 AMEGIC++/DipoleSubtraction
+%% 4706fc8cf71fe6c622b0411effa0d002 AMEGIC++/Main
+%% 9e9e1cc956bde293acf4ab67565d9538 AMEGIC++/Phasespace
+%% d883ce5691ced54af560fd95f5a21905 AMEGIC++/String
+%% e0ba558a8eaae3efffb9a198e64a5e9d AMISIC++/Main
+%% e4a27b8efa13c70b84b1378b043b07e7 AMISIC++/Model
+%% 545e7742a5ae0e0f0fedffac221aa7c7 AMISIC++/Tools
+%% 2be4ed43d9dc24f11042c8b7e383e664 ATOOLS/Math
+%% dc4df66c6ed1281e97f778d3cb810608 ATOOLS/Org
+%% 40387766f0c26473b21d92dfb5056cfb ATOOLS/Phys
+%% 459219471ecfcb23c5d7d310c54f1375 BEAM/Main
+%% 5c1b12a4fe1d0504c821f8350cffe766 COMIX/Amplitude
+%% beb62928329fc288e67240983983e997 COMIX/Cluster
+%% 2ae80298e1a888e996c3dc165615db45 COMIX/Currents
+%% efac61e375872fa3d827ce87f4beca78 COMIX/Main
+%% 06114a246ba4328b8bca52afeb4acb5e COMIX/Models
+%% 3a9a6824c268719a9e5fe1b57a3b63db COMIX/Phasespace
+%% 574d40cf836bf3c8ca3d1dfedf8f2a1c COMIX/Vertices
+%% 5a8fc4193452a4091daa4487f5a015c9 CSSHOWER++/Calculators
+%% 50fd092e46c93c77525cd3346e93edce CSSHOWER++/Main
+%% 929bc15f81d38f4ae8e89b190d3de6ff CSSHOWER++/Showers
+%% c7ec1d4bd7889c6062930fdd821a0e9e CSSHOWER++/Tools
+%% bbc25793bfcdf0187cfc9a3ab48d2b6f EXTRA_XS/Cluster
+%% 45aecfd66bcd1ed20975c5de7c573bf7 EXTRA_XS/Main
+%% 8bbbfbb14c02284d5cdd5abfcbf58774 EXTRA_XS/NLO
+%% 8e0ff30d4022f8c67daece9ff5b2f9b5 EXTRA_XS/Two2Two
+%% 7e8af1cc26ccaec096f4b5c002771d8b HADRONS++/Current_Library
+%% 274db50ee1f3ec2d33479edce0478959 HADRONS++/ME_Library
+%% dea71fea5c75a98fe7afb03ea2dbcb9d HADRONS++/Main
+%% 3a1f831bf2711308dbecf024a2566de8 HADRONS++/PS_Library
+%% 26e3bb9d9d5f26b4e462ff397b9a448a METOOLS/Loops
+%% bc886455ac0dfb7bbedbed77b8a68559 METOOLS/Main
+%% 328974c78c5efcf85f7b4a829e35c796 MODEL/Interaction_Models
+%% 480b89e69f4d2f9816db458095059b01 MODEL/Main
+%% 9f775006893913bed4bddb6cf8dc46a0 PDF/CTEQ
+%% 87e1cf7fb522b5f97acfce5f1671d6df PDF/Main
+%% 3abfcc1499013c128eecdf9f0e31a4db PDF/Remnant
+%% 5296bc2423434c80c1db3966fe820141 PHASIC++/Channels
+%% c590b998d61efa088e1a25ac5050610b PHASIC++/Main
+%% 14ab43f5b840d1a7e88f8228488f5e41 PHASIC++/Process
+%% e71c0c5ead8c1582c2779b0d0be662a8 PHASIC++/Scales
+%% b0c7e4d37e20bc02f1dbb1106d35d92f PHASIC++/Selectors
+%% 8eccb4a1dacda8eed631c2b95bfe59ff PHOTONS++/MEs
+%% d71587a97f62b6caecb2965db4745703 PHOTONS++/Main
+%% eb1ff869cd6e53320dc9a312e12aa32e PHOTONS++/PhaseSpace
+%% 3a02e58c75117ed7accc06e91ea6615c PHOTONS++/Tools
+%% 693b82c25407cadef123e0b4dcc78841 POWHEG/Calculators
+%% 57ef7e74253974168836ed1101474caa POWHEG/Main
+%% f19ff7dbcd56311ff99b6adc4414dda3 POWHEG/Showers
+%% 1e2572c597fbfdbe20bf83827aabcef7 POWHEG/Tools
+%% 3841111cd692b3c4f0fb8598a399fa62 SHERPA/Initialization
+%% 49df660a70414df7adf2edd7fcafaa0b SHERPA/LundTools
+%% f264c096cfd06265fb2589e3bc897555 SHERPA/Main
+%% 4a1063f3735b737c39b2aa23d6cb3cce SHERPA/PerturbativePhysics
+%% 66ab42dd87dc012abe6ac94636711fa5 SHERPA/Single_Events
+%% fb37af043a833a37b5657242d08cde44 SHERPA/SoftPhysics
+%% ed954f841c640326701a26a5e317d086 SHERPA/Tools
+%% }
+
Added: schools/2011-Kyoto/setups/day3/shared/Wposbbjj.lhe.bz2
==============================================================================
Binary file. No diff available.
Added: schools/2011-Kyoto/setups/day3/sherpa/ttbar/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/sherpa/ttbar/Run.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,46 @@
+(run){
+ EVENTS = 10000
+ ANALYSIS Rivet
+ QCUT:=20.
+ NJET:=0
+ EVENT_GENERATION_MODE = Weighted
+ MASSIVE[15]=1
+ STABLE[15] = 1 ! Set tau lepton (PDG-Id 15) stable
+}(run)
+
+(beam){
+ BEAM_1 2212; BEAM_ENERGY_1 3500;
+ BEAM_2 2212; BEAM_ENERGY_2 3500;
+}(beam)
+
+(processes){
+ Process 93 93 -> 6[a] -6[b] 93{NJET};
+ Decay 6[a] -> 5 24[c];
+ Decay 24[c] -> 90 91;
+ Decay -6[b] -> -5 -24[d];
+ Decay -24[d] -> 94 94;
+ CKKW sqr(QCUT/E_CMS);
+ Integration_Error 0.05 {7};
+ Integration_Error 0.1 {8};
+ End process;
+
+ Process 93 93 -> 6[a] -6[b] 93{NJET};
+ Decay 6[a] -> 5 24[c];
+ Decay 24[c] -> 94 94;
+ Decay -6[b] -> -5 -24[d];
+ Decay -24[d] -> 90 91;
+ CKKW sqr(QCUT/E_CMS);
+ Integration_Error 0.05 {7};
+ Integration_Error 0.1 {8};
+ End process;
+}(processes)
+
+(mi){
+ MI_HANDLER = Amisic # None or Amisic
+}(mi)
+
+(analysis){
+ BEGIN_RIVET {
+ -a MC_TOP
+ } END_RIVET
+}(analysis)
Added: schools/2011-Kyoto/setups/day3/sherpa/wbbjets/Run.dat
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ schools/2011-Kyoto/setups/day3/sherpa/wbbjets/Run.dat Tue Jul 17 14:18:39 2012 (r3830)
@@ -0,0 +1,45 @@
+(run){
+ EVENTS 10000
+ ANALYSIS Rivet
+ PARTICLE_CONTAINER 98 bquarks 5 -5;
+ MASSIVE[15] 1; ACTIVE[6] 0;
+ STABLE[15] = 1 ! Set tau lepton (PDG-Id 15) stable
+ EVENT_GENERATION_MODE Weighted;
+ SCALES LOOSE_METS;
+}(run)
+
+(beam){
+ BEAM_1 = 2212; BEAM_ENERGY_1 = 3500;
+ BEAM_2 = 2212; BEAM_ENERGY_2 = 3500;
+}(beam)
+
+(processes){
+ Process 93 93 -> 90 91 98 98 93{2}
+ Order_EW 2
+ CKKW sqr(20/E_CMS)
+ Cut_Core 1;
+ Print_Graphs Graphs;
+ Enhance_Factor 0.1 {4}
+ Enhance_Factor 0.1 {5}
+ Integration_Error 0.02 {6}
+ Integration_Error 0.05 {7}
+ End process;
+}(processes)
+
+(selector){
+ Mass -11 12 10. E_CMS
+ Mass 11 -12 10. E_CMS
+ Mass -13 14 10. E_CMS
+ Mass 13 -14 10. E_CMS
+}(selector)
+
+(mi){
+ MI_HANDLER Amisic
+}(mi)
+
+(analysis){
+ BEGIN_RIVET {
+ -a MC_TOP
+ } END_RIVET
+}(analysis)
+
More information about the Rivet-svn
mailing list
|