[u/mrichter/AliRoot.git] / PHOS / AliPHOSClusterizerv1.h

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

/* $Id$ */

//_________________________________________________________________________
//  Implementation version 1 of the clusterization algorithm                     
//
//*-- Author: Yves Schutz (SUBATECH)

// --- ROOT system ---

// --- Standard library ---

// --- AliRoot header files ---

#include "AliPHOSClusterizer.h"
#include "AliPHOSDigit.h" 


class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
  
public:
  
  AliPHOSClusterizerv1() ;             // ctor            
  virtual ~AliPHOSClusterizerv1(){} ;  // dtor
  
  Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2) ; // Checks if digits are in neighbour cells 
  Float_t Calibrate(Int_t Amp){ return (fA + fB * Amp) ;}     // Tranforms Amp to energy 
  void FillandSort(const DigitsList * dl, TObjArray * tl) ;   // Sorts the list according to increasing id
  virtual void GetNumberOfClustersFound(Int_t * numb) ; 
 
  virtual void GetCalibrationParameters(Float_t & A, Float_t &B) { A = fA; B = fB; } 
  virtual Float_t GetEmcClusteringThreshold() { return fEmcClusteringThreshold;}
  virtual Float_t GetEmcEnergyThreshold() { return fEmcEnergyThreshold; }  
  virtual Float_t GetLocalMaxCut() { return fLocMaxCut;} 
  virtual Float_t GetLogWeightCut() { return fW0;}  
  virtual Float_t GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold;  } 
  virtual Float_t GetPpsdEnergyThreshold() { return  fPpsdEnergyThreshold;  }

  virtual Bool_t IsInEmc(AliPHOSDigit * digit) ;                      // Tells if id digit is in EMC
  virtual void MakeClusters(const DigitsList * dl, 
			    AliPHOSRecPoint::RecPointsList * emcl, 
			    AliPHOSRecPoint::RecPointsList * ppsdl) ; // does the job 
  virtual void PrintParameters() ;  
  virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;} 
  virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
  virtual void SetEmcEnergyThreshold(Float_t enth) { fEmcEnergyThreshold = enth ; } 
  virtual void SetLocalMaxCut(Float_t cut) { fLocMaxCut = cut ; }
  virtual void SetLogWeightCut(Float_t w) { fW0 = w ; }
  virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; }
  virtual void SetPpsdEnergyThreshold(Float_t enth) { fPpsdEnergyThreshold = enth ; } 
  
private:
  
  Float_t fA ;                       // offset of the energy calibration
  Float_t fB ;                       // gain of the energy calibration
  Float_t fEmcClusteringThreshold ;  // minimum energy to include a EMC digit in a cluster
  Float_t fEmcEnergyThreshold ;      // minimum energy of EMC digit to be considered
  Float_t fLocMaxCut ;               // minimum energy difference to distinguish local maxima in a cluster
  Int_t   fNumberOfEmcClusters ;     // number of EMC clusters found 
  Int_t   fNumberOfPpsdClusters ;    // number of PPSD clusters found
  Float_t fPpsdClusteringThreshold ; // minimum energy to include a PPSD digit in a cluster
  Float_t fPpsdEnergyThreshold ;     // minimum energy of PPSD digit to be considered
  Float_t fW0 ;                      // logarithmic weight for the cluster center of gravity calculation
    
  ClassDef(AliPHOSClusterizerv1,1)  // Clusterizer implementation version 1

};

#endif // AliPHOSCLUSTERIZERV1_H
Commit	Line	Data
d15a28e7	1	#ifndef ALIPHOSCLUSTERIZERV1_H
	2	#define ALIPHOSCLUSTERIZERV1_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
6ad0bfa0	6	/* $Id$ */
6ad0bfa0	7
b2a60966	8	//_________________________________________________________________________
	9	// Implementation version 1 of the clusterization algorithm
	10	//
	11	//*-- Author: Yves Schutz (SUBATECH)
d15a28e7	12
	13	// --- ROOT system ---
	14
	15	// --- Standard library ---
	16
	17	// --- AliRoot header files ---
	18
	19	#include "AliPHOSClusterizer.h"
	20	#include "AliPHOSDigit.h"
	21
	22
	23
	24	class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
	25
	26	public:
	27
	28	AliPHOSClusterizerv1() ; // ctor
	29	virtual ~AliPHOSClusterizerv1(){} ; // dtor
	30
9f616d61	31	Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2) ; // Checks if digits are in neighbour cells
	32	Float_t Calibrate(Int_t Amp){ return (fA + fB * Amp) ;} // Tranforms Amp to energy
	33	void FillandSort(const DigitsList * dl, TObjArray * tl) ; // Sorts the list according to increasing id
4a646085	34	virtual void GetNumberOfClustersFound(Int_t * numb) ;
	35
	36	virtual void GetCalibrationParameters(Float_t & A, Float_t &B) { A = fA; B = fB; }
416692f4	37	virtual Float_t GetEmcClusteringThreshold() { return fEmcClusteringThreshold;}
	38	virtual Float_t GetEmcEnergyThreshold() { return fEmcEnergyThreshold; }
	39	virtual Float_t GetLocalMaxCut() { return fLocMaxCut;}
	40	virtual Float_t GetLogWeightCut() { return fW0;}
	41	virtual Float_t GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold; }
	42	virtual Float_t GetPpsdEnergyThreshold() { return fPpsdEnergyThreshold; }
4a646085	43
92862013	44	virtual Bool_t IsInEmc(AliPHOSDigit * digit) ; // Tells if id digit is in EMC
88714635	45	virtual void MakeClusters(const DigitsList * dl,
	46	AliPHOSRecPoint::RecPointsList * emcl,
	47	AliPHOSRecPoint::RecPointsList * ppsdl) ; // does the job
6ad0bfa0	48	virtual void PrintParameters() ;
	49	virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;}
	50	virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
	51	virtual void SetEmcEnergyThreshold(Float_t enth) { fEmcEnergyThreshold = enth ; }
	52	virtual void SetLocalMaxCut(Float_t cut) { fLocMaxCut = cut ; }
	53	virtual void SetLogWeightCut(Float_t w) { fW0 = w ; }
	54	virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; }
	55	virtual void SetPpsdEnergyThreshold(Float_t enth) { fPpsdEnergyThreshold = enth ; }
d15a28e7	56
	57	private:
	58
b2a60966	59	Float_t fA ; // offset of the energy calibration
	60	Float_t fB ; // gain of the energy calibration
	61	Float_t fEmcClusteringThreshold ; // minimum energy to include a EMC digit in a cluster
	62	Float_t fEmcEnergyThreshold ; // minimum energy of EMC digit to be considered
	63	Float_t fLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster
	64	Int_t fNumberOfEmcClusters ; // number of EMC clusters found
	65	Int_t fNumberOfPpsdClusters ; // number of PPSD clusters found
	66	Float_t fPpsdClusteringThreshold ; // minimum energy to include a PPSD digit in a cluster
	67	Float_t fPpsdEnergyThreshold ; // minimum energy of PPSD digit to be considered
	68	Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
92862013	69
b2a60966	70	ClassDef(AliPHOSClusterizerv1,1) // Clusterizer implementation version 1
d15a28e7	71
	72	};
	73
	74	#endif // AliPHOSCLUSTERIZERV1_H