User Tools

Site Tools


gauss:gauss

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Next revision
Previous revision
Last revisionBoth sides next revision
gauss:gauss [2011/10/29 15:22] – created decianmgauss:gauss [2012/04/05 11:59] – [How to write out MCParticles in Gauss] che
Line 2: Line 2:
 ===== How to write out MCParticles in Gauss ===== ===== How to write out MCParticles in Gauss =====
 It can be very useful to let Gauss write out MCParticles in the Generation directly (without the detector simulation), which then can be analysed. It can be very useful to let Gauss write out MCParticles in the Generation directly (without the detector simulation), which then can be analysed.
-This can be achieved the following way:+This can be achieved the following way (write this at the end of your version of Gauss-Job.py):
 <code python> <code python>
 gtos = GenerationToSimulation() gtos = GenerationToSimulation()
Line 15: Line 15:
 The root file can then be run over with an MCTupleTool in DaVinci (with certain restrictions, i.e. only tools which only access 4-vectors will run...) The root file can then be run over with an MCTupleTool in DaVinci (with certain restrictions, i.e. only tools which only access 4-vectors will run...)
  
 +To avoid error messages like: 
 +
 +MainEventGaussSim         FATAL SimInit:: Exception throw: put():: could not register 'LHCb::MCHeader' at address '/Event/MC/Header' StatusCode=2
 +
 +you have to make sure that you have included
 +<code python>
 +Gauss().Phases = ["Generator"]
 +Gauss().Output = 'NONE'
 +Gauss().Histograms = 'NONE'
 +</code>
 +in your Gauss-Job.py file.
 +
 +===== Running EvtGen standalone =====
 +You can run EvtGen standalone, letting the particles decay at rest. This is very fast, as no underlying event has to be generated.
 +A prototype file could look like (Gauss v41r1):
 +<code python>
 +#Event Type: 11114001
 +from Configurables import Generation
 +Generation().EventType = 11114001
 +from Configurables import StandAloneDecayTool
 +Generation().SampleGenerationTool = "StandAloneDecayTool";
 +Generation().addTool( StandAloneDecayTool )
 +Generation().StandAloneDecayTool.ProductionTool = "PythiaProduction";
 +
 +from Configurables import ToolSvc
 +from Configurables import EvtGenDecay
 +
 +ToolSvc().addTool( EvtGenDecay )
 +ToolSvc().EvtGenDecay.UserDecayFile ="$DECFILESROOT/dkfiles/Bd_Kstmumu=DecProdCut.dec"
 +Generation().StandAloneDecayTool.SignalPIDList = [ 511,-511 ]
 +</code>
 +
 +These both things (running EvtGen standalone and writing out MCParticles) can of course be combined.
gauss/gauss.txt · Last modified: 2012/04/10 16:16 by che