track matching macros from Alberto

[u/mrichter/AliRoot.git] / ANALYSIS / AliClusterMap.cxx

diff --git a/ANALYSIS/AliClusterMap.cxx b/ANALYSIS/AliClusterMap.cxx
index 94e501862a06754ad83005440ce5c30f15eccb46..711c1ba75a4275778ae134711b4667cfbd1bce09 100644 (file)
--- a/ANALYSIS/AliClusterMap.cxx
+++ b/ANALYSIS/AliClusterMap.cxx
@@ -1,56 +1,64 @@
-#include "AliClusterMap.h"
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/* $Id$ */
 
-//_________________________________________________
-///////////////////////////////////////////////////
+///////////////////////////////////////////////////////////////////////////
 //
 // class AliClusterMap
-//
 // class that describes cluster occupation at TPC
 // Each padraw has a corresponding bit in fPadRawMap
-// 
-//
 // more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
 // Piotr.Skowronski@cern.ch
 //
 ///////////////////////////////////////////////////////////////////////////
 
+#include <TString.h>
 
+#include "AliClusterMap.h"
 #include "AliESDtrack.h"
+#include "AliLog.h"
 #include "AliTPCtrack.h"
-#include "AliVAODParticle.h"
-#include <TString.h>
-const Int_t AliClusterMap::fNPadRows = 159;
+
+const Int_t AliClusterMap::fgkNPadRows = 159;
+
+ClassImp(AliClusterMap)
 
 AliClusterMap::AliClusterMap():
- fPadRawMap(fNPadRows)
+ fPadRawMap(fgkNPadRows)
 {
 //ctor
 
 }
 /***********************************************************************/
 AliClusterMap::AliClusterMap(AliESDtrack* track):
- fPadRawMap( (track)?track->GetTPCClusterMap():fNPadRows )
+ fPadRawMap( (track)?track->GetTPCClusterMap():fgkNPadRows )
 {
  //ctor
  
- if (AliVAODParticle::GetDebug() > 2)
-  { 
-    Info("AliClusterMap(AliESDtrack*)","");
-    Print();
-  } 
+ StdoutToAliDebug(3,Print());
+
 } 
 /***********************************************************************/
 
 AliClusterMap::AliClusterMap(AliTPCtrack* track):
- fPadRawMap(fNPadRows)
+ fPadRawMap(fgkNPadRows)
 {
  //ctor
  
- //Does not work since indeces in the claster index array 
- //in the TPC track does not correspond to the padraw segmatation
- if (AliVAODParticle::GetDebug() > 9) 
-   Info("AliClusterMap",
-      "#####################################################################"); 
+ AliDebug(10,"#####################################################################"); 
  if (track == 0x0)
   {
     Error("AliClusterMap","Pointer to TPC track is NULL");
@@ -64,8 +72,7 @@ AliClusterMap::AliClusterMap(AliTPCtrack* track):
     Int_t sect = (idx&0xff000000)>>24;
     Int_t row = (idx&0x00ff0000)>>16;
     if (sect > 18) row +=63; //if it is outer sector, add number of inner sectors
-    if (AliVAODParticle::GetDebug() > 9)  
-      Info("AliClusterMap","Cl.idx is %d, sect %d, row %d",idx,sect,row);
+    AliDebug(9,Form("Cl.idx is %d, sect %d, row %d",idx,sect,row));
       
     fPadRawMap.SetBitNumber(row,kTRUE);
     
@@ -96,11 +103,8 @@ AliClusterMap::AliClusterMap(AliTPCtrack* track):
      }
   }
   
- if (AliVAODParticle::GetDebug() > 2)
-  { 
-    Info("AliClusterMap(AliTPCtrack*)","");
-    Print();
-  } 
+ StdoutToAliDebug(3,Print());
+
 }
 /***********************************************************************/
 
@@ -108,7 +112,7 @@ void AliClusterMap::Print() const
 {
 //Prints the bit map 
   TString msg;
-  for ( Int_t i = 0; i < fNPadRows; i++)
+  for ( Int_t i = 0; i < fgkNPadRows; i++)
    {
      if ( fPadRawMap.TestBitNumber(i) )
       {
@@ -135,19 +139,14 @@ Float_t AliClusterMap::GetOverlapFactor(const AliClusterMap& clmap) const
   //  -0.5 (low probability that these tracks are a split track)
   //  and
   //   1.0 (high probability that these tracks are a split track)
-  TString msg1;
-  TString msg2;
   
   Int_t nh = 0;
   Int_t an = 0;
-  for ( Int_t i = 0; i < fNPadRows; i++)
+  for ( Int_t i = 0; i < fgkNPadRows; i++)
    {
      Bool_t x = HasClAtPadRow(i);
      Bool_t y = clmap.HasClAtPadRow(i);
      
-     if (x) msg1+="1";else msg1+="0";
-     if (y) msg2+="1";else msg2+="0";
-     
      if (x && y)//both have clasters
       {
        an--;
@@ -169,24 +168,5 @@ Float_t AliClusterMap::GetOverlapFactor(const AliClusterMap& clmap) const
   if (nh > 0) retval = ((Float_t)an)/((Float_t)nh);
   else Warning("GetOverlapFactor","Number of counted cluters is 0.");
   
-  if (AliVAODParticle::GetDebug() > 2)
-   {
-     Info("GetOverlapFactor","Splitting Quality Factor is %f. SumAn = %d, SumClusters %d",retval,an,nh); 
-     if (retval == 1.0) 
-      { 
-        Print();
-        Info("AliClusterMap","BitMap is\n  %s\n",msg1.Data());
-        clmap.Print(); 
-        Info("AliClusterMap","BitMap is\n  %s\n\n\n\n",msg2.Data());
-      }
-     if (retval == -.5) 
-      { 
-        Print();
-        Info("AliClusterMap","BitMap is\n  %s\n",msg1.Data());
-        clmap.Print(); 
-        Info("AliClusterMap","BitMap is\n  %s\n\n\n\n",msg2.Data());
-      }
-   } 
- 
   return retval;
 }