CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EgHLTMonElemFuncs.cc
Go to the documentation of this file.
1 #include "DQMOffline/Trigger/interface/EgHLTMonElemFuncs.h"
2 
3 #include "DQMOffline/Trigger/interface/EgHLTMonElemMgrEBEE.h"
4 #include "DQMOffline/Trigger/interface/EgHLTEgCutCodes.h"
5 #include "DQMOffline/Trigger/interface/EgHLTDQMCut.h"
6 #include "DQMOffline/Trigger/interface/EgHLTCutMasks.h"
7 
8 
9 using namespace egHLT;
10 
11 void MonElemFuncs::initStdEleHists(std::vector<MonElemManagerBase<OffEle>*>& histVec,const std::string& filterName,const std::string& baseName,const BinData& bins)
12 {
13  addStdHist<OffEle,float>(histVec,baseName+"_et",baseName+" E_{T};E_{T} (GeV)",bins.et,&OffEle::et);
14  addStdHist<OffEle,float>(histVec,baseName+"_etHigh",baseName+" E_{T};E_{T} (GeV)",bins.etHigh,&OffEle::et);
15  addStdHist<OffEle,float>(histVec,baseName+"_etSC",baseName+" E^{SC}_{T};E^{SC}_{T} (GeV)",bins.et,&OffEle::etSC);
16  addStdHist<OffEle,float>(histVec,baseName+"_eta",baseName+" #eta;#eta",bins.eta,&OffEle::detEta);
17  addStdHist<OffEle,float>(histVec,baseName+"_phi",baseName+" #phi;#phi (rad)",bins.phi,&OffEle::phi);
18  // addStdHist<OffEle,int>(histVec,baseName+"_charge",baseName+" Charge; charge",bins.charge,&OffEle::charge);
19 
20  addStdHist<OffEle,float>(histVec,baseName+"_hOverE",baseName+" H/E; H/E",bins.hOverE,&OffEle::hOverE);
21  //----Morse
22  addStdHist<OffEle,float>(histVec,baseName+"_maxr9",baseName+" MAXR9 ; MAXR9",bins.maxr9,&OffEle::r9);
23  //-------
24  addStdHist<OffEle,float>(histVec,baseName+"_dPhiIn",baseName+" #Delta #phi_{in}; #Delta #phi_{in}",bins.dPhiIn,&OffEle::dPhiIn);
25  addStdHist<OffEle,float>(histVec,baseName+"_dEtaIn",baseName+" #Delta #eta_{in}; #Delta #eta_{in}",bins.dEtaIn,&OffEle::dEtaIn);
26  addStdHist<OffEle,float>(histVec,baseName+"_sigmaIEtaIEta",baseName+"#sigma_{i#etai#eta}; #sigma_{i#etai#eta}",bins.sigEtaEta,&OffEle::sigmaIEtaIEta);
27  addStdHist<OffEle,float>(histVec,baseName+"_epIn",baseName+"E/p_{in}; E/p_{in}",bins.eOverP,&OffEle::epIn);
28  addStdHist<OffEle,float>(histVec,baseName+"_epOut",baseName+"E/p_{out}; E/p_{out}",bins.eOverP,&OffEle::epOut);
29  addStdHist<OffEle,float>(histVec,baseName+"_invEInvP",baseName+"1/E -1/p; 1/E - 1/p",bins.invEInvP,&OffEle::invEInvP);
30 
31  addStdHist<OffEle,float>(histVec,baseName+"_e2x5Over5x5",baseName+"E^{2x5}/E^{5x5}; E^{2x5}/E^{5x5}",bins.e2x5,&OffEle::e2x5MaxOver5x5);
32  addStdHist<OffEle,float>(histVec,baseName+"_e1x5Over5x5",baseName+"E^{1x5}/E^{5x5}; E^{1x5}/E^{5x5}",bins.e1x5,&OffEle::e1x5Over5x5);
33  // addStdHist<OffEle,float>(histVec,baseName+"_isolEM",baseName+"Isol EM; Isol EM (GeV)",bins.isolEm,&OffEle::isolEm);
34  //addStdHist<OffEle,float>(histVec,baseName+"_isolHad",baseName+"Isol Had; Isol Had (GeV)",bins.isolHad,&OffEle::isolHad);
35  //addStdHist<OffEle,float>(histVec,baseName+"_isolPtTrks",baseName+"Isol Pt Trks; Isol Pt Tracks (GeV/c)",bins.isolPtTrks,&OffEle::isolPtTrks);
36  addStdHist<OffEle,float>(histVec,baseName+"_hltIsolTrksEle",baseName+"HLT Ele Isol Trks; HLT Ele Iso Tracks (GeV/c)",bins.isolPtTrks,&OffEle::hltIsolTrksEle);
37  //addStdHist<OffEle,float>(histVec,baseName+"_hltIsolTrksPho",baseName+"HLT Pho Isol Trks; HLT Pho Iso Tracks (GeV/c)",bins.isolPtTrks,&OffEle::hltIsolTrksPho);
38  addStdHist<OffEle,float>(histVec,baseName+"_hltIsolHad",baseName+"HLT Isol Had; HLT Isol Had (GeV)",bins.isolHad,&OffEle::hltIsolHad);
39  addStdHist<OffEle,float>(histVec,baseName+"_hltIsolEm",baseName+"HLT Isol Em; HLT Isol Em (GeV)",bins.isolEm,&OffEle::hltIsolEm);
40 
41  histVec.push_back(new MonElemManager2D<OffEle,float,float>(baseName+"_etaVsPhi",
42  baseName+" #eta vs #phi;#eta;#phi (rad)",
43  bins.etaVsPhi.nrX,bins.etaVsPhi.xMin,bins.etaVsPhi.xMax,
44  bins.etaVsPhi.nrY,bins.etaVsPhi.yMin,bins.etaVsPhi.yMax,
45  &OffEle::detEta,&OffEle::phi));
46 
47 }
48 
49 void MonElemFuncs::initStdPhoHists(std::vector<MonElemManagerBase<OffPho>*>& histVec,const std::string& filterName,const std::string& baseName,const BinData& bins)
50 {
51  addStdHist<OffPho,float>(histVec,baseName+"_et",baseName+" E_{T};E_{T} (GeV)",bins.et,&OffPho::et);
52  addStdHist<OffPho,float>(histVec,baseName+"_etHigh",baseName+" E_{T};E_{T} (GeV)",bins.etHigh,&OffPho::et);
53  addStdHist<OffPho,float>(histVec,baseName+"_etSC",baseName+" E^{SC}_{T};E^{SC}_{T} (GeV)",bins.et,&OffPho::etSC);
54  addStdHist<OffPho,float>(histVec,baseName+"_eta",baseName+" #eta;#eta",bins.eta,&OffPho::detEta);
55  addStdHist<OffPho,float>(histVec,baseName+"_phi",baseName+" #phi;#phi (rad)",bins.phi,&OffPho::phi);
56 
57  addStdHist<OffPho,float>(histVec,baseName+"_hOverE",baseName+" H/E; H/E",bins.hOverE,&OffPho::hOverE);
58  //----Morse
59  //addStdHist<OffPho,float>(histVec,baseName+"_r9",baseName+" R9 ; R9",bins.r9,&OffPho::r9);
60  //addStdHist<OffPho,float>(histVec,baseName+"_minr9",baseName+" MINR9 ; MINR9",bins.minr9,&OffPho::r9);
61  addStdHist<OffPho,float>(histVec,baseName+"_maxr9",baseName+" MAXR9 ; MAXR9",bins.maxr9,&OffPho::r9);
62  //-------
63  addStdHist<OffPho,float>(histVec,baseName+"_sigmaIEtaIEta",baseName+"#sigma_{i#etai#eta}; #sigma_{i#etai#eta}",bins.sigEtaEta,&OffPho::sigmaIEtaIEta);
64  addStdHist<OffPho,float>(histVec,baseName+"_e2x5Over5x5",baseName+"E^{2x5}/E^{5x5}; E^{2x5}/E^{5x5}",bins.e2x5,&OffPho::e2x5MaxOver5x5);
65  addStdHist<OffPho,float>(histVec,baseName+"_e1x5Over5x5",baseName+"E^{1x5}/E^{5x5}; E^{1x5}/E^{5x5}",bins.e1x5,&OffPho::e1x5Over5x5);
66  addStdHist<OffPho,float>(histVec,baseName+"_isolEM",baseName+"Isol EM; Isol EM (GeV)",bins.isolEm,&OffPho::isolEm);
67  addStdHist<OffPho,float>(histVec,baseName+"_isolHad",baseName+"Isol Had; Isol Had (GeV)",bins.isolHad,&OffPho::isolHad);
68  addStdHist<OffPho,float>(histVec,baseName+"_isolPtTrks",baseName+"Isol Pt Trks; Isol Pt Tracks (GeV/c)",bins.isolPtTrks,&OffPho::isolPtTrks);
69  addStdHist<OffPho,int>(histVec,baseName+"_isolNrTrks",baseName+"Isol Nr Trks; Isol Nr Tracks",bins.isolNrTrks,&OffPho::isolNrTrks);
70  // addStdHist<OffPho,float>(histVec,baseName+"_hltIsolTrks",baseName+"HLT Isol Trks; HLT Iso Tracks (GeV/c)",bins.isolPtTrks,&OffPho::hltIsolTrks);
71  //addStdHist<OffPho,float>(histVec,baseName+"_hltIsolHad",baseName+"HLT Isol Had; HLT Isol Had (GeV)",bins.isolPtTrks,&OffPho::hltIsolHad);
72 
73  histVec.push_back(new MonElemManager2D<OffPho,float,float>(baseName+"_etaVsPhi",
74  baseName+" #eta vs #phi;#eta;#phi (rad)",
75  bins.etaVsPhi.nrX,bins.etaVsPhi.xMin,bins.etaVsPhi.xMax,
76  bins.etaVsPhi.nrY,bins.etaVsPhi.yMin,bins.etaVsPhi.yMax,
77  &OffPho::detEta,&OffPho::phi));
78 }
79 
80 void MonElemFuncs::initStdEffHists(std::vector<MonElemWithCutBase<OffEle>*>& histVec,const std::string& filterName,const std::string& baseName,const BinData::Data1D& bins,float (OffEle::*vsVarFunc)()const,const CutMasks& masks)
81 {
82  initStdEffHists(histVec,filterName,baseName,bins.nr,bins.min,bins.max,vsVarFunc,masks);
83 }
84 
85 void MonElemFuncs::initStdEffHists(std::vector<MonElemWithCutBase<OffPho>*>& histVec,const std::string& filterName,const std::string& baseName,const BinData::Data1D& bins,float (OffPho::*vsVarFunc)()const,const CutMasks& masks)
86 {
87  initStdEffHists(histVec,filterName,baseName,bins.nr,bins.min,bins.max,vsVarFunc,masks);
88 }
89 
90 void MonElemFuncs::initStdEffHists(std::vector<MonElemWithCutBase<OffEle>*>& histVec,const std::string& filterName,const std::string& baseName,int nrBins,double xMin,double xMax,float (OffEle::*vsVarFunc)()const,const CutMasks& masks)
91 {
92  //some convience typedefs, I hate typedefs but atleast here where they are defined is obvious
93  typedef EgHLTDQMVarCut<OffEle> VarCut;
94  typedef MonElemWithCutEBEE<OffEle,float> MonElemFloat;
95  int stdCutCode = masks.stdEle;
96 
97  //first do the zero and all cuts histograms
98  histVec.push_back(new MonElemFloat(baseName+"_noCuts",baseName+" NoCuts",nrBins,xMin,xMax,vsVarFunc));
99  histVec.push_back(new MonElemFloat(baseName+"_allCuts",baseName+" All Cuts",nrBins,xMin,xMax,vsVarFunc,
100  new VarCut(stdCutCode,&OffEle::cutCode)));
101 
102  //now for the n-1
103  histVec.push_back(new MonElemFloat(baseName+"_n1_dEtaIn",baseName+" N1 #Delta#eta_{in}",nrBins,xMin,xMax,vsVarFunc,
104  new VarCut(~EgCutCodes::DETAIN&stdCutCode,&OffEle::cutCode)));
105  histVec.push_back(new MonElemFloat(baseName+"_n1_dPhiIn",baseName+" N1 #Delta#phi_{in}",nrBins,xMin,xMax,vsVarFunc,
106  new VarCut(~EgCutCodes::DPHIIN&stdCutCode,&OffEle::cutCode)));
107  histVec.push_back(new MonElemFloat(baseName+"_n1_sigmaIEtaIEta",baseName+" N1 #sigma_{i#etai#eta}",nrBins,xMin,xMax,vsVarFunc,
108  new VarCut(~EgCutCodes::SIGMAIETAIETA&stdCutCode,&OffEle::cutCode)));
109  histVec.push_back(new MonElemFloat(baseName+"_n1_hOverE",baseName+" N1 H/E",nrBins,xMin,xMax,vsVarFunc,
110  new VarCut(~EgCutCodes::HADEM&stdCutCode,&OffEle::cutCode)));
111  /* histVec.push_back(new MonElemFloat(baseName+"_n1_isolEm",baseName+" N1 Isol Em",nrBins,xMin,xMax,vsVarFunc,
112  new VarCut(~EgCutCodes::ISOLEM&stdCutCode,&OffEle::cutCode)));
113  histVec.push_back(new MonElemFloat(baseName+"_n1_isolHad",baseName+" N1 Isol Had",nrBins,xMin,xMax,vsVarFunc,
114  new VarCut(~EgCutCodes::ISOLHAD&stdCutCode,&OffEle::cutCode)));
115  histVec.push_back(new MonElemFloat(baseName+"_n1_isolPtTrks",baseName+" N1 Isol Tracks",nrBins,xMin,xMax,vsVarFunc,
116  new VarCut(~EgCutCodes::ISOLPTTRKS&stdCutCode,&OffEle::cutCode)));*/
117  histVec.push_back(new MonElemFloat(baseName+"_n1_hltIsolEm",baseName+" N1 HLT Isol Em",nrBins,xMin,xMax,vsVarFunc,
118  new VarCut(~EgCutCodes::HLTISOLEM&stdCutCode,&OffEle::cutCode)));
119  histVec.push_back(new MonElemFloat(baseName+"_n1_hltIsolHad",baseName+" N1 HLT Isol Had",nrBins,xMin,xMax,vsVarFunc,
120  new VarCut(~EgCutCodes::HLTISOLHAD&stdCutCode,&OffEle::cutCode)));
121  histVec.push_back(new MonElemFloat(baseName+"_n1_hltIsolTrksEle",baseName+" N1 HLT Isol Tracks Ele ",nrBins,xMin,xMax,vsVarFunc,
122  new VarCut(~EgCutCodes::HLTISOLTRKSELE&stdCutCode,&OffEle::cutCode)));
123  histVec.push_back(new MonElemFloat(baseName+"_single_dEtaIn",baseName+" Single #Delta#eta_{in}",nrBins,xMin,xMax,vsVarFunc,
124  new VarCut(EgCutCodes::DETAIN,&OffEle::cutCode)));
125  histVec.push_back(new MonElemFloat(baseName+"_single_dPhiIn",baseName+" Single #Delta#phi_{in}",nrBins,xMin,xMax,vsVarFunc,
126  new VarCut(EgCutCodes::DPHIIN,&OffEle::cutCode)));
127  histVec.push_back(new MonElemFloat(baseName+"_single_sigmaIEtaIEta",baseName+" Single #sigma_{i#etai#eta}",nrBins,xMin,xMax,vsVarFunc,
128  new VarCut(EgCutCodes::SIGMAIETAIETA,&OffEle::cutCode)));
129  histVec.push_back(new MonElemFloat(baseName+"_single_hOverE",baseName+" Single H/E",nrBins,xMin,xMax,vsVarFunc,
130  new VarCut(EgCutCodes::HADEM,&OffEle::cutCode)));
131  /* histVec.push_back(new MonElemFloat(baseName+"_single_isolEm",baseName+" Single Isol Em",nrBins,xMin,xMax,vsVarFunc,
132  new VarCut(EgCutCodes::ISOLEM,&OffEle::cutCode)));
133  histVec.push_back(new MonElemFloat(baseName+"_single_isolHad",baseName+" Single Isol Had",nrBins,xMin,xMax,vsVarFunc,
134  new VarCut(EgCutCodes::ISOLHAD,&OffEle::cutCode)));
135  histVec.push_back(new MonElemFloat(baseName+"_single_isolPtTrks",baseName+" Single Isol Tracks",nrBins,xMin,xMax,vsVarFunc,
136  new VarCut(EgCutCodes::ISOLPTTRKS,&OffEle::cutCode)));*/
137  histVec.push_back(new MonElemFloat(baseName+"_single_hltIsolEm",baseName+" Single HLT Isol Em",nrBins,xMin,xMax,vsVarFunc,
138  new VarCut(EgCutCodes::HLTISOLEM,&OffEle::cutCode)));
139  histVec.push_back(new MonElemFloat(baseName+"_single_hltIsolHad",baseName+" Single HLT Isol Had",nrBins,xMin,xMax,vsVarFunc,
140  new VarCut(EgCutCodes::HLTISOLHAD,&OffEle::cutCode)));
141  histVec.push_back(new MonElemFloat(baseName+"_single_hltIsolTrksEle",baseName+" Single HLT Isol Tracks Ele ",nrBins,xMin,xMax,vsVarFunc,
142  new VarCut(EgCutCodes::HLTISOLTRKSELE,&OffEle::cutCode)));
143  /*histVec.push_back(new MonElemFloat(baseName+"_single_hltIsolTrksPho",baseName+" Single HLT Isol Tracks Pho ",nrBins,xMin,xMax,vsVarFunc,
144  new VarCut(EgCutCodes::HLTISOLTRKSPHO,&OffEle::cutCode)));*/
145 
146 }
147 
148 void MonElemFuncs::initStdEffHists(std::vector<MonElemWithCutBase<OffPho>*>& histVec,const std::string& filterName,const std::string& baseName,int nrBins,double xMin,double xMax,float (OffPho::*vsVarFunc)()const,const CutMasks& masks)
149 {
150  //some convenience typedefs, I hate typedefs but atleast here where they are defined is obvious
151  typedef EgHLTDQMVarCut<OffPho> VarCut;
152  typedef MonElemWithCutEBEE<OffPho,float> MonElemFloat;
153  int stdCutCode = masks.stdPho;
154 
155  histVec.push_back(new MonElemFloat(baseName+"_noCuts",baseName+" NoCuts",nrBins,xMin,xMax,vsVarFunc));
156  histVec.push_back(new MonElemFloat(baseName+"_allCuts",baseName+" All Cuts",nrBins,xMin,xMax,vsVarFunc,
157  new VarCut(stdCutCode,&OffPho::cutCode)));
158 
159 
160  histVec.push_back(new MonElemFloat(baseName+"_n1_sigmaIEtaIEta",baseName+" N1 #sigma_{i#etai#eta}",nrBins,xMin,xMax,vsVarFunc,
161  new VarCut(~EgCutCodes::SIGMAIETAIETA&stdCutCode,&OffPho::cutCode)));
162  //-----Morse------
163  histVec.push_back(new MonElemFloat(baseName+"_n1_hOverE",baseName+" N1 H/E",nrBins,xMin,xMax,vsVarFunc,
164  new VarCut(~EgCutCodes::HADEM&stdCutCode,&OffPho::cutCode)));//---BUG FIX!!--
165  /*histVec.push_back(new MonElemFloat(baseName+"_n1_minr9",baseName+" N1 MINR9",nrBins,xMin,xMax,vsVarFunc,
166  new VarCut(~EgCutCodes::MINR9&stdCutCode,&OffPho::cutCode)));
167  histVec.push_back(new MonElemFloat(baseName+"_n1_maxr9",baseName+" N1 MAXR9",nrBins,xMin,xMax,vsVarFunc,
168  new VarCut(~EgCutCodes::MAXR9&stdCutCode,&OffPho::cutCode)));*/
169  //----------------
170  histVec.push_back(new MonElemFloat(baseName+"_n1_isolEm",baseName+" N1 Isol Em",nrBins,xMin,xMax,vsVarFunc,
171  new VarCut(~EgCutCodes::ISOLEM&stdCutCode,&OffPho::cutCode)));
172  histVec.push_back(new MonElemFloat(baseName+"_n1_isolHad",baseName+" N1 Isol Had",nrBins,xMin,xMax,vsVarFunc,
173  new VarCut(~EgCutCodes::ISOLHAD&stdCutCode,&OffPho::cutCode)));
174  histVec.push_back(new MonElemFloat(baseName+"_n1_isolPtTrks",baseName+" N1 Pt Isol Tracks",nrBins,xMin,xMax,vsVarFunc,
175  new VarCut(~EgCutCodes::ISOLPTTRKS&stdCutCode,&OffPho::cutCode)));
176  histVec.push_back(new MonElemFloat(baseName+"_n1_isolNrTrks",baseName+" N1 Nr Isol Tracks",nrBins,xMin,xMax,vsVarFunc,
177  new VarCut(~EgCutCodes::ISOLNRTRKS&stdCutCode,&OffPho::cutCode)));
178 
179  histVec.push_back(new MonElemFloat(baseName+"_single_hOverE",baseName+" Single H/E",nrBins,xMin,xMax,vsVarFunc,
180  new VarCut(EgCutCodes::HADEM,&OffPho::cutCode)));
181  histVec.push_back(new MonElemFloat(baseName+"_single_sigmaIEtaIEta",baseName+" Single #sigma_{i#etai#eta}",nrBins,xMin,xMax,vsVarFunc,
182  new VarCut(EgCutCodes::SIGMAIETAIETA,&OffPho::cutCode)));
183  histVec.push_back(new MonElemFloat(baseName+"_single_isolEm",baseName+" Single Isol Em",nrBins,xMin,xMax,vsVarFunc,
184  new VarCut(~EgCutCodes::ISOLEM,&OffPho::cutCode)));
185  histVec.push_back(new MonElemFloat(baseName+"_single_isolHad",baseName+" Single Isol Had",nrBins,xMin,xMax,vsVarFunc,
186  new VarCut(~EgCutCodes::ISOLHAD,&OffPho::cutCode)));
187  histVec.push_back(new MonElemFloat(baseName+"_single_isolPtTrks",baseName+" Single Pt Isol Tracks",nrBins,xMin,xMax,vsVarFunc,
188  new VarCut(~EgCutCodes::ISOLPTTRKS,&OffPho::cutCode)));
189  /*histVec.push_back(new MonElemFloat(baseName+"_single_hltIsolHad",baseName+" Single HLT Isol Had",nrBins,xMin,xMax,vsVarFunc,
190  new VarCut(EgCutCodes::HLTISOLHAD,&OffPho::cutCode)));
191  histVec.push_back(new MonElemFloat(baseName+"_single_hltIsolTrksPho",baseName+" Single HLT Isol Tracks Pho ",nrBins,xMin,xMax,vsVarFunc,
192  new VarCut(EgCutCodes::HLTISOLTRKSPHO,&OffPho::cutCode)));*/
193 
194 
195 }
196 
197 
198 //we own the passed in cut, so we give it to the first mon element and then clone it after that
199 //only currently used for trigger tag and probe
200 void MonElemFuncs::initStdEleCutHists(std::vector<MonElemWithCutBase<OffEle>*>& histVec,const std::string& filterName,const std::string& baseName,const BinData& bins,EgHLTDQMCut<OffEle>* cut)
201 {
202  histVec.push_back(new MonElemWithCutEBEE<OffEle,float>(baseName+"_et",
203  baseName+" E_{T};E_{T} (GeV)",
204  bins.et.nr,bins.et.min,bins.et.max,&OffEle::et,cut));
205  histVec.push_back(new MonElemWithCutEBEE<OffEle,float>(baseName+"_eta",
206  baseName+" #eta;#eta",
207  bins.eta.nr,bins.eta.min,bins.eta.max,
208  &OffEle::detEta,cut ? cut->clone(): NULL));
209  histVec.push_back(new MonElemWithCutEBEE<OffEle,float>(baseName+"_phi",
210  baseName+" #phi;#phi (rad)",
211  bins.phi.nr,bins.phi.min,bins.phi.max,
212  &OffEle::phi,cut ? cut->clone():NULL));
213  histVec.push_back(new MonElemWithCutEBEE<OffEle,int>(baseName+"_nVertex",
214  baseName+" nVertex;nVertex",
215  bins.nVertex.nr,bins.nVertex.min,bins.nVertex.max,
216  &OffEle::NVertex,cut ? cut->clone():NULL));
217  /* histVec.push_back(new MonElemWithCutEBEE<OffEle,int>(baseName+"_charge",
218  baseName+" Charge; charge",
219  bins.charge.nr,bins.charge.min,bins.charge.max,
220  &OffEle::charge,cut ? cut->clone():NULL)); */
221 }
222 
223 
224 void MonElemFuncs::initStdPhoCutHists(std::vector<MonElemWithCutBase<OffPho>*>& histVec,const std::string& filterName,const std::string& baseName,const BinData& bins,EgHLTDQMCut<OffPho>* cut)
225 {
226  histVec.push_back(new MonElemWithCutEBEE<OffPho,float>(baseName+"_et",
227  baseName+" E_{T};E_{T} (GeV)",
228  bins.et.nr,bins.et.min,bins.et.max,&OffPho::et,cut));
229  histVec.push_back(new MonElemWithCutEBEE<OffPho,float>(baseName+"_eta",
230  baseName+" #eta;#eta",
231  bins.eta.nr,bins.eta.min,bins.eta.max,
232  &OffPho::detEta,cut ? cut->clone(): NULL));
233  histVec.push_back(new MonElemWithCutEBEE<OffPho,float>(baseName+"_phi",
234  baseName+" #phi;#phi (rad)",
235  bins.phi.nr,bins.phi.min,bins.phi.max,
236  &OffPho::phi,cut ? cut->clone():NULL));
237  /* histVec.push_back(new MonElemWithCutEBEE<OffPho,int>(baseName+"_charge",
238  baseName+" Charge; charge",
239  bins.charge.nr,bins.charge.min,bins.charge.max,
240  &OffPho::charge,cut ? cut->clone():NULL)); */
241 }
242 
243 //we transfer ownership of eleCut to addTrigLooseTrigHist which transfers it to the eleMonElems
244 void MonElemFuncs::initTightLooseTrigHists(std::vector<MonElemContainer<OffEle>*>& eleMonElems,const std::vector<std::string>& tightLooseTrigs,const BinData& bins,EgHLTDQMCut<OffEle>* eleCut)
245 {
246  for(size_t trigNr=0;trigNr<tightLooseTrigs.size();trigNr++){
247  std::vector<std::string> splitString;
248  boost::split(splitString,tightLooseTrigs[trigNr],boost::is_any_of(std::string(":")));
249  if(splitString.size()!=2) continue; //format incorrect
250  const std::string& tightTrig = splitString[0];
251  const std::string& looseTrig = splitString[1];
252  //this step is necessary as we want to transfer ownership of eleCut to the addTrigLooseTrigHist func on the last iteration
253  //but clone it before that
254  //perhaps my object ownership rules need to be re-evalulated
255  if(trigNr!=tightLooseTrigs.size()-2) addTightLooseTrigHist(eleMonElems,tightTrig,looseTrig,eleCut->clone(),"gsfEle",bins);
256  else addTightLooseTrigHist(eleMonElems,tightTrig,looseTrig,eleCut,"gsfEle",bins);
257  }
258 }
259 
260 
261 void MonElemFuncs::initTightLooseTrigHists(std::vector<MonElemContainer<OffPho>*>& phoMonElems,const std::vector<std::string>& tightLooseTrigs,const BinData& bins,EgHLTDQMCut<OffPho>* phoCut)
262 {
263  for(size_t trigNr=0;trigNr<tightLooseTrigs.size();trigNr++){
264  std::vector<std::string> splitString;
265  boost::split(splitString,tightLooseTrigs[trigNr],boost::is_any_of(std::string(":")));
266  if(splitString.size()!=2) continue; //format incorrect
267  const std::string& tightTrig = splitString[0];
268  const std::string& looseTrig = splitString[1];
269 
270  //this step is necessary as we want to transfer ownership of phoCut to the addTrigLooseTrigHist func on the last iteration
271  //but clone it before that
272  //perhaps my object ownership rules need to be re-evalulated
273  if(trigNr!=tightLooseTrigs.size()-2) addTightLooseTrigHist(phoMonElems,tightTrig,looseTrig,phoCut->clone(),"pho",bins);
274  else addTightLooseTrigHist(phoMonElems,tightTrig,looseTrig,phoCut,"pho",bins);
275  }
276 }
277 
278 
279 
280 
281 //there is nothing electron specific here, will template at some point
282 void MonElemFuncs::addTightLooseTrigHist(std::vector<MonElemContainer<OffEle>*>& eleMonElems,
283  const std::string& tightTrig,const std::string& looseTrig,
284  EgHLTDQMCut<OffEle>* eleCut,
285  const std::string& histId,const BinData& bins)
286 {
287  MonElemContainer<OffEle>* passMonElem = NULL;
288  passMonElem = new MonElemContainer<OffEle>(tightTrig+"_"+looseTrig+"_"+histId+"_passTrig","",
289  &(*(new EgMultiCut<OffEle>) <<
290  new EgObjTrigCut<OffEle>(TrigCodes::getCode(tightTrig+":"+looseTrig),EgObjTrigCut<OffEle>::AND) <<
291  eleCut->clone()));
292 
293 
294  MonElemContainer<OffEle>* failMonElem = NULL;
295  failMonElem = new MonElemContainer<OffEle>(tightTrig+"_"+looseTrig+"_"+histId+"_failTrig","",
296  &(*(new EgMultiCut<OffEle>) <<
297  new EgObjTrigCut<OffEle>(TrigCodes::getCode(looseTrig),EgObjTrigCut<OffEle>::AND,TrigCodes::getCode(tightTrig)) <<
298  eleCut));
299 
300  MonElemFuncs::initStdEleHists(passMonElem->monElems(),tightTrig+"_"+looseTrig,passMonElem->name(),bins);
301  MonElemFuncs::initStdEleHists(failMonElem->monElems(),tightTrig+"_"+looseTrig,failMonElem->name(),bins);
302  eleMonElems.push_back(passMonElem);
303  eleMonElems.push_back(failMonElem);
304 }
305 
306 
307 //there is nothing photon specific here, will template at some point
308 void MonElemFuncs::addTightLooseTrigHist(std::vector<MonElemContainer<OffPho>*>& phoMonElems,
309  const std::string& tightTrig,const std::string& looseTrig,
310  EgHLTDQMCut<OffPho>* phoCut,
311  const std::string& histId,const BinData& bins)
312 {
313  MonElemContainer<OffPho>* passMonElem = NULL;
314  passMonElem = new MonElemContainer<OffPho>(tightTrig+"_"+looseTrig+"_"+histId+"_passTrig","",
315  &(*(new EgMultiCut<OffPho>) <<
316  new EgObjTrigCut<OffPho>(TrigCodes::getCode(tightTrig+":"+looseTrig),EgObjTrigCut<OffPho>::AND) <<
317  phoCut->clone()));
318 
319 
320  MonElemContainer<OffPho>* failMonElem = NULL;
321  failMonElem = new MonElemContainer<OffPho>(tightTrig+"_"+looseTrig+"_"+histId+"_failTrig","",
322  &(*(new EgMultiCut<OffPho>) <<
323  new EgObjTrigCut<OffPho>(TrigCodes::getCode(looseTrig),EgObjTrigCut<OffPho>::AND,TrigCodes::getCode(tightTrig)) <<
324  phoCut));
325 
326  MonElemFuncs::initStdPhoHists(passMonElem->monElems(),tightTrig+"_"+looseTrig,passMonElem->name(),bins);
327  MonElemFuncs::initStdPhoHists(failMonElem->monElems(),tightTrig+"_"+looseTrig,failMonElem->name(),bins);
328  phoMonElems.push_back(passMonElem);
329  phoMonElems.push_back(failMonElem);
330 }
331 
332 
333 
334 //we transfer ownership of eleCut to the monitor elements
335 void MonElemFuncs::initTightLooseTrigHistsTrigCuts(std::vector<MonElemContainer<OffEle>*>& eleMonElems,const std::vector<std::string>& tightLooseTrigs,const BinData& bins)
336 {
337  for(size_t trigNr=0;trigNr<tightLooseTrigs.size();trigNr++){
338  std::vector<std::string> splitString;
339  boost::split(splitString,tightLooseTrigs[trigNr],boost::is_any_of(std::string(":")));
340  if(splitString.size()!=2) continue; //format incorrect
341  const std::string& tightTrig = splitString[0];
342  const std::string& looseTrig = splitString[1];
343  EgHLTDQMCut<OffEle>* eleCut = new EgHLTDQMUserVarCut<OffEle,TrigCodes::TrigBitSet>(&OffEle::trigCutsCutCode,TrigCodes::getCode(tightTrig));
344  addTightLooseTrigHist(eleMonElems,tightTrig,looseTrig,eleCut,"gsfEle_trigCuts",bins);
345  }
346 }
347 
348 //we transfer ownership of phoCut to the monitor elements
349 void MonElemFuncs::initTightLooseTrigHistsTrigCuts(std::vector<MonElemContainer<OffPho>*>& phoMonElems,const std::vector<std::string>& tightLooseTrigs,const BinData& bins)
350 {
351  for(size_t trigNr=0;trigNr<tightLooseTrigs.size();trigNr++){
352  std::vector<std::string> splitString;
353  boost::split(splitString,tightLooseTrigs[trigNr],boost::is_any_of(std::string(":")));
354  if(splitString.size()!=2) continue; //format incorrect
355  const std::string& tightTrig = splitString[0];
356  const std::string& looseTrig = splitString[1];
357  EgHLTDQMCut<OffPho>* phoCut = new EgHLTDQMUserVarCut<OffPho,TrigCodes::TrigBitSet>(&OffPho::trigCutsCutCode,TrigCodes::getCode(tightTrig));
358  addTightLooseTrigHist(phoMonElems,tightTrig,looseTrig,phoCut,"pho_trigCuts",bins);
359  }
360 }
361 
362 
363 //we transfer ownership of eleCut to the monitor elements
364 void MonElemFuncs::initTightLooseDiObjTrigHistsTrigCuts(std::vector<MonElemContainer<OffEle>*>& eleMonElems,const std::vector<std::string>& tightLooseTrigs,const BinData& bins)
365 {
366  for(size_t trigNr=0;trigNr<tightLooseTrigs.size();trigNr++){
367  std::vector<std::string> splitString;
368  boost::split(splitString,tightLooseTrigs[trigNr],boost::is_any_of(std::string(":")));
369  if(splitString.size()!=2) continue; //format incorrect
370  const std::string& tightTrig = splitString[0];
371  const std::string& looseTrig = splitString[1];
372  EgHLTDQMCut<OffEle>* eleCut = new EgDiEleUserCut<TrigCodes::TrigBitSet>(&OffEle::trigCutsCutCode,TrigCodes::getCode(tightTrig));
373  addTightLooseTrigHist(eleMonElems,tightTrig,looseTrig,eleCut,"gsfEle_trigCuts",bins);
374  }
375 }
376 
377 
378 //we transfer ownership of phoCut to the monitor elements
379 void MonElemFuncs::initTightLooseDiObjTrigHistsTrigCuts(std::vector<MonElemContainer<OffPho>*>& phoMonElems,const std::vector<std::string>& tightLooseTrigs,const BinData& bins)
380 {
381  for(size_t trigNr=0;trigNr<tightLooseTrigs.size();trigNr++){
382  std::vector<std::string> splitString;
383  boost::split(splitString,tightLooseTrigs[trigNr],boost::is_any_of(std::string(":")));
384  if(splitString.size()!=2) continue; //format incorrect
385  const std::string& tightTrig = splitString[0];
386  const std::string& looseTrig = splitString[1];
387  EgHLTDQMCut<OffPho>* phoCut = new EgDiPhoUserCut<TrigCodes::TrigBitSet>(&OffPho::trigCutsCutCode,TrigCodes::getCode(tightTrig));
388  addTightLooseTrigHist(phoMonElems,tightTrig,looseTrig,phoCut,"pho_trigCuts",bins);
389  }
390 }
391 
392 
393 
394 //tag and probe trigger efficiencies
395 //this is to measure the trigger efficiency with respect to a fully selected offline electron
396 //using a tag and probe technique (note: this will be different to the trigger efficiency normally calculated)
397 void MonElemFuncs::initTrigTagProbeHists(std::vector<MonElemContainer<OffEle>*>& eleMonElems,const std::vector<std::string> filterNames,int cutMask,const BinData& bins)
398 {
399  for(size_t filterNr=0;filterNr<filterNames.size();filterNr++){
400 
401  std::string trigName(filterNames[filterNr]);
402  // float etCutValue = trigTools::getSecondEtThresFromName(trigName);
403  float etCutValue = 0.;
404  //std::cout<<"TrigName= "<<trigName<<" etCutValue= "<<etCutValue<<std::endl;
405  MonElemContainer<OffEle>* monElemCont = new MonElemContainer<OffEle>("trigTagProbe","Trigger Tag and Probe",new EgTrigTagProbeCut_New(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFiltesr"),cutMask,&OffEle::cutCode));
406  //this is all that pass trigtagprobecut
407  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_all",bins,new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC));
408  //this is all that pass trigtagprobecut and the probe passes the test trigger
409  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_pass",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND)));
410  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is NOT a tag
411  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_passNotTag",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND,TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffEle>::AND)));
412  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is also a tag
413  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_passTagTag",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND) << new EgObjTrigCut<OffEle>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffEle>::AND) ));
414  //this is all that pass trigtagprobecut and the probe fails the test trigger
415  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_fail",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),EgObjTrigCut<OffEle>::AND,TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND)));
416  /*
417  monElemCont->monElems().push_back(new MonElemMgrEBEE<OffEle,float>(trigName+"_"+monElemCont->name()+"_gsfEle_all_etUnCut",monElemCont->name()+"_gsfEle_all E_{T} (Uncut);E_{T} (GeV)",
418  bins.et.nr,bins.et.min,bins.et.max,&OffEle::et));
419  monElemCont->cutMonElems().push_back(new MonElemWithCutEBEE<OffEle,float>(trigName+"_"+monElemCont->name()+"_gsfEle_pass_etUnCut",monElemCont->name()+"_gsfEle_pass E_{T} (Uncut);E_{T} (GeV)",
420  bins.et.nr,bins.et.min,bins.et.max,&OffEle::et,new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND)));*/
421  eleMonElems.push_back(monElemCont);
422  } //end filter names loop
423 
424 }
425 
426 //Only one at a time so I can set the folder
427 void MonElemFuncs::initTrigTagProbeHist(std::vector<MonElemContainer<OffEle>*>& eleMonElems,const std::string filterName,int cutMask,const BinData& bins)
428 {
429  std::string trigName(filterName);
430  //float etCutValue = 1.1*trigTools::getSecondEtThresFromName(filterName);
431  float etCutValue = 0.;
432  //std::cout<<"TrigName= "<<trigName<<" etCutValue= "<<etCutValue<<std::endl;
433  MonElemContainer<OffEle>* monElemCont = new MonElemContainer<OffEle>("trigTagProbe","Trigger Tag and Probe",new EgTrigTagProbeCut_New(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),cutMask,&OffEle::cutCode));
434  //this is all that pass trigtagprobecut
435  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_all",bins,new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC));
436  //this is all that pass trigtagprobecut and the probe passes the test trigger
437  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_pass",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND)));
438  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is NOT a tag
439  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_passNotTag",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND,TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffEle>::AND)));
440  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is also a tag
441  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_passTagTag",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND) << new EgObjTrigCut<OffEle>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffEle>::AND) ));
442  //this is all that pass trigtagprobecut and the probe fails the test trigger
443  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_gsfEle_fail",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),EgObjTrigCut<OffEle>::AND,TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND)));
444  /*
445  monElemCont->monElems().push_back(new MonElemMgrEBEE<OffEle,float>(trigName+"_"+monElemCont->name()+"_gsfEle_all_etUnCut",monElemCont->name()+"_gsfEle_all E_{T} (Uncut);E_{T} (GeV)",
446  bins.et.nr,bins.et.min,bins.et.max,&OffEle::et));
447  monElemCont->cutMonElems().push_back(new MonElemWithCutEBEE<OffEle,float>(trigName+"_"+monElemCont->name()+"_gsfEle_pass_etUnCut",monElemCont->name()+"_gsfEle_pass E_{T} (Uncut);E_{T} (GeV)",
448  bins.et.nr,bins.et.min,bins.et.max,&OffEle::et,new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigName),EgObjTrigCut<OffEle>::AND)));*/
449  eleMonElems.push_back(monElemCont);
450 }
451 
452 
453 void MonElemFuncs::initTrigTagProbeHist_2Leg(std::vector<MonElemContainer<OffEle>*>& eleMonElems,const std::string filterName,int cutMask,const BinData& bins)
454 {
455 
456  std::string trigNameLeg1 = filterName.substr(0,filterName.find("::"));
457  std::string trigNameLeg2 = filterName.substr(filterName.find("::")+2);
458 
459  float etCutValue = 0.;
460  MonElemContainer<OffEle>* monElemCont = new MonElemContainer<OffEle>("trigTagProbe","Trigger Tag and Probe",new EgTrigTagProbeCut_New(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),cutMask,&OffEle::cutCode));
461  //this is all that pass trigtagprobecut
462  //MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigNameLeg2,trigNameLeg2+"_"+monElemCont->name()+"_gsfEle_allEt20",bins,new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC));
463  //this is all that pass trigtagprobecut and the probe passes the first trigger
464  //MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigNameLeg2,trigNameLeg2+"_"+monElemCont->name()+"_gsfEle_passEt20",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigNameLeg1),EgObjTrigCut<OffEle>::AND)));
465  //this is all that pass trigtagprobecut and the probe passes the second trigger and fails the first trigger
466  MonElemFuncs::initStdEleCutHists(monElemCont->cutMonElems(),trigNameLeg2,trigNameLeg2+"_"+monElemCont->name()+"_gsfEle_passLeg2failLeg1",bins,&(*(new EgMultiCut<OffEle>) << new EgGreaterCut<OffEle,float>(etCutValue,&OffEle::etSC) << new EgObjTrigCut<OffEle>(TrigCodes::getCode(trigNameLeg2),EgObjTrigCut<OffEle>::AND,TrigCodes::getCode(trigNameLeg1),EgObjTrigCut<OffEle>::AND)));
467 
468 }
469 
470 
471 //Now same for photons
472 void MonElemFuncs::initTrigTagProbeHists(std::vector<MonElemContainer<OffPho>*>& phoMonElems,const std::vector<std::string> filterNames,int cutMask,const BinData& bins)
473 {
474  for(size_t filterNr=0;filterNr<filterNames.size();filterNr++){
475 
476  std::string trigName(filterNames[filterNr]);
477  //float etCutValue = trigTools::getSecondEtThresFromName(trigName);
478  float etCutValue = 0.;
479  //std::cout<<"TrigName= "<<trigName<<" etCutValue= "<<etCutValue<<std::endl;
480  MonElemContainer<OffPho>* monElemCont = new MonElemContainer<OffPho>("trigTagProbe","Trigger Tag and Probe",new EgTrigTagProbeCut_NewPho(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),cutMask,&OffPho::cutCode));
481  //this is all that pass trigtagprobecut
482  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_all",bins,new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC));
483  //this is all that pass trigtagprobecut and the probe passes the test trigger
484  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_pass",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND)));
485  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is NOT a tag
486  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_passNotTag",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND,TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffPho>::AND)));
487  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is also a tag
488  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_passTagTag",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND) << new EgObjTrigCut<OffPho>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffPho>::AND) ));
489  //this is all that pass trigtagprobecut and the probe fails the test trigger
490  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_fail",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),EgObjTrigCut<OffPho>::AND,TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND)));
491  /*
492  monElemCont->monElems().push_back(new MonElemMgrEBEE<OffPho,float>(trigName+"_"+monElemCont->name()+"_pho_all_etUnCut",monElemCont->name()+"_gsfEle_all E_{T} (Uncut);E_{T} (GeV)",
493  bins.et.nr,bins.et.min,bins.et.max,&OffPho::et));
494  monElemCont->cutMonElems().push_back(new MonElemWithCutEBEE<OffPho,float>(trigName+"_"+monElemCont->name()+"_pho_pass_etUnCut",monElemCont->name()+"_gsfEle_pass E_{T} (Uncut);E_{T} (GeV)",
495  bins.et.nr,bins.et.min,bins.et.max,&OffPho::et,new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND)));*/
496  phoMonElems.push_back(monElemCont);
497  } //end filter names loop
498 
499 }
500 
501 void MonElemFuncs::initTrigTagProbeHist(std::vector<MonElemContainer<OffPho>*>& phoMonElems,const std::string filterName,int cutMask,const BinData& bins)
502 {
503  std::string trigName(filterName);
504  //float etCutValue = 1.1*trigTools::getSecondEtThresFromName(trigName);
505  float etCutValue = 0.;
506  //std::cout<<"TrigName= "<<trigName<<" etCutValue= "<<etCutValue<<std::endl;
507  MonElemContainer<OffPho>* monElemCont = new MonElemContainer<OffPho>("trigTagProbe","Trigger Tag and Probe",new EgTrigTagProbeCut_NewPho(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),cutMask,&OffPho::cutCode));
508  //this is all that pass trigtagprobecut
509  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_all",bins,new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC));
510  //this is all that pass trigtagprobecut and the probe passes the test trigger
511  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_pass",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND)));
512  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is NOT a tag
513  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_passNotTag",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND,TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffPho>::AND)));
514  //this is all that pass trigtagprobecut and the probe passes the test trigger and the probe is also a tag
515  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_passTagTag",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND) << new EgObjTrigCut<OffPho>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17TrackIsolFilter"),EgObjTrigCut<OffPho>::AND) ));
516  //this is all that pass trigtagprobecut and the probe fails the test trigger
517  MonElemFuncs::initStdPhoCutHists(monElemCont->cutMonElems(),trigName,trigName+"_"+monElemCont->name()+"_pho_fail",bins,&(*(new EgMultiCut<OffPho>) << new EgGreaterCut<OffPho,float>(etCutValue,&OffPho::etSC) << new EgObjTrigCut<OffPho>(TrigCodes::getCode("hltEle32CaloIdTCaloIsoTTrkIdTTrkIsoTSC17HEDoubleFilter"),EgObjTrigCut<OffPho>::AND,TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND)));
518  /*
519  monElemCont->monElems().push_back(new MonElemMgrEBEE<OffPho,float>(trigName+"_"+monElemCont->name()+"_pho_all_etUnCut",monElemCont->name()+"_gsfEle_all E_{T} (Uncut);E_{T} (GeV)",
520  bins.et.nr,bins.et.min,bins.et.max,&OffPho::et));
521  monElemCont->cutMonElems().push_back(new MonElemWithCutEBEE<OffPho,float>(trigName+"_"+monElemCont->name()+"_pho_pass_etUnCut",monElemCont->name()+"_gsfEle_pass E_{T} (Uncut);E_{T} (GeV)",
522  bins.et.nr,bins.et.min,bins.et.max,&OffPho::et,new EgObjTrigCut<OffPho>(TrigCodes::getCode(trigName),EgObjTrigCut<OffPho>::AND)));*/
523  phoMonElems.push_back(monElemCont);
524 }
525 
526 
527 
528 
egHLT::OffPho::e1x5Over5x5
float e1x5Over5x5() const
Definition: EgHLTOffPho.h:116
egHLT::OffEle::e2x5MaxOver5x5
float e2x5MaxOver5x5() const
Definition: EgHLTOffEle.h:153
egHLT::OffEle::detEta
float detEta() const
Definition: EgHLTOffEle.h:121
egHLT::TrigCodes::getCode
static TrigBitSet getCode(const std::string &descript)
Definition: EgHLTTrigCodes.h:66
egHLT::BinData::isolNrTrks
Data1D isolNrTrks
Definition: EgHLTBinData.h:55
egHLT::OffEle::epIn
float epIn() const
Definition: EgHLTOffEle.h:144
egHLT::OffPho::et
float et() const
Definition: EgHLTOffPho.h:91
egHLT::OffPho::isolHad
float isolHad() const
Definition: EgHLTOffPho.h:121
egHLT::MonElemContainer::name
const std::string & name() const
Definition: EgHLTMonElemContainer.h:60
egHLT::OffEle::r9
float r9() const
Definition: EgHLTOffEle.h:156
egHLT::OffPho::isolPtTrks
float isolPtTrks() const
Definition: EgHLTOffPho.h:123
egHLT::MonElemFuncs::initTightLooseTrigHists
void initTightLooseTrigHists(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::vector< std::string > &tightLooseTrigs, const BinData &bins, EgHLTDQMCut< OffEle > *eleCut)
Definition: EgHLTMonElemFuncs.cc:244
egHLT::OffEle::invEInvP
float invEInvP() const
Definition: EgHLTOffEle.h:159
reco_skim_cfg_mod.filterName
tuple filterName
Definition: reco_skim_cfg_mod.py:159
NULL
#define NULL
Definition: scimark2.h:8
egHLT::MonElemFuncs::initStdEleCutHists
void initStdEleCutHists(std::vector< MonElemWithCutBase< OffEle > * > &histVec, const std::string &filterName, const std::string &baseName, const BinData &bins, EgHLTDQMCut< OffEle > *cut=NULL)
Definition: EgHLTMonElemFuncs.cc:200
egHLT::OffEle::et
float et() const
Definition: EgHLTOffEle.h:113
egHLT::MonElemFuncs::initTightLooseDiObjTrigHistsTrigCuts
void initTightLooseDiObjTrigHistsTrigCuts(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::vector< std::string > &tightLooseTrigs, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:364
egHLT::OffEle::dPhiIn
float dPhiIn() const
Definition: EgHLTOffEle.h:141
egHLT::OffPho::trigCutsCutCode
int trigCutsCutCode(const TrigCodes::TrigBitSet &trigger) const
Definition: EgHLTOffPho.cc:14
egHLT::MonElemFuncs::initStdEleHists
void initStdEleHists(std::vector< MonElemManagerBase< OffEle > * > &histVec, const std::string &filterName, const std::string &baseName, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:11
splitString
std::vector< std::string > splitString(const std::string &fLine)
Definition: HcalDbASCIIIO.cc:55
egHLT::MonElemFuncs::initStdEffHists
void initStdEffHists(std::vector< MonElemWithCutBase< OffEle > * > &histVec, const std::string &filterName, const std::string &baseName, int nrBins, double xMin, double xMax, float(OffEle::*vsVarFunc)() const, const CutMasks &masks)
Definition: EgHLTMonElemFuncs.cc:90
egHLT::BinData::isolPtTrks
Data1D isolPtTrks
Definition: EgHLTBinData.h:54
egHLT::EgTrigTagProbeCut_NewPho
Definition: EgHLTDQMCut.h:485
egHLT::OffPho::detEta
float detEta() const
Definition: EgHLTOffPho.h:98
egHLT::MonElemFuncs::initStdPhoCutHists
void initStdPhoCutHists(std::vector< MonElemWithCutBase< OffPho > * > &histVec, const std::string &filterName, const std::string &baseName, const BinData &bins, EgHLTDQMCut< OffPho > *cut=NULL)
Definition: EgHLTMonElemFuncs.cc:224
egHLT::OffPho::hOverE
float hOverE() const
Definition: EgHLTOffPho.h:106
egHLT::OffEle::dEtaIn
float dEtaIn() const
Definition: EgHLTOffEle.h:140
egHLT::OffEle::phi
float phi() const
Definition: EgHLTOffEle.h:117
egHLT::OffPho::r9
float r9() const
Definition: EgHLTOffPho.h:117
egHLT::EgTrigTagProbeCut_New
Definition: EgHLTDQMCut.h:468
egHLT::MonElemFuncs::initTrigTagProbeHists
void initTrigTagProbeHists(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::vector< std::string > filterNames, int cutMask, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:397
egHLT::OffEle::sigmaIEtaIEta
float sigmaIEtaIEta() const
Definition: EgHLTOffEle.h:150
egHLT::EgHLTDQMCut::clone
virtual EgHLTDQMCut< T > * clone() const =0
egHLT::OffPho::e2x5MaxOver5x5
float e2x5MaxOver5x5() const
Definition: EgHLTOffPho.h:115
egHLT::OffEle::hltIsolHad
float hltIsolHad() const
Definition: EgHLTOffEle.h:171
align_tpl.cut
tuple cut
Definition: align_tpl.py:88
egHLT::MonElemFuncs::initTrigTagProbeHist_2Leg
void initTrigTagProbeHist_2Leg(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::string filterName, int cutMask, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:453
egHLT::OffEle::cutCode
int cutCode() const
Definition: EgHLTOffEle.h:200
egHLT::OffEle::epOut
float epOut() const
Definition: EgHLTOffEle.h:145
egHLT::OffPho::etSC
float etSC() const
Definition: EgHLTOffPho.h:96
egHLT::OffPho::sigmaIEtaIEta
float sigmaIEtaIEta() const
Definition: EgHLTOffPho.h:112
egHLT::OffEle::etSC
float etSC() const
Definition: EgHLTOffEle.h:118
egHLT::BinData::sigEtaEta
Data1D sigEtaEta
Definition: EgHLTBinData.h:43
egHLT::OffEle::hltIsolTrksEle
float hltIsolTrksEle() const
Definition: EgHLTOffEle.h:169
egHLT::OffPho::phi
float phi() const
Definition: EgHLTOffPho.h:95
egHLT::OffEle::trigCutsCutCode
int trigCutsCutCode(const TrigCodes::TrigBitSet &trigger) const
Definition: EgHLTOffEle.cc:18
egHLT::MonElemContainer::monElems
std::vector< MonElemManagerBase< T > * > & monElems()
Definition: EgHLTMonElemContainer.h:56
egHLT::OffPho::isolNrTrks
int isolNrTrks() const
Definition: EgHLTOffPho.h:122
egHLT::MonElemFuncs::addTightLooseTrigHist
void addTightLooseTrigHist(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::string &tightTrig, const std::string &looseTrig, EgHLTDQMCut< OffEle > *eleCut, const std::string &histId, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:282
egHLT::OffPho::cutCode
int cutCode() const
Definition: EgHLTOffPho.h:129
egHLT::OffEle::e1x5Over5x5
float e1x5Over5x5() const
Definition: EgHLTOffEle.h:154
egHLT::OffEle::NVertex
int NVertex() const
Definition: EgHLTOffEle.h:103
egHLT::MonElemFuncs::initTrigTagProbeHist
void initTrigTagProbeHist(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::string filterName, int cutMask, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:427
egHLT::MonElemFuncs::initStdPhoHists
void initStdPhoHists(std::vector< MonElemManagerBase< OffPho > * > &histVec, const std::string &filterName, const std::string &baseName, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:49
egHLT::OffEle::hOverE
float hOverE() const
Definition: EgHLTOffEle.h:139
split
double split
Definition: MVATrainer.cc:139
egHLT::OffPho::isolEm
float isolEm() const
Definition: EgHLTOffPho.h:120
egHLT::MonElemFuncs::initTightLooseTrigHistsTrigCuts
void initTightLooseTrigHistsTrigCuts(std::vector< MonElemContainer< OffEle > * > &eleMonElems, const std::vector< std::string > &tightLooseTrigs, const BinData &bins)
Definition: EgHLTMonElemFuncs.cc:335
egHLT::MonElemContainer::cutMonElems
std::vector< MonElemWithCutBase< T > * > & cutMonElems()
Definition: EgHLTMonElemContainer.h:54
egHLT::OffEle::hltIsolEm
float hltIsolEm() const
Definition: EgHLTOffEle.h:172
masks
static const uint32_t masks[]
Definition: CaloRecHit.cc:14
Generated for CMSSW Reference Manual by doxygen 1.8.5