TWiki> LHCb Web>LHCbComputing>LHCbSimulation (revision 19)EditAttachPDF

Simulation in LHCb

Introduction

Gauss is the LHCb Simulation software, built on the Gaudi Framework.
It consist of a first phase where the events are generated (e.g. pp collisions at 14 TeV) and a second phase where the particles are propagated trough the LHCb detectors.

Coordinators

The Gauss release manager is Gloria Corti. The list of topics and coordinators/contacts is here.

Gauss Web

The main Gauss web page containing release information and links to user guides and tutorials, is here.

Documentation

FAQ

Gauss Tutorial

There are 3 tutorials for self-teaching. See here.

Gauss User Guide

Ongoing tasks/projects, open issues

See:

Survey geometry for simulation

Last update: 8th July 2008

It is agreed that, ultimately, the 2008 simulation should simulate a realistic detector based on the surveyed detector, not the ideal ("Optimisation TDR") detector simulated in DC06
See here for the various issues

Status

Detector SQLDDDB tags Gauss version Contents and Comments Responsible
VELO velo-20080611 Gauss v33r2 (changing options) Marco Cl. received some recently from Marco G., it should be final one: is it perfect in DDDB and survey in LHCBCOND, stepping motor status in ONLINE not yet ready David/Marco G/Silvia
TT     Ideal final geometry in DDDB, Survey cannot be used, overlaps ok for LHCb and G4 Jeroen/Cristophe
IT it-20080613, it-20080715 v34r0, v34r1 Ideal geometry in DDDB. Survey cannot be used, overlaps ok for LHCb, to be checked for g4, some fixes in latest it-20080715 Vincent/Matt
OT     Restructuring with multiple LV for stations and layers, overlaps ok in LHCb (jan), to do for G4 Jan
RICH     Still small work, ideal and survey. The conditions db being filled with measured properties of various things like hpds. Also in dddb, the properties are being added for aerogel, rich2 mirror surface. For the rich1 and rich2 hpd boxes location survey info not available yet. We may add them in dddb or condb depending upon the situation. Sajan/Antonis
Calorimeters calo-20080708 v34r1 calo split in 2 halves done, real position for plugs and z ?, overlaps ok in LHCb. to be done in G4 Vladimir/Olivier D.
Muon     survey of chamber to bein DDDB, half stations final positions in SIMCOND, M1 open to be in DDDB(all chambers there), off in SIMCOND Simone/Alessia
BCM head-20080603 not used Ideal in DDDB, cables? Magnus
Magnet and Beam Pipe     To do as installed in DDDB: Bakeout in RICH removed, collars in magnet, bellow in VELO window, and clear overlaps. Daria/Gloria
Magnetic field map     Current and polarity in CondDB to be in SIMCOND first implementation, to be taken in MagneticFieldSvc. Measured maps for one polarity in FieldMap v5r0, other maps in preparation Adlene
RMS     togheter with IT2 already in Kiev(/Vincent/Matt)

Conditions for simulation

Last update: 4th July 2008

The simulation (and reconstruction of simulated data) will get a copy of the necessary online conditions from a SIMCOND database.

See here for the various issues

Marco Cl. is going to provide a first version of SIMCOND containing the existing ONLINE data structure, and the conditions included in the latest LHCBCOND tag. There will also be a version with all misalignements set to zeros. The first version SIMCOND will contain the magnet variables from online. Propose to have all directory in ONLINE copied to SIMCOND: will need to decide what values for which tag!

CondDB release notes

Samples for validation

Last update: 7th July 2008

A set of samples are necessary for validating a new production version, new versions of Geant4 and of generators. The production of these samples should become automated and the results published on a dedicated web page. This samples are produced by the production system, so that the all chain is exercised. It should also be investigated if the Data Quality/Monitoring infrastructure is suitable for this. Finally a new facility called QMTests is available within the Gaudi context, its suitability for this should also be investigated.

Two different macro-sets are produced depending on what need to be validated:

  • particle-gun and min bias samples (only sim files or sim+digi)
  • physics samples (dst files)

The list of sample to be produced is below:

