p2->cd();
p2->Draw();
TH1F* his2 = new TH1F("Ptresolution","Ptresolution",40,-5.,5.);
- fTree->Draw("100.*(abs(1./fTPCTrack.Get1Pt())-MC.fTrackRef.Pt())/MC.fTrackRef.Pt()>>Ptresolution","MC.fRowsWithDigits>100&&RC.fTPCTrack.fN>50&&RC.fMultiple==1"+cteta1+cpos1+cprim);
+ fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt()>>Ptresolution","MC.fRowsWithDigits>100&&RC.fTPCTrack.fN>50&&RC.fMultiple==1"+cteta1+cpos1+cprim);
AliLabelAxes(his2, "#Delta p_{t} / p_{t} [%]", "entries");
his2->Fit("gaus");
his2->Draw();
Double_t maxRes = 10.;
TH2F* hRes2 = new TH2F("hRes2", "residuals", nBins, bins, nBinsRes, -maxRes, maxRes);
- fTree->Draw("100.*(abs(1./fTPCTrack.Get1Pt())-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
+ fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
TH1F* hMean=0;
TH1F* hRes = CreateResHisto(hRes2, &hMean);
Double_t maxRes = 10.;
TH2F* hRes2 = new TH2F("hRes2", "residuals", nBins, bins, nBinsRes, -maxRes, maxRes);
- fTree->Draw("100.*(1./fTPCTrack.Get1Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
+ fTree->Draw("100.*(fTPCTrack.Pt()-MC.fTrackRef.Pt())/MC.fTrackRef.Pt():MC.fTrackRef.Pt()>>hRes2", selection, "groff");
TH1F* hMean=0;
TH1F* hRes = CreateResHisto(hRes2, &hMean);
if (drawBinFits) canBinFits->cd(bin + dBin);
TH1D* hBin = hRes2->ProjectionY("hBin", bin, bin);
//
- if (hBin->GetEntries() > 10) {
+ if (hBin->GetEntries() > 5) {
fitFunc->SetParameters(hBin->GetMaximum(),hBin->GetMean(),hBin->GetRMS());
hBin->Fit(fitFunc,"s");
Double_t sigma = TMath::Abs(fitFunc->GetParameter(2));