#include "BeAnalysis.hh"
#include "BeRunAction.hh"

#include "G4Run.hh"
#include "G4RunManager.hh"
#include "G4UnitsTable.hh"
#include "G4SystemOfUnits.hh"

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

BeRunAction::BeRunAction()
 : G4UserRunAction()
{
  G4RunManager::GetRunManager()->SetPrintProgress(1);
  auto analysisManager = G4AnalysisManager::Instance();
  analysisManager->SetVerboseLevel(4);
  analysisManager->SetFirstHistoId(1);
//  analysisManager->SetFileName("Det1.csv");
  // Creating 1D histograms
  analysisManager->CreateH1("Energy","Detector 1 Edep", 100, 0., 1*MeV); // h1 Id = 0
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
BeRunAction::~BeRunAction()
{
  delete G4AnalysisManager::Instance();
}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void BeRunAction::BeginOfRunAction(const G4Run*)
{
  auto analysisManager = G4AnalysisManager::Instance();
  analysisManager->SetActivation(true);

  analysisManager->OpenFile("test");
  //  analysisManager->CreateNtuple("DET1", "Edep, Cell ID and Position");
//  analysisManager->CreateNtupleDColumn("Edep");
//  analysisManager->CreateNtupleDColumn("CellID");
//  analysisManager->CreateNtupleDColumn("Position");
}


//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

void BeRunAction::EndOfRunAction(const G4Run*)
{
  // save histograms & ntuple
  auto analysisManager = G4AnalysisManager::Instance();
//  analysisManager->FillH1(0, 1*MeV);
//  G4cout << " Edep : mean = " << G4BestUnit(analysisManager->GetH1(1)->mean(), "Energy");
  analysisManager->Write();
  analysisManager->CloseFile();
}