From b67517ef4bc8437d4c784d9fe28f057e61af9b72 Mon Sep 17 00:00:00 2001
From: belikov <belikov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Mon, 30 Oct 2006 13:39:16 +0000
Subject: [PATCH] A new way of representing the RICH info (K.Shileev)

---
 STEER/AliESDtrack.cxx | 42 +++++++++++++++---------------
 STEER/AliESDtrack.h   | 60 ++++++++++++++++++++-----------------------
 2 files changed, 49 insertions(+), 53 deletions(-)

diff --git a/STEER/AliESDtrack.cxx b/STEER/AliESDtrack.cxx
index 02c393f5e7c..172b0c62fcc 100644
--- a/STEER/AliESDtrack.cxx
+++ b/STEER/AliESDtrack.cxx
@@ -96,13 +96,13 @@ AliESDtrack::AliESDtrack() :
   fTOFsignal(-1),
   fTOFsignalToT(0),
   fRICHchi2(1e10),
-  fRICHncls(0),
-  fRICHindex(0),
+  fRICHqn(-1),
+  fRICHcluIdx(-1),
   fRICHsignal(-1),
-  fRICHtheta(-1),
-  fRICHphi(-1),
-  fRICHdx(-1),
-  fRICHdy(-1),
+  fRICHtrkTheta(-1),
+  fRICHtrkPhi(-1),
+  fRICHtrkX(-1),
+  fRICHtrkY(-1),
   fRICHmipX(-1),
   fRICHmipY(-1),
   fEMCALindex(kEMCALNoMatch),
@@ -175,13 +175,13 @@ AliESDtrack::AliESDtrack(const AliESDtrack& track):
   fTOFsignal(track.fTOFsignal),
   fTOFsignalToT(track.fTOFsignalToT),
   fRICHchi2(track.fRICHchi2),
-  fRICHncls(track.fRICHncls),
-  fRICHindex(track.fRICHindex),
+  fRICHqn(track.fRICHqn),
+  fRICHcluIdx(track.fRICHcluIdx),
   fRICHsignal(track.fRICHsignal),
-  fRICHtheta(track.fRICHtheta),
-  fRICHphi(track.fRICHphi),
-  fRICHdx(track.fRICHdx),
-  fRICHdy(track.fRICHdy),
+  fRICHtrkTheta(track.fRICHtrkTheta),
+  fRICHtrkPhi(track.fRICHtrkPhi),
+  fRICHtrkX(track.fRICHtrkX),
+  fRICHtrkY(track.fRICHtrkY),
   fRICHmipX(track.fRICHmipX),
   fRICHmipY(track.fRICHmipY),
   fEMCALindex(track.fEMCALindex),
