pmom=fPMin+random[1]*(fPMax-fPMin);
theta=fThetaMin+random[2]*(fThetaMax-fThetaMin);
phi=fPhiMin+random[3]*(fPhiMax-fPhiMin);
p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
p[2] = pmom*TMath::Cos(theta);
pmom=fPMin+random[1]*(fPMax-fPMin);
theta=fThetaMin+random[2]*(fThetaMax-fThetaMin);
phi=fPhiMin+random[3]*(fPhiMax-fPhiMin);
p[0] = pmom*TMath::Cos(phi)*TMath::Sin(theta);
p[1] = pmom*TMath::Sin(phi)*TMath::Sin(theta);
p[2] = pmom*TMath::Cos(theta);