[u/mrichter/AliRoot.git] / EMCAL / AliEMCALRecParam.h

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

/* $Id$ */

//-----------------------------------------------------------------------------
// Container of EMCAL reconstruction parameters
// The purpose of this object is to store it to OCDB
// and retrieve it in AliEMCALClusterizerv1, AliEMCALPID,
// AliEMCALTracker and use it to configure AliEMCALRawUtils
// 
// 
// Author: Yuri Kharlov
//-----------------------------------------------------------------------------

// --- ROOT system ---

#include "AliDetectorRecoParam.h" 
#include "AliLog.h"

class AliEMCALRecParam : public AliDetectorRecoParam
{
 public:
  
  AliEMCALRecParam() ;
  AliEMCALRecParam(const AliEMCALRecParam& recParam);
  AliEMCALRecParam& operator = (const AliEMCALRecParam& recParam);
  virtual ~AliEMCALRecParam() {}
  
  //Clustering (Unfolding : Cynthia)
  Float_t GetClusteringThreshold() const     {return fClusteringThreshold ;}
  Float_t GetW0                 () const     {return fW0                  ;}
  Float_t GetMinECut            () const     {return fMinECut             ;}
  Float_t GetLocMaxCut          () const     {return fLocMaxCut           ;}
  Float_t GetTimeCut            () const     {return fTimeCut             ;}
  Bool_t  GetUnfold             () const     {return fUnfold              ;}
  void SetClusteringThreshold(Float_t thrsh)     {fClusteringThreshold = thrsh;}
  void SetW0                 (Float_t w0)        {fW0 = w0                ;}
  void SetMinECut            (Float_t minEcut)   {fMinECut = minEcut      ;}
  void SetLocMaxCut          (Float_t locMaxCut) {fLocMaxCut = locMaxCut  ;}
  void SetTimeCut            (Float_t timeCut)   {fTimeCut = timeCut  ;}
  void SetUnfold             (Bool_t unfold)     {fUnfold = unfold ; if(fUnfold) AliWarning("Cluster Unfolding ON. Implementing only for eta=0 case!!!");}
  
  //PID (Guenole)
  Double_t GetGamma(Int_t i, Int_t j) const       {return fGamma[i][j];} 
  Double_t GetGammaEnergyProb(Int_t i) const      {return fGammaEnergyProb[i];} 
  Double_t GetGamma1to10(Int_t i, Int_t j) const  {return fGamma1to10[i][j];}   // not used
  Double_t GetHadron(Int_t i, Int_t j) const      {return fHadron[i][j];}
  Double_t GetHadron1to10(Int_t i, Int_t j) const {return fHadron1to10[i][j];}   // not used
  Double_t GetHadronEnergyProb(Int_t i) const     {return fHadronEnergyProb[i];}
  Double_t GetPiZero(Int_t i, Int_t j) const      {return fPiZero[i][j];}
  Double_t GetPiZeroEnergyProb(Int_t i) const     {return fPiZeroEnergyProb[i];}
  
  void SetGamma(Int_t i, Int_t j,Double_t param )       {fGamma[i][j]=param;}
  void SetGammaEnergyProb(Int_t i, Double_t param )     {fGammaEnergyProb[i]=param;}
  void SetGamma1to10(Int_t i, Int_t j,Double_t param )  {fGamma1to10[i][j]=param;}
  void SetHadron(Int_t i, Int_t j,Double_t param )      {fHadron[i][j]=param;}
  void SetHadron1to10(Int_t i, Int_t j,Double_t param ) {fHadron1to10[i][j]=param;}
  void SetHadronEnergyProb(Int_t i,Double_t param )     {fHadronEnergyProb[i]=param;}
  void SetPiZero(Int_t i, Int_t j,Double_t param)       {fPiZero[i][j]=param;}
  void SetPiZeroEnergyProb(Int_t i,Double_t param)      {fPiZeroEnergyProb[i]=param;}
  
