/* $Id$ */
+/* History of cvs commits:
+ *
+ * $Log$
+ * Revision 1.50 2007/03/28 19:18:15 kharlov
+ * RecPoints recalculation in TSM removed
+ *
+ * Revision 1.49 2007/03/06 06:51:27 kharlov
+ * Calculation of cluster properties dep. on vertex posponed to TrackSegmentMaker
+ *
+ * Revision 1.48 2006/08/30 16:12:52 kharlov
+ * Reconstruction of raw data from beam test 2006 (B.Polichtchouk)
+ *
+ * Revision 1.47 2006/08/25 16:56:30 kharlov
+ * Compliance with Effective C++
+ *
+ * Revision 1.46 2006/08/01 12:20:17 cvetan
+ * 1. Adding a possibility to read and reconstruct an old rcu formatted raw data. This is controlled by an option of AliReconstruction and AliPHOSReconstructor. 2. In case of raw data processing (without galice.root) create the default AliPHOSGeometry object. Most likely this should be moved to the CDB
+ *
+ * Revision 1.45 2006/04/29 20:26:46 hristov
+ * Separate EMC and CPV calibration (Yu.Kharlov)
+ *
+ * Revision 1.44 2005/09/02 14:32:07 kharlov
+ * Calibration of raw data
+ *
+ * Revision 1.43 2005/05/28 14:19:04 schutz
+ * Compilation warnings fixed by T.P.
+ *
+ */
+
//_________________________________________________________________________
// 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 ---
-
+class TClonesArray ;
+class TVector3 ;
// --- Standard library ---
// --- AliRoot header files ---
#include "AliPHOSClusterizer.h"
-#include "AliPHOSDigit.h"
-#include "AliPHOSGeometry.h"
-
-
+class AliPHOSEmcRecPoint ;
+class AliPHOSDigit ;
+class AliPHOSDigitizer ;
+class AliPHOSGeometry ;
+class AliPHOSCalibData ;
class AliPHOSClusterizerv1 : public AliPHOSClusterizer {
public:
- AliPHOSClusterizerv1() ; // ctor
- virtual ~AliPHOSClusterizerv1(){} ; // dtor
+ AliPHOSClusterizerv1() ;
+ AliPHOSClusterizerv1(const TString alirunFileNameFile, const TString eventFolderName = AliConfig::GetDefaultEventFolderName());
+ AliPHOSClusterizerv1(const AliPHOSClusterizerv1 & clu) ;
+ virtual ~AliPHOSClusterizerv1() ;
- 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 Int_t AreNeighbours(AliPHOSDigit * d1, AliPHOSDigit * d2)const ;
+ // Checks if digits are in neighbour cells
+
+ virtual Float_t CalibrateCPV(Int_t amp, Int_t absId) ; // Tranforms CPV Amp to energy
+ virtual Float_t CalibrateEMC(Float_t amp, Int_t absId) ; // Tranforms EMC 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 Int_t GetRecPointsInRun() const {return fRecPointsInRun ;}
+
+ virtual void Exec(Option_t *option); // Does the job
+
+ void Print(const Option_t * = "")const ;
+
+ void SetEmcMinE(Float_t e){fEmcMinE = e ;}
+ void SetCpvMinE(Float_t e){fCpvMinE = e ;}
+ 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 SetUnfolding(Bool_t toUnfold = kTRUE ) { fToUnfold = toUnfold ;}
+ //Switch to "on flyght" mode, without writing to TreeR and file
+ void SetWriting(Bool_t toWrite = kFALSE){fWrite = toWrite;}
+ static Double_t ShowerShape(Double_t x, Double_t z) ; // 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 passed to MINUIT
+ void Unload() ;
+ virtual const char * Version() const { return "clu-v1"; }
+
+ virtual void SetOldRCUFormat(Bool_t rcuFormat = kFALSE)
+ { fIsOldRCUFormat = rcuFormat; };
+
+protected:
+
+ void WriteRecPoints() ;
+ 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
+ void CleanDigits(TClonesArray * digits) ;
+ void GetCalibrationParameters(void);
+ void SetDistancesToBadChannels();
+
+ AliPHOSClusterizerv1 & operator = (const AliPHOSClusterizerv1 & obj);
- 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:
+
+ const TString BranchName() const ;
- 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
- 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
+ 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)
+ Int_t fEmcCrystals ; // number of EMC cristals in PHOS
+
+ Bool_t fToUnfold ; // To perform unfolding
+ Bool_t fWrite ; // Write RecPoints to TreeR
+
+ Int_t fNumberOfEmcClusters ; // number of EMC clusters found
+ Int_t fNumberOfCpvClusters ; // number of CPV clusters found
+
+ //Calibration parameters
+ AliPHOSCalibData * fCalibData ; //! Calibration database if aval.
+ 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 start EMC cluster
+ Float_t fCpvClusteringThreshold ; // minimum energy to start CPV cluster
+ Float_t fEmcMinE ; // minimum energy of digit to be included into cluster
+ Float_t fCpvMinE ; // minimum energy of digit to be included into 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
-
- ClassDef(AliPHOSClusterizerv1,1) // Clusterizer implementation version 1
+ 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
+
+ Bool_t fIsOldRCUFormat; // assume old RCU raw data format
+
+ ClassDef(AliPHOSClusterizerv1,4) // Clusterizer implementation version 1
};
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff