EvtGen

Home Computing DAQ E-mail Notes Meetings Subsystems Search

 

This page is obsolete and no longer maintained.

The new address for the Gauss/EvtGen web site is:

http://cern.ch/LHCb-release-area/DOC/gauss/generator/evtgen.php

 

EvtGen is the package used in Gauss to generate B decays. It was originally developed and it is maintained by David Lange and Anders Ryd from BaBar and CLEO.

Documentation

A detailed documentation about the conventions used in EvtGen or about the decay models available and the way to use them can be found on the EvtGen web page or at this location in pdf. Instructions on how to setup EvtGen in Gauss and how to write decay files are also given in the Gauss manual.

A documentation of the decay model PVVCPLH is also available here.

Brief description of the generation sequences in Gauss

There are 1 algorithm (Generation) and several tools used in the generation of signal events in Gauss:

And there are 3 methods used to generate signal events (corresponding to the last 3 tools):

The type of method used is decided by the 1 job options: Generation.SampleGenerationTool which has to be set to the name of the tool to use to generate events.

How to run stand alone EvtGen

To perform checks or studies with the generator, it is possible to use Gauss to run only the generation part of the simulation program. This first part may consist :

To d, it is possible to use the Gauss executable with the option file GenStandAlone.opts available in the Gauss package. Using for example Gauss version v25r4 (which is the most recent DC06 version), you have to type on a lxplus machine:

  1. GaussEnv v25r4
  2. cd ~/cmtuser
  3. getpack Sim/Gauss v25r4
  4. getpack Gen/DecFiles head
  5. cd ~/cmtuser/Gen/DecFiles/v13r0/cmt
  6. source setup.csh
  7. cd ~/cmtuser/Sim/Gauss/v25r4/cmt
  8. edit requirement file to "use DecFiles v13r0 Gen"
  9. source setup.csh
  10. gmake
  11. cd ~/cmtuser/Sim/Gauss/v25r4/options
  12. edit "Gauss.opts" file and:
    1. Comment or remove ApplicationMgr.TopAlg += { "Sequencer/Simulation" } ;
    2. Comment or remove #include "$GAUSSOPTS/Simulation.opts"
    3. Comment or remove ApplicationMgr.TopAlg += { "Sequencer/Monitor" } ;
    4. Add #include "$DECFILESROOT/options/XXXXXXXX.opts" at the end of the file.
    5. Uncomment #include "$GAUSSOPTS/GenStandAlone.opts"
    6. In file GenStandAlone.opts comment or remove #include "$GAUSSOPTS/GenStandAlone.opts"
    7. In file GenStandAlone.opts add the lines:

    GaussTape.ItemList = {
    "/Event#1",
    "/Event/Gen#1",
    "/Event/Gen/Header#1",
    "/Event/Gen/Collisions#1",
    "/Event/Gen/HepMCEvents#1"
    };

  13. cd ~/cmtuser/Sim/Gauss/v25r4/job
  14. ../slc3_ia32_gcc323/Gauss.exe ../options/Gauss.opts

You will obtain a file named Gauss.sim in POOL/ROOT format containing the decay trees generated by EvtGen and optionally Pythia. The behavior of the job is controlled by the option files in the options directory of Gauss. The main properties that you may change are:

The output file contains the result of the generation in HepMC format (HepMCEvents containers) in three different locations:

By default, only /Event/Gen/HepMCEvents is saved in the output file. You can choose what to save and add things to save in the option file GenStandAlone.opts. The location "/Event/Gen/SignalDecayTree" is only available in case of events generated with EvtGen only (see at the end of the page) or when generating clean events. The generation of clean events can be activated by job options in the file GenStandAlone.opts.

Please note that as indicated above, not all locations are available depending on the options activated. When trying to save a location that is not available, the program may crash. In this case one should check the consistency between the different options.

By default, EvtGen will be used in conjunction with Pythia and B hadrons will be produced with the correct momentum distribution. If you need only the decay kinematics in the B rest frame, you can use EvtGen alone (it is much faster but generates only B at rest). For this you need to modify the following options:

 By default, Gauss generates a mixture of b and b hadrons (with 50-50 proportions if the signal decay mode is not a CP mode and taking into account the CP asymmetries if it is a CP mode decayed according to a CP decay model). If you want to generate only one type of B hadron, you can specify the PDG Id of the hadron you want to generate in Gauss.opts, at the very end of the file. For B0 for example:

The random seed is controlled by the run number and the event number. To generate independent events,  you have to change for each job the run number with the option GaussInit.RunNumber in Gauss.opts.

Some conventions

In EvtGen, the units are the following:

In HepMC format, the mixing is represented this way: if a B meson is produced as B0 for example and decays as B0, the HepMC record will contain the decay chain B0B0. The first B meson is the type at production (tag) and the second B is the type at decay. If there is only one B, both types are the same and the B has not oscillated. If the decay mode of the B is a CP mode, the flavor of the second B is meaningless and the tag is given by the type of the first B. Since the second B has 0 lifetime, both B's will decay at the same position but the production position of the B is the production vertex of the first B.

Decay files available

The following decay modes are predefined in the package DecFiles. To generate a sample according to the decay mode, change the number of the included option file at the end of GenStandAlone.opts. To have more details about what is generated, look at the corresponding options files and decay files in DecFiles/options and DecFiles/dkfiles.

At the beginning of the job, the branching fraction in the generic decay table of the decay mode under study is printed. If this branching fraction is 0 or if there is a message saying that the decay mode is not in DECAY.DEC, please tell me.

When testing decay files, one should use the latest version of the DecFiles package (typing "getpack Gen/DecFiles head" in ~/cmtuser and updating the version of the package in cmt/requirements of Gauss.)

Event Type, Nickname and Decay Mode follow some conventions described in the following LHCb internal note: LHCb-2005-034.

