added more histograms to fill. Got rid of the mylar energy loss

Author: William Quinn

SensitivityModule.cpp

@@ -230,15 +230,16 @@ SensitivityModule::process(datatools::things& workItem) {
                 else if (first_bool[pos] && !second_bool[pos] && enter_vol_bool) {
                     second_bool[pos] = true;
                 }
-                else if (first_bool[pos] && second_bool[pos] && !third_bool[pos] && exit_vol_bool) {
-                    third_bool[pos] = true;
-                }
+                // This should be uncommented to look at the effect of the mylar
+                //else if (first_bool[pos] && second_bool[pos] && !third_bool[pos] && exit_vol_bool) {
+                //    third_bool[pos] = true;
+                //}
 
-                // If the e has exited the foil, entered the mylar and existed the mylar
+                // If the e has exited the foil
                 // This is a candidate
-                if (first_bool[pos] && second_bool[pos] && third_bool[pos]){
+                if (first_bool[pos] && second_bool[pos]){
                     //escaped_foil[pos] = true;
-                    simdata_.escaped_foil_ = true;
+                    if (simdata_.escaped_foil_){simdata_.escaped_foil_ = true;}else{simdata_.escaped_foil_ = false;}
                     continue;
                 }
 
@@ -362,6 +363,6 @@ void SensitivityModule::ResetVars()
     simdata_.energy_loss_.clear();
     simdata_.displacement_.clear();
     simdata_.track_ids_.clear();
-    simdata_.escaped_foil_ = false;
+    simdata_.escaped_foil_ = true;
 
 }

SensitivityModule.h

@@ -78,7 +78,7 @@ typedef struct SimDataStorage{
     std::vector<double> energy_loss_;
     std::vector<double> displacement_;
     std::vector<int> track_ids_;
-    bool escaped_foil_;
+    bool escaped_foil_ = true;
 
 
 }simdatastorage;

read_file.py

@@ -54,6 +54,14 @@ def main():
                                               n_bins, min_energy, max_energy)
     summed_sim_energy_foil_hist = ROOT.TH1D("summed_sim_energy_foil_hist", "summed_sim_energy_foil_hist",
                                             n_bins, min_energy, max_energy)
+    summed_sim_energy_foil_hist_same = ROOT.TH1D("summed_sim_energy_foil_hist_same", "summed_sim_energy_foil_hist_same",
+                                                 n_bins, min_energy, max_energy)
+    summed_sim_energy_foil_hist_opposite = ROOT.TH1D("summed_sim_energy_foil_hist_opposite",
+                                                     "summed_sim_energy_foil_hist_opposite",
+                                                     n_bins, min_energy, max_energy)
+    summed_sim_energy_foil_hist_rejected = ROOT.TH1D("summed_sim_energy_foil_hist_rejected",
+                                                     "summed_sim_energy_foil_hist_rejected",
+                                                     n_bins, min_energy, max_energy)
     summed_sim_energy_foil_hist_eff = ROOT.TH1D("summed_sim_energy_foil_hist_eff", "summed_sim_energy_foil_hist_eff",
                                                 n_bins, min_energy, max_energy)
     displacement_hist = ROOT.TH1D("displacement_hist", "displacement_hist",
@@ -82,14 +90,15 @@ def main():
             previous_time = t
 
         escape_event = True
+        sides_event = True
         summed_vertex_energy = 0
         summed_foil_energy = 0
 
         for index in range(len(event.og_energy)):
             vertex_energy = event.og_energy[index]
             energy_loss = event.energy_loss[index]
             foil_energy = vertex_energy - energy_loss
-            if foil_energy == 0.0 or foil_energy < 0.0001:
+            if foil_energy == 0.0:
                 escape_event = False
                 eff_energy_bin = int(vertex_energy/energy_bin_width)
                 val = efficiency_hist.GetBinContent(eff_energy_bin)
@@ -108,9 +117,29 @@ def main():
             summed_vertex_energy += vertex_energy
             summed_foil_energy += foil_energy
 
+        index0 = np.where(np.array(list(event.step_track_id)) == event.track_ids[0])[0][-1]
+        index1 = np.where(np.array(list(event.step_track_id)) == event.track_ids[1])[0][-1]
+        last_x0 = list(event.step_stop_vertex_x)[index0]
+        last_x1 = list(event.step_stop_vertex_x)[index1]
+
+        t = thickness/1000
+        if abs(last_x0) < t / 2 and abs(last_x1) < t / 2:
+            sides_event = False
+
+        if (last_x0 < 0 and last_x1 < 0) or (last_x0 > 0 and last_x1 > 0):
+            case = 0
+        else:
+            case = 1
+
         summed_sim_energy_vertex_hist.Fill(summed_vertex_energy)
-        if escape_event:
+        if escape_event and sides_event:
             summed_sim_energy_foil_hist.Fill(summed_foil_energy)
+            if case == 0:
+                summed_sim_energy_foil_hist_same.Fill(summed_foil_energy)
+            if case == 1:
+                summed_sim_energy_foil_hist_opposite.Fill(summed_foil_energy)
+        else:
+            summed_sim_energy_foil_hist_rejected.Fill(summed_foil_energy)
         summed_sim_energy_foil_hist_eff.Fill(summed_foil_energy)
 
         i_event += 1