/**************************************************************************
 * Copyright(c) 1998-2003, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
////////////////////////////////////////////////////
//  Stand alone track class                       //
//  Origin:  Elisabetta Crescio                   //
//  e-mail:  crescio@to.infn.it                   //
//  it is a V2 track with a possible number       //
//  of cluster equal to fgkMaxNumberOfClusters    //
////////////////////////////////////////////////////

#include "AliRun.h"
#include "AliITS.h"
#include "AliITSgeom.h"
#include "AliITStrackSA.h"


ClassImp(AliITStrackSA)

//_____________________________________
AliITStrackSA:: AliITStrackSA() : AliITStrackV2(){
// Default constructor  
  SetNumberOfClusters(0);
  SetNumberOfClustersSA(0);
  ResetIndexSA();
}


//___________________________________________________
AliITStrackSA::AliITStrackSA(const AliITStrackV2& t) : 
AliITStrackV2(t){
// Copy a V2 track into a SA track
  SetNumberOfClustersSA(0);
  ResetIndexSA();
}
//___________________________________________________
AliITStrackSA::AliITStrackSA(const AliITStrackSA& t) : 
AliITStrackV2(t){
// Copy constructor

  ResetIndexSA();
  Int_t number = t.GetNumberOfClustersSA();
  SetNumberOfClustersSA(number);
  for(Int_t i=0;i<number;i++){
    fSain[i]=t.fSain[i];
  }

}
//____________________________________________________
AliITStrackSA::AliITStrackSA(Int_t layer, Int_t ladder, Int_t detector, Double_t Ycoor, Double_t Zcoor, Double_t phi, Double_t tanlambda, Double_t curv, Int_t lab ) {
  // standard constructor. Used for ITS standalone tracking
  AliITS* iTS = (AliITS*)gAlice->GetDetector("ITS");
  if(!iTS){
    Fatal("AliITStrackSA","ITS object not found - Abort\n");
    return;
  }
  AliITSgeom *geom = iTS->GetITSgeom();
  if(!geom){
    Fatal("AliITStrackSA","ITS geometry not found - Abort\n");
    return;
  }
  // get the azimuthal angle of the detector containing the innermost
  // cluster of this track (data member fAlpha)
  Float_t rotmatr[9];
  geom->GetRotMatrix(layer,ladder,detector,rotmatr);
  fAlpha=TMath::ATan2(rotmatr[1],rotmatr[0])+TMath::Pi();
  fAlpha+=TMath::Pi()/2.;
  if(layer==1) fAlpha+=TMath::Pi();


  // get the radius of this detector. Procedure taken from the 
  // AliITStrackerV2 constructor
  Float_t x=0,y=0,z=0;
  geom->GetTrans(layer,ladder,detector,x,y,z);

  Double_t fi=TMath::ATan2(rotmatr[1],rotmatr[0])+TMath::Pi();
  fi+=TMath::Pi()/2;
  if (layer==1) fi+=TMath::Pi();
  Double_t cp=TMath::Cos(fi), sp=TMath::Sin(fi);
  fX=x*cp+y*sp;


  fdEdx = 0;

  fC00 = 0.000009; // 0.000009
  fC10 = 0.;
  fC11 = 0.000003; //0.000030
  fC20 = 0.;
  fC21 = 0.;
  fC22 = 0.000001; //0.000001
  fC30 = 0.;
  fC31 = 0.;
  fC32 = 0.;
  fC33 = 0.000002; //0.000002
  fC40 = 0.;
  fC41 = 0.;
  fC42 = 0.;
  fC43 = 0.;
  fC44 = 0.000001; //0.0000001

  fP0 = Ycoor;
  fP1 = Zcoor;
  
  fP2 = TMath::Sin(phi-fAlpha);
  fP3 = tanlambda;
  fP4 = curv;
  for(Int_t i=0; i<kMaxLayer; i++) fIndex[i] = 0;  // to be set explicitely

  for(Int_t i=0; i<4; i++) fdEdxSample[i] = 0; 

  SetNumberOfClusters(0);
  SetNumberOfClustersSA(0);
  ResetIndexSA();
  SetChi2(0);
  SetMass(0.139);    // pion mass
  SetLabel(lab); 
  
}

//____________________________________________________________
void AliITStrackSA::AddClusterSA(Int_t layer, Int_t clnumb) {
  // add one clusters to the list (maximum number=fgkMaxNumberOfClusters)
  Int_t presnum = GetNumberOfClustersSA();
  if(presnum>=fgkMaxNumberOfClusters){
    Warning("AddClusterSA","Maximum number of clusters already reached. Nothing is done\n");
    return;
  }

  fSain[presnum] = (layer<<28)+clnumb;  
  presnum++;
  SetNumberOfClustersSA(presnum);
}

//____________________________________________________________
void AliITStrackSA::AddClusterV2(Int_t layer,Int_t clnumb) {
  // add one clusters to the list (maximum number=6)
  Int_t presnum = GetNumberOfClusters();
  if(presnum>=kMaxLayer){
    Warning("AddClusterV2","Maximum number of clusters already reached. Nothing is done\n");
    return;
   }    

  fIndex[presnum] = (layer<<28)+clnumb;  
  presnum++;
  SetNumberOfClusters(presnum);
}