AliPhysics  ff07904 (ff07904)
AliEmcalTriggerSelectionCuts.cxx
Go to the documentation of this file.
1 /************************************************************************************
2  * Copyright (C) 2017, Copyright Holders of the ALICE Collaboration *
3  * All rights reserved. *
4  * *
5  * Redistribution and use in source and binary forms, with or without *
6  * modification, are permitted provided that the following conditions are met: *
7  * * Redistributions of source code must retain the above copyright *
8  * notice, this list of conditions and the following disclaimer. *
9  * * Redistributions in binary form must reproduce the above copyright *
10  * notice, this list of conditions and the following disclaimer in the *
11  * documentation and/or other materials provided with the distribution. *
12  * * Neither the name of the <organization> nor the *
13  * names of its contributors may be used to endorse or promote products *
14  * derived from this software without specific prior written permission. *
15  * *
16  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND *
17  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED *
18  * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE *
19  * DISCLAIMED. IN NO EVENT SHALL ALICE COLLABORATION BE LIABLE FOR ANY *
20  * DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES *
21  * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; *
22  * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND *
23  * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT *
24  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
25  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
26  ************************************************************************************/
27 #include <iostream>
28 #include <string>
29 #include <map>
31 #include "AliEMCALTriggerPatchInfo.h"
32 
34 
35 namespace PWG {
36 namespace EMCAL {
37 
38 
39 AliEmcalTriggerSelectionCuts::AliEmcalTriggerSelectionCuts() :
40  TObject(),
41  fSelectionMethod(kADC),
42  fPatchType(kAnyPatch),
43  fAcceptanceType(kEMCALDCALAcceptance),
44  fThreshold(0),
45  fUseSimpleOffline(kFALSE),
46  fUseRecalc(kFALSE)
47 {
48 }
49 
50 Bool_t AliEmcalTriggerSelectionCuts::IsSelected(const AliEMCALTriggerPatchInfo * const patch) const {
51  //if(fUseSimpleOffline && !patch->IsOfflineSimple()) return kFALSE;
52  //else if(!fUseSimpleOffline && patch->IsOfflineSimple()) return kFALSE;
53  if(!SelectAcceptance(patch)) return kFALSE;
54  if(!SelectPatchType(patch)) return kFALSE;
55  if(GetCutPrimitive(patch) < fThreshold) return kFALSE;
56  return kTRUE;
57 }
58 
59 Int_t AliEmcalTriggerSelectionCuts::CompareTriggerPatches(const AliEMCALTriggerPatchInfo *first, const AliEMCALTriggerPatchInfo *second) const {
60  Double_t valfirst = GetCutPrimitive(first), valsecond = GetCutPrimitive(second);
61  if(valfirst == valsecond) return 0;
62  if(valfirst > valsecond) return 1;
63  return -1;
64 }
65 
66 Double_t AliEmcalTriggerSelectionCuts::GetCutPrimitive(const AliEMCALTriggerPatchInfo * const patch) const{
67  double energy(0);
68  switch(fSelectionMethod){
69  case kADC: energy = static_cast<Double_t>(patch->GetADCAmp()); break;
70  case kEnergyRough: energy = patch->GetADCAmpGeVRough(); break;
71  case kEnergyOfflineSmeared: energy = patch->GetSmearedEnergy(); break;
72  case kEnergyOffline: energy = patch->GetPatchE(); break;
73  default: energy = -1.;
74  };
75  return energy;
76 }
77 
78 Bool_t AliEmcalTriggerSelectionCuts::SelectPatchType(const AliEMCALTriggerPatchInfo * const patch) const{
79  if(fPatchType == kAnyPatch) return kTRUE;
81  if(patch->IsJetLowSimple() && ((fPatchType == kL1JetPatch) || (fPatchType == kL1JetLowPatch))) return kTRUE;
82  if(patch->IsJetHighSimple() && ((fPatchType == kL1JetPatch) || (fPatchType == kL1JetHighPatch))) return kTRUE;
83  if(patch->IsGammaLowSimple() && ((fPatchType == kL1GammaPatch) || (fPatchType == kL1GammaLowPatch))) return kTRUE;
84  if(patch->IsGammaHighSimple() && ((fPatchType == kL1GammaPatch) || (fPatchType == kL1GammaHighPatch))) return kTRUE;
85  } else if(fUseRecalc) {
86  if(patch->IsJetLowRecalc() && ((fPatchType == kL1JetPatch) || (fPatchType == kL1JetLowPatch))) return kTRUE;
87  if(patch->IsJetHighRecalc() && ((fPatchType == kL1JetPatch) || (fPatchType == kL1JetHighPatch))) return kTRUE;
88  if(patch->IsGammaLowRecalc() && ((fPatchType == kL1GammaPatch) || (fPatchType == kL1GammaLowPatch))) return kTRUE;
89  if(patch->IsGammaHighRecalc() && ((fPatchType == kL1GammaPatch) || (fPatchType == kL1GammaHighPatch))) return kTRUE;
90  } else {
91  if(patch->IsJetLow() && ((fPatchType == kL1JetPatch) || (fPatchType == kL1JetLowPatch))) return kTRUE;
92  if(patch->IsJetHigh() && ((fPatchType == kL1JetPatch) || (fPatchType == kL1JetHighPatch))) return kTRUE;
93  if(patch->IsGammaLow() && ((fPatchType == kL1GammaPatch) || (fPatchType == kL1GammaLowPatch))) return kTRUE;
94  if(patch->IsGammaHigh() && ((fPatchType == kL1GammaPatch) || (fPatchType == kL1GammaHighPatch))) return kTRUE;
95  if(patch->IsLevel0() && fPatchType == kL0Patch) return kTRUE;
96  }
97  return kFALSE;
98 }
99 
100 Bool_t AliEmcalTriggerSelectionCuts::SelectAcceptance(const AliEMCALTriggerPatchInfo * const patch) const {
101  Bool_t selected(false);
102  switch(fAcceptanceType){
103  case kEMCALAcceptance: selected = patch->IsEMCal(); break;
104  case kDCALAcceptance: selected = patch->IsDCalPHOS(); break;
105  case kEMCALDCALAcceptance: selected = patch->IsEMCal() || patch->IsDCalPHOS(); break;
106  default: selected = false;
107  };
108  return selected;
109 }
110 
111 void AliEmcalTriggerSelectionCuts::PrintStream(std::ostream &stream) const {
112  std::map<AcceptanceType_t, std::string> acceptancetext = {{kEMCALAcceptance, "EMCAL"},
113  {kDCALAcceptance, "DCAL"},
114  {kEMCALDCALAcceptance, "EMCAL-DCAL"}};
115  std::map<PatchType_t, std::string> patchtypetext = {{kL0Patch, "Level0"},
116  {kL1GammaPatch, "L1-gamma"},
117  {kL1GammaHighPatch, "L1-gamma, high threshold"},
118  {kL1GammaLowPatch, "L1-gamma, low threshold"},
119  {kL1JetPatch, "L1-jet"},
120  {kL1JetHighPatch, "L1-jet, high threshold"},
121  {kL1JetLowPatch, "L1-jet, low threshold"}};
122  std::map<SelectionMethod_t, std::string> selmodetext = {{kADC, "FastOR ADC"},
123  {kEnergyRough, "FastOR Energy"},
124  {kEnergyOffline, "FEE Energy"},
125  {kEnergyOfflineSmeared, "FEE Energy, decalibrated"}};
126  stream << " Cut settings:" << std::endl;
127  stream << " acceptance: " << acceptancetext.find(fAcceptanceType)->second << std::endl;
128  stream << " patchtype: " << patchtypetext.find(fPatchType)->second << std::endl;
129  stream << " sel mode: " << selmodetext.find(fSelectionMethod)->second << std::endl;
130  stream << " Offline Patches: " << (fUseSimpleOffline ? "yes" : "no") << std::endl;
131  stream << " Recalc Patches: " << (fUseRecalc ? "yes" : "no") << std::endl;
132  stream << " Threshold: " << fThreshold << std::endl;
133 }
134 
135 }
136 }
137 
138 std::ostream &operator<<(std::ostream &stream, const PWG::EMCAL::AliEmcalTriggerSelectionCuts &cuts){
139  cuts.PrintStream(stream);
140  return stream;
141 }
Class for the selection of trigger patches in the EMCAL triggered event selection.
Bool_t SelectPatchType(const AliEMCALTriggerPatchInfo *const patch) const
Select type of the patch according the definitions in the header file.
double Double_t
Definition: External.C:58
Bool_t SelectAcceptance(const AliEMCALTriggerPatchInfo *const patch) const
Select detector acceptance.
energy
Definition: HFPtSpectrum.C:44
Bool_t fUseSimpleOffline
Request simple offline patches.
Namespace for EMCAL framework classes and task.
int Int_t
Definition: External.C:63
Bool_t IsSelected(const AliEMCALTriggerPatchInfo *const patch) const
Apply selection of the given trigger patch according to the selections described in the object...
Namespace for PWG framework classes.
void PrintStream(std::ostream &stream) const
Helper function for output stream operator.
Double_t GetCutPrimitive(const AliEMCALTriggerPatchInfo *const patch) const
Return (energy) measure we cut on, depending on the selection method specified.
friend std::ostream & operator<<(std::ostream &stream, const AliEmcalTriggerSelectionCuts &cuts)
Output stream operator.
SelectionMethod_t fSelectionMethod
Variable to cut on.
Int_t CompareTriggerPatches(const AliEMCALTriggerPatchInfo *first, const AliEMCALTriggerPatchInfo *second) const
Compare two patches according to the energy measure specified in the cut object.
AcceptanceType_t fAcceptanceType
Acceptance type (EMCAL or DCAL acceptance)
bool Bool_t
Definition: External.C:53
PatchType_t fPatchType
Type of the patch to be selected.