Selaa lähdekoodia

Adds plot for number of good seeds

Caleb 6 vuotta sitten
vanhempi
commit
86ba1bfac8
3 muutettua tiedostoa jossa 12 lisäystä ja 6 poistoa
  1. 8 2
      looper/analysis/tracking_eff.cpp
  2. 1 1
      looper/filval
  3. 3 3
      looper/run_all.condor

+ 8 - 2
looper/analysis/tracking_eff.cpp

@@ -58,8 +58,8 @@ bool in_lum_region(const SimVertex& vertex) {
 bool is_good_sim(const SimTrack& sim_track) {
     const auto& vertex = sim_vertices[sim_track.parentVtxIdx()];
     return abs(sim_track.pdgId()) == 11 and  // electrons
-           in_lum_region(vertex) and  // from luminous region
-           abs(sim_track.eta())<3.0;  // inside ECAL acceptance
+           in_lum_region(vertex) and         // from luminous region
+           abs(sim_track.eta())<3.0;         // inside ECAL acceptance
 };
 
 //bool is_good_fake(const SimTrack& sim_track) {
@@ -341,6 +341,7 @@ void run(){
     auto& ecal_energy_resolution = *tds.register_container<ContainerTH1<float>>("ecal_energy_resolution ", THParams::lookup("ecal_energy_resolution"));
 
     auto& n_seeds = *tds.register_container<ContainerTH1<size_t>>("n_seeds", THParams::lookup("n_seeds"));
+    auto& n_good_seeds = *tds.register_container<ContainerTH1<size_t>>("n_good_seeds", THParams::lookup("n_seeds"));
     auto& n_good_sim = *tds.register_container<ContainerTH1<int>>("n_good_sim", THParams::lookup("n_seeds"));
     auto& n_gsf_tracks = *tds.register_container<ContainerTH1<int>>("n_gsf_track", THParams::lookup("n_seeds"));
 
@@ -375,6 +376,11 @@ void run(){
     while (tds.next()) {
 
         n_seeds.fill(seeds.size());
+        size_t _n_good_seeds = 0;
+        for (const auto& seed : seeds)
+            if (is_good_seed(seed, hoe_cut)) _n_good_seeds++;
+        n_good_seeds.fill(_n_good_seeds);
+
         auto match_stats =  get_match_stats(false);
         n_good_sim.fill(std::get<0>(match_stats));
         n_gsf_tracks.fill(std::get<1>(match_stats));

+ 1 - 1
looper/filval

@@ -1 +1 @@
-Subproject commit 0df63c2446483e4a12db1fcb3c67fa2ba61c4c1e
+Subproject commit a83921cf4b2033e4278b43fa69f7059cd76a0bf1

+ 3 - 3
looper/run_all.condor

@@ -10,10 +10,10 @@ output     = condor_logs/out.$(Process)
 log        = condor_logs/foo.log
 
 queue arguments from (
-    zee-new-default
-    zee-new-wide
-    zee-old-default
     tt-new-default
     tt-new-wide
     tt-old-default
+    zee-new-default
+    zee-new-wide
+    zee-old-default
 )