monosbb.in

read snippets/PPCollider.in
read FRModel.model

set /Herwig/Partons/RemnantDecayer:AllowTop Yes

### me_pdf_commands
## Configure LO PDF set for the hard process
# set /Herwig/Partons/HardLOPDF:PDFName NNPDF23_lo_as_0130_qed
# set /Herwig/Partons/HardLOPDF:Member 0
# set /Herwig/Partons/HardLOPDF:MaxFlav 5
# set /Herwig/Partons/HardLOPDF:RemnantHandler /Herwig/Partons/HadronRemnants
# set /Herwig/Particles/p+:PDF /Herwig/Partons/HardLOPDF
# set /Herwig/Particles/pbar-:PDF /Herwig/Partons/HardLOPDF
# set /Herwig/Partons/PPExtractor:FirstPDF  /Herwig/Partons/HardLOPDF
# set /Herwig/Partons/PPExtractor:SecondPDF /Herwig/Partons/HardLOPDF


## ----------------------------
## Read in Events from LHE File
## ----------------------------
## Create the Handler and Reader
library LesHouches.so
create ThePEG::LesHouchesFileReader /Herwig/EventHandlers/LHEReader
create ThePEG::LesHouchesEventHandler /Herwig/EventHandlers/LHEHandler

## Set LHE filename
set /Herwig/EventHandlers/LHEReader:FileName unweighted_events.lhe

## Setup LHE event handler and LHE reader
set /Herwig/Generators/EventGenerator:EventHandler /Herwig/EventHandlers/LHEHandler
insert /Herwig/EventHandlers/LHEHandler:LesHouchesReaders 0 /Herwig/EventHandlers/LHEReader
set /Herwig/EventHandlers/LHEHandler:PartonExtractor /Herwig/Partons/PPExtractor
set /Herwig/EventHandlers/LHEHandler:CascadeHandler /Herwig/Shower/ShowerHandler
set /Herwig/EventHandlers/LHEHandler:HadronizationHandler /Herwig/Hadronization/ClusterHadHandler
set /Herwig/EventHandlers/LHEHandler:DecayHandler /Herwig/Decays/DecayHandler
set /Herwig/EventHandlers/LHEHandler:WeightNormalization CrossSection
set /Herwig/EventHandlers/LHEHandler:WeightOption VarNegWeight
set /Herwig/EventHandlers/LHEReader:IncludeSpin Yes
set /Herwig/EventHandlers/LHEReader:MomentumTreatment RescaleEnergy
set /Herwig/EventHandlers/LHEReader:AllowedToReOpen No

## Parton shower settings
set /Herwig/Shower/ShowerHandler:MaxPtIsMuF Yes
set /Herwig/Shower/ShowerHandler:RestrictPhasespace Yes
# treatment of wide angle radiation
set /Herwig/Shower/PartnerFinder:PartnerMethod Random
set /Herwig/Shower/PartnerFinder:ScaleChoice Partner

## Don't use any cuts on LHE files
create ThePEG::Cuts /Herwig/Cuts/NoCuts
set /Herwig/EventHandlers/LHEReader:Cuts /Herwig/Cuts/NoCuts


### lhef_mg5amc_commands
## Commands specific to showering of events produced with MG5_aMC@NLO
set /Herwig/Shower/KinematicsReconstructor:ReconstructionOption General
set /Herwig/Shower/KinematicsReconstructor:InitialInitialBoostOption LongTransBoost
set /Herwig/Shower/KinematicsReconstructor:InitialStateReconOption Rapidity
# set /Herwig/Shower/KinematicsReconstructor:FinalStateReconOption Default
set /Herwig/Shower/ShowerHandler:SpinCorrelations No


