Home » PANDA » PandaRoot » Analysis » [FIXED] PndAnalysis::McTruthMatch always returns false, again ....
[FIXED] PndAnalysis::McTruthMatch always returns false, again .... [message #15003] |
Tue, 23 July 2013 16:25 |
Simone Esch
Messages: 37 Registered: August 2010 Location: Forschungszentrum Jülich...
|
continuous participant |
From: *ext.kfa-juelich.de
|
|
Hi all,
i have again the problem that McTruthMatch always returns false, but at another level.
Now I try to match the whole decaytree of my Lambda LambdaBar decay.
What I do:
RhoCandList ppbarsystem;
ppbarsystem.Combine(lambda, lambdabar);
ppbarsystem.SetType(88888);
for(int b=0;b<ppbarsystem.GetLength();b++)
{
if(theAnalysis->McTruthMatch(ppbarsystem[b],2,true))
{
// fill some histogram
}
}
Output created:
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 88888 and 2 daughters
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 3122 and 2 daughters
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 2212 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=2212)
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: -211 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=-211)
Info in <PndMcTruthMatch::MctMatch>: Particle 3122 reached level of matching, leaves, tree topology and intermediate particle types are matched
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: -3122 and 2 daughters
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: -2212 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=-2212)
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 211 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=211)
Info in <PndMcTruthMatch::MctMatch>: Particle -3122 reached level of matching, leaves, tree topology and intermediate particle types are matched
Info in <PndMcTruthMatch::MctMatch>: rejected by not existing MC truth of daughter zero
This appears for all my llbar combinations.
Seems that the MCTruth tree is not build. Should this work and something is going wrong, or am I doing something wrong?
Best regards
Simone
[Updated on: Fri, 23 August 2013 17:42] by Moderator Report message to a moderator
|
|
|
|
Re: PndAnalysis::McTruthMatch always returns false, again .... [message #15008 is a reply to message #15007] |
Wed, 24 July 2013 09:28 |
Simone Esch
Messages: 37 Registered: August 2010 Location: Forschungszentrum Jülich...
|
continuous participant |
From: *ext.kfa-juelich.de
|
|
Hello Klaus,
yes, I can confirm this. Would it help to insert the PDG code before in the files?
############################### Begin of Event Analysis Exec
-W- CreateMcCandidate: strange PDG code:88888
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2487c0) to candidate (uid=1)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2494d0) to candidate (uid=2)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a24b3a0) to candidate (uid=3)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a24bd90) to candidate (uid=4)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a24c650) to candidate (uid=5)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2ba0a0) to candidate (uid=6)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2ba7c0) to candidate (uid=7)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2baa20) to candidate (uid=8)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2bbf80) to candidate (uid=9)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2bc440) to candidate (uid=10)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a248c50) to candidate (uid=11)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2b6300) to candidate (uid=12)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2b6530) to candidate (uid=13)
Info in <PndAnalysis::BuildMcCands()>: Now setting truth (0x7f9c2a2b7710) to candidate (uid=14)
Info in <PndAnalysis::GetEvent()>: Finished loading event fEvtCount=4.
smallest chi2 263.672
filling histos
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 88888 and 2 daughters
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 3122 and 2 daughters
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 2212 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=2212)
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: -211 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=-211)
Info in <PndMcTruthMatch::MctMatch>: Particle 3122 reached level of matching, leaves, tree topology and intermediate particle types are matched
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: -3122 and 2 daughters
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: -2212 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=-2212)
Info in <PndMcTruthMatch::MctMatch>: Investigating in particle with pdg: 211 and 0 daughters
Info in <PndMcTruthMatch::MctMatch>: accepted final state by PDG code (pdg=211)
Info in <PndMcTruthMatch::MctMatch>: Particle -3122 reached level of matching, leaves, tree topology and intermediate particle types are matched
Info in <PndMcTruthMatch::MctMatch>: rejected by not existing MC truth of daughter zero
In addition I have a new problem.
In the next line my analysis task crashes with the following output :
*** Break *** segmentation violation
sed: RE error: illegal byte sequence
===========================================================
There was a crash.
This is the entire stack trace of all threads:
===========================================================
Thread 1 (process 72203):
#0 0x00007fff8bb906ac in wait4 ()
#1 0x00007fff8c34a03a in system ()
#2 0x000000010dad5f0f in TUnixSystem::StackTrace ()
#3 0x000000010dad3d67 in TUnixSystem::DispatchSignals ()
#4 <signal handler called>
#5 RhoCandidate::NDaughters (this=0x0) at /Users/simone/fairroot/pandaroot_newest/rho/RhoBase/RhoCandidate.cxx:969
#6 0x000000011377929d in PndAnalysis::MctMatch (this=0x7f9c2a1178b0, c=0x0, mct=
0x7f9c2a117a48, level=2, verbose=true) at /Users/simone/fairroot/pandaroot_newest/PndTools/AnalysisTools/PndAnalysis.cxx:783
#7 0x0000000113779680 in PndAnalysis::McTruthMatch (this=0x7f9c2a1178b0, list=
0x7fff521d7f40, level=2, verbose=<value temporarily unavailable, due to optimizations>) at /Users/simone/fairroot/pandaroot_newest/PndTools/AnalysisTools/PndAnalysis.cxx:774
===========================================================
The lines below might hint at the cause of the crash.
If they do not help you then please submit a bug report at
http://root.cern.ch/bugs. Please post the ENTIRE stack trace
from above as an attachment in addition to anything else
that might help us fixing this issue.
===========================================================
#5 RhoCandidate::NDaughters (this=0x0) at /Users/simone/fairroot/pandaroot_newest/rho/RhoBase/RhoCandidate.cxx:969
#6 0x000000011377929d in PndAnalysis::MctMatch (this=0x7f9c2a1178b0, c=0x0, mct=
0x7f9c2a117a48, level=2, verbose=true) at /Users/simone/fairroot/pandaroot_newest/PndTools/AnalysisTools/PndAnalysis.cxx:783
#7 0x0000000113779680 in PndAnalysis::McTruthMatch (this=0x7f9c2a1178b0, list=
0x7fff521d7f40, level=2, verbose=<value temporarily unavailable, due to optimizations>) at /Users/simone/fairroot/pandaroot_newest/PndTools/AnalysisTools/PndAnalysis.cxx:774
===========================================================
This is not correlated with the verbose level, I had this before. I cant see a correlation with my code. Sometimes the codes runs without a crash, sometimes not
Best Regards
Simone
|
|
|
|
|
|
|
|
|
Re: PndAnalysis::McTruthMatch always returns false, again .... [message #15046 is a reply to message #15007] |
Tue, 30 July 2013 15:04 |
Simone Esch
Messages: 37 Registered: August 2010 Location: Forschungszentrum Jülich...
|
continuous participant |
From: *ext.kfa-juelich.de
|
|
Hi Klaus, Hi all
Quote: | I believe this it due to a mismatch of the PDG codes used by EvtGen and those in TDatabasePDG.
|
Can't I insert somehow the 88888 pdg code in the database manualy so that the match works?
I saw in this post:
https://forum.gsi.de/index.php?t=tree&goto=14689&rid=1562&S= ab06a719b66b530ae5b7a8f7a9f6cdd3&srch=TDatabasePDG#msg_14689
that one can do it, but I dont know what to fill in for some the variables.
virtual TParticlePDG* AddParticle(const char* Name,
const char* Title,
Double_t Mass,
Bool_t Stable,
Double_t DecayWidth,
Double_t Charge,
const char* ParticleClass,
Int_t PdgCode,
Int_t Anti=-1,
Int_t TrackingCode=0);
Does someone know?
best regards
simone
|
|
|
Re: PndAnalysis::McTruthMatch always returns false, again .... [message #15050 is a reply to message #15046] |
Tue, 30 July 2013 15:53 |
Klaus Götzen
Messages: 293 Registered: June 2006 Location: GSI
|
first-grade participant |
From: *gsi.de
|
|
Hi Simone,
you could try something like
TDatabasePDG::Instance()->AddParticle("pbarpSystem","pbarpSystem",1.9,kFALSE,0.1,0,"",88888)
I think, for the truth match only the pdg code is important anyways.
Best regards,
Klaus
|
|
|
|
|
|
|
|
Goto Forum:
Current Time: Tue Oct 08 20:06:25 CEST 2024
Total time taken to generate the page: 0.00923 seconds
|