GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Bugs, Fixes, Releases » MC truth (cont.)
MC truth (cont.) [message #16039] Tue, 25 March 2014 12:10 Go to next message
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 #16040 is a reply to message #16039] Tue, 25 March 2014 12:26 Go to previous messageGo to next message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: 91.255.27*
Are these primary pions or this number includes also secondaries?
Re: MC truth (cont.) [message #16043 is a reply to message #16040] Tue, 25 March 2014 16:02 Go to previous messageGo to next message
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 #16051 is a reply to message #16039] Tue, 25 March 2014 17:43 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *gsi.de
Hi Karin,


I think simply counting the number of pi+ is not sufficient. My last observation was btw, that the MC truth doublers are still present, but reduced since the track doublers have been removed.

However, checking how often more than one reco track points to the same MC truth object (i.e. GetMcTruth() delivers the same pointer) would be a more useful number.


Best,
Klaus
Re: MC truth (cont.) [message #16071 is a reply to message #16043] Wed, 26 March 2014 17:28 Go to previous messageGo to next message
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 Go to previous message
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.
Previous Topic: [OBSOLETE] run_sim_tpcorsttcombi.C
Next Topic: [CLOSED] pandaroot/macro/run error
Goto Forum:
  


Current Time: Thu Dec 02 11:14:02 CET 2021

Total time taken to generate the page: 0.02065 seconds