@@ -325,16 +325,16 @@ void AliESDtrack::MakeMiniESDtrack(){
 
   // Reset RICH related track information
   fRICHchi2 = 0;     
-  fRICHncls = 0;     
-  fRICHindex = 0;     
+  fRICHqn = -1;     
+  fRICHcluIdx = -1;     
   fRICHsignal = 0;     
   for (Int_t i=0;i<AliPID::kSPECIES;i++) fRICHr[i] = 0;
-  fRICHtheta = 0;     
-  fRICHphi = 0;      
-  fRICHdx = 0;     
-  fRICHdy = 0;      
-  fRICHmipX = 0;
-  fRICHmipY = 0;
+  fRICHtrkTheta = -1;     
+  fRICHtrkPhi = -1;      
+  fRICHtrkX = -1;     
+  fRICHtrkY = -1;      
+  fRICHmipX = -1;
+  fRICHmipY = -1;
   fEMCALindex = kEMCALNoMatch;
 
   delete fFriendTrack; fFriendTrack = 0;
@@ -873,7 +873,7 @@ Bool_t AliESDtrack::RelateToVertex
   alpha=GetAlpha(); sn=TMath::Sin(alpha); cs=TMath::Cos(alpha);
   Double_t s2ylocvtx = cov[0]*sn*sn + cov[2]*cs*cs - 2.*cov[1]*cs*sn;
   fCdd = GetCovariance()[0] + s2ylocvtx;   // neglecting correlations
-  fCdz = GetCovariance()[1];               // between (x,y) and z    
+  fCdz = GetCovariance()[1];               // between (x,y) and z
   fCzz = GetCovariance()[2] + cov[5];      // in vertex's covariance matrix
   //*****
 
diff --git a/STEER/AliESDtrack.h b/STEER/AliESDtrack.h
index 70a8f54e84f..2ca361f0c03 100644
--- a/STEER/AliESDtrack.h
+++ b/STEER/AliESDtrack.h
@@ -195,37 +195,33 @@ public:
   Int_t   GetTOFcluster() const {return fTOFindex;}
   void    SetTOFcluster(Int_t index) {fTOFindex=index;}
   void    SetTOFCalChannel(Int_t index) {fTOFCalChannel=index;}
-  
-  void    SetRICHsignal(Double_t beta) {fRICHsignal=beta;}
+
+// RICH methodes +++++++++++++++++++++++++++++++++ (kir)
+  void    SetRICHsignal(Double_t theta) {fRICHsignal=theta;}
   Float_t GetRICHsignal() const {return fRICHsignal;}
   void    SetRICHpid(const Double_t *p);
-  void    GetRICHpid(Double_t *p) const;
+  void    GetRICHpid(Double_t *p) const;  
   void    SetRICHchi2(Double_t chi2) {fRICHchi2=chi2;}
   Float_t GetRICHchi2() const {return fRICHchi2;}
-  void    SetRICHcluster(Int_t index) {fRICHindex=index;}
-  Int_t   GetRICHcluster() const {return fRICHindex;}
-  void    SetRICHnclusters(Int_t n) {fRICHncls=n;}
-  Int_t   GetRICHnclusters() const {return fRICHncls;}
-  void    SetRICHthetaPhi(Float_t theta, Float_t phi) {
-    fRICHtheta=theta; fRICHphi=phi;
-  }
-  void    GetRICHthetaPhi(Float_t &theta, Float_t &phi) const {
-    theta=fRICHtheta; phi=fRICHphi;
+  void    SetRICHcluster(Int_t index) {fRICHcluIdx=index;}
+  Int_t   GetRICHcluster() const {return fRICHcluIdx;}
+  void    SetRICHcluIdx(Int_t ch,Int_t idx) {fRICHcluIdx=ch*1000000+idx;}
+  Int_t   GetRICHcluIdx() const {return fRICHcluIdx;}
+  void    SetRICHtrk(Float_t  x, Float_t  y, Float_t  th, Float_t  ph) {
+     fRICHtrkX=x; fRICHtrkY=y; fRICHtrkTheta=th; fRICHtrkPhi=ph;
   }
-  void    SetRICHdxdy(Float_t dx, Float_t dy) {
-    fRICHdx=dx;  fRICHdy=dy;
+  void    GetRICHtrk(Float_t &x, Float_t &y, Float_t &th, Float_t &ph) const {
+     x=fRICHtrkX; y=fRICHtrkY; th=fRICHtrkTheta; ph=fRICHtrkPhi;
   }
-  void    GetRICHdxdy(Float_t &dx, Float_t &dy) const {
-    dx=fRICHdx;  dy=fRICHdy;
+  void    SetRICHmip(Float_t  x, Float_t  y, Int_t q, Int_t nph=0) {
+     fRICHmipX=x; fRICHmipY=y; fRICHqn=100000*q+nph;
   }
-  void    SetRICHmipXY(Float_t x, Float_t y) {
-    fRICHmipX=x; fRICHmipY=y;
-  } 
-  void    GetRICHmipXY(Float_t &x, Float_t &y) const {
-    x=fRICHmipX; y=fRICHmipY;
+  void    GetRICHmip(Float_t &x,Float_t &y,Int_t &q,Int_t &nph) const {
+     x=fRICHmipX; y=fRICHmipY; q=fRICHqn/1000000; nph=fRICHqn%1000000;
   }
-  Bool_t IsRICH()  const {return fFlags&kRICHpid;}
-  
+  Bool_t  IsRICH() const {return fFlags&kRICHpid;}
+
+
   Int_t GetEMCALcluster() {return fEMCALindex;}
   void SetEMCALcluster(Int_t index) {fEMCALindex=index;}
   Bool_t IsEMCAL() const {return fFlags&kEMCALmatch;}
@@ -329,16 +325,16 @@ protected:
   Int_t   fTOFLabel[3];    // TOF label 
   Float_t fTOFInfo[10];    //! TOF informations
 
-  // HMPID related track information
+  // HMPID related track information                 (kir)
   Float_t fRICHchi2;       // chi2 in the RICH
-  Int_t   fRICHncls;       // number of photon clusters
-  Int_t   fRICHindex;      // index of the assigned MIP cluster
-  Float_t fRICHsignal;     // RICH PID signal
+  Int_t   fRICHqn;         // 1000000*QDC + number of photon clusters
+  Int_t   fRICHcluIdx;     // 1000000*chamber id + cluster idx of the assigned MIP cluster
+  Float_t fRICHsignal;     // RICH PID signal (Theta ckov, rad)
   Float_t fRICHr[AliPID::kSPECIES];// "detector response probabilities" (for the PID)
-  Float_t fRICHtheta;      // theta of the track extrapolated to the RICH
-  Float_t fRICHphi;        // phi of the track extrapolated to the RICH
-  Float_t fRICHdx;         // x of the track impact minus x of the MIP
-  Float_t fRICHdy;         // y of the track impact minus y of the MIP
+  Float_t fRICHtrkTheta;   // theta of the track extrapolated to the RICH, LORS
+  Float_t fRICHtrkPhi;     // phi of the track extrapolated to the RICH, LORS
+  Float_t fRICHtrkX;       // x of the track impact, LORS 
+  Float_t fRICHtrkY;       // y of the track impact, LORS 
   Float_t fRICHmipX;       // x of the MIP in LORS
   Float_t fRICHmipY;       // y of the MIP in LORS
   
@@ -353,7 +349,7 @@ protected:
 
   AliESDtrack & operator=(const AliESDtrack & ) {return *this;}
 
-  ClassDef(AliESDtrack,31)  //ESDtrack 
+  ClassDef(AliESDtrack,32)  //ESDtrack 
 };
 
 #endif 
-- 
2.43.0