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 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 ResetDigits();
61 virtual void SetSens(Int_t sens);
62 virtual void SetSide(Float_t side);
65 virtual void StepManager()=0;
66 virtual void DrawDetector() const {}
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 MakeBranch(Option_t *opt=" "); //
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 FinishPrimary();
84 virtual void RemapTrackHitIDs(Int_t *map);
85 void SetHitType(Int_t type){fHitType =type;} //set type of hit container
86 void SetDigitsSwitch(Int_t sw){fDigitsSwitch = sw;}
87 void SetDefSwitch(Int_t def){fDefaults = def;}
88 Float_t GetNoise(); //get Current noise
89 void GenerNoise(Int_t tablasize); // make noise table
90 Bool_t IsSectorActive(Int_t sec) const; // check if the sector is active
91 void SetActiveSectors(Int_t * sectors, Int_t n); //set active sectors
92 Int_t GetHitType() const {return fHitType;}
93 void SetActiveSectors(Int_t flag=1); //loop over al hits and set active only hitted sectors
94 Bool_t TrackInVolume(Int_t id,Int_t track); //return true if current track is in volume
95 void SetPrimaryIonisation(Bool_t flag = kTRUE) {fPrimaryIonisation = flag;}
96 void SetGainFactor(Float_t gain){fGainFactor=gain;} //gas gain scaling factor
97 Float_t GetGainFactor()const {return fGainFactor;}//gas gain scaling factor
99 static AliTPCParam* LoadTPCParam(TFile *file);
101 Int_t fDefaults; // defaults switch
102 Int_t fSens; // ISENS
103 Int_t fNsectors; // Number of sectors in TPC
105 AliTPCDigitsArray * fDigitsArray; //! detector digit object
106 AliTPCParam *fTPCParam; // pointer to TPC parameters
107 AliTPCTrackHitsV2 *fTrackHits; //! hits for given track M.I.
108 // AliTPCTrackHits *fTrackHitsOld; //! hits for given track M.I. MIold -
110 Int_t fHitType; // if fNewHit = 1 old data structure if 2 new hits if 4 old MI stucture
112 Int_t fDigitsSwitch; // digits type, 0->normal, 1->summable
116 Float_t fSide; // selects left(-1), right(+1), or both(0) sides of the TPC
117 Bool_t fPrimaryIonisation; //switch between Fluka(true) and geant3(false)
119 AliTPC(const AliTPC& t);
120 AliTPC &operator = (const AliTPC & param);
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
138 Float_t fGainFactor; // scaling factor
139 TTreeSRedirector *fDebugStreamer; //!debug streamer
140 ClassDef(AliTPC,13) // 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,3) // Time Projection Chamber hits