runaction.cc

#include "runaction.hh"


// define the constructor
runaction::runaction(){
  // add new units for dose
  //
  const G4double milligray = 1.e-3*gray;
  const G4double microgray = 1.e-6*gray;
  const G4double nanogray  = 1.e-9*gray;
  const G4double picogray  = 1.e-12*gray;

  new G4UnitDefinition("milligray", "milliGy" , "Dose", milligray);
  new G4UnitDefinition("microgray", "microGy" , "Dose", microgray);
  new G4UnitDefinition("nanogray" , "nanoGy"  , "Dose", nanogray);
  new G4UnitDefinition("picogray" , "picoGy"  , "Dose", picogray);

  // Register accumulable to the accumulable manager
  G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
  accumulableManager->RegisterAccumulable(fEvEdep1);
  accumulableManager->RegisterAccumulable(fEvEdep2);
  accumulableManager->RegisterAccumulable(fEvEdep3);
  accumulableManager->RegisterAccumulable(fEvEdep4);

  G4AnalysisManager*man = G4AnalysisManager::Instance();
  //man->SetVerboseLevel(1);
  man->SetNtupleMerging(true,0);  
}
// define the destructor
runaction::~runaction(){
    // delete G4AnalysisManager::Instance();
}

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

    // inform the runManager to save random number seed
    G4RunManager::GetRunManager()->SetRandomNumberStore(false);

    // reset accumulables to their initial values
    G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
    accumulableManager->Reset();


    // initializate the run manager to store the hits
    G4AnalysisManager*man = G4AnalysisManager::Instance();  
    //man->SetHistoDictoryName("histograms");
    // create a file for each run

    G4int RunID = run->GetRunID();
    std::stringstream strRunID;
    strRunID<<RunID;

    man->OpenFile("output"+strRunID.str()+".root");   


    man->CreateNtuple("Detector_1","Edep");
    man->CreateNtupleDColumn("fedep1");
    man->FinishNtuple(0);   

    man->CreateNtuple("Detector_2","Edep");
    man->CreateNtupleDColumn("fedep2");
    man->FinishNtuple(1);

    man->CreateNtuple("Detector_3","Edep");
    man->CreateNtupleDColumn("fedep3");
    man->FinishNtuple(2);

    man->CreateNtuple("Detector_4","Edep");
    man->CreateNtupleDColumn("fedep4");
    man->FinishNtuple(3);

    man->CreateNtuple("TOF_Det1","TOF");
    man->CreateNtupleDColumn("TOF1");
    man->FinishNtuple(4);

    man->CreateNtuple("TOF_Det2","TOF");
    man->CreateNtupleDColumn("TOF2");
    man->FinishNtuple(5);

    man->CreateNtuple("TOF_Det3","TOF");
    man->CreateNtupleDColumn("TOF3");
    man->FinishNtuple(6);

    man->CreateNtuple("TOF_Det4","TOF");
    man->CreateNtupleDColumn("TOF4");
    man->FinishNtuple(7);

}
    

void runaction::EndOfRunAction(const G4Run*){


    // Merge accumulables
    G4AccumulableManager* accumulableManager = G4AccumulableManager::Instance();
    accumulableManager->Merge();

    // Compute dose = total energy deposit in a run and its variance
    //
    G4double edep1  = fEvEdep1.GetValue();
    G4double edep2  = fEvEdep2.GetValue();
    G4double edep3  = fEvEdep3.GetValue();
    G4double edep4  = fEvEdep4.GetValue();

    const detectorconstruction* detectorConstruction = static_cast<const detectorconstruction*>
        (G4RunManager::GetRunManager()->GetUserDetectorConstruction());

    // get mass of detectors

    G4double mass1 = detectorConstruction->GetScoringVolume1()->GetMass();
    G4double mass2 = detectorConstruction->GetScoringVolume2()->GetMass();
    G4double mass3 = detectorConstruction->GetScoringVolume3()->GetMass();
    G4double mass4 = detectorConstruction->GetScoringVolume4()->GetMass();

    //mass=mass/kg;
    //edep=edep/joule;
    G4double dose1=edep1/mass1;
    G4double dose2=edep2/mass2;
    G4double dose3=edep3/mass3;
    G4double dose4=edep4/mass4;


    // initializate the run manager to store the hits
    G4AnalysisManager*man = G4AnalysisManager::Instance();
    // write all the edep in the root file
    man->Write();
    // close root file
    man->CloseFile();
    #ifndef G4MULTITHREADED
        G4cout<<"Radioactive Particle Tracking Application"<<G4endl;
        G4cout<<"This example simulates a radioactive source inside a pipe"<<G4endl;
        G4cout<<"Code made by: Sebastian Sarasti Zambonino"<<G4endl;
        G4cout<<"Nuclear Science Department"<<G4endl;
        G4cout<<"Escuela Politecnica Nacional"<<G4endl;
    #endif
    G4cout<<"Dose deposited in detector 1 is: "<< G4BestUnit(dose1,"Dose")<<G4endl;
    G4cout<<"Dose deposited in detector 2 is: "<< G4BestUnit(dose2,"Dose")<<G4endl;
    G4cout<<"Dose deposited in detector 3 is: "<< G4BestUnit(dose3,"Dose")<<G4endl;
    G4cout<<"Dose deposited in detector 4 is: "<< G4BestUnit(dose4,"Dose")<<G4endl;


}