description of the new python script used to create job options from dkfile is in pdf

Event Type Nickname Decay Mode Total BR Tested By Date
13144001 Bs_Jpsiphi,mm=CPV Bs->Jpsi(mm) phi (with CP Violation) 3.92 * 10^-5 YES Tristan Du Pree 16 Oct 2006
13144002 Bs_Jpsiphi,mm=CPV,DecProdCut Bs->Jpsi(mm) phi (with CP Violation and decay products in acceptance) 3.92 * 10^-5 YES Tristan Du Pree 16 Oct 2006
13154001 Bs_Jpsiphi,ee=CPV Bs->Jpsi(ee) phi (with CP Violation) 3.92 * 10^-5 YES Tristan Du Pree 16 Oct 2006
13154002 Bs_Jpsiphi,ee=CPV,DecProdCut Bs->Jpsi(ee) phi (with CP Violation and decay products in acceptance) 3.92 * 10^-5 YES Tristan Du Pree 16 Oct 2006
13142400 Bs_Jpsieta,mm,gg Bs->Jpsi(mm) eta(gg) 7.41 * 10^-6 YES Benjamin Carron 5 May 2004
13152400 Bs_Jpsieta,ee,gg Bs->Jpsi(ee) eta(gg) 7.41 * 10^-6 YES Benjamin Carron 5 May 2004
13144400 Bs_Jpsieta,mm,pipipi Bs->Jpsi(mm) eta(3pi) 4.38 * 10^-6 YES Arnaud Robert 26 Apr 2004
13144200 Bs_Jpsietap,mm,rhog Bs->Jpsi(mm) eta'(rho g) 1.14 * 10^-5 YES Matthias Moritz 12 Jul 2004
13144410 Bs_Jpsietap,mm,etapipi Bs->Jpsi(mm) eta'(eta 2pi) 6.48 * 10^-6 YES Sergio Jimenez-Otero 7 May 2004
13442000 Bs_JpsiX,mm Bs->Jpsi(mm) X   YES Gerhard Raven 9 Oct 2006
13442001 Bs_JpsiX,mm=JpsiInAcc Bs->Jpsi(mm) X (mu from J/psi in acceptance)   YES Gerhard Raven 9 Oct 2006
13336000 Bs_etacphi Bs->etac phi 1.61 * 10^-5 YES Luis Fernandez 20 Apr 2004
13264000 Bs_Dspi Bs->Ds pi 1.27 * 10^-4 YES Jeroen van Hunen 20 Apr 2004
13264001 Bs_Dspi=DecProdCut Bs->Ds pi (Decay products in acceptance) 1.27 * 10^-4 YES Jeroen van Hunen 20 Apr 2004
13264400 Bs_Dsrho Bs->Ds rho 3.51 * 10^-4 YES Stephane Monteil 11 Aug 2005
13266000 Bs_Dspipipi Bs->Ds pi pi pi 2.91 * 10^-4 YES Gerhard Raven 19 Sep 2005
13264010 Bs_DsK Bs->Ds K 7.30 * 10^-6 YES Jeroen van Hunen 20 Apr 2004
13264200 Bs_Dsstpi Bs->Ds* pi 1.31 * 10^-4 YES Jeroen van Hunen 6 May 2004
13264210 Bs_DsstK Bs->Ds* K 7.30 * 10^-6 YES Jeroen van Hunen 6 May 2004
13296000 Bs_DsDs Bs->Ds Ds 2.04 * 10^-5 YES Stephane Monteil 20 Apr 2004
13491400 Bs_DsDs,cocktail Bs -> Ds (cocktail) Ds (cocktail) 7.29 * 10^-4 YES Alessio Sarti 4 Nov 2004
13296200 Bs_DsstDs Bs->Ds* Ds 2.13 * 10^-5 YES David Borras 20 Apr 2004
13102010 Bs_pi+K- Bs->K- pi+ 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
13102011 Bs_pi+K-=DecProdCut Bs->K- pi+ (Decay products in acceptance) 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
13102012 Bs_pi+K-=CPV Bs->K- pi+ (Direct CP Violation) 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
13102013 Bs_pi+K-=CPV,DecProdCut Bs->K- pi+ (Direct CP violation and decay products in acceptance) 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
13102000 Bs_K+K- Bs->K+ K- 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
13102001 Bs_K+K-=DecProdCut Bs->K+ K- (Decay products in acceptance) 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
13102002 Bs_K+K-=CPV Bs->K+ K- (With CP Violation) 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
13102003 Bs_K+K-=CPV,DecProdCut Bs->K+ K- (With CP Violation and Decay products in acceptance) 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
13102020 Bs_pi+pi- Bs->pi+ pi- 1.00 * 10^-8 YES Vincenzo Vagnoni 14 Apr 2004
13102030 Bs_pp Bs->p pbar 1.00 * 10^-9      
13104100 Bs_KSKS Bs->KS KS 4.71 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
13104110 Bs_KstKS Bs->K* KS 9.13 * 10^-7 YES Vincenzo Vagnoni 14 Apr 2004
13104000 Bs_Kst0Kst0 Bs->K* K* 1.77 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
13104010 Bs_phiphi Bs->phi phi 1.94 * 10^-6 YES Ivan Belyaev 18 May 2004
13114000 Bs_phimumu Bs->phi mu mu 1.13 * 10^-6 YES Kirill Voronchev 4 Jul 2004
13124000 Bs_phiee Bs->phi e e 1.13 * 10^-6 YES Kirill Voronchev 4 Jul 2004
13102200 Bs_phigamma Bs->phi gamma 2.21 * 10^-5 YES Galina Pakhlova 14 Apr 2004
13112000 Bs_mumu Bs->mu mu 3.50 * 10^-9 YES Franciole da Cunha Marinho 30 Apr 2004
13100200 Bs_gammagamma Bs->gamma gamma 5.00 * 10^-7 YES Thomas Schietinger 14 Apr 2004
13674000 Bs_Dsmunu Bs->Ds mu nu 0.10 % YES Olivier Leroy 30 Oct 2005
13410000 Bs_DsX=cocktail Bs->Ds X   YES Gerhard Raven 9 Oct 2006
13410001 Bs_DsX=cocktail,DsInAcc Bs->Ds X (Decay products from Ds in acceptance)   YES Gerhard Raven 9 Oct 2006
13774000 Bs_Dsmunu=cocktail Bs->Ds(*(*)) (tau->)mu nu 0.4 %      
13774001 Bs_Dsmunu=cocktail,DsmuInAcc Bs->Ds(*(*)) (tau->)mu nu (mu and Ds decay products in acceptance) 0.4 %      
13136200 Bs_psi2Sphi,chic0g,KKpipi Bs->psi(2S)( -> gamma chi_c0 (-> KK pipi)) phi 6.73 * 10^-7 YES Roberta SantaCesaria 22 Feb 2006
13136210 Bs_psi2Sphi,chic0g,pipipipi Bs->psi(2S)( -> gamma chi_c0 (-> pipipipi)) phi 7.65 * 10^-7 YES Roberta SantaCesaria 22 Feb 2006
13136220 Bs_psi2Sphi,chic0g,rhopipi Bs->psi(2S)( -> gamma chi_c0 (-> rho0pipi)) phi 4.88 * 10^-7 YES Roberta SantaCesaria 22 Feb 2006
13144420 Bs_psi2Sphi,Jpsipi0pi0,mm Bs->psi(2S)( -> pi0 pi0 J/psi (-> mu mu)) phi 6.34 * 10^-6 YES Roberta SantaCesaria 22 Feb 2006
13146000 Bs_psi2Sphi,Jpsipipi,mm Bs->psi(2S)( -> pi pi J/psi (-> mu mu)) phi 1.06 * 10^-5 YES Roberta SantaCesaria 22 Feb 2006
13144010 Bs_psi2Sphi,mm Bs->psi(2S)( ->  mu mu) phi 4.13 * 10^-6 YES Roberta SantaCesaria 22 Feb 2006
13112010 Bs_emu Bs -> e mu   YES Ulrik Egede 24 Feb 2006
13122000 Bs_ee Bs -> e+ e-        
13122001 Bs_ee=DecProdCut Bs -> e+ e- (decay products in acceptance)        
13806000 Bs_tautau Bs->tau tau 0      
11144100 Bd_JpsiKS,mm B0->Jpsi(mm) KS 1.82 * 10^-5 YES Herv Terrier 29 Apr 2004
11144101 Bd_JpsiKS,mm=DecProdCut B0->Jpsi(mm) KS (Decay products in acceptance) 1.82 * 10^-5 YES Herv Terrier 29 Apr 2004
11144102 Bd_JpsiKS,mm=CPV B0->Jpsi(mm) KS (with CP violation) 1.82 * 10^-5      
11144103 Bd_JpsiKS,mm=CPV,DecProdCut B0->Jpsi(mm) KS (with CP violation and decay products in acceptance) 1.82 * 10^-5      
11154100 Bd_JpsiKS,ee B0->Jpsi(ee) KS 1.82 * 10^-5 YES Herv Terrier 29 Apr 2004
11154101 Bd_JpsiKS,ee=DecProdCut B0->Jpsi(ee) KS (Decay products in acceptance) 1.82 * 10^-5 YES Herv Terrier 29 Apr 2004
11154102 Bd_JpsiKS,ee=CPV B0->Jpsi(ee) KS (with CP violation) 1.82 * 10^-5      
11154103 Bd_JpsiKS,ee=CPV,DecProdCut B0->Jpsi(ee) KS (with CP violation and decay products in acceptance) 1.82 * 10^-5      
11144000 Bd_JpsiKst,mm B0->Jpsi(mm) K* 5.11 * 10^-5 YES Gerhard Raven 29 Mar 2004
11144001 Bd_JpsiKst,mm=DecProdCut B0->Jpsi(mm) K* (Decay products in acceptance) 5.11 * 10^-5 YES Gerhard Raven 29 Mar 2004
11154000 Bd_JpsiKst,ee B0->Jpsi(ee) K* 5.11 * 10^-5 YES Patrick Robbe 21 Jan 2005
11442000 Bd_JpsiX,mm B0->Jpsi(mm) X   YES Gerhard Raven 9 Oct 2006
11442001 Bd_JpsiX,mm=JpsiInAcc B0->Jpsi(mm) X (mu from J/psi in acceptance)   YES Gerhard Raven 9 Oct 2006
1264000 Bd_D-pi+ B0->D-pi+ 2.70 * 10^-4 YES Vladimir Gligorov 25 Aug 2004
11962000 Bd_Dst-pi+,D0pi,incl B0->D*-(D0~ pi) pi 0.18 % YES Jonas Rademacker 27 Jun 2004
11164000 Bd_Dst-pi+ B0->D*-(D0~(Kpi)pi)pi 7.06 * 10^-5 YES Ulrik Egede 5 May 2004
11166000 Bd_Dst-a1+ B0->D*-(D0~(Kpi)pi)a1 1.54 * 10^-4 YES Agnieszka Oblakowska-Mucha 16 Apr 2004
11164010 Bd_D0Kst,Kpi B0->D0(Kpi) K* 1.22 * 10^-6 YES Miriam Gandelman 6 May 2004
11164020 Bd_D0Kst,KK B0->D0(KK) K* 1.37 * 10^-7 YES Miriam Gandelman 6 May 2004
11164030 Bd_D0Kst,pipi B0->D0(pi pi) K* 4.79 * 10^-8 YES Miriam Gandelman 6 May 2004
11166100 Bd_D0Kst,KSpipi B0->D0(KS pi pi) K* 5.92 * 10^-7 YES Kirill Voronchev 21 Jun 2004
11166110 Bd_D0Kst,KSKK B0->D0(Ks K K) K*   YES Patrick Robbe 11 Oct 2004
11296000 Bd_D-D+ B0->D+D- 2.19 * 10^-6 YES Cristina Lazzeroni 25 Aug 2004
11296001 Bd_D-D+=DecProdCut B0->D+D- (Decay products in acceptance) 2.19 * 10^-6 YES Cristina Lazzeroni 25 Aug 2004
11296010 Bd_DsD B0->Ds+ D- 3.25 * 10^-5 YES David Borras 10 Dec 2005
11102010 Bd_pi+pi- B0->pi+ pi- 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
11102011 Bd_pi+pi-=DecProdCut B0->pi+ pi- (Decay products in acceptance) 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
11102012 Bd_pi+pi-=CPV B0->pi+ pi- (With CP Violation) 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
11102013 Bd_pi+pi-=CPV,DecProdCut B0->pi+ pi- (With CP Violation and Decay products in acceptance) 4.80 * 10^-6 YES Vincenzo Vagnoni 14 Apr 2004
11102000 Bd_K+pi- B0->K+ pi- 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
11102001 Bd_K+pi-=DecProdCut B0->K+ pi- (Decay products in acceptance) 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
11102002 Bd_K+pi-=CPV B0->K+ pi- (Direct CP violation) 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
11102003 Bd_K+pi-=CPV,DecProdCut B0->K+ pi- (Direct CP violation and decay products in acceptance) 1.85 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
11102020 Bd_K+K- B0->K+ K- 3.00 * 10^-7 YES Vincenzo Vagnoni 14 Apr 2004
11202400 Bd_K+pi-pi0 B0->K+ pi- pi0 3.88 * 10^-5 YES Olivier Deschamps 7 Nov 2005
11202401 Bd_K+pi-pi0=DecProdCut B0->K+ pi- pi0 (Decay products in acceptance) 3.88 * 10^-5 YES Olivier Deschamps 7 Nov 2005
11202410 Bd_K+pi-pi0pi0 B0->K+ pi- pi0 pi0 3.36 * 10^-5 YES Olivier Deschamps 7 Nov 2005
11204000 Bd_K+pi-pi+pi- B0->K+ pi- pi+ pi- 3.33 * 10^-5 YES Olivier Deschamps 7 Nov 2005
11202420 Bd_pi+pi-pi0pi0 B0->pi+ pi- pi0 pi0 5.95 * 10^-5 YES Olivier Deschamps 7 Nov 2005
11204010 Bd_pi+pi-pi+pi- B0->pi+ pi- pi+ pi- 1.5 * 10^-5 YES Olivier Deschamps 7 Nov 2005
11102030 Bd_pp B0->p pbar 1.00 * 10^-7      
11104100 Bd_phiKS B0->phi KS 1.69 * 10^-6 YES Yuehong Xie 7 May 2004
11104110 Bd_K+K-KS B0->K+ K- KS 3.43 * 10^-7 YES Yuehong Xie 11 May 2004
11104120 Bd_KSpi+pi- B0->KS pi+ pi-   YES Gabriel Guerrer 14 Sep 2006
11104300 Bd_etapKS,rhog B0->eta'(rho g) KS 6.22 * 10^-6 YES Matthias Moritz 30 Apr 2004
11104500 Bd_etapKS,etapipi B0->eta'(eta pi pi) KS 3.53 * 10^-6 YES Yuehong Xie 13 Feb 2006
11102400 Bd_pipipi=CPV B0->pi+pi-pi0 (CP) 2.17 * 10^-5 YES Arnaud Robert 4 May 2004
11102401 Bd_pipipi=phsp B0->pi+pi-pi0 (PhSp) 2.17 * 10^-5 YES Arnaud Robert 4 May 2004
11102410 Bd_rho+rho- B0->rho+ rho- 2.93 * 10^-5 YES Patrick Robbe 5 May 2004
11104010 Bd_rho0rho0 B0->rho0 rho0 2.00 * 10^-6 YES Patrick Robbe 5 May 2004
11114100 Bd_KSmumu B0->KS (-> pi+ pi-) mu+ mu-   YES Marko Zupan 20 Feb 2006
11114000 Bd_Kstmumu B0->K* mu mu 1.26 * 10^-6 YES Jose Helder Lopes 4 May 2004
11124000 Bd_Kstee B0->K* e e 1.53 * 10^-6 YES Jose Helder Lopes 30 Jul 2004
11102200 Bd_Kstgamma B0->K* gamma 3.00 * 10^-5 YES Galina Pakhlova 01 Apr 2004
11102600 Bd_omegagamma B0->omega gamma 4.30 * 10^-7 YES Galina Pakhlova 14 Apr 2004
11142400 Bd_Jpsipi,mm B0->Jpsi(mm) pi0 1.17 * 10^-6 YES Herve Terrier 9 Jul 2004
11574000 Bd_Dstmunu,Kpi B0->D*- (D0 (K pi) pi-) mu+ nu_mu 1.46 * 10^-3 YES Marta Calvi 23 Sep 2005
11874000 Bd_Dstmunu,Kpi=cocktail B0->D*-(D0 (K pi) pi-) mu+ nu_mu 0.17 % YES Marta Calvi 14 Jul 2006
11874001 Bd_Dstmunu,Kpi=cocktail,D0muInAcc B0->D*-(D0 (K pi) pi-) mu+ nu_mu (mu and D0 decay products in acceptance) 0.17 % YES Marta Calvi 14 Jul 2006
11100400 Bd_pi0pi0 B0->pi0 pi0 2.00 * 10^-6 YES Patrick Robbe 6 Oct 2004
11112010 Bd_emu B0-> e mu   YES Ulrik Egede 24 Feb 2006
12143000 Bu_JpsiK,mm B+->Jpsi(mm) K+ 5.31 * 10^-5 YES Patrick Robbe 7 May 2004
12143001 Bu_JpsiK,mm=DecProdCut B+->Jpsi(mm) K+ (Decay products in acceptance) 5.31 * 10^-5 YES Patrick Robbe 7 May 2004
12153000 Bu_JpsiK,ee B+->Jpsi(ee) K+ 5.31 * 10^-5 YES Patrick Robbe 7 May 2004
12442000 Bu_JpsiX,mm B+->Jpsi(mm) X   YES Gerhard Raven 9 Oct 2006
12442001 Bu_JpsiX,mm=JpsiInAcc B+->Jpsi(mm) X (mu from J/psi in acceptance)   YES Gerhard Raven 9 Oct 2006
12573000 Bu_D0munu,Kpi B+->D0bar(Kpi) mu nu 8.04 * 10^-4 YES Julien Babel 8 Dec 2005
12873000 Bu_D0munu,Kpi=cocktail B+->D(*(*))0bar(Kpi) mu nu 0.396 % YES Julien Babel 4 May 2006
12873001 Bu_D0munu,Kpi=cocktail,D0muInAcc B+->D(*(*))0bar(Kpi) mu nu (mu and D0 decay products in acceptance) 0.396 % YES Julien Babel 4 May 2006
12463010 Bu_D0h,Kpi=cocktail B+->D(*(*))0bar(Kpi) h 0.339 % YES Olivier Leroy 23 Oct 2006
12463011 Bu_D0h,Kpi=cocktail,D0InAcc B+->D(*(*))0bar(Kpi) h (D0 decay products in acceptance) 0.339 % YES Olivier Leroy 23Oct 2006
12165100 Bu_D0K,KSpipi B+->D0(KS pi pi) K+ 8.15 * 10^-6 YES Ulrik Egede 09 Mar 2004
12165101 Bu_D0K,KSpipi=CPV B+->D0(KS pi pi) K+ (with CPV) 8.15 * 10^-6      
12163000 Bu_D0pi,Kpi B+->D0bar(Kpi) pi+ 2.03 * 10^-4 YES Julien Babel 8 Dec 2005
12163001 Bu_D0pi,Kpi=DecProdCut B+->D0bar(Kpi) pi+ (Decay products in acceptance) 2.03 * 10^-4 YES Julien Babel 8 Dec 2005
12265010 Bu_D0pi,Kpipipi B+ -> D0(Kpipipi) pi+   YES James Libby 22 Feb 2006
12163010 Bu_D0K,Kpi B+->D0(Kpi) K+ 1.69 * 10^-5 YES Agnieszka Oblakowska-Mucha 16 Apr 2004
12163020 Bu_D0K,KK B+->D0(KK) K+ 1.89 * 10^-6 YES Agnieszka Oblakowska-Mucha 16 Apr 2004
12163030 Bu_D0K,pipi B+->D0(pipi) K+ 6.60 * 10^-7 YES Agnieszka Oblakowska-Mucha 16 Apr 2004
12265000 Bu_D0K,Kpipipi B+->D0(Kpipipi) K+ 3.46 * 10^-5 YES Vladimir Gligorov 21 Sep 2005
12165000 Bu_D0K,KKpipi B+->D0(KKpipi) K+ 1.1 * 10^-6 YES James Libby 2 Nov 2005
12265100 Bu_D0K,KSKK B+->D0(Ks K K) K+   YES James Libby 22 Feb 2006
12465500 Bu_Dst0K,KSpipi B+->D*0(Ks pi pi) K+ 6.67 * 10^-6 YES James Libby 2 Nov 2005
12463410 Bu_Dst0K,Kpi B+->D*0(K pi) K+ 1.38 * 10^-5 YES James Libby 2 Nov 2005
12103100 Bu_KSpi B+->KS pi+ 6.17 * 10^-6 YES Andy Buckley 4 May 2004
12115100 Bu_Kstmumu B+->K*+ mu+ mu-   YES Marko Zupan 20 Feb 2006
12113010 Bu_pimumu B+->pi+ mu+ mu-        
12103000 Bu_pi+pi-pi+=CPV B+->pi+ pi- pi+ (with resonances) 1.1 * 10^-5 YES Gabriel Guerrer 20 Sep 2005
12103001 Bu_pi+pi-pi+=phsp B+->pi+ pi- pi+ (phase space) 1.1 * 10^-5 YES Gabriel Guerrer 20 Sep 2005
12103011 Bu_K+K-K+=phsp B+ -> K+ K- K+ (phase space) 2.5 * 10^-5 YES James Libby 22 Feb 2006
12103021 Bu_pi+pi-K+=phsp B+ -> pi+ pi- K+ (phase space) 1.0 * 10^-5 YES Pascal David 30 Jan 2006
12203020 Bu_pi+pi-K+=res B+ -> pi+ pi- K+ (with resonances) 1.0 * 10^-5 YES Pascal David 10 Feb 2006
12103031 Bu_pi+K-K+=phsp B+ -> pi+ K- K+ (phase space) 2.0 * 10^-6 YES James Libby 22 Feb 2006
12203410 Bu_pi+pi-pi+pi0 B+->pi+ pi- pi+ pi0 5.43 * 10^-5 YES Olivier Deschamps 7 Nov 2005
12203400 Bu_K+pi-pi+pi0 B+->K+ pi- pi+ pi0 5.43 * 10^-5 YES Olivier Deschamps 7 Nov 2005
12103400 Bu_rho+rho0 B+->rho+ rho0 1.48 * 10^-5 YES Patrick Robbe 5 May 2004
14143000 Bc_Jpsipi,mm Bc->Jpsi(mm) pi 7.67 * 10^-5 YES Oleg Iouchtchenko 10 May 2004
14543000 Bc_Jpsimunu,mm Bc->Jpsi(mm) mu nu 1.12 * 10^-3      
14103000 Bc_rho0pi+ Bc->rho pi 2.00 * 10^-5 YES Oleg Iouchtchenko 10 May 2004
15164100 Lb_D0Lambda Lb->Lambda D0~ 8.35 * 10^-5 YES James Storey 7 Jul 2004
15144100 Lb_JpsiLambda,mm=phsp Lb->Lambda Jpsi(mm) (Phase space) 7.55 * 10^-5 YES Eric Conte 27 May 2004
15144101 Lb_JpsiLambda,mm Lb->Lambda Jpsi(mm) 7.55 * 10^-5 YES Eric Conte 23 Apr 2006
15442000 Lb_JpsiX,mm Lb->X Jpsi(mm)   YES Gerhard Raven 9 Oct 2006
15442001 Lb_JpsiX,mm=JpsiInAcc Lb->X Jpsi(mm) (mu from J/psi in acceptance)   YES Gerhard Raven 9 Oct 2006
15104100 Lb_rhoLambda=phsp Lb->Lambda rho0 (Phase space) 1.36 * 10^-5 YES Eric Conte 2 Jul 2004
15104101 Lb_rhoLambda Lb->Lambda rho0 1.36 * 10^-5 YES Eric Conte 23 Apr 2006
15104102 Lb_rhoLambda=omegaMix Lb->Lambda rho0 (Omega mixing) 1.36 * 10^-5 YES Eric Conte 23 Apr 2006
15104103 Lb_rhoLambda=DecProdCut Lb->Lambda rho0 (Decay products in acceptance) 1.36 * 10^-5 YES Eric Conte 23 Apr 2006
15104110 Lb_phiLambda Lb->Lambda phi 2.67 * 10^-5 YES Eric Conte 2 Jul 2004
15102300 Lb_gammaLambda=pol Lb->Lambda gamma (Pol. Lambda_b) 4.15 * 10^-5 YES Federica Legger 3 Dec 2004
15102301 Lb_gammaLambda=phsp Lb->Lambda gamma (Phase space) 4.15 * 10^-5 YES Federica Legger 9 Mar 2005
15102302 Lb_gammaLambda=pol,onlyLb Lb->Lambda gamma (Only Lambda_b) 4.15 * 10^-5 YES Federica Legger 9 Mar 2005
15102303 Lb_gammaLambda=trpol Lb->Lambda gamma (Transversally Pol. Lambda_b) 4.15 * 10^-5 YES Federica Legger 17 Oct 2005
15102201 Lb_gammaLambda1520 Lb->Lambda(1520) gamma 1.37 * 10^-5 YES Federica Legger 9 Mar 2005
15102210 Lb_gammaLambda1600=pol Lb->Lambda(1600) gamma (Pol. Lambda_b) 1.00 * 10^-5 YES Federica Legger 9 Mar 2005
15102211 Lb_gammaLambda1600=phsp Lb->Lambda(1600) gamma (Phase space) 1.00 * 10^-5 YES Federica Legger 9 Mar 2005
15102213 Lb_gammaLambda1600=trpol Lb->Lambda(1600) gamma (Transversally Pol. Lambda_b) 1.00 * 10^-5 YES Federica Legger 17 Oct 2005
15102220 Lb_gammaLambda1670=pol Lb->Lambda(1670) gamma (Pol. Lambda_b) 5.60 * 10^-6 YES Federica Legger 9 Mar 2005
15102221 Lb_gammaLambda1670=phsp Lb->Lambda(1670) gamma (Phase space) 5.60 * 10^-6 YES Federica Legger 9 Mar 2005
15102223 Lb_gammaLambda1670=trpol Lb->Lambda(1670) gamma (Transversally Pol. Lambda_b) 5.60 * 10^-6 YES Federica Legger 17 Oct 2005
15114100 Lb_Lambdamumu=phsp Lb->Lambda mu mu (Phase space)   YES Ulrik Egede 30 Jan 2006
15104500 Lb_omegaLambda Lb->Lambda omega   YES Eric Conte 24 Feb 2006
15104121 Lb_omegaLambda,pipi Lb->Lambda omega(pi pi)   YES Eric Conte 23 Apr 2006
15102000 Lb_pK Lb->p K 7.88 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
15102010 Lb_ppi Lb->p pi 2.05 * 10^-5 YES Vincenzo Vagnoni 14 Apr 2004
11912000 Bd_smumu B0->s mu mu 4 * 10^-6 YES Marko Zupan 5 May 2004
12912000 Bu_smumu B+->s mu mu 4 * 10^-6 YES Marko Zupan 5 May 2004
23573000 Ds_taunu,mumumu Ds+ -> tau (mu mu mu) nu_tau   YES Ulrik Egede 24 Feb 2006
22112000 D0_mumu D0->mu mu 1.0 * 10^-12 YES Davide Pinci 6 May 2004
30000000 minbias minimum bias        
10000000 incl_b inclusive bb        
10000010 incl_b=Biased biased inclusive bb        
24142000 incl_Jpsi,mm inclusive Jpsi(mm) 0.06 YES Patrick Robbe 10 May 2004
24152000 incl_Jpsi,ee inclusive Jpsi(ee) 0.06 YES Patrick Robbe 10 May 2004
23263000 incl_Ds inclusive Ds(KKpi) 4.87 % YES Patrick Robbe 1 Jul 2005
11900000 Bd_charmless B0 charmless cocktail 0.58 % YES Olivier Deschamps 7 Nov 2005
12900000 Bu_charmless B+ charmless cocktail 0.61 % YES Olivier Deschamps 7 Nov 2005
19912000 bb_mumu b->mu, b->mu   YES Franciole da Cunha Marinho 5 May 2004
19912010 bb_e+mu- b->e+, bbar -> mu-   YES Nicola Serra 25 Feb 2006
19912020 bb_e-mu+ b->e-, bbar -> mu+   YES Nicola Serra 25 Feb 2006
11902000 Bd_phiX Bd->phi X 1.63 % YES Stephanie Moreau 7 Jul 2004
12902000 Bu_phiX Bu->phi X 1.58 % YES Stephanie Moreau 7 Jul 2004
11911000 Bd_muX Bd->mu X 10.61 % YES Patrick Robbe 17 Sep 2004
12911000 Bu_muX Bu->mu X 10.61 % YES Patrick Robbe 17 Sep 2004
13911000 Bs_muX Bs->mu X 9.30 % YES Patrick Robbe 17 Sep 2004
15911000 Lb_muX Lb->mu X 7.70 % YES Patrick Robbe 17 Sep 2004
11973000 Bd_tauX,mumumu Bd -> tau(mu+ mu- mu+) X   YES Ulrik Egede 23 Feb 2006
11972000 Bd_cmu,mu B0->c(mu)mu 0.9894 % YES Franciole da Cunha Marinho 23 Jul 2004
12972000 Bu_cmu,mu B+->c(mu)mu 0.7328 % YES Franciole da Cunha Marinho 23 Jul 2004
13972000 Bs_cmu,mu Bs->c(mu)mu 0.6122 % YES Franciole da Cunha Marinho 23 Jul 2004
15972000 Lb_cmu,mu Lb->c(mu)mu 0.3464 % YES Franciole da Cunha Marinho 23 Jul 2004
11463000 Bd_Dst+X,2body Bd->D*+/-(->D0(Kpi,KK,pipi,piK)pi) X 0.8 % YES Patrick Robbe 13 Feb 2006
12463000 Bu_Dst+X,2body B+ -> D*+/-(->D0(Kpi,KK,pipi,piK)pi) X 0.22 % YES Patrick Robbe 18 Feb 2006
13463000 Bs_Dst+X,2body Bs -> D*+/-(->D0(Kpi,KK,pipi,piK)pi) X 0.13 % YES Patrick Robbe 19 Feb 2006
11444000 Bd_psi2SX,Jpsipipi,mm Bd->psi(2S)(->Jpsi (mm) pi+ pi-) X   YES Alessia Satta 22 Feb 2006
12444000 Bu_psi2SX,Jpsipipi,mm B+ -> psi(2S)(->Jpsi (mm) pi+ pi-) X   YES Alessia Satta 22 Feb 2006
13444000 Bs_psi2SX,Jpsipipi,mm Bs -> psi(2S)(->Jpsi (mm) pi+ pi-) X   YES Alessia Satta 22 Feb 2006
15444000 Lb_psi2SX,Jpsipipi,mm Lambda_b -> psi(2S)(->Jpsi (mm) pi+ pi-) X   YES Alessia Satta 22 Feb 2006
40911000 Higgs_bb=mH115GeV pp->H(bb)W/Z0 (mH = 115 GeV, both b quarks and 1 lepton in acceptance)   YES Patrick Robbe 5 Jan 2005
40911001 Higgs_bb=mH115GeV,onlyl pp->H(bb)W/Z0 (mH = 115 GeV, only 1 lepton in acceptance)   YES Patrick Robbe 7 Jul 2005
40911002 Higgs_bb=mH120GeV pp->H(bb)W/Z0 (mH = 120 GeV, both b quarks and 1 lepton in acceptance)   YES Patrick Robbe 7 Jul 2005
40911003 Higgs_bb=mH120GeV,onlyl pp->H(bb)W/Z0 (mH = 120 GeV, only 1 lepton in acceptance)   YES Patrick Robbe 7 Jul 2005
40911006 Higgs_bb=mH120GeV,onlyl,10GeV pp->H(bb)W/Z0 (mH = 120 GeV, only 1 lepton in acceptance with pT > 10 GeV)   YES Victor Coco 18 Feb 2006
40911004 Higgs_bb=mH140GeV pp->H(bb)W/Z0 (mH = 140 GeV, both b quarks and 1 lepton in acceptance)   YES Patrick Robbe 7 Jul 2005
40911005 Higgs_bb=mH140GeV,onlyl pp->H(bb)W/Z0 (mH = 140 GeV, only 1 lepton in acceptance)   YES Patrick Robbe 7 Jul 2005

