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