#include "event.hh"
#include "G4UnitsTable.hh"
#include "G4SystemOfUnits.hh"

MyEventAction::MyEventAction(MyRunAction*)
{
    fEdep = 0.;
    fEdep1 = 0.;
    fEdep_2 = 0.;
    fEdep_3 = 0.;
    fEdep_4 = 0.;
    fEdep_5 = 0.;
    fEdep_6 = 0.;
    fEdep_7 = 0.;
    fEdep_8 = 0.;
    fEdep_9 = 0.;
    fEdep_10 = 0.;

    thickness_count = 0.;


}

MyEventAction::~MyEventAction()
{}

void MyEventAction::BeginOfEventAction(const G4Event*)
{
    fEdep = 0.;
    fEdep1 = 0.;
    fEdep_2 = 0.;
    fEdep_3 = 0.;
    fEdep_4 = 0.;
    fEdep_5 = 0.;
    fEdep_6 = 0.;
    fEdep_7 = 0.;
    fEdep_8 = 0.;
    fEdep_9 = 0.;
    fEdep_10 = 0.;

    thickness_count = 0.;
}

void MyEventAction::EndOfEventAction(const G4Event*)
{
    //G4cout << "Energy deposition: " << fEdep*MeV << G4endl;
   // G4cout << "  Detector positions: " << fPosition.getX()*cm << " " << fPosition.getY()*cm << " " << fPosition.getZ()*cm << G4endl;
       G4int evt = G4RunManager::GetRunManager()->GetCurrentEvent()->GetEventID();
       G4int numberOfEvents = G4RunManager::GetRunManager()->GetCurrentRun()->GetNumberOfEventToBeProcessed();
       G4int count = 0;

       G4AnalysisManager *man = G4AnalysisManager::Instance();

       if (fEdep > 0.0) {
         // G4cout << "Number of events: " << evt << G4endl;
          //G4cout << "Energy deposition: " << fEdep*MeV << G4endl;
          man->FillNtupleDColumn(1, 0, fEdep*MeV);
          man->AddNtupleRow(1);
          man->FillNtupleDColumn(0, 0, evt);
          man->FillNtupleDColumn(0, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(0, 2, fPosition.getY()*cm);
          man->FillNtupleDColumn(0, 3, fPosition.getZ()*cm);
          man->FillH2(0, fPosition.getX()*cm, fPosition.getY()*cm);
  
          man->AddNtupleRow(0);
          theCollID.push_back(evt);
          G4int Ncol = theCollID.size();
          //G4cout << "total number events OPPAC_1:"<< Ncol << G4endl;


       }

       if (fEdep_2 > 0.0){
          man->FillNtupleDColumn(2,0,fEdep_2*MeV);
          man->FillNtupleDColumn(2, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(2, 2, fPosition.getY()*cm);
          man->FillH2(1, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(2);

       }

       if (fEdep_3 > 0.0){
          man->FillNtupleDColumn(3,0,fEdep_3*MeV);
          man->FillNtupleDColumn(3, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(3, 2, fPosition.getY()*cm);
          man->FillH2(2, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(3);

       }

       if (fEdep_4 > 0.0){
          man->FillNtupleDColumn(4,0,fEdep_4*MeV);
          man->FillNtupleDColumn(4, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(4, 2, fPosition.getY()*cm);
          man->FillH2(3, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(4);

       }

       if (fEdep_5 > 0.0){
          man->FillNtupleDColumn(5,0,fEdep_5*MeV);
          man->FillNtupleDColumn(5, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(5, 2, fPosition.getY()*cm);
          man->FillH2(4, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(5);


       }

       if (fEdep_6 > 0.0){
          man->FillNtupleDColumn(6,0,fEdep_6*MeV);
          man->FillNtupleDColumn(6, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(6, 2, fPosition.getY()*cm);
          man->FillH2(5, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(6);

       }

       if (fEdep_7 > 0.0){
          man->FillNtupleDColumn(7,0,fEdep_7*MeV);
          man->FillNtupleDColumn(7, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(7, 2, fPosition.getY()*cm);
          man->FillH2(6, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(7);


       }

       if (fEdep_8 > 0.0){
          man->FillNtupleDColumn(8,0,fEdep_8*MeV);
          man->FillNtupleDColumn(8, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(8, 2, fPosition.getY()*cm);
          man->FillH2(7, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(8);


       }

       if (fEdep_9 > 0.0){
          man->FillNtupleDColumn(9,0,fEdep_9*MeV);
          man->FillNtupleDColumn(9, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(9, 2, fPosition.getY()*cm);
          man->FillH2(8, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(9);


       }

       if (fEdep_10 > 0.0){
          man->FillNtupleDColumn(10,0,fEdep_10*MeV);
          man->FillNtupleDColumn(10, 1, fPosition.getX()*cm);
          man->FillNtupleDColumn(10, 2, fPosition.getY()*cm);
          man->FillH2(9, fPosition.getX()*cm, fPosition.getY()*cm);
          man->AddNtupleRow(10);

       }

       if (fEdep1 > 0.0){
          man->FillNtupleDColumn(11, 0, fEdep1*MeV);
          man->AddNtupleRow(11);
       }
           // count thickness histogram
//       man->FillNtupleDColumn(0, thickness);
  //     man->AddNtupleRow(3);
}