Remove obselete method Hit2Digit for digitalization

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

/**************************************************************************
 * 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$ */

////////////////////////////////////
//
// MUON cluster reconstructor for MUON
//
// Should implement a virtual class ClusterFinder to chose between VS and AZ method
////////////////////////////////////

#include "AliMUONClusterReconstructor.h"
#include "AliRun.h" // for gAlice
#include "AliRunLoader.h"
#include "AliLoader.h"

#include "AliMUON.h"
#include "AliMUONDigit.h"
#include "AliMUONConstants.h"
#include "AliMUONData.h"
#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)
{
  // Standard Constructor

  // initialize loader's
  fLoader = loader;

  // initialize container
  fMUONData  = new AliMUONData(fLoader,"MUON","MUON");

  // reconstruction model
  fRecModel = new AliMUONClusterFinderVS();
  //fRecModel = new AliMUONClusterFinderAZ();

}

//__________________________________________________________________________
AliMUONClusterReconstructor::AliMUONClusterReconstructor()
  : TObject(),
    fMUONData(0),
    fPrintLevel(fgkDefaultPrintLevel),
    fDebug(0),
    fLoader(0)
{
  // Default Constructor
}

//_______________________________________________________________________
AliMUONClusterReconstructor::AliMUONClusterReconstructor (const AliMUONClusterReconstructor& rhs)
  : TObject(rhs)
{
// Protected copy constructor

  AliFatal("Not implemented.");
}

//_______________________________________________________________________
AliMUONClusterReconstructor & 
AliMUONClusterReconstructor::operator=(const AliMUONClusterReconstructor& rhs)
{
// Protected assignement operator

  if (this == &rhs) return *this;

  AliFatal("Not implemented.");
    
  return *this;  
}

