]>
Commit | Line | Data |
---|---|---|
a38a7850 | 1 | // @(#) $Id$ |
4aa41877 | 2 | // Original: AliHLTClustFinderNew.h,v 1.13 2004/06/18 10:55:26 loizides |
a38a7850 | 3 | |
01f43166 | 4 | #ifndef AliHLTTPC_CLUSTERFINDER |
5 | #define AliHLTTPC_CLUSTERFINDER | |
6 | /* This file is property of and copyright by the ALICE HLT Project * | |
7 | * ALICE Experiment at CERN, All rights reserved. * | |
8 | * See cxx source for full Copyright notice */ | |
a38a7850 | 9 | |
01f43166 | 10 | /** @file AliHLTTPCClusterFinder.h |
11 | @author Anders Vestbo, Constantin Loizides, Jochen Thaeder | |
12 | Kenneth Aamodt kenneth.aamodt@student.uib.no | |
13 | @date | |
14 | @brief Cluster Finder for the TPC | |
15 | */ | |
16 | ||
17 | // see below for class documentation | |
18 | // or | |
19 | // refer to README to build package | |
20 | // or | |
21 | // visit http://web.ift.uib.no/~kjeks/doc/alice-hlt | |
46b33a24 | 22 | |
d8b2f74a | 23 | |
01f43166 | 24 | #include "AliHLTLogging.h" |
25 | #include "AliHLTTPCPadArray.h" | |
a38a7850 | 26 | class AliHLTTPCSpacePointData; |
27 | class AliHLTTPCDigitReader; | |
28 | ||
46b33a24 | 29 | class AliHLTTPCClusterFinder : public AliHLTLogging { |
a38a7850 | 30 | |
31 | public: | |
32 | struct AliClusterData | |
33 | { | |
34 | UInt_t fTotalCharge; //tot charge of cluster | |
35 | UInt_t fPad; //pad value | |
36 | UInt_t fTime; //time value | |
37 | ULong64_t fPad2; //for error in XY direction | |
38 | ULong64_t fTime2; //for error in Z direction | |
39 | UInt_t fMean; //mean in time | |
40 | UInt_t fFlags; //different flags | |
41 | UInt_t fChargeFalling; //for deconvolution | |
42 | UInt_t fLastCharge; //for deconvolution | |
43 | UInt_t fLastMergedPad; //dont merge twice per pad | |
d8b2f74a | 44 | Int_t fRow; //row value |
a38a7850 | 45 | }; |
46 | typedef struct AliClusterData AliClusterData; //! | |
47 | ||
8f8bf0af | 48 | /** standard constructor */ |
49 | AliHLTTPCClusterFinder(); | |
50 | /** destructor */ | |
51 | virtual ~AliHLTTPCClusterFinder(); | |
52 | ||
53 | void Read(void* ptr,unsigned long size); | |
54 | ||
55 | void InitSlice(Int_t slice,Int_t patch,Int_t firstrow, Int_t lastrow,Int_t maxpoints); | |
56 | void InitSlice(Int_t slice,Int_t patch,Int_t maxpoints); | |
57 | void ProcessDigits(); | |
58 | ||
59 | void SetOutputArray(AliHLTTPCSpacePointData *pt); | |
60 | void WriteClusters(Int_t n_clusters,AliClusterData *list); | |
61 | ||
62 | void SetXYError(Float_t f) {fXYErr=f;} | |
63 | void SetZError(Float_t f) {fZErr=f;} | |
64 | void SetDeconv(Bool_t f) {fDeconvPad=f; fDeconvTime=f;} | |
65 | void SetThreshold(UInt_t i) {fThreshold=i;} | |
66 | void SetOccupancyLimit(Float_t f) {fOccupancyLimit=f;} | |
67 | void SetSignalThreshold(Int_t i) {fSignalThreshold=i;} | |
b1c46961 | 68 | void SetNSigmaThreshold(Double_t d) {fNSigmaThreshold=d;} |
8f8bf0af | 69 | void SetMatchWidth(UInt_t i) {fMatch=i;} |
70 | void SetSTDOutput(Bool_t f=kFALSE) {fStdout=f;} | |
71 | void SetCalcErr(Bool_t f=kTRUE) {fCalcerr=f;} | |
72 | void SetRawSP(Bool_t f=kFALSE) {fRawSP=f;} | |
73 | void SetReader(AliHLTTPCDigitReader* f){fDigitReader = f;} | |
74 | Int_t GetNumberOfClusters() const {return fNClusters;} | |
75 | ||
76 | //----------------------------------Methods for the new unsorted way of reading data ---------- | |
77 | void SetPadArray(AliHLTTPCPadArray *padArray); | |
78 | void ReadDataUnsorted(void* ptr,unsigned long size); | |
79 | void FindClusters(); | |
b1c46961 | 80 | Int_t GetActivePads(AliHLTTPCPadArray::AliHLTTPCActivePads* activePads,Int_t maxActivePads); |
8f8bf0af | 81 | void WriteClusters(Int_t nclusters,AliHLTTPCClusters *list); |
aff6e981 | 82 | void SetUnsorted(Int_t unsorted){fUnsorted=unsorted;} |
83 | ||
a38a7850 | 84 | private: |
8f8bf0af | 85 | /** copy constructor prohibited */ |
86 | AliHLTTPCClusterFinder(const AliHLTTPCClusterFinder&); | |
87 | /** assignment operator prohibited */ | |
88 | AliHLTTPCClusterFinder& operator=(const AliHLTTPCClusterFinder&); | |
89 | ||
2a083ac4 | 90 | AliHLTTPCSpacePointData *fSpacePointData; //! array of space points |
91 | AliHLTTPCDigitReader *fDigitReader; //! reader instance | |
a38a7850 | 92 | |
2a083ac4 | 93 | UChar_t* fPtr; //! pointer to packed block |
a38a7850 | 94 | unsigned long fSize; //packed block size |
95 | Bool_t fDeconvTime; //deconv in time direction | |
96 | Bool_t fDeconvPad; //deconv in pad direction | |
97 | Bool_t fStdout; //have print out in write clusters | |
98 | Bool_t fCalcerr; //calculate centroid sigmas | |
99 | Bool_t fRawSP; //store centroids in raw system | |
100 | ||
101 | ||
102 | Int_t fFirstRow; //first row | |
103 | Int_t fLastRow; //last row | |
104 | Int_t fCurrentRow; //current active row | |
105 | Int_t fCurrentSlice; //current slice | |
106 | Int_t fCurrentPatch; //current patch | |
107 | Int_t fMatch; //size of match | |
108 | UInt_t fThreshold; //threshold for clusters | |
84645eb0 | 109 | /** threshold for zero suppression (applied per bin) */ |
b1c46961 | 110 | Int_t fSignalThreshold;// see above |
111 | /** threshold for zero suppression 2007 December run */ | |
112 | Double_t fNSigmaThreshold; // see above | |
a38a7850 | 113 | Int_t fNClusters; //number of found clusters |
114 | Int_t fMaxNClusters; //max. number of clusters | |
115 | Float_t fXYErr; //fixed error in XY | |
116 | Float_t fZErr; //fixed error in Z | |
5235c3e9 | 117 | |
118 | Float_t fOccupancyLimit; // Occupancy Limit | |
119 | ||
01f43166 | 120 | AliHLTTPCPadArray * fPadArray; //! transient |
121 | ||
aff6e981 | 122 | Int_t fUnsorted; // enable for processing of unsorted digit data |
a38a7850 | 123 | |
64defa03 | 124 | /** list of active pads if PadArray is not used */ |
125 | vector<AliHLTTPCPadArray::AliHLTTPCActivePads> fActivePads; //!transient | |
126 | ||
a38a7850 | 127 | #ifdef do_mc |
128 | void GetTrackID(Int_t pad,Int_t time,Int_t *trackID); | |
129 | #endif | |
130 | ||
64defa03 | 131 | ClassDef(AliHLTTPCClusterFinder,3) //Fast cluster finder |
a38a7850 | 132 | }; |
133 | #endif |