1 #ifndef ALIITSCLUSTERFINDERSDD_H
2 #define ALIITSCLUSTERFINDERSDD_H
4 ////////////////////////////////////////////////
5 // ITS Cluster Finder Class //
6 ////////////////////////////////////////////////
13 #include "AliITSClusterFinder.h"
14 #include "AliITSDetTypeRec.h"
17 class AliITSCalibration;
18 class AliITSCalibrationSDD;
19 class AliITSsegmentation;
20 class AliITSsegmentationSDD;
23 class AliITSClusterFinderSDD : public AliITSClusterFinder{
25 AliITSClusterFinderSDD
26 (AliITSDetTypeRec* dettyp,
27 TClonesArray *digits,TClonesArray *recpoints);
28 AliITSClusterFinderSDD();
29 virtual ~AliITSClusterFinderSDD(){};
31 virtual void SetCutAmplitude(Int_t mod,Double_t nsigma=4);
32 virtual Int_t CutAmplitude(Int_t anode) const {// get cut amplitude
33 return fCutAmplitude[anode];}
34 virtual void SetDAnode(Double_t danode=4.2) {// setDAnode
36 virtual Double_t DAnode() const {// get DAnode
38 virtual void SetDTime(Double_t dtime=75) {// SetDTime
40 virtual Double_t DTime() const {// get DTime
42 virtual void SetMinPeak(Int_t minpeak=10) {// SetMinPeak
44 virtual Int_t MinPeak() const {// get MinPeak
46 virtual void SetMinCharge(Int_t mincharge=30) {// SetMinCharge
47 fMinCharge=mincharge;}
48 virtual Int_t MinCharge() const {// get MinCharge
50 virtual void SetMinNCells(Int_t minc=3) {// setNCells
52 virtual Int_t MinNCells() const {// get MinNCells
54 virtual void SetMaxNCells(Int_t maxc=10) {// setNCells
56 virtual Int_t MaxNCells() const {// get MaxNCells
58 virtual void SetTimeCorr(Double_t timec=19.3) {// setNCells
60 virtual Double_t TimeCorr() const{// get Time Correction (ns)
63 // Search for clusters
64 virtual void FindRawClusters(Int_t mod=0);
65 void Find1DClusters();
66 void Find1DClustersE();
68 void SelectClusters();
69 void GetRecPoints(AliITSCalibrationSDD* cal);
70 void ResolveClusters(); // Ernesto Lopez Torres
71 Int_t SearchPeak(Double_t *spect,Int_t xdim,Int_t zdim,Int_t *peakX,
72 Int_t *peakZ,Double_t *peakAmp,Double_t minpeak);//Ernesto
73 Int_t NoLinearFit( Int_t xdim, Int_t zdim, Double_t *param, Double_t *spe,
74 Int_t *niter, Double_t *chir );
75 void Minim(Int_t xdim,Int_t zdim,Double_t *param,Double_t *prm0,
76 Double_t *steprm,Double_t *chisqr,Double_t *spe,Double_t *speFit);
77 Double_t ChiSqr(Int_t xdim,Int_t zdim,Double_t *spe,Double_t *speFit)const;
78 void PeakFunc( Int_t xdim, Int_t zdim, Double_t *par, Double_t *spe,
79 Double_t *Integral=0 );
80 void PrintStatus() const;
81 void CorrectPosition(Double_t &z, Double_t&y);
83 virtual AliITSCalibrationSDD* GetResp(Int_t mod)const{
84 return (AliITSCalibrationSDD*) fDetTypeRec->GetCalibrationModel(mod);}//Return Response
85 //Returns fSegmentation
86 virtual AliITSsegmentationSDD* GetSeg()const{
87 return (AliITSsegmentationSDD*)fDetTypeRec->GetSegmentationModel(1);}
88 AliITSClusterFinderSDD(const AliITSClusterFinderSDD &source); // copy ctor
89 AliITSClusterFinderSDD& operator=(const AliITSClusterFinderSDD &source);
91 Int_t fNclusters; //! num of clusters
92 Double_t fDAnode; //! fDanode
93 Double_t fDTime; //! fDtime
94 Double_t fTimeCorr; //! Correction factor along time coord
95 TArrayI fCutAmplitude; //! cut amplitude
96 Int_t fMinPeak; //! min peak
97 Int_t fMinCharge; //! min charge
98 Int_t fMinNCells; //! min num of cells
99 Int_t fMaxNCells; //! max num of cells
101 ClassDef(AliITSClusterFinderSDD,3) // SDD clustering - Piergiorgio C. algo