| [Solved]Problem with PndEvtGenDirect [message #12216] |
Wed, 06 July 2011 10:00  |
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  |
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.
|
|
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
