Shadowing of variabled corrected.

[u/mrichter/AliRoot.git] / TRD / AliTRDtrackerV1.h

diff --git a/TRD/AliTRDtrackerV1.h b/TRD/AliTRDtrackerV1.h
index aa7e361926a6abe272ab4007d43ddc4be56630a4..95120f9ddb2629a2995a0d92e4bb62b490cfd608 100644 (file)
--- a/TRD/AliTRDtrackerV1.h
+++ b/TRD/AliTRDtrackerV1.h
@@ -67,8 +67,10 @@ public:
   
        Int_t          Clusters2Tracks(AliESDEvent *esd);
   static TTreeSRedirector* DebugStreamer() {return fgDebugStreamer;}
+  static void SetNTimeBins(Int_t nTimeBins){fgNTimeBins = nTimeBins; }
   AliCluster*    GetCluster(Int_t index) const;
        static void    GetExtrapolationConfig(Int_t iconfig, Int_t planes[2]);
+       static const Int_t   GetNTimeBins() {return fgNTimeBins;}
        static void    GetSeedingConfig(Int_t iconfig, Int_t planes[4]);
        static TLinearFitter *GetTiltedRiemanFitter();
        static TLinearFitter *GetTiltedRiemanFitterConstraint();
@@ -88,38 +90,58 @@ public:
   
   static Int_t   Freq(Int_t n, const Int_t *inlist, Int_t *outlist, Bool_t down); // to be removed 
 
+       class AliTRDLeastSquare{
+       public:
+                                       AliTRDLeastSquare();
+                                       ~AliTRDLeastSquare(){};
+               
+               void            AddPoint(Double_t *x, Double_t y, Double_t sigmaY);
+               void            RemovePoint(Double_t *x, Double_t y, Double_t sigmaY);
+               void            Eval();
+               
+               Double_t        GetFunctionParameter(Int_t ParNumber) const {return fParams[ParNumber];}
+               Double_t        GetFunctionValue(Double_t *xpos) const;
+               void            GetCovarianceMatrix(Double_t *storage) const;
+       private:
+                                       AliTRDLeastSquare(const AliTRDLeastSquare &);
+               AliTRDLeastSquare& operator=(const AliTRDLeastSquare &);
+               Double_t        fParams[2];                                             // Fitparameter 
+               Double_t        fCovarianceMatrix[3];                   // Covariance Matrix
+               Double_t        fSums[6];                                               // Sums
+       };
+
 protected:
   Bool_t         AdjustSector(AliTRDtrackV1 *track); 
        Double_t       BuildSeedingConfigs(AliTRDtrackingChamber **stack, Int_t *configs);
-       static Float_t CalculateChi2Z(AliTRDseedV1 *tracklets, Double_t offset, Double_t slope);
+       static Float_t CalculateChi2Z(AliTRDseedV1 *tracklets, Double_t offset, Double_t slope, Double_t xref);
        Int_t          Clusters2TracksSM(Int_t sector, AliESDEvent *esd);
        Int_t          Clusters2TracksStack(AliTRDtrackingChamber **stack, TClonesArray *esdTrackList);
-       void           CookLabel(AliKalmanTrack *pt, Float_t wrong) const;
-       AliTRDcluster* FindSeedingCluster(AliTRDtrackingChamber *chamber, AliTRDseedV1 *sfit) const;
-       Int_t          GetSeedingLayers(AliTRDtrackingChamber *chamber, Double_t *params);
-       void           GetMeanCLStack(AliTRDtrackingChamber *chamber, Int_t *planes, Double_t *params);
        AliTRDseedV1*  GetTracklet(AliTRDtrackV1 *trk, Int_t plane, Int_t &idx);
        Bool_t         GetTrackPoint(Int_t index, AliTrackPoint &p) const;      
        Int_t          MakeSeeds(AliTRDtrackingChamber **stack, AliTRDseedV1 *sseed, Int_t *ipar);
        AliTRDtrackV1* MakeTrack(AliTRDseedV1 *seeds, Double_t *params);
   Int_t          PropagateToX(AliTRDtrackV1 &t, Double_t xToGo, Double_t maxStep);
        Int_t          SetTracklet(AliTRDseedV1 *tracklet);
