AliMUONClusterFinderPeakFit.h

Go to the documentation of this file.

#ifndef ALIMUONCLUSTERFINDERPEAKFIT_H
#define ALIMUONCLUSTERFINDERPEAKFIT_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

/* $Id$ */

//  Author Alexander Zinchenko, JINR Dubna; Laurent Aphecetche, SUBATECH
//

#include "AliMUONVClusterFinder.h"

#ifndef ROOT_TObjArray
#  include "TObjArray.h"
#endif
#ifndef ROOT_TVector2
#  include "TVector2.h"
#endif

class AliMUONMathieson;
//class AliMUONPad;

class TH2D;

class AliMUONClusterFinderPeakFit : public AliMUONVClusterFinder
{
public:
  AliMUONClusterFinderPeakFit(Bool_t plot, AliMUONVClusterFinder* clusterFinder); // Constructor
  virtual ~AliMUONClusterFinderPeakFit(); // Destructor

  virtual Bool_t NeedSegmentation() const { return kTRUE; }

  using AliMUONVClusterFinder::Prepare;

  virtual Bool_t Prepare(Int_t detElemId, TObjArray* pads[2],
                         const AliMpArea& area, const AliMpVSegmentation* seg[2]);

  virtual AliMUONCluster* NextCluster();
  
  Int_t GetNMax() const { return fNMax; }
  
  virtual void Print(Option_t* opt="") const;

private:
  AliMUONClusterFinderPeakFit(const AliMUONClusterFinderPeakFit& rhs);
  AliMUONClusterFinderPeakFit& operator=(const AliMUONClusterFinderPeakFit& rhs);

  Bool_t WorkOnPreCluster();

  AliMUONCluster* CheckPrecluster(const AliMUONCluster& cluster); 
  AliMUONCluster* CheckPreclusterTwoCathodes(AliMUONCluster* cluster); 
  
  Bool_t Overlap(const AliMUONPad& pad, const AliMUONPad& pixel); 
  
  void BuildPixArray(AliMUONCluster& cluster); 
  void BuildPixArrayOneCathode(AliMUONCluster& cluster); 
  void PadOverHist(Int_t idir, Int_t ix0, Int_t iy0, AliMUONPad *pad, TH2D *h1, TH2D *h2);

  void RemovePixel(Int_t i);
  
  AliMUONPad* Pixel(Int_t i) const;
  
  Int_t FindLocalMaxima(TObjArray *pixArray, Int_t *localMax, Double_t *maxVal); // find local maxima 
  void  FlagLocalMax(TH2D *hist, Int_t i, Int_t j, Int_t *isLocalMax); // flag local max
  void  FindClusterCOG(AliMUONCluster& cluster, const Int_t *localMax, Int_t iMax); // find cluster around local max with COG
  void  FindClusterFit(AliMUONCluster& cluster, const Int_t *localMax, const Int_t *maxPos, Int_t nMax); // find cluster around local max with fitting
  void  PadsInXandY(AliMUONCluster& cluster, Int_t &nInX, Int_t &nInY) const; // get number of pads in X and Y

  void CheckOverlaps();

private:
  // Status flags for pads
  static const Int_t fgkZero; 
  static const Int_t fgkMustKeep; 
  static const Int_t fgkUseForFit; 
  static const Int_t fgkOver; 
  static const Int_t fgkModified; 
  static const Int_t fgkCoupled; 
    
  // Some constants
  static const Double_t fgkZeroSuppression; 
  static const Double_t fgkDistancePrecision; 
  static const TVector2 fgkIncreaseSize; 
  static const TVector2 fgkDecreaseSize; 
  
  AliMUONVClusterFinder* fPreClusterFinder; 
  AliMUONCluster* fPreCluster; 
  TObjArray fClusterList; 
  AliMUONMathieson* fMathieson; 
  
  Int_t fEventNumber; 
  Int_t fDetElemId; 
  Int_t fClusterNumber; 
  Int_t fNMax; 
  TH2D *fHistAnode; 
  
  const AliMpVSegmentation *fkSegmentation[2]; 
  
  TObjArray* fPixArray; 
  Int_t fDebug; 
  Bool_t fPlot; 
  
  Int_t fNClusters; 
  Int_t fNAddVirtualPads; 
  
  ClassDef(AliMUONClusterFinderPeakFit,0) // cluster finder in MUON arm of ALICE
};

#endif