Page MenuHomeHEPForge
Authored By
aidin
May 31 2023, 4:32 PM
Size
3 KB
Referenced Files
None
Subscribers
None

LesHouches.in

# -*- ThePEG-repository -*-
##################################################
# Usage: Herwig read LesHouches.in
##################################################
# set up the collider
# we launch an electron-electron collider and
# modify it to handle photon-photon collisions
read snippets/EECollider.in
##################################################
# Set up photon-photon collisions
##################################################
# use Herwig alpha_em
cd /Herwig
create Herwig::AlphaEM AlphaEM2
set Model:EW/RunningAlphaEM AlphaEM2
# Create GammaGammaHandler with photon beams
cd /Herwig/EventHandlers
set Luminosity:Energy 1000.
set EventHandler:BeamA /Herwig/Particles/gamma
set EventHandler:BeamB /Herwig/Particles/gamma
set EventHandler:Sampler /Herwig/ACDCSampler
# the cuts
cd /Herwig/Cuts
insert Cuts:OneCuts[0] JetKtCut
insert Cuts:OneCuts[1] PhotonKtCut
insert Cuts:OneCuts[2] LeptonKtCut
insert Cuts:OneCuts[3] TopKtCut
insert Cuts:OneCuts[4] WBosonKtCut
insert Cuts:OneCuts[5] ZBosonKtCut
set JetKtCut:MinKT 10.
set LeptonKtCut:MinKT 0.
set WBosonKtCut:MinKT 0.
##################################################
# Technical parameters for this run
##################################################
cd /Herwig/Generators
set EventGenerator:NumberOfEvents 100000000
set EventGenerator:RandomNumberGenerator:Seed 31122001
set EventGenerator:DebugLevel 1
set EventGenerator:PrintEvent 0
set EventGenerator:MaxErrors 10000
set EventGenerator:EventHandler:StatLevel Full
set EventGenerator:EventHandler:CascadeHandler:MPIHandler NULL
cd /Herwig/Generators
##################################################
# Setup the Les Houches Reader(s)
##################################################
library LesHouches.so
cd /Herwig/EventHandlers
create ThePEG::LesHouchesFileReader MyLesHouchesReader
set MyLesHouchesReader:FileName yyll/Events/run_01/unweighted_events.lhe.gz
create ThePEG::Cuts /Herwig/Cuts/NoCuts
set MyLesHouchesReader:Cuts /Herwig/Cuts/NoCuts
set MyLesHouchesReader:InitPDFs 0
set MyLesHouchesReader:PDFA /Herwig/Partons/NoPDF
set MyLesHouchesReader:PDFB /Herwig/Partons/NoPDF
create ThePEG::LesHouchesEventHandler myLesHouchesHandler
set myLesHouchesHandler:CascadeHandler /Herwig/Shower/ShowerHandler
set myLesHouchesHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler
set myLesHouchesHandler:DecayHandler /Herwig/Decays/DecayHandler
set myLesHouchesHandler:PartonExtractor /Herwig/Partons/EEExtractor
set myLesHouchesHandler:WeightOption VarNegWeight
set MyLesHouchesReader:WeightWarnings false
set /Herwig/Particles/gamma:PDF /Herwig/Partons/NoPDF
insert myLesHouchesHandler:LesHouchesReaders 0 MyLesHouchesReader
###################################################
# Turn handlers on and off
##################################################
cd /Herwig/Generators
set EventGenerator:EventHandler /Herwig/EventHandlers/myLesHouchesHandler
# set EventGenerator:EventHandler:CascadeHandler NULL
# set EventGenerator:EventHandler:HadronizationHandler NULL
# set EventGenerator:EventHandler:DecayHandler NULL
set /Herwig/Shower/ShowerHandler:Interactions ALL # QCD QED QEDQCD EWOnly ALL
###################################################
# Rivet analysis
##################################################
create ThePEG::RivetAnalysis /Herwig/Analysis/Rivet RivetAnalysis.so
insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/Rivet
insert /Herwig/Analysis/Rivet:Analyses 0 RivetAnalysis_yy
###################################################
# Save run for later usage with 'Herwig run'
##################################################
saverun LesHouches EventGenerator

File Metadata

Mime Type
text/plain
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
1283815
Default Alt Text
LesHouches.in (3 KB)

Event Timeline