3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 ////////////////////////////////////////////////
9 // Manager class for TPC //
10 ////////////////////////////////////////////////
15 #include <TMatrixFfwd.h>
18 class AliTPCDigitsArray;
21 class AliTPCTrackHits; // M.I. -MI4 old hits - to be removed later
22 class AliTPCTrackHitsV2; // M.I.
24 #include "AliDetector.h"
29 class AliTPC : public AliDetector {
33 AliTPC(const char *name, const char *title);
35 virtual AliLoader* MakeLoader(const char* topfoldername);
38 virtual void AddHit(Int_t a1, Int_t *a2, Float_t *a3);
39 virtual void BuildGeometry();
40 virtual void CreateGeometry() {}
41 virtual void CreateMaterials();
42 virtual void AddAlignableVolumes() const {}
44 virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* manager) const;
45 virtual void SDigits2Digits(){;} //MI change -cycling to production
46 virtual void SDigits2Digits2(Int_t eventnumber=0);
48 virtual void Hits2SDigits(); // MI - cycling around
49 virtual void Hits2SDigits2(Int_t eventnumber=0);
51 virtual void Hits2Digits();
52 virtual void Hits2Digits(Int_t eventnumber); //MI change
53 virtual void Hits2DigitsSector(Int_t isec); //MI change
55 virtual Int_t IsVersion() const =0;
56 virtual void Digits2Raw();
57 virtual Bool_t Raw2SDigits(AliRawReader* rawReader);
58 Int_t GetNsectors() const {return fNsectors;}
59 virtual void MakeBranch(Option_t *opt=" ");
60 virtual void ResetDigits();
61 virtual void SetSens(Int_t sens);
62 virtual void SetSide(Float_t side);
65 virtual void StepManager()=0;
66 virtual void DrawDetector() {}
67 AliTPCDigitsArray* GetDigitsArray() {return fDigitsArray;} //MI change
68 AliTPCParam *GetParam(){return fTPCParam;} // M.K, M.I changes
69 void SetParam(AliTPCParam *param){fTPCParam=param;} // M.K, M.I changes
70 void SetDigitsArray(AliTPCDigitsArray* param) {fDigitsArray=param;} //MI change
72 // additional function neccesary for the new hits
73 virtual void MakeBranch2(Option_t *opt=" ", const char *file=0); //
74 virtual void SetTreeAddress();
75 virtual void SetTreeAddress2();
76 virtual void AddHit2(Int_t a1, Int_t *a2, Float_t *a3); //
77 virtual void ResetHits();
78 virtual void ResetHits2();
79 virtual AliHit* FirstHit(Int_t track);
80 virtual AliHit* NextHit();
81 virtual AliHit* FirstHit2(Int_t track);
82 virtual AliHit* NextHit2();
83 virtual void LoadPoints(Int_t dummy);
84 virtual void LoadPoints2(Int_t dummy);
85 virtual void LoadPoints3(Int_t dumy);
86 virtual void FinishPrimary();
87 virtual void RemapTrackHitIDs(Int_t *map);
88 void SetHitType(Int_t type){fHitType =type;} //set type of hit container
89 void SetDigitsSwitch(Int_t sw){fDigitsSwitch = sw;}
90 void SetDefSwitch(Int_t def){fDefaults = def;}
91 Float_t GetNoise(); //get Current noise
92 void GenerNoise(Int_t tablasize); // make noise table
93 Bool_t IsSectorActive(Int_t sec) const; // check if the sector is active
94 void SetActiveSectors(Int_t * sectors, Int_t n); //set active sectors
95 Int_t GetHitType() const {return fHitType;}
96 void SetActiveSectors(Int_t flag=0); //loop over al hits and set active only hitted sectors
97 Bool_t TrackInVolume(Int_t id,Int_t track); //return true if current track is in volume
100 static AliTPCParam* LoadTPCParam(TFile *file);
102 Int_t fDefaults; // defaults switch
103 Int_t fSens; // ISENS
104 Int_t fNsectors; // Number of sectors in TPC
106 AliTPCDigitsArray * fDigitsArray; //! detector digit object
107 AliTPCParam *fTPCParam; // pointer to TPC parameters
108 AliTPCTrackHitsV2 *fTrackHits; //! hits for given track M.I.
109 // AliTPCTrackHits *fTrackHitsOld; //! hits for given track M.I. MIold -
111 Int_t fHitType; // if fNewHit = 1 old data structure if 2 new hits if 4 old MI stucture
113 Int_t fDigitsSwitch; // digits type, 0->normal, 1->summable
117 Float_t fSide; // selects left(-1), right(+1), or both(0) sides of the TPC
120 AliTPC(const AliTPC& t);
121 AliTPC &operator = (const AliTPC & param);
124 void DigitizeRow(Int_t irow,Int_t isec,TObjArray **rowTriplet);
125 Float_t GetSignal(TObjArray *p1, Int_t ntr, TMatrixF *m1,
126 TMatrixF *m2,Int_t *IndexRange);
127 void GetList (Float_t label,Int_t np,TMatrixF *m,Int_t *IndexRange,
129 void MakeSector(Int_t isec,Int_t nrows,TTree *TH,Stat_t ntracks,TObjArray **row);
130 void TransportElectron(Float_t *xyz, Int_t *index);
131 Int_t fCurrentIndex[4];// index[0] indicates coordinate system,
132 // index[1] sector number,
133 // index[2] pad row number
134 // index[3] pad row number for which signal is calculated
135 Int_t fNoiseDepth; //! noise table
136 Float_t * fNoiseTable; //! table with noise
137 Int_t fCurrentNoise; //! index of the noise in the noise table
138 Bool_t* fActiveSectors; //! bool indicating which sectors are active
140 ClassDef(AliTPC,12) // Time Projection Chamber class
144 //_____________________________________________________________________________
146 class AliTPChit : public AliHit {
148 Int_t fSector; //sector number
149 Int_t fPadRow; //Pad Row number
150 Float_t fQ ; //charge
151 Float_t fTime; //hit time
155 AliTPChit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits);
156 virtual ~AliTPChit() {}
157 void SetX(Float_t x){fX = x;}
158 void SetY(Float_t y){fY = y;}
159 void SetZ(Float_t z){fZ = z;}
161 Float_t Time() const {return fTime;}
163 ClassDef(AliTPChit,2) // Time Projection Chamber hits