Fix Coverity report

[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemtoUser / AliFemtoShareQualityCorrFctn.cxx

diff --git a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
index 045c9fb740f3ae44cb8da2a94f91a985d7f54455..47115efab799db52f0ce2f90da45d0835335a9a2 100644 (file)
--- a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
@@ -26,22 +26,22 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(char* title, const in
 {
   // set up numerator
   //  title = "Num Qinv (MeV/c)";
-  char tTitNum[100] = "NumShare";
-  strcat(tTitNum,title);
+  char tTitNum[101] = "NumShare";
+  strncat(tTitNum,title, 100);
   fShareNumerator = new TH2D(tTitNum,title,nbins,QinvLo,QinvHi,100,0.0,1.00001);
   // set up denominator
   //title = "Den Qinv (MeV/c)";
-  char tTitDen[100] = "DenShare";
-  strcat(tTitDen,title);
+  char tTitDen[101] = "DenShare";
+  strncat(tTitDen,title, 100);
   fShareDenominator = new TH2D(tTitDen,title,nbins,QinvLo,QinvHi,100,0.0,1.00001);
 
-  char tTit2Num[100] = "NumQuality";
-  strcat(tTit2Num,title);
+  char tTit2Num[101] = "NumQuality";
+  strncat(tTit2Num,title, 100);
   fQualityNumerator = new TH2D(tTit2Num,title,nbins,QinvLo,QinvHi,150,-0.500001,1.000001);
   // set up denominator
   //title = "Den Qinv (MeV/c)";
-  char tTit2Den[100] = "DenQuality";
-  strcat(tTit2Den,title);
+  char tTit2Den[101] = "DenQuality";
+  strncat(tTit2Den,title, 100);
   fQualityDenominator = new TH2D(tTit2Den,title,nbins,QinvLo,QinvHi,150,-0.500001,1.000001);
   // set up ratio
   //title = "Ratio Qinv (MeV/c)";
@@ -51,13 +51,13 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(char* title, const in
   //mShareDenominator->SetDirectory(0);
   //mRatio->SetDirectory(0);
 
-  char tTit3Num[100] = "NumTPCSep";
-  strcat(tTit3Num,title);
+  char tTit3Num[101] = "NumTPCSep";
+  strncat(tTit3Num,title, 100);
   fTPCSepNumerator = new TH2D(tTit3Num,title,nbins,QinvLo,QinvHi,150,0.0,100.0);
   // set up denominator
   //title = "Den Qinv (MeV/c)";
-  char tTit3Den[100] = "DenTPCSep";
-  strcat(tTit3Den,title);
+  char tTit3Den[101] = "DenTPCSep";
+  strncat(tTit3Den,title, 100);
   fTPCSepDenominator = new TH2D(tTit3Den,title,nbins,QinvLo,QinvHi,150,0.0,100.0);
 
   // to enable error bar calculation...
@@ -77,7 +77,9 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(const AliFemtoShareQu
   fShareNumerator(0),
   fShareDenominator(0),
   fQualityNumerator(0),
-  fQualityDenominator(0)
+  fQualityDenominator(0),
+  fTPCSepNumerator(0),
+  fTPCSepDenominator(0)
 {
   // copy constructor
   if (aCorrFctn.fShareNumerator)
@@ -154,9 +156,9 @@ AliFemtoString AliFemtoShareQualityCorrFctn::Report(){
   // create report
   string stemp = "Qinv Correlation Function Report:\n";
   char ctemp[100];
-  sprintf(ctemp,"Number of entries in numerator:\t%E\n",fShareNumerator->GetEntries());
+  snprintf(ctemp , 100, "Number of entries in numerator:\t%E\n",fShareNumerator->GetEntries());
   stemp += ctemp;
-  sprintf(ctemp,"Number of entries in denominator:\t%E\n",fShareDenominator->GetEntries());
+  snprintf(ctemp , 100, "Number of entries in denominator:\t%E\n",fShareDenominator->GetEntries());
   stemp += ctemp;
   //  stemp += mCoulombWeight->Report();
   AliFemtoString returnThis = stemp;
@@ -165,6 +167,9 @@ AliFemtoString AliFemtoShareQualityCorrFctn::Report(){
 //____________________________
 void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
   // add real (effect) pair
+  if (fPairCut)
+    if (!fPairCut->Pass(pair)) return;
+
   double tQinv = fabs(pair->QInv());   // note - qInv() will be negative for identical pairs...
   Int_t nh = 0;
   Int_t an = 0;
@@ -289,6 +294,9 @@ void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
 //____________________________
 void AliFemtoShareQualityCorrFctn::AddMixedPair( AliFemtoPair* pair){
   // add mixed (background) pair
+  if (fPairCut)
+    if (!fPairCut->Pass(pair)) return;
+
   double weight = 1.0;
   double tQinv = fabs(pair->QInv());   // note - qInv() will be negative for identical pairs...
   Int_t nh = 0;