/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
-
-#include <TArrayS.h>
-#include <TRandom.h>
-#include <TMath.h>
-
+//#include <TMath.h>
#include "AliITSClusterFinder.h"
+//#include "AliITSsegmentationSSD.h"
+#include "AliITSDetTypeRec.h"
+class TArrayI;
class AliITSclusterSSD;
class AliITSpackageSSD;
+class AliITSsegmentation;
+class AliITSsegmentationSSD;
+class AliITSresponse;
+class AliITSresponseSSD;
+class AliITSCalibration;
+class AliITSCalibrationSSD;
+class AliITSClusterFinderSSD: public AliITSClusterFinder{
-class AliITSClusterFinderSSD: public AliITSClusterFinder
-{
-
-public:
-
- AliITSClusterFinderSSD(AliITSsegmentation *seg, TClonesArray *digits, TClonesArray *recp);
-
- virtual ~AliITSClusterFinderSSD();
-
-
- void FindRawClusters();
-
-
- void SetAlpha1(Float_t a) {falpha1 =a;}
- void SetAlpha2(Float_t a) {falpha2 =a;}
- void SetAlpha3(Float_t a) {falpha3 =a;}
-
-
+ public:
+ AliITSClusterFinderSSD();
+ AliITSClusterFinderSSD(AliITSDetTypeRec* dettyp);
+ AliITSClusterFinderSSD(AliITSDetTypeRec* dettyp, TClonesArray *digits);
+ // AliITSClusterFinderSSD(AliITSsegmentation *seg, TClonesArray *digits);
+ virtual ~AliITSClusterFinderSSD();
+ void FindRawClusters(Int_t module);
+
protected:
-
+
+ AliITSClusterFinderSSD(const AliITSClusterFinderSSD &source); // Copy constructor
+ AliITSClusterFinderSSD& operator=(const AliITSClusterFinderSSD &source);// = operator
+
+ virtual AliITSsegmentationSSD* GetSeg()const{
+ return (AliITSsegmentationSSD*)fDetTypeRec->GetSegmentationModel(2);}
void InitReconstruction();
Bool_t CreateNewRecPoint(Float_t P, Float_t dP, Float_t N, Float_t dN,
- Float_t Sig,Float_t dSig,
- AliITSclusterSSD *clusterP, AliITSclusterSSD *clusterN,
- Stat_t prob);
- Bool_t CreateNewRecPoint(AliITSclusterSSD *clusterP, AliITSclusterSSD *clusterN, Stat_t prob);
-
+ Float_t Sig,Float_t dSig,
+ AliITSclusterSSD *clusterP,
+ AliITSclusterSSD *clusterN,Stat_t prob);
AliITSclusterSSD* GetPSideCluster(Int_t idx);
AliITSclusterSSD* GetNSideCluster(Int_t idx);
AliITSclusterSSD* GetCluster(Int_t idx, Bool_t side);
-
-
- void FindNeighbouringDigits();
+ void FindNeighbouringDigits(Int_t module);
void SeparateOverlappedClusters();
- void SplitCluster(TArrayI *list, Int_t nsplits, Int_t index, Bool_t side);
+ void SplitCluster(TArrayI *list,Int_t nsplits,Int_t indx,Bool_t side);
Int_t SortDigitsP(Int_t start, Int_t end);
Int_t SortDigitsN(Int_t start, Int_t end);
void FillDigitsIndex();
void SortDigits();
- void FillClIndexArrays(Int_t* arrayP, Int_t *arrayN);
+ void FillClIndexArrays(Int_t* arrayP, Int_t *arrayN) const;
void SortClusters(Int_t* arrayP, Int_t *arrayN);
Int_t SortClustersP(Int_t start, Int_t end,Int_t *array);
Int_t SortClustersN(Int_t start, Int_t end,Int_t *array);
- void ConsumeClusters();
void ClustersToPackages();
- void PackagesToPoints();
- void ReconstructNotConsumedClusters();
- Bool_t Strip2Local( Float_t stripP, Float_t stripN, Float_t &Z,Float_t &X);
- Float_t GetClusterZ(AliITSclusterSSD* clust);
- Bool_t IsCrossing(AliITSclusterSSD* p, AliITSclusterSSD* n);
- //returns index of best combination in "comb"
- Int_t GetBestComb(Int_t** comb,Int_t Ncomb, Int_t Ncl, AliITSpackageSSD * pkg);
-
- //get point that have best signal ratio
- void GetBestMatchingPoint(Int_t & ip, Int_t & in, AliITSpackageSSD* pkg );
+ Int_t GetDiff(Float_t */*retx*/, Float_t */*rety*/) const {return 0;}
+ void CalcStepFactor(Float_t Psteo, Float_t Nsteo );
+ Bool_t GetCrossing(Float_t &x, Float_t &z); //x, y of strips crossing
+ void GetCrossingError(Float_t& dp, Float_t& dn);//x, y of strips crossing err.
+ virtual AliITSCalibrationSSD* GetResp(Int_t mod)const{
+ return (AliITSCalibrationSSD*) fDetTypeRec->GetCalibrationModel(mod);}//Return Response
- //calculates Distance To Perfect Matching Line
- Float_t DistToPML(Float_t psig, Float_t nsig){ return (TMath::Abs( (7.0*nsig - 8.0*psig )/10.630146) );}
+ // Data memebers
+ // AliITS *fITS; //!Pointer to AliITS object
+ TClonesArray *fClusterP; //!
+ Int_t fNClusterP; //!Number of P side clusters in the array
+ TClonesArray *fClusterN; //!Number of N side clusters in the array
+ Int_t fNClusterN; //!
+ TClonesArray *fPackages; //!packages
+ Int_t fNPackages; //!
+ TArrayI *fDigitsIndexP; //!Digits on P side
+ Int_t fNDigitsP; //!Number of Digits on P side
+ TArrayI *fDigitsIndexN; //!Digits on N side
+ Int_t fNDigitsN; //!Number of Digits on N side
+ Float_t fPitch; //!Strip pitch
+ Float_t fTanP; //!Pside stereo angle tangent
+ Float_t fTanN; //!Nside stereo angle tangent
- Int_t GetDiff(Float_t *retx, Float_t *rety) {return 0;}
+ //Float_t **fNoise;
- void CalcStepFactor(Float_t Psteo, Float_t Nsteo );
+ /*************************************************/
+ /** parameters for reconstruction ****/
+ /** to be tune when slow simulation raliable ****/
+ /*************************************************/
-/*************************************************/
-/** methods for resolving packages ****/
-/*************************************************/
-//names may not be meaningful for all, see implementations for descriptions
-
- void ResolveSimplePackage(AliITSpackageSSD *pkg);
- void ResolvePackageWithOnePSideCluster(AliITSpackageSSD *pkg);
- void ResolvePackageWithOneNSideCluster(AliITSpackageSSD *pkg);
- void ResolveTwoForTwoPackage(AliITSpackageSSD *pkg);
-
- void ResolveClusterWithOneCross(AliITSpackageSSD *pkg,
- Int_t clusterIndex, Bool_t clusterSide);
-
- void ResolvePClusterWithOneCross(AliITSpackageSSD *pkg, Int_t clusterIndex);
- void ResolveNClusterWithOneCross(AliITSpackageSSD *pkg, Int_t clusterIndex);
- Bool_t ResolvePackageBestCombin(AliITSpackageSSD *pkg);
- void ResolveOneBestMatchingPoint(AliITSpackageSSD *pkg);
-
- Bool_t GetCrossing(Float_t &x, Float_t &z); //x, y of strips crossing
- void GetCrossingError(Float_t&, Float_t&); //x, y of strips crossing errors
-
- // Data memebers
-
- AliITS *fITS; //Pointer to AliITS object
- TClonesArray *fDigits; //Pointer to TClonesArray of digits
-
- TClonesArray *fRecPoints; //Pointer to TClonesArray of rec points
+ //Float_t fAlpha1; //!
+ //Float_t fAlpha2; //!
+ //Float_t fAlpha3; //!
+ Float_t fPNsignalRatio; //!
+
+ static const Bool_t fgkSIDEP; //!
+ static const Bool_t fgkSIDEN; //!
+
+ static const Int_t fgkNoiseThreshold; // at least one strip in the cluster has to have a signal > fgkNoiseThresold*noise
+
+ Int_t fSFF; //!forward stepping factor
+ Int_t fSFB; //!backward stepping factor
-
- TClonesArray *fClusterP; //
- Int_t fNClusterP; //Number of P side clusters in the array
-
- TClonesArray *fClusterN; //Number of N side clusters in the array
- Int_t fNClusterN;
+ ClassDef(AliITSClusterFinderSSD, 1) //Class for clustering and reconstruction of space points in SSDs
- TClonesArray *fPackages; //packages
- Int_t fNPackages;
+ };
- TArrayI *fDigitsIndexP; //Digits on P side
- Int_t fNDigitsP; //Number of Digits on P side
-
- TArrayI *fDigitsIndexN; //Digits on N side
- Int_t fNDigitsN; //Number of Digits on N side
-
-
- Float_t fPitch; //Strip pitch
- Float_t fTanP; //Pside stereo angle tangent
- Float_t fTanN; //Nside stereo angle tangent
-
-/*************************************************/
-/** parameters for reconstruction ****/
-/** to be tune when slow simulation raliable ****/
-/*************************************************/
-
- Float_t falpha1;
- Float_t falpha2;
- Float_t falpha3;
- Float_t fPNsignalRatio;
-
- static const Bool_t fgkSIDEP;
- static const Bool_t fgkSIDEN;
-
- Int_t fSFF; //forward stepping factor
- Int_t fSFB; //backward stepping factor
-
-public:
- ClassDef(AliITSClusterFinderSSD, 1) //Class for clustering and reconstruction of space points in SSDs
-
-};
-
-
#endif
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff