Line data    Source code

       1             : /**************************************************************************
       2             :  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
       3             :  *                                                                        *
       4             :  * Author: The ALICE Off-line Project.                                    *
       5             :  * Contributors are mentioned in the code where appropriate.              *
       6             :  *                                                                        *
       7             :  * Permission to use, copy, modify and distribute this software and its   *
       8             :  * documentation strictly for non-commercial purposes is hereby granted   *
       9             :  * without fee, provided that the above copyright notice appears in all   *
      10             :  * copies and that both the copyright notice and this permission notice   *
      11             :  * appear in the supporting documentation. The authors make no claims     *
      12             :  * about the suitability of this software for any purpose. It is          *
      13             :  * provided "as is" without express or implied warranty.                  *
      14             :  **************************************************************************/
      15             : 
      16             : //-----------------------------------------------------------------
      17             : //           Implementation of the ESD HLT track class
      18             : //   ESD = Event Summary Data
      19             : //   HLT = High Level Trigger
      20             : //   This is the class to deal with during the phisical analysis of data
      21             : //-----------------------------------------------------------------
      22             : 
      23             : #include "TMath.h"
      24             : #include "AliESDHLTtrack.h"
      25             : #include "AliKalmanTrack.h"
      26             : 
      27         172 : ClassImp(AliESDHLTtrack)
      28             : 
      29             : AliESDHLTtrack::AliESDHLTtrack() :
      30           0 :   TObject(),
      31           0 :   fNHits(0),
      32           0 :   fMCid(0),
      33           0 :   fWeight(0),
      34           0 :   fFromMainVertex(kFALSE),
      35           0 :   fSector(0),
      36           0 :   fQ(0),
      37           0 :   fTanl(0),
      38           0 :   fPsi(0),
      39           0 :   fPt(0),
      40           0 :   fPterr(0),
      41           0 :   fPsierr(0),
      42           0 :   fTanlerr(0),
      43           0 :   fBinX(0),
      44           0 :   fBinY(0),
      45           0 :   fSizeX(0),
      46           0 :   fSizeY(0),
      47           0 :   fPID(0)
      48           0 : {
      49             :   // Default constructor
      50           0 :   fRowRange[0] = fRowRange[1] = 0;
      51           0 :   fFirstPoint[0] = fFirstPoint[1] = fFirstPoint[2] = 0;
      52           0 :   fLastPoint[0] = fLastPoint[1] = fLastPoint[2] = 0;
      53           0 : }
      54             : 
      55             : Double_t AliESDHLTtrack::GetP() const
      56             : {
      57             :   // Returns total momentum.  
      58           0 :   return TMath::Abs(GetPt())*sqrt(1. + GetTgl()*GetTgl());
      59             : }
      60             : 
      61             : Double_t AliESDHLTtrack::GetPseudoRapidity() const
      62             : {
      63             :   // Calculates the pseudorapidity
      64           0 :   return 0.5 * TMath::Log((GetP() + GetPz()) / (GetP() - GetPz()));
      65             : }
      66             : 
      67             : Bool_t AliESDHLTtrack::UpdateTrackParams(const AliKalmanTrack *t)
      68             : {
      69             :   // Updates the track parameters
      70             : 
      71           0 :   fNHits = t->GetNumberOfClusters();
      72           0 :   fMCid = t->GetLabel();
      73             : 
      74           0 :   Double_t alpha = t->GetAlpha();
      75           0 :   fSector = (UShort_t)(alpha/(2*TMath::Pi()/18));
      76           0 :   Double_t x,p[5]; t->GetExternalParameters(x,p);
      77           0 :   if(p[4]<=0)
      78           0 :     fQ = 1;
      79             :   else
      80           0 :     fQ = -1;
      81           0 :   fPt = TMath::Abs(1./p[4]);
      82           0 :   fTanl = p[3];
      83           0 :   fPsi = alpha + TMath::ASin(p[2]);
      84             : 
      85           0 :   fFirstPoint[0] = x*TMath::Cos(alpha) - p[0]*TMath::Sin(alpha);
      86           0 :   fFirstPoint[1] = x*TMath::Sin(alpha) + p[0]*TMath::Cos(alpha);
      87           0 :   fFirstPoint[2] = p[1];
      88             : 
      89           0 :   return kTRUE;
      90           0 : }