GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Monte Carlo Engines » [Solved]Problem with PndEvtGenDirect
[Solved]Problem with PndEvtGenDirect [message #12216] Wed, 06 July 2011 10:00 Go to next message
Tobias Weber is currently offline  Tobias Weber
Messages: 9
Registered: November 2010
occasional visitor
From: *kph.uni-mainz.de
Hi,

I'm simulating the decay X(3872) to mu+mu- pi+pi- in the final state. But I observe a strange behaviour for EvtGenDirect and Geant3. Instead of 4 particles in the final state Geant tells me that I gets 6 particles from EvtGen. These particles seem to be additional muons with p=7Gev/c and Th=Phi=0.
If I use EvtGen to create a output.evt and use this for the simulation everything looks fine. I attaced the decay file, simulation macro and some plots.
I use PandaRoot 11145 and Fairsoft may.

Cheers,
Tobias

[Updated on: Wed, 06 July 2011 11:17]

Report message to a moderator

Re: Problem with PndEvtGenDirect [message #12217 is a reply to message #12216] Wed, 06 July 2011 10:35 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
Hi,
first of all, your macro is "very" old, muon and dirc detectors have changed the initialization. You could take a look into macro/run/tdrct macros.

  PndMdt *Muo = new PndMdt("MDT",kTRUE);
  Muo->SetBarrel("fast");
  Muo->SetEndcap("fast");
  Muo->SetMuonFilter("fast");
  Muo->SetMdtMagnet(kTRUE);
  Muo->SetMdtMFIron(kTRUE);
  fRun->AddModule(Muo);

  PndDrc *Drc = new PndDrc("DIRC", kTRUE);
  Drc->SetGeometryFileName("dirc_l0_p0.root");
  Drc->SetRunCherenkov(kFALSE); // for fast sim Cherenkov -> kFALSE
  fRun->AddModule(Drc);


Second. The problem is when you are storing the original EvtGen tree. Are you seeing many messages such as:
 -W FairPrimaryGenerator: PDG code 88888 not found in database. This warning can be savely ignored.


?
This means that the StoreTree is on, maybe because you are using an old version of the code as your run_sim shows.
The command:

EvtGen->SetStoreTree(kFALSE);


should prevent those entries.
However, those tracks are not reconstructed, you will see them in MCTrack but never in the detectors, they should be quite harmless.
Re: Problem with PndEvtGenDirect [message #12218 is a reply to message #12216] Wed, 06 July 2011 11:16 Go to previous message
Tobias Weber is currently offline  Tobias Weber
Messages: 9
Registered: November 2010
occasional visitor
From: *kph.uni-mainz.de
Thanks Stefano. That solved the problem.
Previous Topic: Group velocity for Cherenkov photon propagation in G3/G4
Next Topic: dipole field in Geant4
Goto Forum:
  


Current Time: Wed Oct 16 02:40:41 CEST 2024

Total time taken to generate the page: 0.00647 seconds