Differences

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

--- gauss:gauss [2011/10/29 15:22] – decianm
+++ gauss:gauss [2012/04/05 11:59] – [How to write out MCParticles in Gauss] che
@@ 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...)
+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.

Physik-Institut

LHCb Wiki Pages

Differences

Physik-Institut

LHCb Wiki Pages

User Tools

Site Tools

Differences

Page Tools