GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Monte Carlo Engines » MC true match with geant4
MC true match with geant4 [message #15626] Thu, 31 October 2013 12:14 Go to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *kph.uni-mainz.de
Dear all,

I am wondering how one can handle MC true match for Geant4 case.
I am using TGeant4 and PndEvtGenDirect for evtgen generator.

If I turn off saving the event tree,
evtGen->SetStoreTree(kFALSE);
then simulation process works fine.

But If I turn on StoreTree to use MCtruth match later on for signal MC analysis,
evtGen->SetStoreTree(kTRUE);
then Geant4 cannot recognize this event tree due to different particle definition.

How can we avoid this problem?

We need a event tree to make a truth matching in analysis level.
It must to be saved!

Do you have any idea or solution for this problem?


Quote:


### Run 0 start.
-W FairPrimaryGenerator: PDG code 88881 not found in database. This warning can be savely ignored.
Warning in <TParticle::SetPdgCode>: PDG code 88881 unknown from TDatabasePDG
[INFO ] FairPrimaryGenerator: (Event 1) 12 primary tracks from vertex (0.000000, 0.000000, 0.000000 ) with beam gradiant (0.000000, 0.000000) Event Time = 0.013075 (ns)

TG4PrimaryGeneratorAction::TransformPrimaries:
G4ParticleTable::FindParticle() failed for XXX pdgEncoding=88881.
*** TG4Exception: Aborting execution ***

WARNING - Attempt to delete the physical volume store while geometry closed !
WARNING - Attempt to delete the logical volume store while geometry closed !
WARNING - Attempt to delete the solid store while geometry closed !
WARNING - Attempt to delete the region store while geometry closed !
Warning in <TStreamerInfo::Build:>: TStreamerBase: base class TStreamerElement has no streamer or dictionary it will not be saved
root [0]

Re: MC true match with geant4 [message #15632 is a reply to message #15626] Fri, 01 November 2013 12:39 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 2.235.190*
This is an old and well known problem, and up to now I was not able to find a solution. Maybe there are some problems in our mcapplication, since the problem does not appear with g3, but up to know nothing was found. If you need to store the tree then you must use g3 for the moment.
Re: MC true match with geant4 [message #15711 is a reply to message #15632] Tue, 12 November 2013 18:27 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
Please try the following:

in your sim macro add the following lines after FairRunSim:

FairParticle *pbarp = new FairParticle(88888,"pbarpSystem",kPTUndefined,1.9,0,0);
fRun->AddNewParticle(pbarp);


You should see some warnings with G4, but the simulation does not crash. Most probably you need to do the same with your analysis macro.
Please tell me if it works and if the results are reasonable.
Re: MC true match with geant4 [message #15721 is a reply to message #15711] Fri, 15 November 2013 15:20 Go to previous messageGo to next message
donghee is currently offline  donghee
Messages: 385
Registered: January 2009
Location: Germnay
first-grade participant
From: *dip0.t-ipconnect.de
Hi Stefano,

I tried to add new particle via FairParticle.
Upto sim->dig->rec->pid it works fine with Geant4,
but I cannot see anything about 88888 in the MCtrack list.

Best wishes,
Donghee
Re: MC true match with geant4 [message #15733 is a reply to message #15721] Fri, 15 November 2013 16:41 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *to.infn.it
It seems there are problems in our MC application for g4, Mohammad will take a look.
Previous Topic: Reco error at Geant4 simulation
Next Topic: Difference of neutral candidate between Geant4 and Geant3
Goto Forum:
  


Current Time: Sat Nov 26 22:53:20 CET 2022

Total time taken to generate the page: 0.00802 seconds