SteppingAction.cc

Go to the documentation of this file.

#include "SteppingAction.hh"


SteppingAction::SteppingAction(DetectorConstruction* det,
                                         EventAction* evt)
:detector(det), eventaction(evt)
{
  
}

SteppingAction::~SteppingAction()
{ }

void SteppingAction::UserSteppingAction(const G4Step* aStep)
{

  G4cout.precision(15); //increase precision of cout.
  
  // get volume of the current step
  G4VPhysicalVolume* volume 
  = aStep->GetPreStepPoint()->GetTouchableHandle()->GetVolume();
  
  // collect energy and track length step by step
  //G4double edep = aStep->GetTotalEnergyDeposit();
  G4double edep = aStep->GetTotalEnergyDeposit();

  // record the time
  G4double time = aStep->GetTrack()->GetGlobalTime();
  
  if (volume == detector->GetScintillator())
  {
    //Determine the scintillator ID.
    G4int scintillatorID = aStep->GetPreStepPoint()->GetTouchable()->GetReplicaNumber(1);
    scintillatorID += (aStep->GetPreStepPoint()->GetTouchable()->GetReplicaNumber(2))*(detector->GetNumberOfScintillatorsPerBox());
    if (aStep->GetPreStepPoint()->GetTouchable()->GetVolume(3) == detector->GetLowerBetaDetector())
      scintillatorID += (detector->GetNumberOfScintillatorBoxes())*(detector->GetNumberOfScintillatorsPerBox());


    eventaction->AddScintillatorTubeEvent(scintillatorID,time, edep);
    std::stringstream ss;
    ss.precision(15);
    ss << time/ns << " " << scintillatorID << " " << edep/MeV << " " << aStep->GetTrack()->GetDefinition()->GetParticleName().c_str() << " " << aStep->GetTrack()->GetTotalEnergy()/MeV << std::endl;
    eventaction->PrintToVerboseFile(ss.str());
  }

  if (volume == detector->GetGermaniumDetector())
  {
    eventaction->AddGeDetectorEvent(time, edep);


    std::stringstream ss;
    ss.precision(15);
    ss << time/ns << " GermaniumDetector " << edep/MeV << " " << aStep->GetTrack()->GetDefinition()->GetParticleName().c_str() << " " << aStep->GetTrack()->GetTotalEnergy()/MeV << std::endl;
    eventaction->PrintToVerboseFile(ss.str());
  }

  if (volume == detector->GetLowerFrontScintillatorPlate())
    {
      eventaction->AddScintillatorPlateEvent(AllAccumulatedStatistics::LowerFront,time/ns,edep);
    }

  
  if (volume == detector->GetUpperFrontScintillatorPlate())
    {
      eventaction->AddScintillatorPlateEvent(AllAccumulatedStatistics::UpperFront,time, edep);
    }
  
  if (volume == detector->GetBackScintillatorPlate())
  {
    eventaction->AddScintillatorPlateEvent(AllAccumulatedStatistics::Back,time, edep);
  }
 
  if (volume == detector->GetUpperFrontScintillatorPlate() || volume == detector->GetBackScintillatorPlate() || volume == detector->GetLowerFrontScintillatorPlate())
  {
    if (edep/eV > 1)
      {
        
        std::stringstream ss;
        ss.precision(15);
        ss << time/ns << " " << volume->GetName().c_str() << " " << edep/MeV << " " << aStep->GetTrack()->GetDefinition()->GetParticleName().c_str() << " " << aStep->GetTrack()->GetTotalEnergy()/MeV << std::endl;
        eventaction->PrintToVerboseFile(ss.str());
      }
  }
  
  //if (aStep->GetTrack()->GetDefinition()->GetParticleName().compare("gamma")==0)
  //    fprintf(file, "GAMMA: Time: %.15f, Energy: %f, Region: %s \n", time, aStep->GetTrack()->GetTotalEnergy()/MeV, volume->GetName().c_str());
                        
 // G4cout << "Time: " << aStep->GetTrack()->GetGlobalTime()/s << G4endl;
  
}