| MC truth (cont.) [message #16039] |
Tue, 25 March 2014 12:10  |
Karin Schönning
Messages: 65
Registered: August 2012
Location: Uppsala University
|
continuous participant |
From: *physics.uu.se
|
|
Hi again,
I asked previously about MC truth PID, then we thought the issue was solved but now it appears that it is not solved after all. Using MC truth PID and matching with a MC truth mother in the following way:
int ipiprec =0;
int ipipmoth=0;
mcidpip.SetType(211);
for (k=0; k<mcidpip.GetLength(); ++k) {
if ( theAnalysis->McTruthMatch(mcidpip[k]) ) {
RhoCandidate *piptr=mcidpip[k]->GetMcTruth();
if(piptr && piptr->PdgCode()==211){
ipiprec++;
RhoCandidate *mother = 0;
mother = piptr->TheMother();
if (mother && mother->PdgCode()==411)
{
ipipmoth++;
}
}
}
}
for Psi'(3770)->D+D-, D+->K-pi+pi+ and D-->K+pi-pi-
gives in about 12% of the events, cases with more than 2 pi+ or pi- or more than 1 K+ or K- in the final state. I have the same problem in the January 2014 release as in the trunk version of pandaroot.
Did anyone else still have this problem?
Best regards,
/Karin
|
|
|
|
|
|
| Re: MC truth (cont.) [message #16043 is a reply to message #16040] |
Tue, 25 March 2014 16:02  |
Karin Schönning
Messages: 65
Registered: August 2012
Location: Uppsala University
|
continuous participant |
From: *physics.uu.se
|
|
I thought that this
mother = piptr->TheMother();
if (mother && mother->PdgCode()==411)
was supposed to make sure that the pi+ comes from a D+ decay?
|
|
|
|
|
|
| Re: MC truth (cont.) [message #16071 is a reply to message #16043] |
Wed, 26 March 2014 17:28 |
Karin Schönning
Messages: 65
Registered: August 2012
Location: Uppsala University
|
continuous participant |
From: *physics.uu.se
|
|
|
There seems to be both reco doublets and different reconstructed tracks that point to the same MC truth object.
|
|
|
|
| Re: MC truth (cont.) [message #16072 is a reply to message #16051] |
Thu, 27 March 2014 11:21 |
Karin Schönning
Messages: 65
Registered: August 2012
Location: Uppsala University
|
continuous participant |
From: *physics.uu.se
|
|
There are no reco doublets (checked my code again and made a correction) but there in about 19% of the Psi(3770)->D+D-->K-2pi+K+2pi- events, it happens at least once that two different reco tracks point to the same MC truth track.
About 1-2% of the kaon MC truth tracks have more than one reco track pointing to it, and for about 4-5% of the pion tracks.
|
|
|
|
GSI Forum
Search
Help
Members
Login
Home
