[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                     
//  Performs clusterization (collects neighbouring active cells) and 
//  unfolding of the clusters with several local maxima.  
//  results are stored in TreeR#, branches PHOSEmcRP (EMC recPoints),
//  PHOSCpvRP (CPV RecPoints) and AliPHOSClusterizer
//
//*-- Author: Yves Schutz (SUBATECH)

// --- ROOT system ---

// --- Standard library ---

// --- AliRoot header files ---

#include "AliPHOSClusterizer.h"
class AliPHOSEmcRecPoint ; 
class AliPHOSDigit ;
class AliPHOSDigitizer ;
class AliPHOSGeometry ;


class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
  
public:
  
  AliPHOSClusterizerv1() ;         
  AliPHOSClusterizerv1(const char * headerFile, const char * name = "Default", const Bool_t toSplit=kFALSE);
  virtual ~AliPHOSClusterizerv1()  ;
  
  virtual Int_t   AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2)const ; 
                               // Checks if digits are in neighbour cells 

  virtual Float_t Calibrate(Int_t amp, Int_t absId)const ;  // Tranforms Amp to energy 

  virtual void    GetNumberOfClustersFound(int * numb )const{  numb[0] = fNumberOfEmcClusters ; 
                                                               numb[1] = fNumberOfCpvClusters ; }

  virtual Float_t GetEmcClusteringThreshold()const{ return fEmcClusteringThreshold;}
  virtual Float_t GetEmcLocalMaxCut()const        { return fEmcLocMaxCut;} 
  virtual Float_t GetEmcLogWeight()const          { return fW0;}  
  virtual Float_t GetEmcTimeGate() const          { return fEmcTimeGate ; }
  virtual Float_t GetCpvClusteringThreshold()const{ return fCpvClusteringThreshold;  } 
  virtual Float_t GetCpvLocalMaxCut()const        { return fCpvLocMaxCut;} 
  virtual Float_t GetCpvLogWeight()const          { return fW0CPV;}  
  virtual const char *  GetRecPointsBranch() const{ return GetName() ;}
  virtual const Int_t GetRecPointsInRun() const   {return fRecPointsInRun ;} 
  //  virtual const char *  GetDigitsBranch() const         { return (char*) fDigitsBranchTitle.Data() ;}

  void    Exec(Option_t *option);                // Does the job

  virtual void Print(Option_t * option)const ;

  virtual void SetEmcClusteringThreshold(Float_t cluth)  { fEmcClusteringThreshold = cluth ; }
  virtual void SetEmcLocalMaxCut(Float_t cut)            { fEmcLocMaxCut = cut ; }
  virtual void SetEmcLogWeight(Float_t w)                { fW0 = w ; }
  virtual void SetEmcTimeGate(Float_t gate)              { fEmcTimeGate = gate ;}
  virtual void SetCpvClusteringThreshold(Float_t cluth)  { fCpvClusteringThreshold = cluth ; }
  virtual void SetCpvLocalMaxCut(Float_t cut)            { fCpvLocMaxCut = cut ; }
  virtual void SetCpvLogWeight(Float_t w)                { fW0CPV = w ; }
  //  virtual void SetDigitsBranch(const char * title) { fDigitsBranchTitle = title  ;}
  //  virtual void SetRecPointsBranch(const char *title){fRecPointsBranchTitle = title; }
  virtual void SetUnfolding(Bool_t toUnfold = kTRUE )    { fToUnfold = toUnfold ;}  
  static Double_t ShowerShape(Double_t r) ; // Shape of EM shower used in unfolding; 
                                            //class member function (not object member function)
  static void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret, Double_t * x, Int_t iflag)  ;
                                            // Chi^2 of the fit. Should be static to be passes to MINUIT
  virtual const char * Version() const { return "clu-v1" ; }  

protected:

  void           WriteRecPoints(Int_t event) ;
  virtual void   MakeClusters( ) ;            
  virtual Bool_t IsInEmc (AliPHOSDigit * digit)const ;     // Tells if id digit is in EMC
  virtual Bool_t IsInCpv (AliPHOSDigit * digit)const ;     // Tells if id digit is in CPV

  
private:

  const   TString BranchName() const ; 
  void    GetCalibrationParameters(void) ;
  
  Bool_t  FindFit(AliPHOSEmcRecPoint * emcRP, AliPHOSDigit ** MaxAt, Float_t * maxAtEnergy, 
		  Int_t NPar, Float_t * FitParametres) const; //Used in UnfoldClusters, calls TMinuit
  void    Init() ;
  void    InitParameters() ;

  virtual void   MakeUnfolding() ;
  void           UnfoldCluster(AliPHOSEmcRecPoint * iniEmc,Int_t Nmax, 
		       AliPHOSDigit ** maxAt,Float_t * maxAtEnergy ) ; //Unfolds cluster using TMinuit package
  void           PrintRecPoints(Option_t * option) ;

