Float_t eta = -TMath::Log(TMath::Tan(theta/2.));
Float_t y = 0.5*TMath::Log((E+Pz+1.e-13)/(E-Pz+1.e-13));
- // if (child1 >= 0) continue;
- if (mother > 0) continue;
+ if (child1 >= 0) continue;
if (DataBase->GetParticle(mpart)->Charge() == 0) continue;
Float_t wgt = 1./(Float_t ((evNumber2-evNumber1)+1.));
thetaH->Fill(theta*180./TMath::Pi(),wgt);
phiH->Fill(phi*180./TMath::Pi(),wgt);
else
phiH1->Fill(phi*180./TMath::Pi(),wgt);
- phiH->Fill(phi*180./TMath::Pi(),wgt);
etaH->Fill(eta,5.*wgt);
eetaH->Fill(eta,E);
yH->Fill(y,5.*wgt);