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
|
|
|