Suppressed warning from fkSegmentation->PadByPosition(...)

[u/mrichter/AliRoot.git] / MUON / AliMUONClusterReconstructor.cxx

diff --git a/MUON/AliMUONClusterReconstructor.cxx b/MUON/AliMUONClusterReconstructor.cxx
index ef444b0f16630a8d6da2819c348b22ec416c558f..c659f6745eac79783a32c51c14ac19f4ac9e1ebe 100644 (file)
--- a/MUON/AliMUONClusterReconstructor.cxx
+++ b/MUON/AliMUONClusterReconstructor.cxx
@@ -34,20 +34,13 @@
 #include "AliMUONClusterFinderVS.h"
 #include "AliMUONClusterInput.h"
 #include "AliMUONRawCluster.h"
-#include "AliRawReader.h" // for raw data
 #include "AliLog.h"
 
-
-const Int_t AliMUONClusterReconstructor::fgkDefaultPrintLevel = 0;
-
 ClassImp(AliMUONClusterReconstructor) // Class implementation in ROOT context
 
 //__________________________________________________________________________
-AliMUONClusterReconstructor::AliMUONClusterReconstructor(AliLoader* loader)
-  : TObject(),
-    fMUONData(0),
-    fPrintLevel(fgkDefaultPrintLevel),
-    fDebug(0)
+  AliMUONClusterReconstructor::AliMUONClusterReconstructor(AliLoader* loader, AliMUONData* data)
+  : TObject()
 {
   // Standard Constructor
 
@@ -55,8 +48,8 @@ AliMUONClusterReconstructor::AliMUONClusterReconstructor(AliLoader* loader)
   fLoader = loader;
 
   // initialize container
-  fMUONData  = new AliMUONData(fLoader,"MUON","MUON");
-
+  fMUONData = data;
+  
   // reconstruction model
   fRecModel = new AliMUONClusterFinderVS();
   //fRecModel = new AliMUONClusterFinderAZ();
@@ -67,8 +60,7 @@ AliMUONClusterReconstructor::AliMUONClusterReconstructor(AliLoader* loader)
 AliMUONClusterReconstructor::AliMUONClusterReconstructor()
   : TObject(),
     fMUONData(0),
-    fPrintLevel(fgkDefaultPrintLevel),
-    fDebug(0),
+    fRecModel(0),
     fLoader(0)
 {
   // Default Constructor
@@ -100,17 +92,17 @@ AliMUONClusterReconstructor::operator=(const AliMUONClusterReconstructor& rhs)
 AliMUONClusterReconstructor::~AliMUONClusterReconstructor(void)
 {
 
-  if (fMUONData)
-    delete fMUONData;
+  if (fRecModel)
+    delete fRecModel;
 
   return;
 }
 //____________________________________________________________________
-void AliMUONClusterReconstructor::Digits2Clusters()
+void AliMUONClusterReconstructor::Digits2Clusters(Int_t chBeg)
 {
 
     TClonesArray *dig1, *dig2, *digAll;
-    Int_t ndig, k, idDE, idDE_prev;
+    Int_t ndig, k, idDE, idDEprev;
     dig1 = new TClonesArray("AliMUONDigit",1000);
     dig2 = new TClonesArray("AliMUONDigit",1000);
     digAll = new TClonesArray("AliMUONDigit",2000);
@@ -124,60 +116,32 @@ void AliMUONClusterReconstructor::Digits2Clusters()
     TClonesArray* muonDigits;
     Int_t n2;
     Int_t n1;
-    Int_t flag = 0;
 
-    for (Int_t ich = 0; ich < AliMUONConstants::NTrackingCh(); ich++) {
+    for (Int_t ich = chBeg; ich < AliMUONConstants::NTrackingCh(); ich++) {
  
       id.Reset();
       n1 = 0;
       n2 = 0;
-      //cathode 0
-      fMUONData->ResetDigits();
-      fMUONData->GetCathode(0);
+
+      //cathode 0 & 1
       muonDigits = fMUONData->Digits(ich); 
       ndig = muonDigits->GetEntriesFast();
       TClonesArray &lDigit = *digAll;
 
-      idDE_prev = 0;
-
+      idDEprev = 0;
+      muonDigits->Sort();
       for (k = 0; k < ndig; k++) {
 
        digit = (AliMUONDigit*) muonDigits->UncheckedAt(k);
        new(lDigit[n1++]) AliMUONDigit(*digit);
        idDE = digit->DetElemId();
-       if (idDE != idDE_prev) {
+       if (idDE != idDEprev) {
          id.AddAt(idDE,n2++);
        }
-       idDE_prev = idDE;
+       idDEprev = idDE;
       }
 
-      //cathode 1
-      fMUONData->ResetDigits();
-      fMUONData->GetCathode(1);
-      muonDigits =  fMUONData->Digits(ich);  
-      ndig = muonDigits->GetEntriesFast();
-
       Int_t idSize = n2;
-    
-      for (k = 0; k < ndig; k++) {
-
-       digit = (AliMUONDigit*) muonDigits->UncheckedAt(k);
-       new(lDigit[n1++]) AliMUONDigit(*digit);
-       idDE = digit->DetElemId();
-       flag = 0;
-
-       // looking for new idDE in cathode 1 (method to be checked CF)
-       for (Int_t n = 0; n < idSize; n++) {
-         if (idDE == id[n]) {
-           flag = 1;
-           break;
-         }
-       }
-       if (flag) continue;
-       id.AddAt(idDE,n2++);
-      }
-
-      idSize = n2;
 
       // loop over id DE
       for (idDE = 0; idDE < idSize; idDE++) {
@@ -198,8 +162,6 @@ void AliMUONClusterReconstructor::Digits2Clusters()
          }
        }
 
-       //      if (id[idDE] < 500 && id[idDE] > 299) continue; // temporary patch til St2 geometry is not yet ok (CF)
-
        // cluster finder
        if (fRecModel) {
          AliMUONClusterInput::Instance()->SetDigits(ich, id[idDE], dig1, dig2);
@@ -221,14 +183,6 @@ void AliMUONClusterReconstructor::Digits2Clusters()
     delete digAll;
 }
 
-//____________________________________________________________________
-void AliMUONClusterReconstructor::Digits2Clusters(AliRawReader* /*rawReader*/)
-{
-
-//  Perform cluster finding form raw data
-
-   AliFatal("clusterization not implemented for raw data input");
-}
 //_______________________________________________________________________
 void AliMUONClusterReconstructor::Trigger2Trigger() 
 {
@@ -237,11 +191,3 @@ void AliMUONClusterReconstructor::Trigger2Trigger()
   fMUONData->SetTreeAddress("GLT");
   fMUONData->GetTriggerD();
 }
-//_______________________________________________________________________
-void AliMUONClusterReconstructor::Trigger2Trigger(AliRawReader* /*rawReader*/) 
-{
-// call the Trigger Algorithm from raw data and fill TreeR 
-
-   AliFatal("Trigger not implemented for raw data input");
-
-}