Instead of using the FitSlices command from root to evaluate the pull, I took the distribution of the pullXY and pullZ for different bins in P : &delta P =0.2 GeV/c
Then for each bin, a double gaussian was better than a simple or a simple gaussian+pol2