#include "run.hh"
#include "G4Run.hh"
#include "G4SDManager.hh"
#include "G4HCofThisEvent.hh"
#include "G4Event.hh"
//#include "G4ios.hh"
#include "G4SystemOfUnits.hh"
#include "TFile.h"
#include "TNtuple.h"

MyRunAction::MyRunAction()
{
    
}

MyRunAction::~MyRunAction()
{}

void MyRunAction::BeginOfRunAction(const G4Run* run)
{

    //G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();

    //int runID = run->GetRunID();

    //d::stringstream strRunID;
    //rRunID << runID;
    //analysisManager->OpenFile("output.root");
   //analysisManager->OpenFile("output" + strRunID.str() + ".root");

    TFile* rootFile = new TFile("output.root", "RECREATE");
    TNtuple* ntuple = new TNtuple("trajectory", "Proton trajectory", "x:y:z");
}
  

void MyRunAction::EndOfRunAction(const G4Run*)
{
    // Close the ROOT file
    TFile* rootFile = TFile::Open("output.root", "UPDATE");
    rootFile->Write();
    rootFile->Close();

}
/*
 void MyRunAction::AddEventInfo(G4int eventID, G4double edep, G4ThreeVector position)
{
    fTotalEdep += edep;
    fPositions.push_back(position);
    fEventIDs.push_back(eventID);
}
*/