Remove opt mode

We'll likely reuse the scan results.

README.md

@@ -57,5 +57,5 @@ get_latex_tables.py -i ~/Entanglement/analysis/SuperKEKB/2023Oct22 \
 # to produce tables for a particular cut in |cos(theta)|
 get_latex_tables.py -i ~/Entanglement/analysis/SuperKEKB/2023Oct22 \
                        ~/Entanglement/analysis/SuperKEKB/2023Oct22_hadHad \
-                    -C absCosTheta0p55 # or "opt" if you have run the analysis with -M opt
+                    -C absCosTheta0p55
 ```

scripts/get_latex_tables.py

@@ -6,13 +6,11 @@
 #
 # In case the JSON files are stored in multiple directories, you can list those directories as argument to
 # -i/--input as long as there are no files with the same name in those directories.
-# If you want to produce Latex tables for a particular cut in |cos(theta)|, you can do that by adding
+# If you want to produce Latex tables for a particular cut in |cos(theta)|, you can do that by adding, e.g.,
 #
-# -c <cut>
-#
-# to the command, where <cut> can be "opt" if you want tables for the optimal cut, or something like
-# "absCosTheta0p55" if you want tables for the cut |cos(theta)| <= 0.55.
+# -c absCosTheta0p55
 #
+# if you want tables for the cut |cos(theta)| <= 0.55.
 # To colorize the terminal output, pipe it to: pygmentize -l latex
 
 import argparse

test/analyzeNtuples.py

@@ -19,7 +19,7 @@
 decayMode_choices = [ "pi_pi", "pi_rho", "pi_a1", "rho_rho", "rho_a1", "a1_a1" ]
 decayMode_choices_all = decayMode_choices + [ "had_had" ]
 spinAnalyzer_choices = [ "by_summation", "by_mlfit", "by_differentialXsec1d", "by_differentialXsec2d", "by_asymmetry" ]
-analysis_choices = [ "inclusive", "scan", "optimal" ]
+analysis_choices = [ "inclusive", "scan" ]
 plot_choices = [ "makeResolutionPlots", "makeControlPlots" ]
 
 parser = argparse.ArgumentParser(formatter_class = argparse.ArgumentDefaultsHelpFormatter)
@@ -72,7 +72,6 @@
     from TauAnalysis.Entanglement.samples import samples_SuperKEKB as samples, PAR_GEN_SUPERKEKB as par_gen, \
       MAX_SUM_PHOTON_EN_SUPERKEKB as maxSumPhotonEn
     absCosTheta_acceptanceCut = 0.92 # obtained with acceptanceCalculator
-    absCosTheta_optimalCut = -1 #TODO TBD
 else:
     assert(False)
 
@@ -182,9 +181,6 @@ def get_all_values(dictionary):
                 absCosTheta_cuts = [ -1 ]
               elif analysis_mode == 'scan':
                 absCosTheta_cuts = absCosTheta_bins
-              elif analysis_mode == 'optimal':
-                assert(absCosTheta_optimalCut > 0)
-                absCosTheta_cuts = [ absCosTheta_optimalCut ]
               else:
                 assert(False)
               for absCosTheta_cut in absCosTheta_cuts:
@@ -193,8 +189,6 @@ def get_all_values(dictionary):
                 suffix = ""
                 if analysis_mode == "scan":
                   suffix = f"absCosTheta{absCosTheta_cut:.2f}".replace(".", "p")
-                elif analysis_mode == "optimal":
-                  suffix = "opt"
                 cfg_baseName = f"analyzeEntanglementNtuple_{sampleName}_{mode}Mode_{hAxis}Axis_{decayMode}DecayMode_{spinAnalyzer}"
                 if suffix:
                   cfg_baseName += f"_{suffix}"