[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" 
#include "AliPHOSGeometry.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 GetLocalMaxCutCPV()          { return fLocMaxCutCPV;} 
  virtual Float_t GetLogWeightCutCPV()         { return fW0CPV;}  
  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 SetLocalMaxCutCPV(Float_t cut)            { fLocMaxCutCPV = cut ; }
  virtual void SetLogWeightCutCPV(Float_t w)             { fW0CPV = 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
  AliPHOSGeometry * fGeom ;          // pointer to geometry
  Float_t fEmcClusteringThreshold ;  // minimum energy to include a EMC digit in a cluster
  Float_t fEmcEnergyThreshold ;      // minimum energy of EMC digit to be considered
  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 fLocMaxCut ;               // minimum energy difference to distinguish local maxima in a cluster
  Float_t fW0 ;                      // logarithmic weight for the cluster center of gravity calculation

  Float_t fLocMaxCutCPV ;            // minimum energy difference to distinguish local maxima in a CPV cluster
  Float_t fW0CPV ;                   // logarithmic weight for the CPV 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"
6b9b768e	21	#include "AliPHOSGeometry.h"
d15a28e7	22
	23
	24
	25	class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
	26
	27	public:
	28
	29	AliPHOSClusterizerv1() ; // ctor
	30	virtual ~AliPHOSClusterizerv1(){} ; // dtor
	31
9f616d61	32	Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2) ; // Checks if digits are in neighbour cells
	33	Float_t Calibrate(Int_t Amp){ return (fA + fB * Amp) ;} // Tranforms Amp to energy
	34	void FillandSort(const DigitsList * dl, TObjArray * tl) ; // Sorts the list according to increasing id
4a646085	35	virtual void GetNumberOfClustersFound(Int_t * numb) ;
	36
	37	virtual void GetCalibrationParameters(Float_t & A, Float_t &B) { A = fA; B = fB; }
d72dfbc3	38	virtual Float_t GetEmcClusteringThreshold() { return fEmcClusteringThreshold;}
	39	virtual Float_t GetEmcEnergyThreshold() { return fEmcEnergyThreshold; }
	40	virtual Float_t GetLocalMaxCut() { return fLocMaxCut;}
	41	virtual Float_t GetLogWeightCut() { return fW0;}
	42	virtual Float_t GetLocalMaxCutCPV() { return fLocMaxCutCPV;}
	43	virtual Float_t GetLogWeightCutCPV() { return fW0CPV;}
416692f4	44	virtual Float_t GetPpsdClusteringThreshold() { return fPpsdClusteringThreshold; }
d72dfbc3	45	virtual Float_t GetPpsdEnergyThreshold() { return fPpsdEnergyThreshold; }
4a646085	46
92862013	47	virtual Bool_t IsInEmc(AliPHOSDigit * digit) ; // Tells if id digit is in EMC
88714635	48	virtual void MakeClusters(const DigitsList * dl,
	49	AliPHOSRecPoint::RecPointsList * emcl,
	50	AliPHOSRecPoint::RecPointsList * ppsdl) ; // does the job
6ad0bfa0	51	virtual void PrintParameters() ;
6ad0bfa0	52	virtual void SetCalibrationParameters(Float_t A,Float_t B){ fA = A ; fB = B;}
d72dfbc3	53	virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
	54	virtual void SetEmcEnergyThreshold(Float_t enth) { fEmcEnergyThreshold = enth ; }
	55	virtual void SetLocalMaxCut(Float_t cut) { fLocMaxCut = cut ; }
	56	virtual void SetLogWeightCut(Float_t w) { fW0 = w ; }
	57	virtual void SetLocalMaxCutCPV(Float_t cut) { fLocMaxCutCPV = cut ; }
	58	virtual void SetLogWeightCutCPV(Float_t w) { fW0CPV = w ; }
6ad0bfa0	59	virtual void SetPpsdClusteringThreshold(Float_t cluth) { fPpsdClusteringThreshold = cluth ; }
d72dfbc3	60	virtual void SetPpsdEnergyThreshold(Float_t enth) { fPpsdEnergyThreshold = enth ; }
d15a28e7	61
	62	private:
	63
b2a60966	64	Float_t fA ; // offset of the energy calibration
b2a60966	65	Float_t fB ; // gain of the energy calibration
6b9b768e	66	AliPHOSGeometry * fGeom ; // pointer to geometry
b2a60966	67	Float_t fEmcClusteringThreshold ; // minimum energy to include a EMC digit in a cluster
b2a60966	68	Float_t fEmcEnergyThreshold ; // minimum energy of EMC digit to be considered
b2a60966	69	Int_t fNumberOfEmcClusters ; // number of EMC clusters found
	70	Int_t fNumberOfPpsdClusters ; // number of PPSD clusters found
	71	Float_t fPpsdClusteringThreshold ; // minimum energy to include a PPSD digit in a cluster
	72	Float_t fPpsdEnergyThreshold ; // minimum energy of PPSD digit to be considered
d72dfbc3	73	Float_t fLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster
b2a60966	74	Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
d72dfbc3	75
	76	Float_t fLocMaxCutCPV ; // minimum energy difference to distinguish local maxima in a CPV cluster
	77	Float_t fW0CPV ; // logarithmic weight for the CPV cluster center of gravity calculation
92862013	78
b2a60966	79	ClassDef(AliPHOSClusterizerv1,1) // Clusterizer implementation version 1
d15a28e7	80
	81	};
	82
	83	#endif // AliPHOSCLUSTERIZERV1_H