//__________________________________________________________________________
AliMUONClusterReconstructor::~AliMUONClusterReconstructor(void)
{

  if (fMUONData)
    delete fMUONData;

  return;
}
//____________________________________________________________________
void AliMUONClusterReconstructor::Digits2Clusters()
{
//
//  Perform cluster finding
//

    AliMUON* pMUON = (AliMUON*) gAlice->GetModule("MUON");
    if (pMUON->WhichSegmentation() == 1)
      Digits2ClustersOld();
    else
      Digits2ClustersNew();

}
//____________________________________________________________________
void AliMUONClusterReconstructor::Digits2ClustersOld()
{

//
//  Perform cluster finding
//
    TClonesArray *dig1, *dig2;
    Int_t ndig, k;
    dig1 = new TClonesArray("AliMUONDigit",1000);
    dig2 = new TClonesArray("AliMUONDigit",1000);
    AliMUONDigit *digit;

// Loop on chambers and on cathode planes     
    TClonesArray * muonDigits;

    for (Int_t ich = 0; ich < 10; ich++) {

        fMUONData->ResetDigits();
        fMUONData->GetCathode(0);
        //TClonesArray *
        muonDigits = fMUONData->Digits(ich); 
        ndig=muonDigits->GetEntriesFast();
        AliDebug(1,Form("1 Found %d digits in %p chamber %d", ndig, (void*)muonDigits,ich));
        TClonesArray &lhits1 = *dig1;
        Int_t n = 0;
        for (k = 0; k < ndig; k++) {
            digit = (AliMUONDigit*) muonDigits->UncheckedAt(k);
	    new(lhits1[n++]) AliMUONDigit(*digit);
        }
        fMUONData->ResetDigits();
        fMUONData->GetCathode(1);
        muonDigits =  fMUONData->Digits(ich);  
        ndig=muonDigits->GetEntriesFast();
        AliDebug(1,Form("2 Found %d digits in %p %d", ndig, (void*)muonDigits, ich));
        TClonesArray &lhits2 = *dig2;
        n=0;
        
        for (k=0; k<ndig; k++) {
            digit= (AliMUONDigit*) muonDigits->UncheckedAt(k);
	    new(lhits2[n++]) AliMUONDigit(*digit);
        }

        if (fRecModel) {         
            AliMUONClusterInput::Instance()->SetDigits(ich, dig1, dig2);
            fRecModel->FindRawClusters();
        }
        // copy into the container
        TClonesArray* tmp = fRecModel->GetRawClusters();
        for (Int_t id = 0; id < tmp->GetEntriesFast(); id++) {
          AliMUONRawCluster* pClus = (AliMUONRawCluster*) tmp->At(id);
          fMUONData->AddRawCluster(ich, *pClus);
        }
        dig1->Delete();
        dig2->Delete();
    } // for ich
    delete dig1;
    delete dig2;
}
//____________________________________________________________________
void AliMUONClusterReconstructor::Digits2ClustersNew()
{

    TClonesArray *dig1, *dig2, *digAll;
    Int_t ndig, k, idDE, idDE_prev;
    dig1 = new TClonesArray("AliMUONDigit",1000);
    dig2 = new TClonesArray("AliMUONDigit",1000);
    digAll = new TClonesArray("AliMUONDigit",2000);

    AliMUONDigit* digit;

    TArrayI id(100); // contains the different IdDE
    id.Reset();
  
// Loop on chambers and on cathode planes     
    TClonesArray* muonDigits;
    Int_t n2;
    Int_t n1;
    Int_t flag = 0;

    for (Int_t ich = 0; ich < AliMUONConstants::NTrackingCh(); ich++) {
      n1 = 0;
      n2 = 0;
      //cathode 0
      fMUONData->ResetDigits();
      fMUONData->GetCathode(0);
      muonDigits = fMUONData->Digits(ich); 
      ndig = muonDigits->GetEntriesFast();
      TClonesArray &lDigit = *digAll;

      idDE_prev = 0;

      for (k = 0; k < ndig; k++) {

	digit = (AliMUONDigit*) muonDigits->UncheckedAt(k);
	new(lDigit[n1++]) AliMUONDigit(*digit);
	idDE = digit->DetElemId();
	if (idDE != idDE_prev)
	  id.AddAt(idDE,n2++);
	idDE_prev = idDE;
      }

      //cathode 1
      fMUONData->ResetDigits();
      fMUONData->GetCathode(1);
      muonDigits =  fMUONData->Digits(ich);  
      ndig = muonDigits->GetEntriesFast();

      for (k = 0; k < ndig; k++) {

	digit = (AliMUONDigit*) muonDigits->UncheckedAt(k);
	new(lDigit[n1++]) AliMUONDigit(*digit);
	idDE = digit->DetElemId();
	flag = 0;

	for (Int_t n = 0; n < id.GetSize(); n++) {
	  if (id[n] == idDE) {
	    flag = 1;
	    break;
	  }
	  if (flag) continue;
	  id[id.GetSize()+1] = idDE;
	}
      }

      // loop over id DE
      for (idDE = 0; idDE < id.GetSize(); idDE++) {

	TClonesArray &lhits1 = *dig1;
	TClonesArray &lhits2 = *dig2;
	n1 = n2 = 0;

	for (k = 0; k < digAll->GetEntriesFast(); k++) {
	  digit = (AliMUONDigit*) digAll->UncheckedAt(k);
	  if (id[idDE] == digit->DetElemId())
	    if (digit->Cathode() == 1)
	      new(lhits1[n1++]) AliMUONDigit(*digit);
	    else 
	      new(lhits2[n2++]) AliMUONDigit(*digit);

	}

	// cluster finder
	if (fRecModel) {
	  AliMUONClusterInput::Instance()->SetDigits(ich, id[idDE], dig1, dig2);
	  fRecModel->FindRawClusters();
	}
	// copy into the container
	TClonesArray* tmp = fRecModel->GetRawClusters();
	for (Int_t id = 0; id < tmp->GetEntriesFast(); id++) {
	  AliMUONRawCluster* pClus = (AliMUONRawCluster*) tmp->At(id);
	  fMUONData->AddRawCluster(ich, *pClus);
	}
	dig1->Delete();
	dig2->Delete();
	digAll->Delete();

      }
    } // for ich
    delete dig1;
    delete dig2;
    delete digAll;
}

//____________________________________________________________________
void AliMUONClusterReconstructor::Digits2Clusters(AliRawReader* /*rawReader*/)
{

//  Perform cluster finding form raw data

   AliFatal("clusterization not implemented for raw data input");
}
//_______________________________________________________________________
void AliMUONClusterReconstructor::Trigger2Trigger() 
{
// copy trigger from TreeD to TreeR

  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");

}