//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
/// \file B1SteppingAction.cc
/// \brief Implementation of the B1SteppingAction class

#include "B1SteppingAction.hh"
#include "B1EventAction.hh"
#include "B1DetectorConstruction.hh"

#include "G4Step.hh"
#include "G4Event.hh"
#include "G4RunManager.hh"
#include "g4analysis.hh"
#include "G4LogicalVolume.hh"

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

B1SteppingAction::B1SteppingAction(B1EventAction* eventAction)
: G4UserSteppingAction(),
  fEventAction(eventAction),
  fScoringVolume(0),
  
  //////////---------------START SETTING FOR FSCORINGVOLUME2------------//////
  
  fScoringVolume2(0)
  //////////---------------END SETTING FOR FSCORINGVOLUME2------------//////
{}

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

B1SteppingAction::~B1SteppingAction()
{}

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

void B1SteppingAction::UserSteppingAction(const G4Step* step)
{
	
	 //-----------START DECLARE ANALYSIS MANAGER-------------------------------//
   auto analysisManager = G4AnalysisManager::Instance();
   //-----------END DECLARE ANALYSIS MANAGER-------------------------------//
  if (!fScoringVolume) { 
    const B1DetectorConstruction* detectorConstruction
      = static_cast<const B1DetectorConstruction*>
        (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
    fScoringVolume = detectorConstruction->GetScoringVolume();   
    fScoringVolume2 = detectorConstruction->GetScoringVolume2();
}

  // get volume of the current step
  G4LogicalVolume* volume 
    = step->GetPreStepPoint()->GetTouchableHandle()
      ->GetVolume()->GetLogicalVolume();
  // check if we are in scoring volume

  // collect energy deposited in this step
//  G4double edepStep = step->GetTotalEnergyDeposit();
//  fEventAction->AddEdep(edepStep);
  
  //////////---------------START FSCORINGVOLUME2------------//////
//   if (!fScoringVolume2) {
//    const B1DetectorConstruction* detectorConstruction2
//      = static_cast<const B1DetectorConstruction*>
//        (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
//    fScoringVolume2 = detectorConstruction2->GetScoringVolume2();
//  }

  // get volume of the current step
//  G4LogicalVolume* volume2
//    = step->GetPreStepPoint()->GetTouchableHandle()
//      ->GetVolume()->GetLogicalVolume();
      
  // check if we are in scoring volume

  
  //--------------START CHECK VOLUME WE ARE IN--------------///////
  //if (volume2 != fScoringVolume2) return;
//  if (volume != fScoringVolume2 && volume != fScoringVolume) return;
  //--------------START CHECK VOLUME WE ARE IN--------------///////

  // collect energy deposited in this step
	if (volume == fScoringVolume) {
  G4double edepStep = step->GetTotalEnergyDeposit();
  fEventAction->AddEdep(edepStep);
	}
	if (volume == fScoringVolume2) {
  G4double edepStep2 = step->GetTotalEnergyDeposit();
  fEventAction->AddEdep2(edepStep2);
	}
	
	//----------START CHECK PARTICLE ID----------///////
		if (volume == fScoringVolume) {
			int pdg = step->GetTrack()->GetParticleDefinition()->GetPDGEncoding();
			 analysisManager->FillNtupleDColumn(2, pdg);
  		 
		}
		if (volume == fScoringVolume2) {
			int pdg2 = step->GetTrack()->GetParticleDefinition()->GetPDGEncoding();
			analysisManager->FillNtupleDColumn(3, pdg2);
		}
	
		//----------END CHECK PARTICLE ID--------///////
		
		

  if (volume != fScoringVolume || volume != fScoringVolume2) return;
  
  

  
  //////////---------------END FSCORINGVOLUME2------------//////
}

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