/*
  The TrackStore class contains all the methods necessary to prepare 
  the tracks to feed the alignment algorithm
*/

#include "TrackStore.h"
#include "Track.h"
#include "Line.h"
#include "constants.h"
#include <cstdlib>
#include <fstream>
using namespace std;

TrackStore::TrackStore()
{  
  cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" << endl; 
  cout << "%                                     %" << endl;
  cout << "%    Entering into the TrackStore     %" << endl;
  cout << "%                                     %" << endl;
  cout << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%" << endl; 
  cout << "" << endl;
}
TrackStore::~TrackStore()
{
}


//=================================================================
// PrepareStore
//=================================================================
//
// read tracks from ASCII file
//
//=================================================================

void TrackStore::PrepareStore(Tracks& my_tracks, int n_event)
{
  int event = 0;

  // while (event<n_event)  
  // {
  //   if (debug_mode) cout << "Read event number " << event << endl; 

    bool sc = TrackStore::ReadEvent(my_tracks,n_event);
  //   if (sc) event++;
  // }
  if (debug_mode) cout << "Number of tracks " << my_tracks.size() << endl; 
}


//=================================================================
// ReadEvent
//=================================================================
//
// Read a set of tracks
//
//=================================================================

bool TrackStore::ReadEvent(Tracks& my_event_tracks,int n_event)
{

  /// read file ASCII file---------------------------------------------------------
  std::ifstream in;
  in.open("/data/FAIRsorf/pandaroot/karavdina/myKnossos/TracksForAligment.txt");
  Double_t xhit,yhit,zhit,errxhit,erryhit,errzhit,sensorID;
  int trkend=0;
  Int_t ntrk = 0;
  // cout<<" ntrk = "<<ntrk<<" "<<n_event<<endl;
  while(ntrk<n_event){
    // cout<<" ntrk = "<<ntrk<<endl;
    Track *new_track = new Track();
    while (trkend!=1) {
      //cout<<"trkend = "<<trkend<<endl;
      in>>xhit>>errxhit>>yhit>>erryhit>>zhit>>errzhit>>sensorID>>trkend;
      zhit+=0.0075;
      new_track->addHit(xhit,errxhit,yhit,erryhit,zhit,errzhit,sensorID);
      if (debug_mode) cout<<xhit<<" "<<errxhit<<" "<<yhit<<" "<<erryhit<<" "<<zhit<<" "<<errzhit<<" "<<sensorID<<endl;
    }
    if (debug_mode)    cout<<" "<<endl;
    new_track->setTracknum(ntrk);
    my_event_tracks.push_back(*new_track);
    trkend=0;
    ntrk++;
  }
  return true;
}