Comment retrouver les particules mères dans la branche GEN

[Xavier Rouby]

But: retrouver les mères et les filles d'une particule dans la branche GEN. Jérémy Andréa, 17/02/2010

Salut Severine,

je suis un postdoc qui decouvre les joies de Delphes. Je suis en train d'essayer de décoder l'info du génerateur et en particulier, pour une TRootC::GenParticle donnée, j'essaye de récupérer le GenParticle parent.

Aurais-tu par hasard un example de code qui ferait ça ??

Merci bien,
Jeremy

[Xavier Rouby]

Voici la réponse:

TSimpleArray<TRootGenParticle> LeptWH(const TClonesArray *GEN)
{

TSimpleArray<TRootGenParticle> array,array2;
TIter itGen((TCollection*)GEN);
TRootGenParticle *gen;
itGen.Reset();
while( (gen = (TRootGenParticle*) itGen.Next()) ) Loop over all particles in event

{

array.Add(gen);

}

itGen.Reset();
TRootGenParticle *gen1;
while( (gen1 = (TRootGenParticle*) itGen.Next()) ) Loop over all particles in event

{

Int_t Lmoth = gen1->M1;

if( (abs(gen1->PID) == pE)

(abs(gen1->PID) == pMU)

(abs(gen1->PID) == pTAU))

{

if(Lmoth < array.GetEntries() && Lmoth > 0 && abs(array[Lmoth]->PID)==24)

{

Int_t Wmoth = array[Lmoth]->M1;
if(Wmoth < array.GetEntries() && Wmoth > 0 && abs(array[Wmoth]->PID)==25)
array2.Add(gen1);

}

}

}

return array2;

}