/*
 *====================================================================
 *
 *  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 "CbmL1StsHit.h"
#include <vector>
using namespace std;

struct CbmL1MCPoint
{ 
  double x, y, z, px, py, pz;
  double p, q, mass;
  int pdg, ID, mother_ID;
  double T[6];
  double TOut[6];
  int iPileUp;
  int UID; // ID*100+iPileUp

  static bool compareIDz( const CbmL1MCPoint &a, const CbmL1MCPoint &b )
  {
    return (a.iPileUp<b.iPileUp) || (a.iPileUp==b.iPileUp)&&( a.ID < b.ID ) || 
      ( (a.iPileUp==b.iPileUp)&&( a.ID == b.ID ) && (a.z < b.z) );
  }

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

struct CbmL1MCTrack
{
  double mass, q, p, x, y, z, px, py, pz;
  int ID, mother_ID, pdg;
  int iPileUp;
  int UID; // ID*100+iPileUp
  double T[6];
  vector<CbmL1MCPoint*> Points;
  vector<CbmL1StsHit*> StsHits;
};


#endif