Page MenuHomeHEPForge
Files F8309046

No OneTemporary
Actions

View Options

diff --git a/Tests/Rivet/LHC-2760-UE.in b/Tests/Rivet/LHC-2760-UE.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-2760-UE.in
@@ -0,0 +1,42 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+# override the beam energy
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 2760.0
+
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+# min bias
+insert SimpleQCD:MatrixElements[0] MEMinBias
+# same as ue process
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+# min ibas cuts
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
+set /Herwig/Cuts/QCDCuts:X1Min 0.01
+set /Herwig/Cuts/QCDCuts:X2Min 0.01
+
+##################################################
+# select the analyses
+##################################################
+
+# underlying event forward rapidities
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_FWD_11_003
+
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+# get the cross section right
+create Herwig::MPIXSecReweighter MPIXSecReweighter
+insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter
+saverun LHC-2760-UE LHCGenerator
diff --git a/Tests/Rivet/LHC-7-GammaGamma.in b/Tests/Rivet/LHC-7-GammaGamma.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-GammaGamma.in
@@ -0,0 +1,60 @@
+##################################################
+# LHC gamma gamma initiated processes
+##################################################
+
+##################################################
+# Technical parameters for this run
+##################################################
+read LHCBase.in
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:Sampler:Ntry 10000
+set /Herwig/Partons/QCDExtractor:FlatSHatY 0
+
+##################################################
+# Change the proton PDFs to those for photon radiation
+##################################################
+set /Herwig/Particles/p+:PDF /Herwig/Partons/BudnevPDF
+set /Herwig/Particles/pbar-:PDF /Herwig/Partons/BudnevPDF
+
+##################################################
+# Cuts
+##################################################
+cd /Herwig/Cuts
+set QCDCuts:ScaleMin 0.0*GeV
+set QCDCuts:X1Min 0
+set QCDCuts:X2Min 0
+set QCDCuts:X1Max 1.
+set QCDCuts:X2Max 1.
+set QCDCuts:MHatMin 1.*GeV
+erase QCDCuts:MultiCuts 0
+set LeptonKtCut:MinKT 3*GeV
+
+##################################################
+# Selected the hard process
+##################################################
+cd /Herwig/MatrixElements
+
+# fermion-antifermion
+insert SimpleQCD:MatrixElements 0 /Herwig/MatrixElements/MEgg2ff
+set /Herwig/MatrixElements/MEgg2ff:Process Muon
+
+##################################################
+# LHC physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 7000.
+set LHCGenerator:EventHandler:CascadeHandler:MPIHandler NULL
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_I954992
+
+
+###################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+saverun LHC-7-GammaGamma LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Heavy-0.in b/Tests/Rivet/LHC-7-Heavy-0.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-Heavy-0.in
@@ -0,0 +1,31 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+cp MEHeavyQuark MEBottom
+set MEBottom:QuarkType Bottom
+set MEBottom:Process Pair
+insert SimpleQCD:MatrixElements[0] MEBottom
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 0.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# forward B production
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2010_I867355
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-7-Heavy-0 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-All-3.in b/Tests/Rivet/LHC-7-Jets-All-3.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-Jets-All-3.in
@@ -0,0 +1,29 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEQCD2to2
+set MEQCD2to2:MaximumFlavour 5
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 90.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# flavour composition of jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1188891
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-7-Jets-All-3 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-All-4.in b/Tests/Rivet/LHC-7-Jets-All-4.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-Jets-All-4.in
@@ -0,0 +1,29 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEQCD2to2
+set MEQCD2to2:MaximumFlavour 5
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 150.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# flavour composition of jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1188891
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-7-Jets-All-4 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-All-5.in b/Tests/Rivet/LHC-7-Jets-All-5.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-Jets-All-5.in
@@ -0,0 +1,29 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEQCD2to2
+set MEQCD2to2:MaximumFlavour 5
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 300.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# flavour composition of jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1188891
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-7-Jets-All-5 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Top-L.in b/Tests/Rivet/LHC-7-Top-L.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-Top-L.in
@@ -0,0 +1,28 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+set MEHeavyQuark:QuarkType Top
+insert SimpleQCD:MatrixElements[0] MEHeavyQuark
+do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# ATLAS top with jet veto
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1094568
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-7-Top-L LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Top-SL.in b/Tests/Rivet/LHC-7-Top-SL.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-7-Top-SL.in
@@ -0,0 +1,31 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+set MEHeavyQuark:QuarkType Top
+insert SimpleQCD:MatrixElements[0] MEHeavyQuark
+set /Herwig/Particles/t:Synchronized Not_synchronized
+set /Herwig/Particles/tbar:Synchronized Not_synchronized
+do /Herwig/Particles/t:SelectDecayModes t->nu_e,e+,b; t->nu_mu,mu+,b;
+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;
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# MC analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_TTBAR
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-7-Top-SL LHCGenerator
diff --git a/Tests/Rivet/LHC-Powheg-VBF.in b/Tests/Rivet/LHC-Powheg-VBF.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-Powheg-VBF.in
@@ -0,0 +1,34 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase-Powheg.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+
+# select h0 -> tau+ tau decays
+do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;
+set /Herwig/Particles/tau-:Stable Stable
+set /Herwig/Particles/h0:NominalMass 120.
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+cd /Herwig/MatrixElements
+##################################################
+# gg -> Higgs
+##################################################
+insert SimpleQCD:MatrixElements[0] PowhegMEPP2HiggsVBF
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_HJETS
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-Powheg-VBF LHCGenerator
diff --git a/Tests/Rivet/LHC-Powheg-WH.in b/Tests/Rivet/LHC-Powheg-WH.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-Powheg-WH.in
@@ -0,0 +1,29 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase-Powheg.in
+
+# select h0 -> b bbar decays
+do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;
+# leptonic W decays
+do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+;
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] PowhegMEPP2WH
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_VH2BB
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-Powheg-WH LHCGenerator
diff --git a/Tests/Rivet/LHC-Powheg-WZ.in b/Tests/Rivet/LHC-Powheg-WZ.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-Powheg-WZ.in
@@ -0,0 +1,36 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase-Powheg.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+# Drell-Yan Z/gamma
+insert SimpleQCD:MatrixElements[0] PowhegMEPP2VV
+set PowhegMEPP2VV:Process WZ
+
+#do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+;
+#do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_ebar,e-; /Herwig/Particles/W-/W-->nu_mubar,mu-;
+#do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+;
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter
+insert /Herwig/EventHandlers/LHCHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter
+
+
+##################################################
+# select the analyses
+##################################################
+
+# general analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I954993
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-Powheg-WZ LHCGenerator
diff --git a/Tests/Rivet/LHC-Powheg-ZH.in b/Tests/Rivet/LHC-Powheg-ZH.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-Powheg-ZH.in
@@ -0,0 +1,29 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase-Powheg.in
+
+# select h0 -> b bbar decays
+do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;
+# leptonic W decays
+do /Herwig/Particles/Z0:SelectDecayModes Z0->e-,e+; Z0->mu-,mu+;
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] PowhegMEPP2ZH
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_VH2BB
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-Powheg-ZH LHCGenerator
diff --git a/Tests/Rivet/LHC-Powheg-ggH.in b/Tests/Rivet/LHC-Powheg-ggH.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-Powheg-ggH.in
@@ -0,0 +1,34 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase-Powheg.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+
+# select h0 -> tau+ tau- decays
+do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;
+set /Herwig/Particles/tau-:Stable Stable
+set /Herwig/Particles/h0:NominalMass 120.
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+cd /Herwig/MatrixElements
+##################################################
+# gg -> Higgs
+##################################################
+insert SimpleQCD:MatrixElements[0] PowhegMEHiggs
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_HJETS
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-Powheg-ggH LHCGenerator
diff --git a/Tests/Rivet/LHC-VBF.in b/Tests/Rivet/LHC-VBF.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-VBF.in
@@ -0,0 +1,34 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+
+# select h0 -> tau+ tau decays
+do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;
+set /Herwig/Particles/tau-:Stable Stable
+set /Herwig/Particles/h0:NominalMass 120.
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# gg -> Higgs
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEPP2HiggsVBF
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_HJETS
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-VBF LHCGenerator
diff --git a/Tests/Rivet/LHC-WH.in b/Tests/Rivet/LHC-WH.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-WH.in
@@ -0,0 +1,29 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+# select h0 -> b bbar decays
+do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;
+# leptonic W decays
+do /Herwig/Particles/W+:SelectDecayModes W+->nu_e,e+; W+->nu_mu,mu+;
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEPP2WH
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_VH2BB
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-WH LHCGenerator
diff --git a/Tests/Rivet/LHC-WZ.in b/Tests/Rivet/LHC-WZ.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-WZ.in
@@ -0,0 +1,37 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+# Drell-Yan Z/gamma
+insert SimpleQCD:MatrixElements[0] MEPP2VV
+set MEPP2VV:Process WZ
+
+do /Herwig/Particles/W+:SelectDecayModes /Herwig/Particles/W+/W+->nu_e,e+; /Herwig/Particles/W+/W+->nu_mu,mu+;
+do /Herwig/Particles/W-:SelectDecayModes /Herwig/Particles/W-/W-->nu_ebar,e-; /Herwig/Particles/W-/W-->nu_mubar,mu-;
+do /Herwig/Particles/Z0:SelectDecayModes /Herwig/Particles/Z0/Z0->e-,e+; /Herwig/Particles/Z0/Z0->mu-,mu+;
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+# reweight for the branching ratios
+create Herwig::BranchingRatioReweighter /Herwig/Generators/BRReweighter
+insert /Herwig/EventHandlers/LHCHandler:PostHadronizationHandlers 0 /Herwig/Generators/BRReweighter
+
+
+##################################################
+# select the analyses
+##################################################
+
+# general analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I954993
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-WZ LHCGenerator
diff --git a/Tests/Rivet/LHC-ZH.in b/Tests/Rivet/LHC-ZH.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-ZH.in
@@ -0,0 +1,29 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+# select h0 -> b bbar decays
+do /Herwig/Particles/h0:SelectDecayModes h0->b,bbar;
+# leptonic W decays
+do /Herwig/Particles/Z0:SelectDecayModes Z0->e-,e+; Z0->mu-,mu+;
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEPP2ZH
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_VH2BB
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-ZH LHCGenerator
diff --git a/Tests/Rivet/LHC-ggH.in b/Tests/Rivet/LHC-ggH.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-ggH.in
@@ -0,0 +1,37 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+
+# select h0 -> tau+ tau decays
+do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;
+set /Herwig/Particles/tau-:Stable Stable
+set /Herwig/Particles/h0:NominalMass 120.
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# gg -> Higgs
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEHiggs
+insert SimpleQCD:MatrixElements[0] MEHiggsJet
+set MEHiggsJet:Process qqbar
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_HJETS
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-ggH LHCGenerator
diff --git a/Tests/Rivet/LHC-ggHJet.in b/Tests/Rivet/LHC-ggHJet.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/LHC-ggHJet.in
@@ -0,0 +1,34 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read LHCBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+
+# select h0 -> tau+ tau decays
+do /Herwig/Particles/h0:SelectDecayModes h0->tau-,tau+;
+set /Herwig/Particles/tau-:Stable Stable
+set /Herwig/Particles/h0:NominalMass 120.
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+cd /Herwig/MatrixElements
+##################################################
+# gg -> Higgs
+##################################################
+insert SimpleQCD:MatrixElements[0] MEHiggsJet
+
+##################################################
+# select the analyses
+##################################################
+
+# General analysis
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 MC_HJETS
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun LHC-ggHJet LHCGenerator
diff --git a/Tests/Rivet/TVT-300-Jets-1.in b/Tests/Rivet/TVT-300-Jets-1.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/TVT-300-Jets-1.in
@@ -0,0 +1,34 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read TVTBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 300.0
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEQCD2to2
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 6.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# run I CDF UE
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CDF_2012_NOTE10874
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-300-Jets-1 LHCGenerator
diff --git a/Tests/Rivet/TVT-300-UE.in b/Tests/Rivet/TVT-300-UE.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/TVT-300-UE.in
@@ -0,0 +1,39 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read TVT-UEBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 300.0
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+# min bias
+insert SimpleQCD:MatrixElements[0] MEMinBias
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+# min ibas cuts
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
+set /Herwig/Cuts/QCDCuts:X1Min 0.01
+set /Herwig/Cuts/QCDCuts:X2Min 0.01
+
+##################################################
+# select the analyses
+##################################################
+
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CDF_2012_NOTE10874
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+# get the cross section right
+create Herwig::MPIXSecReweighter MPIXSecReweighter
+insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter
+saverun TVT-300-UE LHCGenerator
diff --git a/Tests/Rivet/TVT-900-Jets-1.in b/Tests/Rivet/TVT-900-Jets-1.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/TVT-900-Jets-1.in
@@ -0,0 +1,34 @@
+###################################################
+# Rivet analyses at the Z pole
+##################################################
+read TVTBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 900.0
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEQCD2to2
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 10.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# run I CDF UE
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CDF_2012_NOTE10874
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-900-Jets-1 LHCGenerator
diff --git a/Tests/Rivet/TVT-900-UE.in b/Tests/Rivet/TVT-900-UE.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/TVT-900-UE.in
@@ -0,0 +1,39 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read TVT-UEBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 900.0
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+# min bias
+insert SimpleQCD:MatrixElements[0] MEMinBias
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+# min ibas cuts
+set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
+set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
+set /Herwig/Cuts/QCDCuts:X1Min 0.01
+set /Herwig/Cuts/QCDCuts:X2Min 0.01
+
+##################################################
+# select the analyses
+##################################################
+
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CDF_2012_NOTE10874
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+# get the cross section right
+create Herwig::MPIXSecReweighter MPIXSecReweighter
+insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter
+saverun TVT-900-UE LHCGenerator
diff --git a/Tests/Rivet/TVT-Run-II-Jets-0.in b/Tests/Rivet/TVT-Run-II-Jets-0.in
new file mode 100644
--- /dev/null
+++ b/Tests/Rivet/TVT-Run-II-Jets-0.in
@@ -0,0 +1,34 @@
+##################################################
+# Rivet analyses at the Z pole
+##################################################
+read TVTBase.in
+
+##################################################
+# LEP physics parameters (override defaults)
+##################################################
+cd /Herwig/Generators
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 1960.0
+
+##################################################
+# Use the q qbar matrix element
+##################################################
+cd /Herwig/MatrixElements
+insert SimpleQCD:MatrixElements[0] MEQCD2to2
+set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
+set /Herwig/Cuts/JetKtCut:MinKT 15.
+
+# No QED radiation
+erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
+
+##################################################
+# select the analyses
+##################################################
+
+# CDF jet shapes
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CDF_2012_NOTE10874
+
+##################################################
+# Save run for later usage with 'Herwig++ run'
+##################################################
+cd /Herwig/Generators
+saverun TVT-Run-II-Jets-0 LHCGenerator
diff --git a/Tests/python/merge-TVT-Energy b/Tests/python/merge-TVT-Energy
new file mode 100755
--- /dev/null
+++ b/Tests/python/merge-TVT-Energy
@@ -0,0 +1,189 @@
+#! /usr/bin/env python
+import lighthisto
+import logging
+import sys
+
+if sys.version_info[:3] < (2,4,0):
+ print "rivet scripts require Python version >= 2.4.0... exiting"
+ sys.exit(1)
+
+import os, copy, re
+from math import sqrt
+
+## Try to load faster but non-standard cElementTree module
+try:
+ import xml.etree.cElementTree as ET
+except ImportError:
+ try:
+ import cElementTree as ET
+ except ImportError:
+ try:
+ import xml.etree.ElementTree as ET
+ except:
+ sys.stderr.write("Can't load the ElementTree XML parser: please install it!\n")
+ sys.exit(1)
+
+# #############################################
+
+def fillAbove(desthisto, sourcehistosbyptmin):
+ for i,b in enumerate(desthisto.getBins()):
+ ## Fill bins with pT-ordered histos (so that 'highest always wins')
+ for ptmin, h in sorted(sourcehistosbyptmin.iteritems()):
+ newb = h.getBin(i)
+ if newb.xlow >= ptmin:
+ b.val = newb.val
+ b.errplus = newb.errplus
+ b.errminus = newb.errminus
+ b._focus= newb._focus
+
+def mergeByPt(hpath, sqrts):
+ global inhistos
+ global outhistos
+ try:
+ fillAbove(outhistos[hpath], inhistos[hpath][float(sqrts)])
+ except:
+ pass
+
+def useOnePt(hpath, sqrts, ptmin):
+ global inhistos
+ global outhistos
+ try:
+ ## Find best pT_min match
+ ptmins = inhistos[hpath][float(sqrts)].keys()
+ closest_ptmin = None
+ for ptm in ptmins:
+ if closest_ptmin is None or \
+ abs(ptm-float(ptmin)) < abs(closest_ptmin-float(ptmin)):
+ closest_ptmin = ptm
+ if closest_ptmin != float(ptmin):
+ logging.warning("Inexact match for requested pTmin=%s: " % ptmin + \
+ "using pTmin=%e instead" % closest_ptmin)
+ outhistos[hpath] = inhistos[hpath][float(sqrts)][closest_ptmin]
+ except:
+ pass
+
+# #######################################
+
+if __name__ == "__main__":
+ import logging
+ from optparse import OptionParser, OptionGroup
+ parser = OptionParser(usage="%prog name")
+ verbgroup = OptionGroup(parser, "Verbosity control")
+ verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
+ default=logging.INFO, help="print debug (very verbose) messages")
+ verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
+ default=logging.INFO, help="be very quiet")
+ parser.add_option_group(verbgroup)
+ (opts, args) = parser.parse_args()
+ logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")
+
+
+ ## Prefix used in dat file headers
+ headerprefix = "# "
+
+
+ ## Check args
+ if len(args) < 1:
+ logging.error("Must specify at least the name of the files")
+ sys.exit(1)
+
+aidafiles=["-Run-II-UE.aida","-Run-II-Jets-0.aida",
+ "-300-UE.aida" ,"-300-Jets-1.aida" ,
+ "-900-UE.aida" ,"-900-Jets-1.aida" ]
+
+## Get histos
+inhistos = {}
+outhistos={}
+weights = {}
+for f in aidafiles:
+ file=args[0]+f
+ if(file.find("Run-II-UE")>0) :
+ sqrts=1960
+ ptmin=0.
+ elif(file.find("Run-II-Jets-0")>0) :
+ sqrts=1960
+ ptmin=15
+ elif(file.find("300-UE")>0) :
+ sqrts=300
+ ptmin=0.
+ elif(file.find("300-Jets-1")>0) :
+ sqrts=300
+ ptmin=6.
+ elif(file.find("900-UE")>0) :
+ sqrts=900
+ ptmin=0.
+ elif(file.find("900-Jets-1")>0) :
+ sqrts=900
+ ptmin=10.
+ if not os.access(file, os.R_OK):
+ logging.error("%s can not be read" % file)
+ break
+ try:
+ tree = ET.parse(file)
+ except:
+ logging.error("%s can not be parsed as XML" % file)
+ break
+ tree = ET.parse(file)
+ ## Get histos from this AIDA file
+ for dps in tree.findall("dataPointSet"):
+ h = lighthisto.Histo.fromDPS(dps)
+ # di-jet decorrelations
+ # jet shapes
+ if(h.fullPath().find("NOTE10874")<0) : continue
+ if not inhistos.has_key(h.fullPath()):
+ inhistos[h.fullPath()] = {}
+ tmpE = inhistos[h.fullPath()]
+ if not tmpE.has_key(sqrts):
+ tmpE[sqrts] = {}
+ tmpP = tmpE[sqrts]
+ if not tmpP.has_key(ptmin):
+ tmpP[ptmin] = h
+ else:
+ raise Exception("A set with ptmin = %s already exists" % ( ptmin))
+
+## Make empty output histos if needed
+for hpath,hsets in inhistos.iteritems():
+ workhisto = copy.deepcopy(hsets.values()[0].values()[0])
+ outhistos[hpath] = workhisto
+ ## Empty the bin set for histos which we're going to merge
+ for b in outhistos[hpath]:
+ b.val = 0
+ b.errplus = 0
+ b.errminus = 0
+ b._focus= None
+
+## Field analysis
+logging.info("Processing CDF_2012_NOTE10874")
+## Angular distributions in different pT bins
+mergeByPt("/CDF_2012_NOTE10874/d01-x01-y01", "300")
+mergeByPt("/CDF_2012_NOTE10874/d01-x01-y02", "900")
+mergeByPt("/CDF_2012_NOTE10874/d01-x01-y03","1960")
+mergeByPt("/CDF_2012_NOTE10874/d02-x01-y01", "300")
+mergeByPt("/CDF_2012_NOTE10874/d02-x01-y02", "900")
+mergeByPt("/CDF_2012_NOTE10874/d02-x01-y03","1960")
+mergeByPt("/CDF_2012_NOTE10874/d03-x01-y01", "300")
+mergeByPt("/CDF_2012_NOTE10874/d03-x01-y02", "900")
+mergeByPt("/CDF_2012_NOTE10874/d03-x01-y03","1960")
+
+# Choose output file
+for i in range(1,4) :
+ print i
+ if(i==1) :
+ name = args[0]+"-300-Energy.aida"
+ elif(i==2) :
+ name = args[0]+"-900-Energy.aida"
+ elif(i==3) :
+ name = args[0]+"-1960-Energy.aida"
+ out = open(name, "w")
+ ## Write out merged histos
+ out.write('<?xml version="1.0" encoding="ISO-8859-1" ?>\n')
+ out.write('<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">\n')
+ out.write('<aida version="3.3">\n')
+ out.write(' <implementation version="1.1" package="FreeHEP"/>\n')
+ for hpath, h in sorted(outhistos.iteritems()):
+ if(hpath.find("d0"+str(i))>0) :
+ logging.debug("hpath = %s" % hpath)
+ out.write(h.asAIDA() + "\n\n")
+ out.write('</aida>\n')
+
+sys.exit(0)

File Metadata

Mime Type
text/x-diff
Expires
Sat, Dec 21, 2:05 PM (12 h, 17 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4023060
Default Alt Text
(38 KB)

Event Timeline