private:

  Bool_t  fDefaultInit;              //! Says if the task was created by defaut ctor (only parameters are initialized)
  //  TString fFrom ;                    // name of Digits 
  //  TString fHeaderFileName ;          // name of the file which contains gAlice, Tree headers etc.
  //  TString fDigitsBranchTitle ;       // name of the file, where digits branch is stored
  //  TString fRecPointsBranchTitle ;    // name of the file, where RecPoints branchs are stored

  Int_t   fEmcCrystals ;             // number of EMC cristalls in PHOS

  Bool_t  fToUnfold ;                // To perform unfolding 

  Int_t   fNumberOfEmcClusters ;     // number of EMC clusters found 
  Int_t   fNumberOfCpvClusters ;     // number of CPV clusters found
 
  //Calibration parameters... to be replaced by database 
  Float_t fADCchanelEmc ;           // width of one ADC channel in GeV
  Float_t fADCpedestalEmc ;         //
  Float_t fADCchanelCpv ;           // width of one ADC channel in CPV 'popugais'
  Float_t fADCpedestalCpv ;         // 

  Float_t fEmcClusteringThreshold ;  // minimum energy to include a EMC digit in a cluster
  Float_t fCpvClusteringThreshold ;  // minimum energy to include a CPV digit in a cluster
  Float_t fEmcLocMaxCut ;            // minimum energy difference to distinguish local maxima in a cluster
  Float_t fW0 ;                      // logarithmic weight for the cluster center of gravity calculation
  Float_t fCpvLocMaxCut ;            // minimum energy difference to distinguish local maxima in a CPV cluster
  Float_t fW0CPV ;                   // logarithmic weight for the CPV cluster center of gravity calculation
  Int_t fRecPointsInRun ;            //! Total number of recpoints in one run
  Float_t fEmcTimeGate ;             // Maximum time difference between the digits in ont EMC cluster
    
  ClassDef(AliPHOSClusterizerv1,2)   // Clusterizer implementation version 1

};