  //Track Matching (Alberto)
  /* track matching cut setters */
  void SetTrkCutX(Double_t value)        {fTrkCutX = value;}
  void SetTrkCutY(Double_t value)        {fTrkCutY = value;}
  void SetTrkCutZ(Double_t value)        {fTrkCutZ = value;}
  void SetTrkCutR(Double_t value)        {fTrkCutR = value;}
  void SetTrkCutAlphaMin(Double_t value) {fTrkCutAlphaMin = value;}
  void SetTrkCutAlphaMax(Double_t value) {fTrkCutAlphaMax = value;}
  void SetTrkCutAngle(Double_t value)    {fTrkCutAngle = value;}
  void SetTrkCutNITS(Double_t value)     {fTrkCutNITS = value;}
  void SetTrkCutNTPC(Double_t value)     {fTrkCutNTPC = value;}
  /* track matching cut getters */
  Double_t GetTrkCutX() const        {return fTrkCutX;}
  Double_t GetTrkCutY() const        {return fTrkCutY;}
  Double_t GetTrkCutZ() const        {return fTrkCutZ;}
  Double_t GetTrkCutR() const        {return fTrkCutR;}
  Double_t GetTrkCutAlphaMin() const {return fTrkCutAlphaMin;}
  Double_t GetTrkCutAlphaMax() const {return fTrkCutAlphaMax;}
  Double_t GetTrkCutAngle() const    {return fTrkCutAngle;}
  Double_t GetTrkCutNITS() const     {return fTrkCutNITS;}
  Double_t GetTrkCutNTPC() const     {return fTrkCutNTPC;}
  
  //Raw signal fitting (Jenn)
  /* raw signal setters */
  void SetHighLowGainFactor(Double_t value) {fHighLowGainFactor = value;}
  void SetOrderParameter(Int_t value)       {fOrderParameter = value;}
  void SetTau(Double_t value)               {fTau = value;}
  void SetNoiseThreshold(Int_t value)       {fNoiseThreshold = value;}
  void SetNPedSamples(Int_t value)          {fNPedSamples = value;} 
  void SetRemoveBadChannels(Bool_t val)     {fRemoveBadChannels=val; }
  void SetFittingAlgorithm(Int_t val)       {fFittingAlgorithm=val; }

  /* raw signal getters */
  Double_t GetHighLowGainFactor() const {return fHighLowGainFactor;}
  Int_t    GetOrderParameter()    const {return fOrderParameter;}
  Double_t GetTau()               const {return fTau;}
  Int_t    GetNoiseThreshold()    const {return fNoiseThreshold;}
  Int_t    GetNPedSamples()       const {return fNPedSamples;}
  Bool_t   GetRemoveBadChannels() const {return fRemoveBadChannels;}
  Int_t    GetFittingAlgorithm()  const {return fFittingAlgorithm; }
	
	
	virtual void Print(Option_t * option="") const ;
  
  static AliEMCALRecParam* GetDefaultParameters();
  static AliEMCALRecParam* GetLowFluxParam();
  static AliEMCALRecParam* GetHighFluxParam();
  static AliEMCALRecParam* GetCalibParam();
  static AliEMCALRecParam* GetCosmicParam();
  
  static const  TObjArray* GetMappings();
  
 private:
  //Clustering
  Float_t fClusteringThreshold ; // Minimum energy to seed a EC digit in a cluster
  Float_t fW0 ;                  // Logarithmic weight for the cluster center of gravity calculation
  Float_t fMinECut;              // Minimum energy for a digit to be a member of a cluster
  Bool_t  fUnfold;               // Flag to perform cluster unfolding
  Float_t fLocMaxCut;            // Minimum energy difference to consider local maxima in a cluster
  Float_t fTimeCut ;             // Maximum time of digits in EMC cluster

  //PID (Guenole)
  Double_t fGamma[6][6];         // Parameter to Compute PID for photons     
  Double_t fGamma1to10[6][6];    // Parameter to Compute PID not used
  Double_t fHadron[6][6]; 	     // Parameter to Compute PID for hadrons  	 
  Double_t fHadron1to10[6][6]; 	 // Parameter to Compute PID for hadrons between 1 and 10 GeV  	 
  Double_t fHadronEnergyProb[6]; // Parameter to Compute PID for energy ponderation for hadrons  	 
  Double_t fPiZeroEnergyProb[6]; // Parameter to Compute PID for energy ponderation for Pi0  	 
  Double_t fGammaEnergyProb[6];  // Parameter to Compute PID for energy ponderation for gamma  	 
  Double_t fPiZero[6][6];        // Parameter to Compute PID for pi0  	 
  
  
  //Track-Matching (Alberto)
  Double_t  fTrkCutX;              // X-difference cut for track matching
  Double_t  fTrkCutY;              // Y-difference cut for track matching
  Double_t  fTrkCutZ;              // Z-difference cut for track matching
  Double_t  fTrkCutR;              // cut on allowed track-cluster distance
  Double_t  fTrkCutAlphaMin;       // cut on 'alpha' parameter for track matching (min)
  Double_t  fTrkCutAlphaMax;       // cut on 'alpha' parameter for track matching (min)
  Double_t  fTrkCutAngle;          // cut on relative angle between different track points for track matching
  Double_t  fTrkCutNITS;           // Number of ITS hits for track matching
  Double_t  fTrkCutNTPC;           // Number of TPC hits for track matching
  
