About particle PID

[Jung Chang]

Dear Delphes Team,

I want to count hard scattering level b for p p > di-b + diphoton process.
So I used : if(particle->Status==3 && abs(particle->PID)==5) nb++;
and cout nb, but its all 0.
Then I used Example5.C provided in Delphes/Examples,

for(Int_t i=0; i < branchParticle->GetEntriesFast(); i++)

{

GenParticle *gen = (GenParticle*) branchParticle->At(i);
cout<<"N: "<<i<<", St: "<<gen->Status<<", PID: "<<gen->PID<<", E: "<<gen->E<<", Px: "<<gen->Px<<", Py: "<<gen->Py<<", Pz: "<<gen->Pz<<", M: "<<gen->Mass<<", M1: "<<gen->M1<<", M2: "<<gen->M2<<", D1: "<<gen->D1<<", D2: "<<gen->D2<<endl;

}

And found gen->Status are some 44, 52,....

I know Status 1 = stabel particle, 2 = unstable particle. Then how can I select for hard scattering level? Thanks a lot.

Best,
Justine

[Jung Chang]

Dear Delphes Team,

I just found this Pythia8 page which described the status code.
​http://home.thep.lu.se/~torbjorn/pythia81html/ParticleProperties.html

Now my question are solved.