How to generate data-like events

It's easy to generate a sample of events which represent data. Several generators can be used (Pythia, Herwig, Isajet) to generate the process of interest. You just need to run the cdfSim executable with an appropriate directives file. Instructions for changing the directives file (here, the Pythia generator is used as an example) are given below.

How to generate events

A sample directives file can be seen in /cdf/disk1/utilities/pythia_zmumu.tcl . This file generates Z->mumu decays. Copy this and change it however you want. The relevant parts that control particle generation are shown below:

mod talk Pythia
  call_qq set t
  output set HEPG
  MSEL set 0
  commonMenu 
    for {set i 1} {$i<100} {incr i} {
      set_msub -index=$i -value=0
    }
    # Produce Zs only
    set_msub -index=15 -value=1
    set_msub -index=19 -value=0
    set_msub -index=30 -value=0
    show_msub
    # Z->mu mu only
    set_mstp -index=43 -value=2
    show_mstp
    for {set i 174} {$i<190} {incr i} {
      set_mdme -channelIndex=$i -decayType=0
    }
    set_mdme -channelIndex=184 -decayType=1
  exit
exit

Here we explicitly set the processes that (1) produce Z bosons, and (2) only allow the Z bosons to decay to muon pairs. A full description of all these switches can be found in the Pythia manual (see Useful Links below). In this example, MSEL (which can be used to set the flavour of the initial quarks) is turned off. We turn on Z production inside the commonMenu submenu, by setting the set_msub flags. The allowed Z decays are specified by the set_mdme flags.

How to run the job

Just type executable_path tcl_path , where executable_path is the full path to your cdfSim executable, and tcl_path is the full path to your directives file. Make sure that you have set the name of your output file and the number of events you want to generate in your directives file:
talk DHOutput
  output create gen_stream pythia_zmumu_3_17_0int2.gen
... you will produce a file called pythia_zmumu_3_17_0int2.gen

begin -nev 5000
.... with 5000 events in.

Useful links