Port of changes from v4-07-Release and additional rule conformance

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

diff --git a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
index c9dc4f730778300a8c18e74eb8ad9149fe8438a7..0effcc86d4e12bbf4c4d65903c3700da0588588b 100644 (file)
--- a/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
+++ b/PWG2/FEMTOSCOPY/AliFemtoUser/AliFemtoShareQualityCorrFctn.cxx
@@ -25,21 +25,21 @@ AliFemtoShareQualityCorrFctn::AliFemtoShareQualityCorrFctn(char* title, const in
   //  title = "Num Qinv (MeV/c)";
   char tTitNum[100] = "NumShare";
   strcat(tTitNum,title);
-  fShareNumerator = new TH2D(tTitNum,title,nbins,QinvLo,QinvHi,50,0.0,1.00001);
+  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);
-  fShareDenominator = new TH2D(tTitDen,title,nbins,QinvLo,QinvHi,50,0.0,1.00001);
+  fShareDenominator = new TH2D(tTitDen,title,nbins,QinvLo,QinvHi,100,0.0,1.00001);
 
   char tTit2Num[100] = "NumQuality";
   strcat(tTit2Num,title);
-  fQualityNumerator = new TH2D(tTit2Num,title,nbins,QinvLo,QinvHi,75,-0.500001,1.000001);
+  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);
-  fQualityDenominator = new TH2D(tTit2Den,title,nbins,QinvLo,QinvHi,75,-0.500001,1.000001);
+  fQualityDenominator = new TH2D(tTit2Den,title,nbins,QinvLo,QinvHi,150,-0.500001,1.000001);
   // set up ratio
   //title = "Ratio Qinv (MeV/c)";
   // this next bit is unfortunately needed so that we can have many histos of same "title"
@@ -135,7 +135,7 @@ AliFemtoString AliFemtoShareQualityCorrFctn::Report(){
 //____________________________
 void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
   // add real (effect) pair
-  double tQinv = fabs(pair->QInv());   // note - qInv() will be negative for identical pairs...
+  //  double tQinv = fabs(pair->QInv());   // note - qInv() will be negative for identical pairs...
   Int_t nh = 0;
   Int_t an = 0;
   Int_t ns = 0;
@@ -148,9 +148,9 @@ void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
       if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap) &&
          pair->Track2()->Track()->TPCsharing().TestBitNumber(imap))
        {
-         if (tQinv < 0.01) {
-           cout << "Shared cluster in row " << imap << endl; 
-         }
+//       if (tQinv < 0.01) {
+//         cout << "Shared cluster in row " << imap << endl; 
+//       }
          an++;
          nh+=2;
          ns+=2;
@@ -169,23 +169,23 @@ void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
       nh++;
     }
   }
-  if (tQinv < 0.01) {
-    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;
-    }
-  }
+//    if (tQinv < 0.01) {
+//     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;
+//     }
+//   }
 
   Float_t hsmval = 0.0;
   Float_t hsfval = 0.0;
@@ -195,12 +195,49 @@ void AliFemtoShareQualityCorrFctn::AddRealPair( AliFemtoPair* pair){
     hsfval = ns*1.0/nh;
   }
 
-  if (tQinv < 0.01) {
-    cout << "Quality  Sharity " << hsmval << " " << hsfval << " " << pair->Track1()->Track() << " " << pair->Track2()->Track() << endl;
-  }
+//   if ((tQinv < 0.01) && (hsmval<-0.46)) {
+//     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) << " "
+//      << 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) << " "
+//      << endl;
+//   }
 
-  fShareNumerator->Fill(tQinv, hsfval);
-  fQualityNumerator->Fill(tQinv, hsmval);
+//   fShareNumerator->Fill(tQinv, hsfval);
+//   fQualityNumerator->Fill(tQinv, hsmval);
   //  cout << "AliFemtoShareQualityCorrFctn::AddRealPair : " << pair->qInv() << " " << tQinv <<
   //" " << pair->Track1().FourMomentum() << " " << pair->Track2().FourMomentum() << endl;
 }
@@ -218,7 +255,7 @@ void AliFemtoShareQualityCorrFctn::AddMixedPair( AliFemtoPair* pair){
     if (pair->Track1()->Track()->TPCclusters().TestBitNumber(imap) && 
        pair->Track2()->Track()->TPCclusters().TestBitNumber(imap)) {
       // Do they share it ?
-      if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap) ||
+      if (pair->Track1()->Track()->TPCsharing().TestBitNumber(imap) &&
          pair->Track2()->Track()->TPCsharing().TestBitNumber(imap))
        {
          //      cout << "A shared cluster !!!" << endl;