fortran loop removed

[u/mrichter/AliRoot.git] / TPC / AliTPC.h

diff --git a/TPC/AliTPC.h b/TPC/AliTPC.h
index 5bf32c7a9dad30bcc19509d6b34b789d2e11549a..03f6e10a711a0d6b529657e6974d5ccf648ec812 100644 (file)
--- a/TPC/AliTPC.h
+++ b/TPC/AliTPC.h
@@ -11,16 +11,24 @@
 #include "AliDetector.h"
 #include "AliHit.h" 
 #include "AliDigit.h" 
+#include "TLorentzVector.h" 
+
+#include "AliTPCLoader.h"
 
 
 class TMatrix;
+class AliTPCFastMatrix;  //MI change
+class AliTPCFastVector;  //MI change
+
 class TTree;
 
 class TFile;
 class AliTPCParam;
 class AliTPCDigitsArray;
 class AliTPCClustersArray;
-class AliTPCTrackHits; // M.I.
+class AliTPCTrackHitsV2; // M.I.
+class AliTPCTrackHits; // M.I.  -MI4 old hits - to be removed later
+
 
 class AliTPC : public AliDetector {
 protected:
@@ -35,9 +43,11 @@ protected:
   AliTPCDigitsArray * fDigitsArray;              //!detector digit object  
   AliTPCClustersArray * fClustersArray; //!detector cluster object
   AliTPCParam *fTPCParam;           // pointer to TPC parameters 
-  AliTPCTrackHits *fTrackHits;      //!hits for given track M.I.
-  Int_t  fHitType; // if fNewHit = 1 old data structure if 2 new hits
-  //  3 both types  
+  AliTPCTrackHitsV2 *fTrackHits;      //!hits for given track M.I.
+  AliTPCTrackHits *fTrackHitsOld;      //!hits for given track M.I. MIold -
+
+  Int_t  fHitType; // if fNewHit = 1 old data structure if 2 new hits  if 4  old MI stucture
+  //  3 both types 
   Int_t fDigitsSwitch; // digits type, 0->normal, 1->summable
 
   //MK changes
@@ -50,13 +60,16 @@ protected:
 public:
   AliTPC();
   AliTPC(const char *name, const char *title);
+  
+  virtual AliLoader* MakeLoader(const char* topfoldername);
+  
   virtual      ~AliTPC();
   virtual void  AddHit(Int_t a1, Int_t *a2, Float_t *a3);
   Int_t         DistancetoPrimitive(Int_t px, Int_t py);
   virtual void  BuildGeometry();
   virtual void  CreateGeometry() {}
   virtual void  CreateMaterials();
-  virtual void  Hits2Clusters(TFile *of, Int_t eventn=0);
+  virtual void  Hits2Clusters(Int_t eventn=0);
   virtual void  Hits2ExactClustersSector(Int_t isec); // MI change calculate "exact" cluster position
   
   virtual void  SDigits2Digits(){;} //MI change -cycling to production
@@ -70,17 +83,17 @@ public:
   virtual void  Hits2DigitsSector(Int_t isec);  //MI change
   virtual void  Init();
   virtual Int_t IsVersion() const =0;
-  virtual void  Digits2Clusters(TFile *of, Int_t eventnumber=0);
-  virtual void  Clusters2Tracks(TFile *of);
+  virtual void  Digits2Clusters(Int_t eventnumber=0);
+  virtual void  Clusters2Tracks();
 
   Int_t         GetNsectors()       {return fNsectors;}
-  virtual void  MakeBranch(Option_t *opt=" ", const char *file=0 );
+  virtual void  MakeBranch(Option_t *opt=" ");
   virtual void  ResetDigits();
   virtual void  SetSecAL(Int_t sec);
   virtual void  SetSecAU(Int_t sec);
   virtual void  SetSecLows(Int_t s1,Int_t s2,Int_t s3,Int_t s4,Int_t s5, Int_t s6);
   virtual void  SetSecUps (Int_t s1,Int_t s2,Int_t s3,Int_t s4,Int_t s5, Int_t s6,
-                          Int_t s7,Int_t s8,Int_t s9,Int_t s10, Int_t s11, Int_t s12);
+         Int_t s7,Int_t s8,Int_t s9,Int_t s10, Int_t s11, Int_t s12);
   virtual void  SetSens(Int_t sens);
 
 
@@ -108,7 +121,6 @@ public:
    virtual AliHit* NextHit();
    virtual AliHit* FirstHit2(Int_t track);
    virtual AliHit* NextHit2();
-   
    virtual void LoadPoints(Int_t dummy);
    virtual void LoadPoints2(Int_t dummy);
    virtual void LoadPoints3(Int_t dumy);
@@ -120,15 +132,24 @@ public:
    void SetHitType(Int_t type){fHitType =type;} //set type of hit container
    void SetDigitsSwitch(Int_t sw){fDigitsSwitch = sw;}
    void SetDefSwitch(Int_t def){fDefaults = def;}
-   virtual void  Merge(TTree * intree, Int_t *mask, Int_t nin, Int_t outid);
-
-private:
+   Float_t GetNoise();  //get Current noise  
+   void    GenerNoise(Int_t tablasize);  // make noise table
+   Bool_t  IsSectorActive(Int_t sec);    // check if the sector is active
+   void    SetActiveSectors(Int_t * sectors, Int_t n);  //set active sectors
+   Int_t GetHitType(){return fHitType;}
+   void    SetActiveSectors(Int_t flag=0); //loop over al hits and set active only hitted sectors
+
+// static functions
+   static AliTPCParam* LoadTPCParam(TFile *file);  
+   
+ private:
   //
+   Bool_t  TrackInVolume(Int_t id,Int_t track);  //return true if current track is in volume
   void SetDefaults();
   void DigitizeRow(Int_t irow,Int_t isec,TObjArray **rowTriplet);
-  Float_t GetSignal(TObjArray *p1, Int_t ntr, TMatrix *m1, TMatrix *m2,
-                    Int_t *IndexRange);
-  void GetList (Float_t label,Int_t np,TMatrix *m,Int_t *IndexRange,
+  Float_t GetSignal(TObjArray *p1, Int_t ntr, AliTPCFastMatrix *m1, 
+                   AliTPCFastMatrix *m2,Int_t *IndexRange);
+  void GetList (Float_t label,Int_t np,AliTPCFastMatrix *m,Int_t *IndexRange,
                 Float_t **pList);
   void MakeSector(Int_t isec,Int_t nrows,TTree *TH,Stat_t ntracks,TObjArray **row);
   void TransportElectron(Float_t *xyz, Int_t *index);
@@ -136,8 +157,12 @@ private:
                          // index[1] sector number, 
                          // index[2] pad row number  
                          // index[3] pad row number for which signal is calculated
-  
-  ClassDef(AliTPC,4)  // Time Projection Chamber class
+  Int_t      fNoiseDepth;  //!noise table
+  Float_t *  fNoiseTable;  //![fNoiseDepth] table with noise
+  Int_t      fCurrentNoise; //!index of the noise in  the noise table 
+  Bool_t*    fActiveSectors; //!bool indicating which sectors are active
+
+  ClassDef(AliTPC,10)  // Time Projection Chamber class
 };