+ //+++++++++++++++++//
+ //p array
+ const Int_t pBins = 36;
+ Double_t nArrayP[pBins+1]={0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, 6.0, 7.0, 8.0, 9.0, 10.0, 15.0, 20.0};
+ //nSigma Array
+ const Int_t nSigmaBins = 250;
+ Double_t nArrayS[nSigmaBins+1];
+ for (Int_t i = 0; i <= nSigmaBins; i++){
+ nArrayS[i]=i-125; //i+1
+ //Printf("nS: %lf - i: %d", nSigmaArray[i], i);
+ }
+
+ fHistNSigmaTPCTOFPbefPID = new TH3D ("fHistNSigmaTPCTOFPbefPID","fHistNSigmaTPCTOFPbefPID;#sigma_{TPC};#sigma_{TOF};p_{T} (GeV/c)", nSigmaBins, nArrayS, nSigmaBins, nArrayS, pBins,nArrayP);
+ fHistListPIDQA->Add(fHistNSigmaTPCTOFPbefPID);
+ //++++++++++++++//
+