refactor: HTCC

Author: c-dilks

qa-detectors/cuts/cuts.txt

@@ -1,11 +1,6 @@
 ltcc ltcc_elec_nphe_sec          11     14     counts sec3
 ltcc ltcc_elec_nphe_sec          11     14     counts sec5
 
-htcc htcc_nphe_sec_npheMean      12.5   NB     counts
-htcc htcc_nphe_sec_normFactor    0.9    1.1    unitless
-htcc htcc_vtimediff_sector_mean  -1     1      ns
-htcc htcc_vtimediff_sector_sigma NB     1      ns
-
 ftof ftof_edep_p1a_midangles     8.75   10.75  MeV
 ftof ftof_edep_p1b_midangles     10.5   12.5   MeV
 ftof ftof_edep_p2                8.75   10.75  MeV

src/main/java/org/jlab/clas/timeline/analysis/htcc/htcc_nphe_sector.groovy

@@ -2,6 +2,7 @@ package org.jlab.clas.timeline.analysis
 import java.util.concurrent.ConcurrentHashMap
 import org.jlab.groot.data.TDirectory
 import org.jlab.groot.data.GraphErrors
+import org.jlab.clas.timeline.util.QA
 
 class htcc_nphe_sector {
   def data = new ConcurrentHashMap()
@@ -31,7 +32,7 @@ class htcc_nphe_sector {
       TDirectory out = new TDirectory()
       out.mkdir('/timelines')
 
-      data.each { sec, runs ->
+      def graphList = data.collect { sec, runs ->
         GraphErrors graph
 
         if (plotType == 'npheMean') {
@@ -57,11 +58,15 @@ class htcc_nphe_sector {
             graph.addPoint(it.run, it.correctionFactor, 0, 0)
           }
         }
-
-        out.cd('/timelines')
-        out.addDataSet(graph)
+        graph
       }
-      out.writeFile("htcc_nphe_sec_${plotType}.hipo")
+
+      out.cd('/timelines')
+      if(plotType == 'npheMean')
+        QA.cutGraphs(*graphList, lb: 12.5, out: out)
+      else if(plotType == 'normFactor')
+        QA.cutGraphs(*graphList, lb: 0.9, ub: 1.1, out: out)
+      out.writeFile("htcc_nphe_sec_${plotType}_QA.hipo")
     }
   }
 

src/main/java/org/jlab/clas/timeline/analysis/htcc/htcc_vtimediff_sector.groovy

@@ -4,6 +4,7 @@ import java.util.concurrent.ConcurrentHashMap
 import org.jlab.groot.data.TDirectory
 import org.jlab.groot.data.GraphErrors
 import org.jlab.clas.timeline.fitter.HTCCFitter
+import org.jlab.clas.timeline.util.QA
 
 class htcc_vtimediff_sector {
 
@@ -32,7 +33,7 @@ class htcc_vtimediff_sector {
       TDirectory out = new TDirectory()
       out.mkdir('/timelines')
 
-      data.sort{it.key}.each{ttl, runs->
+      def grtlList = data.sort{it.key}.collect{ttl, runs->
 
         def grtl = new GraphErrors(ttl)
         grtl.setTitle("HTCC vtime - STT, electrons")
@@ -47,11 +48,12 @@ class htcc_vtimediff_sector {
 
           grtl.addPoint(it.run, it[name], 0, 0)
         }
-
-        out.cd('/timelines')
-        out.addDataSet(grtl)
+        grtl
       }
-      out.writeFile("htcc_vtimediff_sector_${name}.hipo")
+
+      out.cd('/timelines')
+      QA.cutGraphsMeanSigma(name, *grtlList, mean_lb: -1, mean_ub: 1, sigma_ub: 1, out: out)
+      out.writeFile("htcc_vtimediff_sector_${name}_QA.hipo")
     }
   }