Sample Type Energy Physics(list and cuts) No. of events !DecFiles/EvtType Status Responsible
ParticleGun pi+/pi- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay off + KineHT 80 K for each energy v15r1, 53230100, 53230200, 53230400, 53210001, 53210005, 53210010, 53210017, 53210034, 53210050, 53210100, 53210120, 53210168 Ask production with Gauss v35r0
e+/e- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay off + KineHT 80 K for each energy v15r1, 51230100, 51230200, 51230400, 51210001, 51210005, 51210010, 51210017, 51210034, 51210050, 51210100, 51210120, 51210168 Ask production with Gauss v35r0
mu+/mu- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay off + KineHT 50 K for each energy v15r1, 52230100, 52230200, 52230400, 52210001, 52210005, 52210010, 52210017, 52210034, 52210050, 52210100, 52210120, 52210168 Ask production with Gauss v35r0
K+/K- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay off + KineHT 50 K for each energy v15r1, 54230100, 54230200, 54230400, 54210001, 54210005, 54210010, 54210017, 54210034, 54210050, 54210100, 54210120, 54210168 Ask production with Gauss v35r0
p/p~ 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay off + KineHT 50 K for each energy v15r1, 55230100, 55230200, 55230400, 55210001, 55210005, 55210010, 55210017, 55210034, 55210050, 55210100, 55210120, 55210168 Ask production with Gauss v35r0
ParticleGun J/Psi 20 GeV < E < 140 GeV LHEP + EM option1 + DRay off + KineHT 50 K    
MinimumBias     LHEP + EM option1 + DRay off + KineHT 5 K    
ParticleGun pi+/pi- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay on+ KineHT 80 K for each energy v15r1, 53230100, 53230200, 53230400, 53210001, 53210005, 53210010, 53210017, 53210034, 53210050, 53210100, 53210120, 53210168 To do Options for physics
e+/e- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay on + KineHT 80 K for each energy v15r1, 51230100, 51230200, 51230400, 51210001, 51210005, 51210010, 51210017, 51210034, 51210050, 51210100, 51210120, 51210168 To do Options for physic
mu+/mu- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay on + KineHT 50 K for each energy v15r1, 52230100, 52230200, 52230400, 52210001, 52210005, 52210010, 52210017, 52210034, 52210050, 52210100, 52210120, 52210168 To do Options for physics
K+/K- 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay on + KineHT 50 K for each energy v15r1, 54230100, 54230200, 54230400, 54210001, 54210005, 54210010, 54210017, 54210034, 54210050, 54210100, 54210120, 54210168 To do Options for physics
p/p~ 0.1 GeV, 0.2 GeV, 0.4 GeV, 1 GeV, 5 GeV, 10 GeV, 16.8 GeV , 33.8 GeV , 50 GeV ,100 GeV, 120 GeV, 168 GeV LHEP + EM option1 + DRay on + KineHT 50 K for each energy v15r1, 55230100, 55230200, 55230400, 55210001, 55210005, 55210010, 55210017, 55210034, 55210050, 55210100, 55210120, 55210168 To do Options for physics
ParticleGun J/Psi 20 GeV < E < 140 GeV LHEP + EM option1 + DRay on + KineHT 50 K    
MinimumBias     LHEP + EM option1 + DRay on + KineHT 5 K    
MinimumBias     QGSP + EM option1 + DRay on + KineHT 5 K    

Legenda for Physics (lists and cuts)

  • LHEP parameterised modeling of hadronic interactions (G4)
  • QGSP modeling using Quark Gluon String model for high energy hadronic interactions of protons, neutrons, pions, and Kaons. The high energy interaction creates an exited nucleus, which is passed to the precompound model modeling the nuclear de-excitation. (G4)
  • EM option1 EM processes using parameters optimised for better CPU performance at the cost slightly less precise physics (G4)
  • EM standard EM processes with default setting using standard parameters tuned for best physics performance (G4)
  • DRay off production cuts as follow: e- 10 km, e+ 5 mm, gamma 10 mm (G4)
  • DRay on production cuts as follow: e- 5 mm, e+ 5 mm, gamma 5 mm (G4)
  • KineHT tracking cuts (kinetic E) as follow: hadrons/mu+/mu- 10 MeV, e+/e-/gamma 1 MeV, others none (Gauss)
  • KineLT tracking cuts (kinetic E) as follow: hadrons/mu+/mu- 1 MeV, e+/e-/gamma 0.1 MeV, others none (Gauss)

Validation of generators

Plots for different generator and Gauss versions can be found here

Migration to HepMC2

Migration done in Gauss v35r0

Recipies on how to migrate user code from HepMC 1.26 to HepMC 2.03 are available

CLHEP::HepLorentzVector and the new class HepMC::FourVector

Next versions

Last update : 24 October 2008

-- GloriaCorti - 04 Jul 2008

Edit | Attach | Watch | Print version | History: r106 | r21 < r20 < r19 < r18 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r19 - 2008-11-14 - GloriaCorti
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    LHCb All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback