Physics "Data" for Atlas Rome meeting

Continuously updated/always incomplete
Last update 12 Feb 2005

Ian Hinchliffe

Link to DC0 and DC1 information
Link to DC2 information

This page is no longer being maintained. For new information status, lists of samples consult the wiki

We expect to generate, simulate and reconstruct approx 15M events. The first production will be a "Background sample of order 10M events. The G-4 simulation will occur with "initial geometry" which is still being validated.
Some of the generated events from DC2 are being resimulated.

The joboptions used for input can be found in atlas CVS. in Generators/GeneratorOptionsRome and Simulation/SimulationOptionsRome and RecExample/RecExCommon/scripts All event generation is being done "privately", but the JobOptions must be send to Ian/Borut for vetting and will then be put in the CVS area (to ensure consistency). The DC2 generation options provide more examples.

Here are the guidlines for making generator events
Here is a link to the jobs that are released for production

Important note regarding Pythia. Note that the events with Pythia (6.226) have tuned underlying event. Do not change the structure functions (CTEQ5L is the default) or any of the following parameters. There is an interplay between these parameters and the structure functions which means that they cannot be varied in a simple way. MSTP(82)=4, PARP(82)=1.8 PARP(84)=0.5. We used mstp(128)=0 in DC2. mstp(128)=2 was used in DC1 and gives a simpler structure. This change was requested by physics users and some people objected later (Ian prefers the DC1 choice). Here is the complete set of lines that should be in all pythia joboptions files.

Pythia.PythiaCommand+= ["pypars mstp 82 4","pydat1 mstj 11 3","pydat1 mstj 22 2","pydat1 parj 54 -0.07"]
Pythia.PythiaCommand+= ["pydat1 parj 55 -0.006","pypars parp 82 1.8", "pypars parp 84 0.5"]

Important note regarding Herwig. All Herwig production should use Jimmy to ensure underlying event consistency. The following parameters should be set. JMUEO=1, PTJIM=3.5, JMRAD(73)=0.57. CTEQ5L structure functions must be used and the following set.

Jimmy.JimmyCommand += ["modpdf 46", "autpdf CTEQ"]
Jimmy.JimmyCommand += ["pltcut 0.0000000000333", "clpow 1.20"]; # Make K0 and lambda stable and fix B meson to baryon ratio
Jimmy.JimmyCommand += ["ptmin  10"]; 
Jimmy.JimmyCommand += ["jmueo   1"]; 
Jimmy.JimmyCommand += ["ptjim 3.5"];                 # 
Jimmy.JimmyCommand += ["prsof 0"];                   # Tell Herwig not to produce the Underlying Event
Jimmy.JimmyCommand += ["jmrad 73 0.57"];

Important note regarding MC@NLO. All this production should use Jimmy to ensure underlying event consistency. However a next ot leading order set of structure funtions is needed and the following set of parameters must be used (note the differant value of ptjim and the structure functions)

Jimmy.JimmyCommand  = ["iproc mcatnlo"]              # signal that hard process is done by MCatNLO
Jimmy.JimmyCommand += ["modpdf 53", "autpdf CTEQ"]   # CTEQ5M1 (see mails to Arthur)
Jimmy.JimmyCommand += ["msflag 1"]                
Jimmy.JimmyCommand += ["jmbug  0"]


#-- Tuning for CTEQ5M from Arthur
Jimmy.JimmyCommand += ["ptmin  10"]; 
Jimmy.JimmyCommand += ["jmueo   1"]; 
Jimmy.JimmyCommand += ["ptjim 2.8"];                 # Still to be checked by Ian (might be 3.5!)
Jimmy.JimmyCommand += ["prsof 0"];                   # Tell Herwig not to produce the Underlying Event
Jimmy.JimmyCommand += ["jmrad 73 0.57"];   
Jimmy.JimmyCommand += ["pltcut 0.0000000000333", "clpow 1.20"]; # Make K0 and lambda stable and fix B meson to baryon ratio

Important note regarding Sherpa. We expect to use Sherpa for some samples. There are some issues to be resolved (Underlying event, tau decays, stability of K_S,Lambda etc). Most of these issues are now resolved. (Jan 19 2005)

Important note regarding Tauola. If you are using Tauola you need this.

Tauola = Algorithm( "Tauola" )
Tauola.TauolaCommand = ["tauola polar 1","tauola radcor 1","tauola phox 0.01","tauola dmode 0",<>"tauola jak1 0","tauola jak2 0"]# radiate only the tau-s; neccessary!!<>
Photos = Algorithm( "Photos" )
Photos.PhotosCommand = ["photos pmode 2","photos xphcut 0.01","photos alpha -1.","photos interf 1","photos isec 1","photos iftop 0"]