Add "acc(MUON...)" cuts and change EE start

These cuts encode the acceptance used in the analysis in which we combine with the muon groups. EE has been changed to starting at |eta| > 1.56 instead of 1.566, as this is the number agreed upon with the muon groups.
jturkewitz · May 7, 2014 · 89f995d · 89f995d
1 parent 2b5ab89
commit 89f995d
Show file tree

Hide file tree

Showing 4 changed files with 78 additions and 36 deletions.
diff --git a/ZFinder/Event/docs/ZDefinition.md b/ZFinder/Event/docs/ZDefinition.md
@@ -17,25 +17,31 @@ ZFinderElectron.CutPassed method.
 Acceptance cuts make selections of electrons based on the geometric acceptance
 of the detector.
 
-| Cut Name      | Explanation                                                 |
-|---------------|-------------------------------------------------------------|
-| acc(ALL)      | Always true (used for counting total events, mainly in MC). |
-| acc(DETECTOR) | In any of the regions defined below, not including ALL.     |
-| acc(EB)       | -1.4442 < eta < 1.4442                                      |
-| acc(EB+)      | 0 < eta < 1.4442                                            |
-| acc(EB-)      | -1.4442 < eta < 0                                           |
-| acc(EE)       | 1.566 < abs(eta) < 2.850                                    |
-| acc(EE+)      | 1.566 < eta < 2.850                                         |
-| acc(EE-)      | -2.850 < eta < -1.566                                       |
-| acc(ET)       | In acc(EB) or acc(EE)                                       |
-| acc(ET+)      | In acc(EB+) or acc(EE+)                                     |
-| acc(ET-)      | In acc(EB-) or acc(EE-)                                     |
-| acc(NT)       | 2.5 < abs(eta) < 2.850                                      |
-| acc(NT+)      | 2.5 < eta < 2.850                                           |
-| acc(NT-)      | -2.850 < eta < -2.5                                         |
-| acc(HF)       | 3.1 < abs(eta) < 4.6                                        |
-| acc(HF+)      | 3.1 < eta < 4.6                                             |
-| acc(HF-)      | -4.6 < eta < -3.1                                           |
+| Cut Name         | Explanation                                                 |
+|------------------|-------------------------------------------------------------|
+| acc(ALL)         | Always true (used for counting total events, mainly in MC). |
+| acc(DETECTOR)    | In any of the regions defined below, not including ALL.     |
+| acc(EB)          | -1.4442 < eta < 1.4442                                      |
+| acc(EB+)         | 0 < eta < 1.4442                                            |
+| acc(EB-)         | -1.4442 < eta < 0                                           |
+| acc(EE)          | 1.566 < abs(eta) < 2.850                                    |
+| acc(EE+)         | 1.566 < eta < 2.850                                         |
+| acc(EE-)         | -2.850 < eta < -1.566                                       |
+| acc(ET)          | In acc(EB) or acc(EE)                                       |
+| acc(ET+)         | In acc(EB+) or acc(EE+)                                     |
+| acc(ET-)         | In acc(EB-) or acc(EE-)                                     |
+| acc(NT)          | 2.5 < abs(eta) < 2.850                                      |
+| acc(NT+)         | 2.5 < eta < 2.850                                           |
+| acc(NT-)         | -2.850 < eta < -2.5                                         |
+| acc(HF)          | 3.1 < abs(eta) < 4.6                                        |
+| acc(HF+)         | 3.1 < eta < 4.6                                             |
+| acc(HF-)         | -4.6 < eta < -3.1                                           |
+| acc(MUON_TIGHT)  | -2.1 < eta < 2.1                                            |
+| acc(MUON_TIGHT+) | 0 < eta < 2.1                                               |
+| acc(MUON_TIGHT-) | -2.1 < eta < 0                                              |
+| acc(MUON_LOOSE)  | -2.4 < eta < 2.4                                            |
+| acc(MUON_LOOSE+) | 0 < eta < 2.4                                               |
+| acc(MUON_LOOSE-) | -2.4 < eta < 0                                              |
 
 ###### Trigger
 