  //Raw signal fitting parameters (Jenn)
  Double_t fHighLowGainFactor;     // gain factor to convert between high and low gain
  Int_t    fOrderParameter;        // order parameter for raw signal fit
  Double_t fTau;                   // decay constant for raw signal fit
  Int_t    fNoiseThreshold;        // threshold to consider signal or noise
  Int_t    fNPedSamples;           // number of time samples to use in pedestal calculation
  Bool_t   fRemoveBadChannels;     // select if bad channels are removed before fitting
  Int_t    fFittingAlgorithm;      // select the fitting algorithm

  static TObjArray* fgkMaps;       // ALTRO mappings for RCU0..RCUX
  
  ClassDef(AliEMCALRecParam,9)     // Reconstruction parameters
    
    } ;

#endif //  ALIEMCALRECPARAM_H
Commit	Line	Data
3a8be91c	1	#ifndef ALIEMCALRECPARAM_H
	2	#define ALIEMCALRECPARAM_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	// Container of EMCAL reconstruction parameters
	10	// The purpose of this object is to store it to OCDB
b4133f05	11	// and retrieve it in AliEMCALClusterizerv1, AliEMCALPID,
	12	// AliEMCALTracker and use it to configure AliEMCALRawUtils
	13	//
	14	//
3a8be91c	15	// Author: Yuri Kharlov
	16	//-----------------------------------------------------------------------------
	17
	18	// --- ROOT system ---
	19