#endif // AliPHOSCLUSTERIZERV1_H
Commit	Line	Data
	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
	6	/* $Id$ */
	7
	8	//_________________________________________________________________________
	9	// Implementation version 1 of the clusterization algorithm
	10	// Performs clusterization (collects neighbouring active cells) and
	11	// unfolding of the clusters with several local maxima.
	12	// results are stored in TreeR#, branches PHOSEmcRP (EMC recPoints),
	13	// PHOSCpvRP (CPV RecPoints) and AliPHOSClusterizer
	14	//
	15	//*-- Author: Yves Schutz (SUBATECH)
	16
	17	// --- ROOT system ---
	18
	19	// --- Standard library ---
	20
	21	// --- AliRoot header files ---
	22
	23	#include "AliPHOSClusterizer.h"
	24	class AliPHOSEmcRecPoint ;
	25	class AliPHOSDigit ;
	26	class AliPHOSDigitizer ;
	27	class AliPHOSGeometry ;
	28
	29
	30	class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
	31
	32	public:
	33
	34	AliPHOSClusterizerv1() ;
	35	AliPHOSClusterizerv1(const char * headerFile, const char * name = "Default", const Bool_t toSplit=kFALSE);
	36	virtual ~AliPHOSClusterizerv1() ;
	37
	38	virtual Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2)const ;
	39	// Checks if digits are in neighbour cells
	40
	41	virtual Float_t Calibrate(Int_t amp, Int_t absId)const ; // Tranforms Amp to energy
	42
	43	virtual void GetNumberOfClustersFound(int * numb )const{ numb[0] = fNumberOfEmcClusters ;
	44	numb[1] = fNumberOfCpvClusters ; }
	45
	46	virtual Float_t GetEmcClusteringThreshold()const{ return fEmcClusteringThreshold;}
	47	virtual Float_t GetEmcLocalMaxCut()const { return fEmcLocMaxCut;}
	48	virtual Float_t GetEmcLogWeight()const { return fW0;}
	49	virtual Float_t GetEmcTimeGate() const { return fEmcTimeGate ; }
	50	virtual Float_t GetCpvClusteringThreshold()const{ return fCpvClusteringThreshold; }
	51	virtual Float_t GetCpvLocalMaxCut()const { return fCpvLocMaxCut;}
	52	virtual Float_t GetCpvLogWeight()const { return fW0CPV;}
	53	virtual const char * GetRecPointsBranch() const{ return GetName() ;}
	54	virtual const Int_t GetRecPointsInRun() const {return fRecPointsInRun ;}
	55	// virtual const char * GetDigitsBranch() const { return (char*) fDigitsBranchTitle.Data() ;}
	56
	57	void Exec(Option_t *option); // Does the job
	58
	59	virtual void Print(Option_t * option)const ;
	60
	61	virtual void SetEmcClusteringThreshold(Float_t cluth) { fEmcClusteringThreshold = cluth ; }
	62	virtual void SetEmcLocalMaxCut(Float_t cut) { fEmcLocMaxCut = cut ; }
	63	virtual void SetEmcLogWeight(Float_t w) { fW0 = w ; }
	64	virtual void SetEmcTimeGate(Float_t gate) { fEmcTimeGate = gate ;}
	65	virtual void SetCpvClusteringThreshold(Float_t cluth) { fCpvClusteringThreshold = cluth ; }
	66	virtual void SetCpvLocalMaxCut(Float_t cut) { fCpvLocMaxCut = cut ; }
	67	virtual void SetCpvLogWeight(Float_t w) { fW0CPV = w ; }
	68	// virtual void SetDigitsBranch(const char * title) { fDigitsBranchTitle = title ;}
	69	// virtual void SetRecPointsBranch(const char *title){fRecPointsBranchTitle = title; }
	70	virtual void SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold ;}
	71	static Double_t ShowerShape(Double_t r) ; // Shape of EM shower used in unfolding;
	72	//class member function (not object member function)
	73	static void UnfoldingChiSquare(Int_t & nPar, Double_t * Grad, Double_t & fret, Double_t * x, Int_t iflag) ;
	74	// Chi^2 of the fit. Should be static to be passes to MINUIT
	75	virtual const char * Version() const { return "clu-v1" ; }
	76
	77	protected:
	78
	79	void WriteRecPoints(Int_t event) ;
	80	virtual void MakeClusters( ) ;
	81	virtual Bool_t IsInEmc (AliPHOSDigit * digit)const ; // Tells if id digit is in EMC
	82	virtual Bool_t IsInCpv (AliPHOSDigit * digit)const ; // Tells if id digit is in CPV
	83
	84
	85	private:
	86
	87	const TString BranchName() const ;
	88	void GetCalibrationParameters(void) ;
	89
	90	Bool_t FindFit(AliPHOSEmcRecPoint * emcRP, AliPHOSDigit ** MaxAt, Float_t * maxAtEnergy,
	91	Int_t NPar, Float_t * FitParametres) const; //Used in UnfoldClusters, calls TMinuit
	92	void Init() ;
	93	void InitParameters() ;
	94
	95	virtual void MakeUnfolding() ;
	96	void UnfoldCluster(AliPHOSEmcRecPoint * iniEmc,Int_t Nmax,
	97	AliPHOSDigit ** maxAt,Float_t * maxAtEnergy ) ; //Unfolds cluster using TMinuit package
	98	void PrintRecPoints(Option_t * option) ;
	99
	100	private:
	101
	102	Bool_t fDefaultInit; //! Says if the task was created by defaut ctor (only parameters are initialized)
	103	// TString fFrom ; // name of Digits
	104	// TString fHeaderFileName ; // name of the file which contains gAlice, Tree headers etc.
	105	// TString fDigitsBranchTitle ; // name of the file, where digits branch is stored
	106	// TString fRecPointsBranchTitle ; // name of the file, where RecPoints branchs are stored
	107
	108	Int_t fEmcCrystals ; // number of EMC cristalls in PHOS
	109
	110	Bool_t fToUnfold ; // To perform unfolding
	111
	112	Int_t fNumberOfEmcClusters ; // number of EMC clusters found
	113	Int_t fNumberOfCpvClusters ; // number of CPV clusters found
	114
	115	//Calibration parameters... to be replaced by database
	116	Float_t fADCchanelEmc ; // width of one ADC channel in GeV
	117	Float_t fADCpedestalEmc ; //
	118	Float_t fADCchanelCpv ; // width of one ADC channel in CPV 'popugais'
	119	Float_t fADCpedestalCpv ; //
	120
	121	Float_t fEmcClusteringThreshold ; // minimum energy to include a EMC digit in a cluster
	122	Float_t fCpvClusteringThreshold ; // minimum energy to include a CPV digit in a cluster
	123	Float_t fEmcLocMaxCut ; // minimum energy difference to distinguish local maxima in a cluster
	124	Float_t fW0 ; // logarithmic weight for the cluster center of gravity calculation
	125	Float_t fCpvLocMaxCut ; // minimum energy difference to distinguish local maxima in a CPV cluster
	126	Float_t fW0CPV ; // logarithmic weight for the CPV cluster center of gravity calculation
	127	Int_t fRecPointsInRun ; //! Total number of recpoints in one run
	128	Float_t fEmcTimeGate ; // Maximum time difference between the digits in ont EMC cluster
	129
	130	ClassDef(AliPHOSClusterizerv1,2) // Clusterizer implementation version 1
	131
	132	};
	133
	134	#endif // AliPHOSCLUSTERIZERV1_H