diff --git a/ZFinder/Event/interface/AcceptanceSetter.h b/ZFinder/Event/interface/AcceptanceSetter.h
@@ -27,7 +27,7 @@ namespace zf {
             static constexpr double EB_MAX_ = 1.4442;
             static constexpr double EB_MIN_ = -EB_MAX_;
             static constexpr double EB_CENT_ = 0;
-            static constexpr double EEP_MIN_ = 1.566;
+            static constexpr double EEP_MIN_ = 1.56;
             static constexpr double EEP_MAX_ = 2.850;
             static constexpr double EEM_MIN_ = -EEP_MAX_;
             static constexpr double EEM_MAX_ = -EEP_MIN_;
@@ -39,6 +39,10 @@ namespace zf {
             static constexpr double HFP_MAX_ = 4.6;
             static constexpr double HFM_MIN_ = -HFP_MAX_;
             static constexpr double HFM_MAX_ = -HFP_MIN_;
+            static constexpr double MUON_TIGHT_MAX_ = 2.1;
+            static constexpr double MUON_TIGHT_MIN_ = -MUON_TIGHT_MAX_;
+            static constexpr double MUON_LOOSE_MAX_ = 2.4;
+            static constexpr double MUON_LOOSE_MIN_ = -MUON_LOOSE_MAX_;
     };
 }  // namespace zf
 #endif  // ZFINDER_ACCEPTANCESETTER_H_
