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