//elisabetta Prencipe 21/5/2014

#include "GeaneMCApplication.h"
#include "TVirtualMC.h"
#include <iostream>
#include "TSystem.h"
#include"assert.h"
#include"FieldManager.h"

GeaneMCApplication::GeaneMCApplication()
{
}

void GeaneMCApplication::ConstructGeometry(){
  gGeoManager->CloseGeometry();
  gMC->SetRootGeometry();  
}

void GeaneMCApplication::InitMC(){
  gMC->Init();
  gMC->BuildPhysics();
}

void GeaneMCApplication::Field(const Double_t* x, Double_t* b) const
{
  //assert(field!=NULL);
  TVector3 pos(x[0],x[1],x[2]);
  //  TVector3 B = FieldManager::getInstance()->getField()->get(pos);
  TVector3 B = FieldManager::getFieldVal(pos);
  //  std::cout << "GeaneMCApplication::Field()" << std::endl;
  //B.Print();
  b[0]=B.X();
  b[1]=B.Y();
  b[2]=B.Z();
}


ClassImp(GeaneMCApplication)