]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
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 07c2fc2c7605fb5ee918af101d234eb81f62b89e..47115efab799db52f0ce2f90da45d0835335a9a2 100644 (file)
--- a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
@@ -20,26 +20,28 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(char* title, const in
   fShareNumerator(0),
   fShareDenominator(0),
   fQualityNumerator(0),
-  fQualityDenominator(0)
+  fQualityDenominator(0),
+  fTPCSepNumerator(0),
+  fTPCSepDenominator(0)
 {
   // 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)";
@@ -49,12 +51,24 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(char* title, const in
   //mShareDenominator->SetDirectory(0);
   //mRatio->SetDirectory(0);
 
+  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[101] = "DenTPCSep";
+  strncat(tTit3Den,title, 100);
+  fTPCSepDenominator = new TH2D(tTit3Den,title,nbins,QinvLo,QinvHi,150,0.0,100.0);
+
   // to enable error bar calculation...
   fShareNumerator->Sumw2();
   fShareDenominator->Sumw2();
 
   fQualityNumerator->Sumw2();
   fQualityDenominator->Sumw2();
+  
+  fTPCSepNumerator->Sumw2();
+  fTPCSepDenominator->Sumw2();
 }
 
 //____________________________
@@ -63,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)
@@ -74,6 +90,10 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(const AliFemtoShareQu
     fQualityNumerator = new TH2D(*aCorrFctn.fQualityNumerator);
   if (aCorrFctn.fQualityDenominator)
     fQualityDenominator = new TH2D(*aCorrFctn.fQualityDenominator);
+  if (aCorrFctn.fTPCSepNumerator)
+    fTPCSepNumerator = new TH2D(*aCorrFctn.fTPCSepNumerator);
+  if (aCorrFctn.fTPCSepDenominator)
+    fTPCSepDenominator = new TH2D(*aCorrFctn.fTPCSepDenominator);
 }
 //____________________________
 AliFemtoShareQualityCorrFctn::~AliFemtoShareQualityCorrFctn(){
@@ -82,6 +102,8 @@ AliFemtoShareQualityCorrFctn::~AliFemtoShareQualityCorrFctn(){
   delete fShareDenominator;
   delete fQualityNumerator;
   delete fQualityDenominator;
+  delete fTPCSepNumerator;
+  delete fTPCSepDenominator;
 }
 //_________________________
 AliFemtoShareQualityCorrFctn& AliFemtoShareQualityCorrFctn::operator=(const AliFemtoShareQualityCorrFctn& aCorrFctn)
@@ -106,6 +128,14 @@ AliFemtoShareQualityCorrFctn& AliFemtoShareQualityCorrFctn::operator=(const AliF
     fQualityDenominator = new TH2D(*aCorrFctn.fQualityDenominator);
   else
     fQualityDenominator = 0;
+  if (aCorrFctn.fTPCSepNumerator)
+    fTPCSepNumerator = new TH2D(*aCorrFctn.fTPCSepNumerator);
+  else
+    fTPCSepNumerator = 0;
+  if (aCorrFctn.fTPCSepDenominator)
+    fTPCSepDenominator = new TH2D(*aCorrFctn.fTPCSepDenominator);
+  else
+    fTPCSepDenominator = 0;
 
   return *this;
 }
@@ -126,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;
@@ -137,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;
@@ -197,63 +230,73 @@ void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
     hsfval = ns*1.0/nh;
   }
 
-  if ((tQinv < 0.005) && (hsmval<-0.0)) {
-    cout << "Quality  Sharity " << hsmval << " " << hsfval << " " << pair->Track1()->Track() << " " << pair->Track2()->Track() << endl;
-    cout << "Qinv of the pair is " << tQinv << endl;
-    cout << "Clusters: " << endl;
-    for (unsigned int imap=0; imap<pair->Track1()->Track()->TPCclusters().GetNbits(); imap++) {
-      cout << imap ;
-      if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
-      else cout << " 0 " ;
-      if (pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
-      else cout << " 0 " ;
-      cout << "     ";
-      if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
-      else cout << " X ";
-      if (pair->Track2()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
-      else cout << " X ";
-      cout << endl;
-    }
-    cout << "Momentum1 " 
-        << pair->Track1()->Track()->P().x() << " " 
-        << pair->Track1()->Track()->P().y() << " "  
-        << pair->Track1()->Track()->P().z() << " "  
-        << pair->Track1()->Track()->Label() << " "  
-        << pair->Track1()->Track()->TrackId() << " "  
-        << pair->Track1()->Track()->Flags() << " "
-        << pair->Track1()->Track()->KinkIndex(0) << " "
-        << pair->Track1()->Track()->KinkIndex(1) << " "
-        << pair->Track1()->Track()->KinkIndex(2) << " "
-        << pair->Track1()->Track()->ITSchi2() << " "
-        << pair->Track1()->Track()->ITSncls() << " "
-        << pair->Track1()->Track()->TPCchi2() << " "
-        << pair->Track1()->Track()->TPCncls() << " "
-        << endl;
-    cout << "Momentum2 " 
-        << pair->Track2()->Track()->P().x() << " "  
-        << pair->Track2()->Track()->P().y() << " "  
-        << pair->Track2()->Track()->P().z() << " "  
-        << pair->Track2()->Track()->Label() << " "  
-        << pair->Track2()->Track()->TrackId() << " "  
-        << pair->Track2()->Track()->Flags() << " " 
-        << pair->Track2()->Track()->KinkIndex(0) << " "
-        << pair->Track2()->Track()->KinkIndex(1) << " "
-        << pair->Track2()->Track()->KinkIndex(2) << " "
-        << pair->Track2()->Track()->ITSchi2() << " "
-        << pair->Track2()->Track()->ITSncls() << " "
-        << pair->Track2()->Track()->TPCchi2() << " "
-        << pair->Track2()->Track()->TPCncls() << " "
-        << endl;
-  }
+//   if ((tQinv < 0.005) && (hsmval<-0.0)) {
+//     cout << "Quality  Sharity " << hsmval << " " << hsfval << " " << pair->Track1()->Track() << " " << pair->Track2()->Track() << endl;
+//     cout << "Qinv of the pair is " << tQinv << endl;
+//     cout << "Clusters: " << endl;
+//     for (unsigned int imap=0; imap<pair->Track1()->Track()->TPCclusters().GetNbits(); imap++) {
+//       cout << imap ;
+//       if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
+//       else cout << " 0 " ;
+//       if (pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) cout << " 1 ";
+//       else cout << " 0 " ;
+//       cout << "     ";
+//       if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
+//       else cout << " X ";
+//       if (pair->Track2()->Track()->TPCsharing().TestBitNumber(imap)) cout << " S ";
+//       else cout << " X ";
+//       cout << endl;
+//     }
+//     cout << "Momentum1 " 
+//      << pair->Track1()->Track()->P().x() << " " 
+//      << pair->Track1()->Track()->P().y() << " "  
+//      << pair->Track1()->Track()->P().z() << " "  
+//      << pair->Track1()->Track()->Label() << " "  
+//      << pair->Track1()->Track()->TrackId() << " "  
+//      << pair->Track1()->Track()->Flags() << " "
+//      << pair->Track1()->Track()->KinkIndex(0) << " "
+//      << pair->Track1()->Track()->KinkIndex(1) << " "
+//      << pair->Track1()->Track()->KinkIndex(2) << " "
+//      << pair->Track1()->Track()->ITSchi2() << " "
+//      << pair->Track1()->Track()->ITSncls() << " "
+//      << pair->Track1()->Track()->TPCchi2() << " "
+//      << pair->Track1()->Track()->TPCncls() << " "
+//      << endl;
+//     cout << "Momentum2 " 
+//      << pair->Track2()->Track()->P().x() << " "  
+//      << pair->Track2()->Track()->P().y() << " "  
+//      << pair->Track2()->Track()->P().z() << " "  
+//      << pair->Track2()->Track()->Label() << " "  
+//      << pair->Track2()->Track()->TrackId() << " "  
+//      << pair->Track2()->Track()->Flags() << " " 
+//      << pair->Track2()->Track()->KinkIndex(0) << " "
+//      << pair->Track2()->Track()->KinkIndex(1) << " "
+//      << pair->Track2()->Track()->KinkIndex(2) << " "
+//      << pair->Track2()->Track()->ITSchi2() << " "
+//      << pair->Track2()->Track()->ITSncls() << " "
+//      << pair->Track2()->Track()->TPCchi2() << " "
+//      << pair->Track2()->Track()->TPCncls() << " "
+//      << endl;
+//   }
+
+  double distx = pair->Track1()->Track()->NominalTpcEntrancePoint().x() - pair->Track2()->Track()->NominalTpcEntrancePoint().x();
+  double disty = pair->Track1()->Track()->NominalTpcEntrancePoint().y() - pair->Track2()->Track()->NominalTpcEntrancePoint().y();
+  double distz = pair->Track1()->Track()->NominalTpcEntrancePoint().z() - pair->Track2()->Track()->NominalTpcEntrancePoint().z();
+  double dist = sqrt(distx*distx + disty*disty + distz*distz);
 
   fShareNumerator->Fill(tQinv, hsfval);
   fQualityNumerator->Fill(tQinv, hsmval);
+  fTPCSepNumerator->Fill(tQinv, dist);
+
   //  cout << "AliFemtoShareQualityCorrFctn::AddRealPair : " << pair->qInv() << " " << tQinv <<
   //" " << pair->Track1().FourMomentum() << " " << pair->Track2().FourMomentum() << endl;
 }
 //____________________________
 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;
@@ -299,8 +342,15 @@ void AliFemtoShareQualityCorrFctn::AddMixedPair( AliFemtoPair* pair){
     hsfval = ns*1.0/nh;
   }
 
+  double distx = pair->Track1()->Track()->NominalTpcEntrancePoint().x() - pair->Track2()->Track()->NominalTpcEntrancePoint().x();
+  double disty = pair->Track1()->Track()->NominalTpcEntrancePoint().y() - pair->Track2()->Track()->NominalTpcEntrancePoint().y();
+  double distz = pair->Track1()->Track()->NominalTpcEntrancePoint().z() - pair->Track2()->Track()->NominalTpcEntrancePoint().z();
+  double dist = sqrt(distx*distx + disty*disty + distz*distz);
+
   fShareDenominator->Fill(tQinv,hsfval,weight);
   fQualityDenominator->Fill(tQinv,hsmval,weight);
+  fTPCSepDenominator->Fill(tQinv, dist, weight);
+
 }
 
 
@@ -311,6 +361,8 @@ void AliFemtoShareQualityCorrFctn::WriteHistos()
   fShareDenominator->Write();
   fQualityNumerator->Write();
   fQualityDenominator->Write();
+  fTPCSepNumerator->Write();
+  fTPCSepDenominator->Write();
   
 }
 //______________________________
@@ -323,6 +375,8 @@ TList* AliFemtoShareQualityCorrFctn::GetOutputList()
   tOutputList->Add(fShareDenominator);  
   tOutputList->Add(fQualityNumerator);  
   tOutputList->Add(fQualityDenominator);  
+  tOutputList->Add(fTPCSepNumerator);  
+  tOutputList->Add(fTPCSepDenominator);  
 
   return tOutputList;
 }