diff --git a/ZFinder/Event/python/ZDefinitions_cfi.py b/ZFinder/Event/python/ZDefinitions_cfi.py
@@ -15,16 +15,16 @@
     # ID cuts
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Gen Cuts"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "gpt>30", "geta<2.1", "geta>-2.1"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "gpt>20", "geta<2.4", "geta>-2.4"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "gpt>30", "geta<2.1", "geta>-2.1"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "gpt>20", "geta<2.4", "geta>-2.4"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(True)
         ),
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Reco Cuts"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>30", "eta<2.1", "eta>-2.1"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>20", "eta<2.4", "eta>-2.4"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(MUON_TIGHT)", "pt>30"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(MUON_LOOSE)", "pt>20"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(False)
@@ -34,16 +34,16 @@
     # trigger (and hence use eg_tight) on one leg
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Double"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "trig(et_et_tight)", "pt>30", "eta<2.1", "eta>-2.1", "eg_medium"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "trig(et_et_loose)", "pt>20", "eta<2.4", "eta>-2.4", "eg_medium"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(MUON_TIGHT)", "trig(et_et_tight)", "pt>30", "eg_medium"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(MUON_LOOSE)", "trig(et_et_loose)", "pt>20", "eg_medium"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(False)
         ),
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Single"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "trig(single_ele)", "pt>30", "eta<2.1", "eta>-2.1", "eg_tight"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "acc(ET)",          "pt>20", "eta<2.4", "eta>-2.4", "eg_medium"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(MUON_TIGHT)", "trig(single_ele)", "pt>30", "eg_tight"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(MUON_LOOSE)", "acc(MUON_LOOSE)",  "pt>20", "eg_medium"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(False)
@@ -90,16 +90,16 @@
     # ID cuts
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Gen Cuts"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "gpt>30", "geta<2.1", "geta>-2.1"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "gpt>20", "geta<2.4", "geta>-2.4"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "gpt>30", "geta<2.1", "geta>-2.1"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "gpt>20", "geta<2.4", "geta>-2.4"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(True)
         ),
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Reco Cuts"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>30", "eta<2.1", "eta>-2.1"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>20", "eta<2.4", "eta>-2.4"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(MUON_TIGHT)", "pt>30"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(MUON_LOOSE)", "pt>20"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(False)
@@ -109,16 +109,16 @@
     # trigger (and hence use eg_tight) on one leg
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Double"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>30", "eta<2.1", "eta>-2.1", "eg_medium"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>20", "eta<2.4", "eta>-2.4", "eg_medium"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(MUON_TIGHT)", "pt>30", "eg_medium"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(MUON_LOOSE)", "pt>20", "eg_medium"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(False)
         ),
     cms.untracked.PSet(
         name = cms.untracked.string("Combined Single"),
-        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>30", "eta<2.1", "eta>-2.1", "eg_tight"),
-        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(ET)", "pt>20", "eta<2.4", "eta>-2.4", "eg_medium"),
+        cuts0 = cms.untracked.vstring("acc(ALL)", "acc(MUON_TIGHT)", "pt>30", "eg_tight"),
+        cuts1 = cms.untracked.vstring("acc(ALL)", "acc(MUON_LOOSE)", "pt>20", "eg_medium"),
         min_mz = cms.untracked.double(MIN_MZ),
         max_mz = cms.untracked.double(MAX_MZ),
         use_truth_mass = cms.untracked.bool(False)

diff --git a/ZFinder/Event/src/AcceptanceSetter.cc b/ZFinder/Event/src/AcceptanceSetter.cc
@@ -42,6 +42,12 @@ namespace zf {
         bool nt = false;
         bool ntm = false;
         bool ntp = false;
+        bool muon_tight = false;
+        bool muon_tight_plus = false;
+        bool muon_tight_minus = false;
+        bool muon_loose = false;
+        bool muon_loose_plus = false;
+        bool muon_loose_minus = false;
         // EB
         if (EB_MIN_ < ETA && ETA < EB_MAX_) {  // In EB
             detector = true;
@@ -86,6 +92,25 @@ namespace zf {
                 hfm = true;
             }
         }
+        // Selection for combination with the muon group, we want to reject in
+        // the gap between EB and EE, so we also require either eb or ee
+        if (FETA < MUON_LOOSE_MAX_ && (eb || ee)) {  // In ML
+            detector = true;
+            muon_loose = true;
+            if (EB_CENT_ < ETA && ETA < MUON_LOOSE_MAX_) {  // In ML+
+                muon_loose_plus = true;
+            } else if (MUON_LOOSE_MIN_ < ETA && ETA < EB_CENT_) {  // In ML-
+                muon_loose_minus = true;
+            }
+            if (FETA < MUON_TIGHT_MAX_) {  // In MT
+                muon_tight = true;
+                if (EB_CENT_ < ETA && ETA < MUON_TIGHT_MAX_) {  // In MT+
+                    muon_tight_plus = true;
+                } else if (MUON_TIGHT_MIN_ < ETA && ETA < EB_CENT_) {  // In MT-
+                    muon_tight_minus = true;
+                }
+            }
+        }
 
         // A special cut to start: ALL is now used to track all events in the
         // sample, which is useful for MC normalization based on the number of
@@ -110,5 +135,12 @@ namespace zf {
         zf_elec->AddCutResult("acc(NT)", nt, WEIGHT);
         zf_elec->AddCutResult("acc(NT+)", ntp, WEIGHT);
         zf_elec->AddCutResult("acc(NT-)", ntm, WEIGHT);
+        zf_elec->AddCutResult("acc(MUON_TIGHT)", muon_tight, WEIGHT);
+        zf_elec->AddCutResult("acc(MUON_TIGHT+)", muon_tight_plus, WEIGHT);
+        zf_elec->AddCutResult("acc(MUON_TIGHT-)", muon_tight_minus, WEIGHT);
+        zf_elec->AddCutResult("acc(MUON_LOOSE)", muon_loose, WEIGHT);
+        zf_elec->AddCutResult("acc(MUON_LOOSE+)", muon_loose_plus, WEIGHT);
+        zf_elec->AddCutResult("acc(MUON_LOOSE-)", muon_loose_minus, WEIGHT);
+
     }
 }  // namespace zf