413e6b81	20	#include "AliDetectorRecoParam.h"
225cd96d	21	#include "AliLog.h"
3a8be91c	22
413e6b81	23	class AliEMCALRecParam : public AliDetectorRecoParam
3a8be91c	24	{
1e7c9b89	25	public:
3a8be91c	26
3a8be91c	27	AliEMCALRecParam() ;
413e6b81	28	AliEMCALRecParam(const AliEMCALRecParam& recParam);
413e6b81	29	AliEMCALRecParam& operator = (const AliEMCALRecParam& recParam);
3a8be91c	30	virtual ~AliEMCALRecParam() {}
1e7c9b89	31
225cd96d	32	//Clustering (Unfolding : Cynthia)
1e7c9b89	33	Float_t GetClusteringThreshold() const {return fClusteringThreshold ;}
	34	Float_t GetW0 () const {return fW0 ;}
	35	Float_t GetMinECut () const {return fMinECut ;}
	36	Float_t GetLocMaxCut () const {return fLocMaxCut ;}
a435f763	37	Float_t GetTimeCut () const {return fTimeCut ;}
1e7c9b89	38	Bool_t GetUnfold () const {return fUnfold ;}
	39	void SetClusteringThreshold(Float_t thrsh) {fClusteringThreshold = thrsh;}
	40	void SetW0 (Float_t w0) {fW0 = w0 ;}
	41	void SetMinECut (Float_t minEcut) {fMinECut = minEcut ;}
	42	void SetLocMaxCut (Float_t locMaxCut) {fLocMaxCut = locMaxCut ;}
a435f763	43	void SetTimeCut (Float_t timeCut) {fTimeCut = timeCut ;}
1e7c9b89	44	void SetUnfold (Bool_t unfold) {fUnfold = unfold ; if(fUnfold) AliWarning("Cluster Unfolding ON. Implementing only for eta=0 case!!!");}
1e7c9b89	45
8ba062b1	46	//PID (Guenole)
1e7c9b89	47	Double_t GetGamma(Int_t i, Int_t j) const {return fGamma[i][j];}
	48	Double_t GetGammaEnergyProb(Int_t i) const {return fGammaEnergyProb[i];}
	49	Double_t GetGamma1to10(Int_t i, Int_t j) const {return fGamma1to10[i][j];} // not used
	50	Double_t GetHadron(Int_t i, Int_t j) const {return fHadron[i][j];}
	51	Double_t GetHadron1to10(Int_t i, Int_t j) const {return fHadron1to10[i][j];} // not used
	52	Double_t GetHadronEnergyProb(Int_t i) const {return fHadronEnergyProb[i];}
	53	Double_t GetPiZero(Int_t i, Int_t j) const {return fPiZero[i][j];}
	54	Double_t GetPiZeroEnergyProb(Int_t i) const {return fPiZeroEnergyProb[i];}
	55
	56	void SetGamma(Int_t i, Int_t j,Double_t param ) {fGamma[i][j]=param;}
	57	void SetGammaEnergyProb(Int_t i, Double_t param ) {fGammaEnergyProb[i]=param;}
	58	void SetGamma1to10(Int_t i, Int_t j,Double_t param ) {fGamma1to10[i][j]=param;}
	59	void SetHadron(Int_t i, Int_t j,Double_t param ) {fHadron[i][j]=param;}
	60	void SetHadron1to10(Int_t i, Int_t j,Double_t param ) {fHadron1to10[i][j]=param;}
	61	void SetHadronEnergyProb(Int_t i,Double_t param ) {fHadronEnergyProb[i]=param;}
	62	void SetPiZero(Int_t i, Int_t j,Double_t param) {fPiZero[i][j]=param;}
	63	void SetPiZeroEnergyProb(Int_t i,Double_t param) {fPiZeroEnergyProb[i]=param;}
	64
8ba062b1	65	//Track Matching (Alberto)
	66	/* track matching cut setters */
	67	void SetTrkCutX(Double_t value) {fTrkCutX = value;}
	68	void SetTrkCutY(Double_t value) {fTrkCutY = value;}
	69	void SetTrkCutZ(Double_t value) {fTrkCutZ = value;}
	70	void SetTrkCutR(Double_t value) {fTrkCutR = value;}
	71	void SetTrkCutAlphaMin(Double_t value) {fTrkCutAlphaMin = value;}
	72	void SetTrkCutAlphaMax(Double_t value) {fTrkCutAlphaMax = value;}
	73	void SetTrkCutAngle(Double_t value) {fTrkCutAngle = value;}
1e7c9b89	74	void SetTrkCutNITS(Double_t value) {fTrkCutNITS = value;}
1e7c9b89	75	void SetTrkCutNTPC(Double_t value) {fTrkCutNTPC = value;}
8ba062b1	76	/* track matching cut getters */
	77	Double_t GetTrkCutX() const {return fTrkCutX;}
	78	Double_t GetTrkCutY() const {return fTrkCutY;}
	79	Double_t GetTrkCutZ() const {return fTrkCutZ;}
	80	Double_t GetTrkCutR() const {return fTrkCutR;}
	81	Double_t GetTrkCutAlphaMin() const {return fTrkCutAlphaMin;}
	82	Double_t GetTrkCutAlphaMax() const {return fTrkCutAlphaMax;}
	83	Double_t GetTrkCutAngle() const {return fTrkCutAngle;}
1e7c9b89	84	Double_t GetTrkCutNITS() const {return fTrkCutNITS;}
	85	Double_t GetTrkCutNTPC() const {return fTrkCutNTPC;}
	86
b4133f05	87	//Raw signal fitting (Jenn)
	88	/* raw signal setters */
	89	void SetHighLowGainFactor(Double_t value) {fHighLowGainFactor = value;}
	90	void SetOrderParameter(Int_t value) {fOrderParameter = value;}
	91	void SetTau(Double_t value) {fTau = value;}
	92	void SetNoiseThreshold(Int_t value) {fNoiseThreshold = value;}
9f467289	93	void SetNPedSamples(Int_t value) {fNPedSamples = value;}
	94	void SetRemoveBadChannels(Bool_t val) {fRemoveBadChannels=val; }
	95	void SetFittingAlgorithm(Int_t val) {fFittingAlgorithm=val; }
	96
b4133f05	97	/* raw signal getters */
	98	Double_t GetHighLowGainFactor() const {return fHighLowGainFactor;}
	99	Int_t GetOrderParameter() const {return fOrderParameter;}
	100	Double_t GetTau() const {return fTau;}
	101	Int_t GetNoiseThreshold() const {return fNoiseThreshold;}
	102	Int_t GetNPedSamples() const {return fNPedSamples;}
9f467289	103	Bool_t GetRemoveBadChannels() const {return fRemoveBadChannels;}
	104	Int_t GetFittingAlgorithm() const {return fFittingAlgorithm; }
	105
	106
	107	virtual void Print(Option_t * option="") const ;
1e7c9b89	108
413e6b81	109	static AliEMCALRecParam* GetDefaultParameters();
f5fc991a	110	static AliEMCALRecParam* GetLowFluxParam();
f5fc991a	111	static AliEMCALRecParam* GetHighFluxParam();
98cf874d	112	static AliEMCALRecParam* GetCalibParam();
98cf874d	113	static AliEMCALRecParam* GetCosmicParam();
1e7c9b89	114
feedcab9	115	static const TObjArray* GetMappings();
1e7c9b89	116
1e7c9b89	117	private:
8ba062b1	118	//Clustering
a435f763	119	Float_t fClusteringThreshold ; // Minimum energy to seed a EC digit in a cluster
a435f763	120	Float_t fW0 ; // Logarithmic weight for the cluster center of gravity calculation
3a8be91c	121	Float_t fMinECut; // Minimum energy for a digit to be a member of a cluster
a435f763	122	Bool_t fUnfold; // Flag to perform cluster unfolding
	123	Float_t fLocMaxCut; // Minimum energy difference to consider local maxima in a cluster
	124	Float_t fTimeCut ; // Maximum time of digits in EMC cluster
	125
8ba062b1	126	//PID (Guenole)
1e7c9b89	127	Double_t fGamma[6][6]; // Parameter to Compute PID for photons
	128	Double_t fGamma1to10[6][6]; // Parameter to Compute PID not used
	129	Double_t fHadron[6][6]; // Parameter to Compute PID for hadrons
	130	Double_t fHadron1to10[6][6]; // Parameter to Compute PID for hadrons between 1 and 10 GeV
	131	Double_t fHadronEnergyProb[6]; // Parameter to Compute PID for energy ponderation for hadrons
	132	Double_t fPiZeroEnergyProb[6]; // Parameter to Compute PID for energy ponderation for Pi0
	133	Double_t fGammaEnergyProb[6]; // Parameter to Compute PID for energy ponderation for gamma
	134	Double_t fPiZero[6][6]; // Parameter to Compute PID for pi0
	135
	136
8ba062b1	137	//Track-Matching (Alberto)
	138	Double_t fTrkCutX; // X-difference cut for track matching
	139	Double_t fTrkCutY; // Y-difference cut for track matching
	140	Double_t fTrkCutZ; // Z-difference cut for track matching
	141	Double_t fTrkCutR; // cut on allowed track-cluster distance
	142	Double_t fTrkCutAlphaMin; // cut on 'alpha' parameter for track matching (min)
	143	Double_t fTrkCutAlphaMax; // cut on 'alpha' parameter for track matching (min)
	144	Double_t fTrkCutAngle; // cut on relative angle between different track points for track matching
1e7c9b89	145	Double_t fTrkCutNITS; // Number of ITS hits for track matching
	146	Double_t fTrkCutNTPC; // Number of TPC hits for track matching
	147
b4133f05	148	//Raw signal fitting parameters (Jenn)
9f467289	149	Double_t fHighLowGainFactor; // gain factor to convert between high and low gain
	150	Int_t fOrderParameter; // order parameter for raw signal fit
	151	Double_t fTau; // decay constant for raw signal fit
	152	Int_t fNoiseThreshold; // threshold to consider signal or noise
	153	Int_t fNPedSamples; // number of time samples to use in pedestal calculation
	154	Bool_t fRemoveBadChannels; // select if bad channels are removed before fitting
	155	Int_t fFittingAlgorithm; // select the fitting algorithm
	156
feedcab9	157	static TObjArray* fgkMaps; // ALTRO mappings for RCU0..RCUX
1e7c9b89	158
9f467289	159	ClassDef(AliEMCALRecParam,9) // Reconstruction parameters
1e7c9b89	160
1e7c9b89	161	} ;
3a8be91c	162
3a8be91c	163	#endif // ALIEMCALRECPARAM_H
8ba062b1	164