3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 ////////////////////////////////////////////////////////////////////////////
10 // The TRD track seed //
12 ////////////////////////////////////////////////////////////////////////////
15 #include "AliTRDseed.h"
19 #include "AliRieman.h"
22 class TTreeSRedirector;
26 class AliTRDstackLayer;
28 class AliTRDrecoParam;
30 class AliTRDseedV1 : public AliTRDseed
35 AliTRDseedV1(Int_t layer = -1, AliTRDrecoParam *p=0x0);
37 AliTRDseedV1(const AliTRDseedV1 &ref, Bool_t owner=kFALSE);
38 AliTRDseedV1& operator=(const AliTRDseedV1 &);
40 Bool_t AttachClustersIter(AliTRDstackLayer *layer, Float_t quality, Bool_t kZcorr = kFALSE, AliTRDcluster *c=0x0);
41 Bool_t AttachClustersProj(AliTRDstackLayer *layer, Float_t quality, Bool_t kZcorr = kFALSE, AliTRDcluster *c=0x0);
42 static Float_t FitRiemanTilt(AliTRDseedV1 * cseed, Bool_t terror);
45 // Bool_t AttachClusters(Double_t *dx, Float_t quality, Bool_t kZcorr=kFALSE, AliTRDcluster *c=0x0);
46 inline Float_t GetChi2Z(const Float_t z = 0.) const;
47 inline Float_t GetChi2Y(const Float_t y = 0.) const;
48 Float_t GetQuality(Bool_t kZcorr) const;
49 Int_t GetLayer() const { return fLayer; }
51 inline void Update(const AliRieman *rieman);
52 void Print(Option_t * /*o*/) const { }
55 void SetLayer(Int_t l) { fLayer = l; }
56 void SetNTimeBins(Int_t nTB) { fTimeBins = nTB; }
57 void SetRecoParam(AliTRDrecoParam *p) { fRecoParam = p; }
61 void Copy(TObject &ref) const;
65 Int_t fLayer; // layer for this seed
66 Int_t fTimeBins; // local copy of the DB info
67 Bool_t fOwner; // owner of the clusters
68 AliTRDrecoParam *fRecoParam; //! local copy of the reco params
70 ClassDef(AliTRDseedV1, 1) // New TRD seed
74 //____________________________________________________________
75 inline Float_t AliTRDseedV1::GetChi2Z(const Float_t z) const
77 Float_t z1 = (z == 0.) ? fMeanz : z;
78 Float_t chi = fZref[0] - z1;
82 //____________________________________________________________
83 inline Float_t AliTRDseedV1::GetChi2Y(const Float_t y) const
85 Float_t y1 = (y == 0.) ? fYfitR[0] : y;
86 Float_t chi = fYref[0] - y1;
90 //____________________________________________________________
91 inline void AliTRDseedV1::Update(const AliRieman *rieman)
93 fZref[0] = rieman->GetZat(fX0);
94 fZref[1] = rieman->GetDZat(fX0);
95 fYref[0] = rieman->GetYat(fX0);
96 fYref[1] = rieman->GetDYat(fX0);