Let's first make one hundred events on the Z pole at LEP:
> cd /home/potter/MG5_aMC_v2_6_1
> emacs lep_eetoztoff.dat &
Put some code like this in there:
import model sm
generate e+ e- > z > e+ e-
add process e+ e- > z > u u~
output mg5amc/lep_eetoztoff
launch
shower=Pythia8
set lpp1 0
set lpp2 0
set mz 91.1876
set ebeam1 45.5988
set ebeam2 45.5988
set nevents 100
set iseed 1
And now let's run, go to the Events directory, unpack the StdHEP file, and move it to the Delphes directory:
> ./bin/mg5_aMC lep_eetoztoff.dat
> cd mg5amc/lep_eetoztoff/Events/run_01/
> gunzip tag_1_pythia8_events.hepmc.gz
> mv tag_1_pythia8_events.hepmc /home/potter/Delphes-3.4.1/lep_eetozff.hepmc
> cd /home/potter/Delphes-3.4.1/
Let's use the CEPC detector card since there's no LEP detector card:
> ./DelphesHepMC cards/delphes_card_CEPC.tcl lep_eetozff.root lep_eetozff.hepmc
Now let's run a sample analysis in Example1 macro on the output. The macro code can be found in the examples/ directory:
> root
root[0] gSystem->Load("libDelphes")
root[0] .X examples/Example1.C("lep_eetozff.root")
Behold!! The Z!
Now let's see how these events look in an Event Display:
> make display
> root -l examples/EventDisplay.C("delphes_card_CEPC.tcl", "lep_eetozff.root")
WOW!