added option AddTRD

[u/mrichter/AliRoot.git] / ITS / UPGRADE / AliITSUTrackerGlo.h

diff --git a/ITS/UPGRADE/AliITSUTrackerGlo.h b/ITS/UPGRADE/AliITSUTrackerGlo.h
index 8d22395ac1efe85a70728cf87c8c6ea5003bd43c..1e73ae23b54166fc9ba3df13091eddb6521ee36e 100644 (file)
--- a/ITS/UPGRADE/AliITSUTrackerGlo.h
+++ b/ITS/UPGRADE/AliITSUTrackerGlo.h
@@ -12,6 +12,7 @@
 #include "AliITSUTrackCond.h"
 #include "AliITSUTrackHyp.h"
 #include "AliITSUAux.h"
+#include "AliITSUMatLUT.h"
 #include <TArrayI.h>
 
 class AliITSUReconstructor;
@@ -21,7 +22,6 @@ class AliESDtrack;
 class AliITSURecoLayer;
 class TTree;
 
-
 //-------------------------------------------------------------------------
 class AliITSUTrackerGlo : public AliTracker {
 
@@ -75,7 +75,9 @@ class AliITSUTrackerGlo : public AliTracker {
   Bool_t                 GoToEntranceToLayer(AliExternalTrackParam* seed, AliITSURecoLayer* lr, Int_t dir, Bool_t check=kFALSE);
   Bool_t                 PropagateSeed(AliITSUSeed *seed, Double_t xToGo, Double_t mass, Double_t maxStep=1.0, Bool_t matCorr=kTRUE);
   Bool_t                 PropagateSeed(AliExternalTrackParam *seed, Double_t xToGo, Double_t mass, Double_t maxStep=1.0, Bool_t matCorr=kTRUE);
-  Double_t               RefitTrack(AliExternalTrackParam* trc, Double_t r, Int_t stopCond=0);
+  Double_t               RefitTrack(AliITSUTrackHyp* trc, Double_t r, Int_t& nclFit, Int_t stopCond=0);
+  Double_t               GetMaterialBudget(const double* pnt0, const double* pnt1, double& x2x0, double& rhol) const;
+
   Int_t                  GetTrackingPhase()                 const {return fTrackPhaseID;}
 
   //
@@ -122,6 +124,7 @@ class AliITSUTrackerGlo : public AliTracker {
  protected:
   AliITSUReconstructor*           fReconstructor;  // ITS global reconstructor 
   AliITSURecoDet*                 fITS;            // interface to ITS, borrowed from reconstructor
+  AliITSUMatLUT*                  fMatLUT;         // material lookup table
   AliESDtrack*                    fCurrESDtrack;   // current esd track in processing
   Int_t                           fCurrESDtrMClb;  // its eventual mc label
   Double_t                        fCurrMass;       // current track mass
@@ -154,6 +157,7 @@ class AliITSUTrackerGlo : public AliTracker {
   AliITSURecoLayer*               fCurrLayer;      //! current layer being processed  (set only when needed, not guaranteed)
   Int_t                           fTrackPhaseID;   //! tracking phase (clusters2tracks, backward, inward)
   Int_t                           fCurrPassID;     //! tracking pass (different tracking conditions)
+  Bool_t                          fUseMatLUT;      //! use material lookup table rather than TGeo
   //
   static const Double_t           fgkToler;        // tracking tolerance
   //
@@ -211,5 +215,23 @@ inline void AliITSUTrackerGlo::AddProlongationHypothesis(AliITSUSeed* seed, Int_
   fCurrHyp->AddSeed(seed,lr);  
 }
 
+//_________________________________________________________________________
+inline Double_t AliITSUTrackerGlo::GetMaterialBudget(const double* pnt0,const double* pnt1, double& x2x0, double& rhol) const
+{
+  double par[7];
+  if (fUseMatLUT && fMatLUT) {
+    double d = fMatLUT->GetMatBudget(pnt0,pnt1,par);
+    x2x0 = par[AliITSUMatLUT::kParX2X0];
+    rhol = par[AliITSUMatLUT::kParRhoL];    
+    return d;
+  }
+  else {
+    MeanMaterialBudget(pnt0,pnt1,par);
+    x2x0 = par[1];
+    rhol = par[0]*par[4];    
+    return par[4];
+  }
+}
+
 #endif