### Herwig71_EvtGen.py
# Quick "fix" to the mismatch between Herwig 7 and EvtGen of the masses below
set /Herwig/Particles/B'_c1+:NominalMass 7.3
set /Herwig/Particles/B'_c1-:NominalMass 7.3
set /Herwig/Particles/B_c1+:NominalMass 7.3
set /Herwig/Particles/B_c1-:NominalMass 7.3
set /Herwig/Particles/B_c2+:NominalMass 7.35
set /Herwig/Particles/B_c2-:NominalMass 7.35
set /Herwig/Particles/B*_c0+:NominalMass 7.25
set /Herwig/Particles/B*_c0-:NominalMass 7.25
set /Herwig/Particles/B_c+:NominalMass 6.277
set /Herwig/Particles/B_c-:NominalMass 6.277
#
# This is a complete kludge.  Herwig7 has switched the mapping of names and pdgids for the D'(s1) and the D(s1)
# That means EvtGen misinterprets which particle it should be decaying, since the PDGID is what is written
# in the record.  So, we'll set things up to change the masses so that the HEPMC record lis OK.  That means
# we want:
#  Mass: 2.45960  ID: 20433 (Pythia8 calls this the D*_1s+ , EvtGen calls it D_s1-.  20433 is what Herwig calls
#                 D'_s1+ and orginally gave it a mass of 2.53535)
#  Mass: 2.53512  ID: 10433 (Pythia8 calls this the D_1s+ , EvtGen calls it D'_s1+.  10433 is what Herwig calls
#                 D_s1+ and orginally gave it a mass of 2.4589
#
#  Since EvtGen will redecay everything anyway, we'll make these particles stable in Herwig
#
set /Herwig/Particles/D'_s1+:NominalMass 2.4595000e+00
set /Herwig/Particles/D'_s1+:Width 0.001
set /Herwig/Particles/D'_s1+:WidthCut 0.01
set /Herwig/Particles/D'_s1+:Width_generator:Initialize Yes
set /Herwig/Particles/D'_s1+:Mass_generator:Initialize Yes
set /Herwig/Particles/D'_s1-:NominalMass 2.4595000e+00
set /Herwig/Particles/D'_s1-:Width 0.001
set /Herwig/Particles/D'_s1-:WidthCut 0.01
set /Herwig/Particles/D'_s1-:Width_generator:Initialize Yes
set /Herwig/Particles/D'_s1-:Mass_generator:Initialize Yes
set /Herwig/Particles/D_s1+:NominalMass 2.5352800e+00
set /Herwig/Particles/D_s1+:Width 0.001
set /Herwig/Particles/D_s1+:WidthCut 0.01
set /Herwig/Particles/D_s1+:Width_generator:Initialize Yes
set /Herwig/Particles/D_s1+:Mass_generator:Initialize Yes
set /Herwig/Particles/D_s1-:NominalMass 2.5352800e+00
set /Herwig/Particles/D_s1+:Width 0.001
set /Herwig/Particles/D_s1+:WidthCut 0.01
set /Herwig/Particles/D_s1+:Width_generator:Initialize Yes
set /Herwig/Particles/D_s1+:Mass_generator:Initialize Yes
set /Herwig/Particles/D'_s1-:Stable Stable
set /Herwig/Particles/D'_s1+:Stable Stable
set /Herwig/Particles/D_s1-:Stable Stable
set /Herwig/Particles/D_s1+:Stable Stable


### HW7 Settings
cd /Herwig/Particles
create /ThePEG/ParticleData hs
setup hs 54  hs MHs WIDTH MAX_WIDTH 0 0 0 1 0
create /ThePEG/ParticleData Zp
setup Zp 55  Zp MZP WIDTH MAX_WIDTH 0 0 0 3 0
create /ThePEG/ParticleData DM
setup DM 1000022  DM MDM WIDTH MAX_WIDTH 0 0 0 2 1
set /Herwig/Particles/DM:Stable Stable

cd /Herwig/Generators
insert EventGenerator:AnalysisHandlers 0 /Herwig/Analysis/HepMCFile

#set EventGenerator:NumberOfEvents 1
#set EventGenerator:RandomNumberGenerator:Seed 31122001
set EventGenerator:DebugLevel 0
set EventGenerator:EventHandler:StatLevel Full
set EventGenerator:PrintEvent 100
set EventGenerator:MaxErrors 1000


set /Herwig/Analysis/HepMCFile:PrintEvent 100000
set /Herwig/Analysis/HepMCFile:Format GenEvent
set /Herwig/Analysis/HepMCFile:Units GeV_mm
set /Herwig/Analysis/HepMCFile:Filename monosbb.hepmc

saverun monosbb EventGenerator