up to output working. missing interface wrapper to python

src/EventAction.cc

@@ -90,6 +90,7 @@ void EventAction::EndOfEventAction(const G4Event* event)
   double tot_energy=0;
   int sensorNumber=0;
   int layerNumber=0;
+  int nactive =0;
   for(const auto& layer: cw->getLayers()){
     for(const auto& sensor: layer.sensors){
       sensorEnergies.push_back(sensor.getEnergy());
@@ -108,6 +109,9 @@ void EventAction::EndOfEventAction(const G4Event* event)
        sensorNumber++;
     }
     layerNumber++;
+    if(layer.isActive){
+        nactive++;
+    }
   }
   //assign the vectors
 
@@ -116,7 +120,8 @@ void EventAction::EndOfEventAction(const G4Event* event)
   analysisManager->FillNtupleDColumn(0, gen->getPartEnergy()); //needs to be true energy DEBUG
   analysisManager->FillNtupleDColumn(1, tot_energy);
   analysisManager->FillNtupleIColumn(2, (int)cw->getLayers().size());
-  analysisManager->FillNtupleIColumn(3, cw->getNSensors());
+  analysisManager->FillNtupleIColumn(3, nactive);
+  analysisManager->FillNtupleIColumn(4, cw->getNSensors());
 
   //vectors are stored automatically
   analysisManager->AddNtupleRow();

src/RunAction.cc

@@ -67,6 +67,7 @@ RunAction::RunAction()
   analysisManager->CreateNtupleDColumn("true_energy");
   analysisManager->CreateNtupleDColumn("total_dep_energy");
   analysisManager->CreateNtupleIColumn("N_layers");
+  analysisManager->CreateNtupleIColumn("N_active_layers");
   analysisManager->CreateNtupleIColumn("N_sensors");
 
   analysisManager->CreateNtupleDColumn("sensor_energy",hitEnergy);
@@ -95,7 +96,7 @@ void RunAction::BeginOfRunAction(const G4Run* /*run*/)
 
   // Open an output file
   //
-  G4String fileName = "B4.root";
+  G4String fileName = "_1234567890_Hits.root";
   // Other supported output types:
   // G4String fileName = "B4.csv";
   // G4String fileName = "B4.hdf5";