X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=blobdiff_plain;f=FASTSIM%2FAliMUONFastTracking.cxx;h=c38b83e4abf3fd9a17911ea99d8bb7fef114a2c0;hp=943e5837d4fca1d35ad81521eae7ef3b0da35a5a;hb=22120c988b36d3f1890be3f6640c19a39d8c17f1;hpb=9edefa045eaccc2bcb513367c96c63f19a0f98e9 diff --git a/FASTSIM/AliMUONFastTracking.cxx b/FASTSIM/AliMUONFastTracking.cxx index 943e5837d4f..c38b83e4abf 100644 --- a/FASTSIM/AliMUONFastTracking.cxx +++ b/FASTSIM/AliMUONFastTracking.cxx @@ -41,12 +41,15 @@ #include "AliMUONFastTracking.h" #include "AliMUONFastTrackingEntry.h" +using std::cout; +using std::endl; ClassImp(AliMUONFastTracking) AliMUONFastTracking* AliMUONFastTracking::fgMUONFastTracking=NULL; static Double_t FitP(Double_t *x, Double_t *par){ +// Fit function Double_t dx = x[0] - par[0]; Double_t dx2 = x[0] - par[4]; Double_t sigma = par[1] * ( 1 + par[2] * dx); @@ -143,8 +146,8 @@ void AliMUONFastTracking::Init(Float_t bkg) } char filename [100]; - if (fClusterFinder==kOld) sprintf (filename,"$(ALICE_ROOT)/FASTSIM/data/MUONtrackLUT.root"); - else sprintf (filename,"$(ALICE_ROOT)/FASTSIM/data/MUONtrackLUT-AZ.root"); + if (fClusterFinder==kOld) snprintf (filename, 100, "$(ALICE_ROOT)/FASTSIM/data/MUONtrackLUT.root"); + else snprintf (filename, 100, "$(ALICE_ROOT)/FASTSIM/data/MUONtrackLUT-AZ.root"); TFile *file = new TFile(filename); ReadLUT(file); @@ -167,13 +170,13 @@ void AliMUONFastTracking::ReadLUT(TFile* file) const Float_t kBkg[4] = {0, 0.5, 1, 2}; for (Int_t ibkg=0; ibkg<4; ibkg++) { - sprintf (tag,"BKG%g",kBkg[ibkg]); + snprintf (tag, 40, "BKG%g",kBkg[ibkg]); file->cd(tag); for (Int_t isplp = 0; isplpGet(tag2); - sprintf (tag2,"hacc[%d][%d]",isplp,ispltheta); + snprintf (tag2, 40, "hacc[%d][%d]",isplp,ispltheta); hacc[isplp][ispltheta] = (TH3F*)gDirectory->Get(tag2); } } @@ -618,17 +621,17 @@ void AliMUONFastTracking::SetSpline(){ ysigmaphi[ip] = fCurrentEntry[ip][itheta][iphi]->GetSigmaphi(); } if (fPrintLevel>3) cout << " creating new spline " << splname << endl; - sprintf (splname,"fSplineEff[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineEff[%d][%d]",ispline,ivar); fSplineEff[ispline][ivar] = new TSpline3(splname,xsp2,yeff,5 * nbins[ivar]); - sprintf (splname,"fSplineAcc[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineAcc[%d][%d]",ispline,ivar); fSplineAcc[ispline][ivar] = new TSpline3(splname,xsp2,yacc,5 * nbins[ivar]); - sprintf (splname,"fSplineSigmap[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmap[%d][%d]",ispline,ivar); fSplineSigmap[ispline][ivar] = new TSpline3(splname,xspl,ysigmap,nbins[ivar]); - sprintf (splname,"fSplineSigma1p[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigma1p[%d][%d]",ispline,ivar); fSplineSigma1p[ispline][ivar] = new TSpline3(splname,xspl,ysigma1p,nbins[ivar]); - sprintf (splname,"fSplineSigmatheta[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmatheta[%d][%d]",ispline,ivar); fSplineSigmatheta[ispline][ivar] = new TSpline3(splname,xspl,ysigmatheta,nbins[ivar]); - sprintf (splname,"fSplineSigmaphi[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmaphi[%d][%d]",ispline,ivar); fSplineSigmaphi[ispline][ivar] = new TSpline3(splname,xspl,ysigmaphi,nbins[ivar]); ispline++; } @@ -647,17 +650,17 @@ void AliMUONFastTracking::SetSpline(){ ysigmaphi[itheta] = fCurrentEntry[ip][itheta][iphi]->GetSigmaphi(); } if (fPrintLevel>3) cout << " creating new spline " << splname << endl; - sprintf (splname,"fSplineEff[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineEff[%d][%d]",ispline,ivar); fSplineEff[ispline][ivar] = new TSpline3(splname,xspl,yeff, nbins[ivar]); - sprintf (splname,"fSplineAcc[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineAcc[%d][%d]",ispline,ivar); fSplineAcc[ispline][ivar] = new TSpline3(splname,xspl,yacc, nbins[ivar]); - sprintf (splname,"fSplineSigmap[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmap[%d][%d]",ispline,ivar); fSplineSigmap[ispline][ivar] = new TSpline3(splname,xspl,ysigmap,nbins[ivar]); - sprintf (splname,"fSplineSigma1p[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigma1p[%d][%d]",ispline,ivar); fSplineSigma1p[ispline][ivar] = new TSpline3(splname,xspl,ysigma1p,nbins[ivar]); - sprintf (splname,"fSplineSigmatheta[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmatheta[%d][%d]",ispline,ivar); fSplineSigmatheta[ispline][ivar] = new TSpline3(splname,xspl,ysigmatheta,nbins[ivar]); - sprintf (splname,"fSplineSigmaphi[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmaphi[%d][%d]",ispline,ivar); fSplineSigmaphi[ispline][ivar] = new TSpline3(splname,xspl,ysigmaphi,nbins[ivar]); ispline++; } @@ -676,17 +679,17 @@ void AliMUONFastTracking::SetSpline(){ ysigmaphi[iphi] = fCurrentEntry[ip][itheta][iphi]->GetSigmaphi(); } if (fPrintLevel>3) cout << " creating new spline " << splname << endl; - sprintf (splname,"fSplineEff[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineEff[%d][%d]",ispline,ivar); fSplineEff[ispline][ivar] = new TSpline3(splname,xspl,yeff, nbins[ivar]); - sprintf (splname,"fSplineAcc[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineAcc[%d][%d]",ispline,ivar); fSplineAcc[ispline][ivar] = new TSpline3(splname,xspl,yacc, nbins[ivar]); - sprintf (splname,"fSplineSigmap[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmap[%d][%d]",ispline,ivar); fSplineSigmap[ispline][ivar] = new TSpline3(splname,xspl,ysigmap,nbins[ivar]); - sprintf (splname,"fSplineSigma1p[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigma1p[%d][%d]",ispline,ivar); fSplineSigma1p[ispline][ivar] = new TSpline3(splname,xspl,ysigma1p,nbins[ivar]); - sprintf (splname,"fSplineSigmatheta[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmatheta[%d][%d]",ispline,ivar); fSplineSigmatheta[ispline][ivar] = new TSpline3(splname,xspl,ysigmatheta,nbins[ivar]); - sprintf (splname,"fSplineSigmaphi[%d][%d]",ispline,ivar); + snprintf (splname, 40, "fSplineSigmaphi[%d][%d]",ispline,ivar); fSplineSigmaphi[ispline][ivar] = new TSpline3(splname,xspl,ysigmaphi,nbins[ivar]); ispline++; } @@ -773,7 +776,7 @@ TF1* AliMUONFastTracking::GetFitP(Int_t ip,Int_t itheta,Int_t iphi) { // gets the correct prec-pgen distribution for a given LUT cell if (!fFitp[ip][itheta][iphi]) { char name[256]; - sprintf(name, "fit_%d_%d_%d", ip, itheta, iphi); + snprintf(name, 256, "fit_%d_%d_%d", ip, itheta, iphi); fFitp[ip][itheta][iphi] = new TF1(name ,FitP,-20.,20.,6); fFitp[ip][itheta][iphi]->SetNpx(500); fFitp[ip][itheta][iphi]->SetParameters(0.,0.,0.,0.,0.,0.);