// -------------------------------------------------------------------------
// -----                      PndMCTrack header file                   -----
// -----                  Created 03/08/04  by V. Friese               -----
// -----                  Adapted to Panda  M. Al-Turany               -----
// -------------------------------------------------------------------------


/** PndMCTrack.h
 *@author V.Friese <v.friese@gsi.de>
 **
 ** Data class for storing Monte Carlo tracks processed by the PndStack.
 ** A MCTrack can be a primary track put into the simulation or a
 ** secondary one produced by the transport through decay or interaction.
 **/


#ifndef CBMMCTRACK_H
#define CBMMCTRACK_H 1

#include "PndDetectorList.h"
#include "TObject.h"
#include "TParticle.h"
#include "TVector3.h"
#include "TLorentzVector.h"

#ifndef ROOT_TParticlePDG
#include "TParticlePDG.h"
#endif
#ifndef ROOT_TDatabasePDG
#include "TDatabasePDG.h"
#endif

class PndMCTrack : public TObject
{

 public:


  /**  Default constructor  **/
  PndMCTrack();


  /**  Standard constructor  **/
////  not used!!
//  PndMCTrack(Int_t pdgCode, Int_t motherID, TVector3 startVvertex,
//	     Double_t startTime, TLorentzVector momentum, Int_t  nPoint=0);


  /**  Copy constructor  **/
  PndMCTrack(const PndMCTrack& track);


  /**  Constructor from TParticle  **/
  PndMCTrack(TParticle* particle);


  /**  Destructor  **/
  virtual ~PndMCTrack();


  /**  Output to screen  **/
  void Print(Int_t iTrack) const;


  /**  Accessors  **/
  Int_t    GetPdgCode()     const { return fPdgCode; }
  Int_t    GetMotherID()    const { return fMotherID; }
  Int_t    GetSecondMotherID()    const { return fSecondMotherID; }
  TVector3 GetStartVertex() const { return TVector3(fStartX, fStartY,fStartZ);}
  Double_t GetStartTime()   const { return fStartT; }
  TVector3 GetMomentum()    const { return TVector3(fPx, fPy, fPz); }
  Double_t GetPt()       const { return TMath::Sqrt(fPx*fPx+fPy*fPy); }
  /*
  Int_t GetStsPoints()  const { return   (fPoints &  15        )        ; }
  */
   
  Bool_t IsGeneratorCreated(void) const { return (fGeneratorFlags&0x1)!=0; }
  Bool_t IsGeneratorDecayed(void) const { return (fGeneratorFlags&0x2)!=0; }
  Bool_t IsGeneratorLast(void) const { return (fGeneratorFlags&0x1)!=0 && (fGeneratorFlags&0x2)==0; }
  void SetGeneratorCreated(void) { fGeneratorFlags|=0x1; }
  void SetGeneratorDecayed(void) { fGeneratorFlags|=0x2; }

  Int_t  GetNPoints(DetectorId detId)  const;  
  /**  Modifiers  **/
  
  void SetMotherID(Int_t id) { fMotherID = id; }
  void SetSecondMotherID(Int_t id) { fSecondMotherID = id; }
  /*
  void SetStsPoints(Int_t np);
  */ 
 
  
  /**  Add one detector point to the fPoint variable  **/
 /*
  void AddStsPoint()  { SetStsPoints(  GetStsPoints()  + 1 ); }
 */ 
 
  
  void SetNPoints(Int_t iDet, Int_t  np);
  
  TLorentzVector Get4Momentum() const;


private:

  /**  PDG particle code  **/
  Int_t  fPdgCode;

  /** Momentum components at production [GeV]  **/
  Double32_t fPx, fPy, fPz, fE;

  /**  Index of mother track. Zero( Minus One???) for primary particles.  **/
  Int_t  fMotherID;
  Int_t  fSecondMotherID;

  /** Flag if particle was created (bit 0) and/or decayed (bit 1) by generator **/
  Int_t fGeneratorFlags;

  /** Coordinates of start vertex [cm, ns]  **/
  Double32_t fStartX, fStartY, fStartZ, fStartT;

  /**  Bitvector representing the number of MCPoints for this track in 
   **  each subdetector. The detectors are represented by
   **  This goes to fPoints:
   **  DRC:  Bit  0 -  1  (2 bit max. value 3)
   **  MDT:  Bit  2 -  3  (2 bit max. value 3)
   **  MVD:  Bit  4 -  5  (2 bit max. value 3)
   **  DCH:  Bit  6 -  7  (2 bit max. value 3)
   **  EMC:  Bit  8 -  9  (2 bit max. value 3)
   **  STT:  Bit 10 - 11  (2 bit max. value 3)
   **  TPC:  Bit 12 - 13  (2 bit max. value 3)
   **  TOF:  Bit 14 - 15  (2 bit max. value 3)
   **  GEM:  Bit 16 - 17  (2 bit max. value 3)
   **  DSK:  Bit 18 - 19  (2 bit max. value 3)
   **  HYP:  Bit 20 - 21  (2 bit max. value 3)
   **  RPC:  Bit 22 - 23  (2 bit max. value 3)
   **  LUMI: Bit 24 - 25  (2 bit max. value 3)
   **  HYPG: Bit 26 - 27  (2 bit max. value 3)
   **  FTS:  Bit 28 - 29  (2 bit max. value 3)
   
   **  The respective point numbers can be accessed and modified **/

  Int_t  fPoints;

  ClassDef(PndMCTrack,1);

};



#endif