ENH: Fit WEPL to angle standard deviation in pctweplfit

[acoussat]

peux-tu (pas trop difficilement) rajouter à pctweplfit un fit sur l'écart type des angles ? Au lieu de regarder la moyenne de l'énergie, tu regardes l'écart type des angles de déviation comme dans pctpaircuts :
const double anglex = std::acos(std::min(1., dInX * dOutX / (dInX.GetNorm() * dOutX.GetNorm())));
const double angley = std::acos(std::min(1., dInY * dOutY / (dInY.GetNorm() * dOutY.GetNorm())));

Would something like this work for you? The code is getting slightly messier but I guess we could still refactor it later on.