/*
 *====================================================================
 *
 *  CBM Level 1 Reconstruction 
 *  
 *  Authors: I.Kisel,  S.Gorbunov
 *
 *  e-mail : ikisel@kip.uni-heidelberg.de 
 *
 *====================================================================
 *
 *  L1 Monte Carlo information
 *
 *====================================================================
 */

#ifndef CbmL1MC_H
#define CbmL1MC_H

#include <vector>

// class L1StsHit;

struct CbmL1MCPoint
{ 
  double x, y, z, px, py, pz;
  double p, q, mass;
  int pdg, ID, mother_ID;
  int iStation;
//   std::vector<int> sdsdf;
  static bool compareIDz( const CbmL1MCPoint &a, const CbmL1MCPoint &b )
  {
    return ( a.ID < b.ID ) || ( ( a.ID == b.ID ) && (a.z < b.z) );
  }

  static bool pcompareIDz( const CbmL1MCPoint *a, const CbmL1MCPoint *b )
  {
    return ( a->ID < b->ID ) || ( ( a->ID == b->ID ) && (a->z < b->z) );
  }
};

struct CbmL1MCTrack
{
  double x, y, z, px, py, pz;	// px,py,pz - Momentum components at start vertex [GeV]  **/
  double p, q, mass;
  int pdg, ID, mother_ID;	//  Index of mother track(what is it??). mother_ID = -1 for primary(=vertex) track and >=0 for secondary
  std::vector<CbmL1MCPoint*> Points;  	// MCpoints in track
  std::vector<int> StsHits;		// hits in track, contain index of this hits in vHitStore
  int nMCContStations;    		// how many station, where it was detect??
  int nHitContStations;

// 	CbmL1MCTrack() {Points.clear(); StsHits.clear();};
// 	~CbmL1MCTrack() {Points.clear(); StsHits.clear();};
};


#endif