# File generated by myLesHouches.sh.
# ==> Shower is switched on
# ==> Hadronization is switched off
# Herwig will perform QCD shower only
# Underlying events will not be generated
#######################################################
# Load libraries needed to read LesHouces Event Files #
#######################################################
library powhegHerwig.so
library LesHouches.so
 
#############################################################
# Create a LesHouches reader object and open the input file #
#############################################################
cd /Herwig/EventHandlers
#create ThePEG::LesHouchesFileReader myReader
create ThePEG::powhegLesHouchesFileReader myReader
set myReader:WeightWarnings 0
set myReader:FileName pwgevents-0098.lhe.gz
set myReader:CacheFileName cache.tmp
set myReader:MaxScan 3

##############################################################
# Create a cuts object		                             #
##############################################################
cd /Herwig/EventHandlers
create ThePEG::Cuts NoCuts
set myReader:Cuts NoCuts

#################################################
#Set up the LHAPDF and set the PDFs             #
#################################################
# Don't try and find PDF index out from the LH file ...
set myReader:InitPDFs 0
# option to ensure momentum conservation is O.K. due rounding errors (recommended)
set myReader:MomentumTreatment RescaleEnergy
cd /Herwig/Partons
create ThePEG::LHAPDF myPDFset ThePEGLHAPDF.so
# Instead set them explicitly here:
set myPDFset:PDFName PDF4LHC21_40
set myPDFset:RemnantHandler HadronRemnants
set /Herwig/Particles/p+:PDF myPDFset
set /Herwig/Partons/PPExtractor:FirstPDF myPDFset
set /Herwig/Partons/PPExtractor:SecondPDF myPDFset
set /Herwig/EventHandlers/myReader:PDFA myPDFset
set /Herwig/EventHandlers/myReader:PDFB myPDFset


#############################################################
# Setup the Les Houches event handler ...                   #
#############################################################
cd /Herwig/EventHandlers

create ThePEG::LesHouchesEventHandler myLesHouchesHandler
insert myLesHouchesHandler:LesHouchesReaders[0] myReader

set myLesHouchesHandler:PartonExtractor /Herwig/Partons/PPExtractor
#set myLesHouchesHandler:WeightOption NegUnitWeight
set myLesHouchesHandler:WeightOption VarNegWeight

#Options: shower + hadr, only shower, no shower and no hard
set myLesHouchesHandler:CascadeHandler /Herwig/Shower/ShowerHandler
set myLesHouchesHandler:DecayHandler          NULL      
set myLesHouchesHandler:HadronizationHandler  NULL     
#########################################################    
# Option to off shower / hadronization / decays / MPI.  #    
#########################################################    
#cd /Herwig/EventHandlers                                         
# The handler for multiple parton interactions          
set /Herwig/Shower/ShowerHandler:MPIHandler    NULL
##################################################
#  Shower parameters
##################################################
# normally, especially for POWHEG, you want
# the scale supplied in the event files (SCALUP)
# to be used as a pT veto scale in the parton shower
set /Herwig/Shower/ShowerHandler:MaxPtIsMuF Yes
set /Herwig/Shower/ShowerHandler:RestrictPhasespace Yes
# Shower parameters
# treatment of wide angle radiation
set /Herwig/Shower/PartnerFinder:PartnerMethod Random
set /Herwig/Shower/PartnerFinder:ScaleChoice Partner


#########################################################  
# Option to off QCD/QED shower                              #   
#########################################################  
#options: QCD QED QCDandQED 
set /Herwig/Shower/ShowerHandler:Interactions QCD 
###################################################################
# Technical parameters for this run and setup the event generator #
###################################################################
cd /Herwig/Generators

# theGenerator represents the Les Houches Handler generator
create ThePEG::EventGenerator myGenerator
set myGenerator:RandomNumberGenerator /Herwig/Random
set myGenerator:StandardModelParameters /Herwig/Model
set myGenerator:EventHandler /Herwig/EventHandlers/myLesHouchesHandler
#set myGenerator:EventHandler:Cuts /Herwig/Cuts/QCDCuts

set myGenerator:NumberOfEvents 500000
set myGenerator:RandomNumberGenerator:Seed 31122001
set myGenerator:DebugLevel 1
set myGenerator:PrintEvent 5
set myGenerator:MaxErrors 200000


######################################################### 
# Set particles properties                              # 
######################################################### 

set /Herwig/Model:EW/Scheme GMuScheme
set /Herwig/Model:EW/FermiConstant 1.16638e-5


set /Herwig/Particles/Z0:NominalMass 91.1876
set /Herwig/Particles/Z0:Width 2.4952
do /Herwig/Particles/Z0:UnsetHardProcessMass
do /Herwig/Particles/Z0:UnsetHardProcessWidth

set /Herwig/Particles/W+:NominalMass 80.379
set /Herwig/Particles/W+:Width 2.085
do /Herwig/Particles/W+:UnsetHardProcessMass
do /Herwig/Particles/W+:UnsetHardProcessWidth

set /Herwig/Particles/h0:NominalMass 125.0
do /Herwig/Particles/h0:UnsetHardProcessMass
set /Herwig/Particles/h0:HardProcessWidth 0.
set /Herwig/Particles/h0:Stable Stable

set /Herwig/MatrixElements/Matchbox/Factory:FixedQEDCouplings Yes


#########################################################   
# Selectively choose to make certain particles stable   #
######################################################### 
#set /Herwig/Particles/[ParticleName]:Stable Stable
# To read the particle list: 
#	> Herwig build
#	> cd /Herwig/Particles
#	> ls


##################################################
# Analyses 					 #
##################################################
cd /Herwig/Analysis
create ThePEG::RivetAnalysis RivetAnalysis RivetAnalysis.so
#create ThePEG::NLORivetAnalysis RivetAnalysis RivetAnalysis.so
insert /Herwig/Generators/myGenerator:AnalysisHandlers 0 RivetAnalysis
insert RivetAnalysis:Analyses 0 MY_ANALYSIS


##################################################
# HepMC						 #
##################################################
cd /Herwig/Generators
create ThePEG::powhegAnalysis powhegAnalysis powhegHerwig.so
insert myGenerator:AnalysisHandlers 1 powhegAnalysis

set powhegAnalysis:Units GeV_mm
set powhegAnalysis:RunNumber 0098


#################################################################
# Save run.							#
#################################################################
cd /Herwig/Generators
saverun HerwigRun-0098 myGenerator