41900000

tt_bb pp->t(b) tbar(bbar)   YES Patrick Robbe 7 Feb 2005
41900002 tt_bb=onlyl,10GeV pp->t(b) tbar(bbar) (only 1 lepton in acceptance, pT > 10 GeV)   YES Victor Coco 18 Feb 2006
41900001 tt_bb=only1b pp->t(b) tbar(bbar) (only 1 b in acceptance)   YES Patrick Robbe 7 Jul 2005
41900010 W_tX pp -> W(td, ts, tb) X   YES Patrick Robbe 7 Jul 2005
41900020 t_b pp-> t(b) X   YES Patrick Robbe 7 Jul 2005
42911000 WZ_munum,bb pp-> W(mu nu_mu) Z0(bb)   YES Patrick Robbe 7 Jul 2005
42921000 WZ_enue,bb pp-> W(e nu_e) Z0(bb)   YES Patrick Robbe 7 Jul 2005
42912000 ZZ_mumu,bb pp-> Z0(mu mu) Z0(bb)   YES Patrick Robbe 7 Jul 2005
42922000 ZZ_ee,bb pp-> Z0(e e) Z0(bb)   YES Patrick Robbe 7 Jul 2005
42900000 Z_bb pp-> Z0(bb) X   YES Patrick Robbe 7 Jul 2005
42311000 W_munumu pp-> W(mu nu_mu) X   YES Patrick Robbe 7 Jul 2005
42321000 W_enue pp-> W(e nu_e) X   YES Patrick Robbe 7 Jul 2005
42112000 Z_mumu pp-> Z0(mu mu) X   YES Karol Hennessy 27 Oct 2005
42122000 Z_ee pp-> Z0(e e) X   YES Karol Hennessy 27 Oct 2005
42900010 Z_ff=ffZ pp->Z0(ff) X (ff -> Z0)   YES Patrick Robbe 23 Feb 2006
42900020 Z_ff=ffgZ pp->Z0(ff) X (ff -> gZ0)   YES Patrick Robbe 23 Feb 2006
42900030 Z_ff=fgfZ pp->Z0(ff) X (fg -> fZ0)   YES Patrick Robbe 23 Feb 2006
42900040 W_ff=ffW pp->W(ff') X (ff'->W)   YES Patrick Robbe 23 Feb 2006
42900050 W_ff=ffgW pp->W(ff') X (ff' -> gW)   YES Patrick Robbe 23 Feb 2006
42900060 W_ff=fgfW pp->W(ff') X (fg -> fW)   YES Patrick Robbe 23 Feb 2006
42711000 pp_qqWbb=AcerMC pp->qq(W bb)   YES Victor Coco 10 Oct 2005
42312000 pp_ggZbb=AcerMC pp->gg(Z bb)   YES Victor Coco 10 Oct 2005
42312010 pp_qqZbb=AcerMC pp->qq(Z bb)   YES Victor Coco 10 Oct 2005

Analysis Ntuple

You can analyze directly the file Gauss.sim with a Gaudi algorithm.

You can also use the monitoring algorithm GeneratorFullMonitor which is already present in Gauss. This algorithm is located in the package Sim/GaussMonitor.

To add the algorithm to the execution sequence when generating events with Gauss,

Then run the job like before.

You will obtain a hbook Ntuple GaussMonitor.hbook with the content of the HepMC container "/Event/Gen/HepMCEvents". You can modify the code directly to change the content of the Ntuple and to adapt it to your needs. By default, the content of the Ntuple will be :

******************************************************************
* Ntuple ID = 1 Entries = 10 MCTruth
******************************************************************
* Var numb * Type * Packing * Range        * Block    * Name *
******************************************************************
* 1        * I*4  *         * [0,2000]     * AUTO_BLK * NPart
* 2        * R*4  *         *              * AUTO_BLK * e(NPart)
* 3        * R*4  *         *              * AUTO_BLK * px(NPart)
* 4        * R*4  *         *              * AUTO_BLK * py(NPart)
* 5        * R*4  *         *              * AUTO_BLK * pz(NPart)
* 6        * R*4  *         *              * AUTO_BLK * vxProd(NPart)
* 7        * R*4  *         *              * AUTO_BLK * vyProd(NPart)
* 8        * R*4  *         *              * AUTO_BLK * vzProd(NPart)
* 9        * R*4  *         *              * AUTO_BLK * vtProd(NPart)
* 10       * R*4  *         *              * AUTO_BLK * vxDecay(NPart)
* 11       * R*4  *         *              * AUTO_BLK * vyDecay(NPart)
* 12       * R*4  *         *              * AUTO_BLK * vzDecay(NPart)
* 13       * R*4  *         *              * AUTO_BLK * vtDecay(NPart)
* 14       * I*4  *         * [-2147483648 * AUTO_BLK * pdgId(NPart)
* 15       * I*4  *         * [-2147483648 * AUTO_BLK * nDau(NPart)
* 16       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdMother(NPart)
* 17       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdDau1(NPart)
* 18       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdDau2(NPart)
* 19       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdDau3(NPart)
* 20       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdDau4(NPart)
* 21       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdDau5(NPart)
* 22       * I*4  *         * [-2147483648 * AUTO_BLK * pdgIdDau6(NPart)
* 23       * I*4  *         * [-2147483648 * AUTO_BLK * indexMother(NPart)
* 24       * I*4  *         * [-2147483648 * AUTO_BLK * indexInter(NPart)
* 25       * I*4  *         * [0,20]       * AUTO_BLK * NInter
* 26       * I*4  *         * [-2147483648 * AUTO_BLK * isBB(NInter)
********************************************************************

The corresponding Fortran include file is :

*********************************************************
* *
* This file was generated by HUWFUN. *
* *
*********************************************************
*
* Ntuple Id: 1
* Ntuple Title: MCTruth
* Creation: 19/03/2004 16.15.50
*
*********************************************************
*
      REAL e,px,py,pz,vxProd,vyProd,vzProd,vtProd,vxDecay,vyDecay
     + ,vzDecay,vtDecay
      INTEGER NPart,pdgId,nDau,pdgIdMother,pdgIdDau1,pdgIdDau2,pdgIdDau3
     + ,pdgIdDau4,pdgIdDau5,pdgIdDau6,indexMother,indexInter,NInter,isBB
      COMMON /AUTO_BLK/ NPart,e(2000),px(2000),py(2000),pz(2000)
     + ,vxProd(2000),vyProd(2000),vzProd(2000),vtProd(2000)
     + ,vxDecay(2000),vyDecay(2000),vzDecay(2000),vtDecay(2000)
     + ,pdgId(2000),nDau(2000),pdgIdMother(2000),pdgIdDau1(2000)
     + ,pdgIdDau2(2000),pdgIdDau3(2000),pdgIdDau4(2000),pdgIdDau5(2000)
     + ,pdgIdDau6(2000),indexMother(2000),indexInter(2000),NInter
     + ,isBB(20)
 

If you prefer, you can also make a ROOT Ntuple instead. For this, replace "HBOOK" by "ROOT" and "HbookCnv" by "RootHistCnv" in Monitor.opts and add "use RootHistCnv v*" to cmt/requirements.

Last update: 13 December 2006

Please send questions or comments to Patrick Robbe