Hi,
I'm trying to rewrite my macros to trunk and need some changes. But I got to a point I'm not sure about. I'm checking the mother candidate after using McTruthMatch(RhoCandidate). I'm checking for primary pions directly from pbarpsystem (McIdx=0 or PdgCode=88888). But I want to use the code for my dpm background too, so (since rel apr13) there is no pdgcode for the first particle(just the McIdx=-1). Therefore I used GetMcMotherIdx() and checked if it's lower than zero, like this:
int motherId = mctrk[piplus[ipiplus].GetMcIdx()].GetMcMotherIdx();
if(motherId<0 || mctrk[motherId].PdgCode()==88888) {
So my question is now: Is mother set or not for primary pions?
theAnalysis->FillList(mctrk, "McTruth");
theAnalysis->FillList(piplus, "PionLoosePlus",pidalgos);
...
for (Int_t ipiplus=0;ipiplus<piplus.GetLength();++ipiplus){
if (theAnalysis->McTruthMatch(piplus[ipiplus])){
pidpiplus.Add(piplus[ipiplus]);
RhoCandidate *mother = piplus[ipiplus]->TheMother();
if(!mother || mother->PdgCode()==88888) { <------------
truepiplus.Add(piplus[ipiplus]);
} else {
falsepiplus.Add(piplus[ipiplus]);
}
}
}
Simon
PS: decay is: pbarpsytem ->pi+ pi- ... -> ...
GSI Forum
Members
Search
Help
Register
Login
Home