+       Bool_t         UpdateTracklet(AliTRDseedV1 *tracklet, Int_t index);
 
 private:
        AliTRDtrackerV1(const AliTRDtrackerV1 &tracker);
        AliTRDtrackerV1 &operator=(const AliTRDtrackerV1 &tracker);
-       Double_t       CookLikelihood(AliTRDseedV1 *cseed, Int_t planes[4], Double_t *chi2);
-       Double_t       CalculateTrackLikelihood(AliTRDseedV1 *tracklets, Double_t *chi2);
-       Int_t          ImproveSeedQuality(AliTRDtrackingChamber **stack, AliTRDseedV1 *tracklet);
+       Double_t        CookLikelihood(AliTRDseedV1 *cseed, Int_t planes[4], Double_t *chi2);
+       Double_t        CalculateTrackLikelihood(AliTRDseedV1 *tracklets, Double_t *chi2);
+       Int_t           ImproveSeedQuality(AliTRDtrackingChamber **stack, AliTRDseedV1 *tracklet);
+       static Float_t  CalculateReferenceX(AliTRDseedV1 *tracklets);
+       
+       Float_t     GetChi2Y(AliTRDseedV1 *tracklets) const;
+       Float_t     GetChi2Z(AliTRDseedV1 *tracklets) const;
 
 private:
   AliTRDgeometry          *fGeom;                          // Pointer to TRD geometry
   AliTRDtrackingSector    fTrSec[kTrackingSectors];       // Array of tracking sectors;    
 
-       TClonesArray        *fClusters;                       // List of clusters
-       TClonesArray        *fTracklets;                      // List of tracklets
-       TClonesArray        *fTracks;                         // List of tracks
-  Int_t                    fTimeBinsPerPlane;              // Timebins per plane in track prolongation 
+  TClonesArray        *fClusters;                       // List of clusters
+  TClonesArray        *fTracklets;                      // List of tracklets
+  TClonesArray        *fTracks;                         // List of tracks
   
   // should go to the recoParam
   static const Double_t    fgkMaxChi2;                     // Max increment in track chi2 
@@ -129,16 +151,18 @@ private:
   static const Double_t    fgkMaxStep;                     // Maximal step for tracking  
        
        // stand alone tracking
-       static Double_t      fgTopologicQA[kNConfigs];         //  Topologic quality
+       static Double_t      fgTopologicQA[kNConfigs];        //  Topologic quality
        Double_t             fTrackQuality[kMaxTracksStack];  //  Track quality 
        Int_t                fSeedLayer[kMaxTracksStack];     //  Seed layer
        Int_t                fSieveSeeding;                   //! Seeding iterator
        
-       static TLinearFitter *fgTiltedRieman;                   //  Fitter for the tilted Rieman fit without vertex constriant
-       static TLinearFitter *fgTiltedRiemanConstrained;        //  Fitter for the tilted Rieman fit with vertex constraint     
-       static AliRieman     *fgRieman;                             //  Fitter for the untilted Rieman fit
-  static TTreeSRedirector *fgDebugStreamer;                 //!Debug streamer
-
+  static Int_t         fgNTimeBins;                     // Timebins per plane in track prolongation 
+       static TLinearFitter *fgTiltedRieman;                 //  Fitter for the tilted Rieman fit without vertex constriant
+       static TLinearFitter *fgTiltedRiemanConstrained;      //  Fitter for the tilted Rieman fit with vertex constraint       
+       static AliRieman     *fgRieman;                       //  Fitter for the untilted Rieman fit
+  
+       static TTreeSRedirector *fgDebugStreamer;             //!Debug streamer
+       
        ClassDef(AliTRDtrackerV1, 2)                          //  TRD tracker development class
 
 };