ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.h
(Generate patch)

Comparing UserCode/OSUT3Analysis/AnaTools/plugins/OSUAnalysis.h (file contents):
Revision 1.1 by lantonel, Mon Jan 21 10:39:36 2013 UTC vs.
Revision 1.27 by biliu, Wed May 1 16:18:54 2013 UTC

# Line 6 | Line 6
6   #include <string>
7   #include <vector>
8  
9 + #include <fstream>
10 + #include <iostream>
11 + #include <algorithm>
12 + #include "TH1.h"
13 + #include "TH2.h"
14   #include "TH1D.h"
15   #include "TH2D.h"
16   #include "TLorentzVector.h"
# Line 18 | Line 23
23   #include "DataFormats/Common/interface/Handle.h"
24   #include "CommonTools/UtilAlgos/interface/TFileService.h"
25   #include "FWCore/Framework/interface/MakerMacros.h"
26 + #include "DataFormats/Math/interface/deltaPhi.h"
27 + #include "DataFormats/Math/interface/deltaR.h"
28  
29   #include "ProductArea/BNcollections/interface/BNbxlumi.h"
30   #include "ProductArea/BNcollections/interface/BNelectron.h"
# Line 35 | Line 42
42   #include "ProductArea/BNcollections/interface/BNtrigobj.h"
43   #include "ProductArea/BNcollections/interface/BNtau.h"
44   #include "ProductArea/BNcollections/interface/BNgenjet.h"
45 <
39 < #include "OSUT3Analysis/AnaTools/interface/CutFlow.h"
45 > #include "ProductArea/BNcollections/interface/BNstop.h"
46  
47  
48 + #include "OSUT3Analysis/AnaTools/interface/CutFlow.h"
49 + #include "OSUT3Analysis/AnaTools/interface/PUWeight.h"
50 + #include "OSUT3Analysis/AnaTools/interface/SFWeight.h"
51  
52 + #ifdef DISPLACED_SUSY
53 + #include "DisplacedSUSY/Configuration/interface/CTauWeight.h"
54 + #endif
55  
56   using namespace std;
57  
# Line 51 | Line 63 | class OSUAnalysis : public edm::EDAnalyz
63  
64        void analyze (const edm::Event &, const edm::EventSetup &);
65        bool evaluateComparison (double, string, double);
66 +      bool evaluateTriggers   (vector<string>,const BNtriggerCollection*);
67 +      double applyFunction(string, double);
68 +
69 +      double valueLookup (const BNjet* object, string variable, string function = "");
70 +      double valueLookup (const BNmuon* object, string variable, string function = "");
71 +      double valueLookup (const BNmuon* object1, const BNmuon* object2, string variable, string function = "");
72 +      double valueLookup (const BNelectron* object, string variable, string function = "");
73 +      double valueLookup (const BNelectron* object1, const BNelectron* object2, string variable, string function = "");
74 +      double valueLookup (const BNelectron* object1, const BNmuon* object2, string variable, string function = "");
75 +      double valueLookup (const BNelectron* object1, const BNjet* object2, string variable, string function = "");
76 +      double valueLookup (const BNelectron* object1, const BNtrack* object2, string variable, string function = "");
77 +      double valueLookup (const BNmuon* object1, const BNjet* object2, string variable, string function = "");
78 +      double valueLookup (const BNmuon* object1, const BNtrack* object2, string variable, string function = "");
79 +      double valueLookup (const BNmuon* object1, const BNtau* object2, string variable, string function = "");
80 +      double valueLookup (const BNtau* object1, const BNtau* object2, string variable, string function = "");
81 +      double valueLookup (const BNtau* object1, const BNtrack* object2, string variable, string function);
82 +      double valueLookup (const BNjet* object1, const BNjet* object2, string variable, string function = "");
83 +      double valueLookup (const BNevent* object, string variable, string function = "");
84 +      double valueLookup (const BNtau* object, string variable, string function = "");
85 +      double valueLookup (const BNmet* object, string variable, string function = "");
86 +      double valueLookup (const BNtrack* object, string variable, string function = "");
87 +      double valueLookup (const BNtrack* object1, const BNevent* object2, string variable, string function = "");
88 +      double valueLookup (const BNgenjet* object, string variable, string function = "");
89 +      double valueLookup (const BNmcparticle* object, string variable, string function = "");
90 +      double valueLookup (const BNprimaryvertex* object, string variable, string function = "");
91 +      double valueLookup (const BNbxlumi* object, string variable, string function = "");
92 +      double valueLookup (const BNphoton* object, string variable, string function = "");
93 +      double valueLookup (const BNsupercluster* object, string variable, string function = "");
94 +      double valueLookup (const BNstop* object, string variable, string function = "");
95 +
96 +      int getTrkIsIso (const BNtrack* track1, const BNtrackCollection* trackColl);
97 +      double getTrkPtTrue (const BNtrack* track1, const BNmcparticleCollection* genPartColl);
98 +      double getTrkPtRes (const BNtrack* track1);
99 +      double getTrkIsIso (const BNtrack* track1);
100 +      double getTrkCaloTotRhoCorr(const BNtrack* track);
101 +      void WriteDeadEcal ();
102 +      int getTrkIsMatchedDeadEcal (const BNtrack* track1);
103  
55 –      double valueLookup (const BNjet*, string);
56 –      double valueLookup (const BNmuon*, string);
57 –      double valueLookup (const BNelectron*, string);
58 –      double valueLookup (const BNevent*, string);
59 –      double valueLookup (const BNtau*, string);
60 –      double valueLookup (const BNmet*, string);
61 –      double valueLookup (const BNtrack*, string);
62 –      double valueLookup (const BNgenjet*, string);
63 –      double valueLookup (const BNmcparticle*, string);
64 –      double valueLookup (const BNprimaryvertex*, string);
65 –      double valueLookup (const BNbxlumi*, string);
66 –      double valueLookup (const BNphoton*, string);
67 –      double valueLookup (const BNsupercluster*, string);
104        //BNskimbits
69 –      //BNtrigger
105        //BNtrigobj
106  
107        vector<string> splitString (string);
108 <
108 >      void getTwoObjs(string tempInputCollection, string& obj1, string& obj2);
109 >      string getObjToGet(string obj);  
110  
111      private:
112  
# Line 80 | Line 116 | class OSUAnalysis : public edm::EDAnalyz
116        //inner vector corresponds to each object in input collection
117        //bool         tells if object passes or fails cuts
118        typedef map<string, vector < vector<bool> > > flagMap;
119 <      
119 >
120        //counterMap:
121        //string holds input collection type
122        //vector corresponds to each cut
# Line 97 | Line 133 | class OSUAnalysis : public edm::EDAnalyz
133        edm::InputTag tracks_;
134        edm::InputTag genjets_;
135        edm::InputTag mcparticles_;
136 +      edm::InputTag stops_;
137        edm::InputTag primaryvertexs_;
138        edm::InputTag bxlumis_;
139        edm::InputTag photons_;
140        edm::InputTag superclusters_;
141 +      edm::InputTag triggers_;
142 +      std::string puFile_;
143 +      std::string deadEcalFile_;
144 +      std::string muonSFFile_;
145 +      std::string dataPU_;
146 +      std::string electronSFID_;
147 +      std::string muonSF_;
148 +      std::string dataset_;
149 +      std::string datasetType_;
150 +      vector<edm::ParameterSet> channels_;
151 +      vector<edm::ParameterSet> histogramSets_;
152 +      bool plotAllObjectsInPassingEvents_;
153 +      bool doPileupReweighting_;
154 +      bool printEventInfo_;  
155 +      bool useTrackCaloRhoCorr_;  // to use the calo-based rho correction for the by-hand calculation of the track isolation energy
156 +      vector<double> stopCTau_;
157 +
158 +      struct DeadEcal {
159 +        double etaEcal;
160 +        double phiEcal;
161 +      };
162  
163 +      std::vector<DeadEcal> DeadEcalVec;
164  
165 <      vector<edm::ParameterSet> channels_;
165 >
166 >      //Collections
167 >      edm::Handle<BNtriggerCollection> triggers;
168 >      edm::Handle<BNjetCollection> jets;
169 >      edm::Handle<BNmuonCollection> muons;
170 >      edm::Handle<BNelectronCollection> electrons;
171 >      edm::Handle<BNeventCollection> events;
172 >      edm::Handle<BNtauCollection> taus;
173 >      edm::Handle<BNmetCollection> mets;
174 >      edm::Handle<BNtrackCollection> tracks;
175 >      edm::Handle<BNgenjetCollection> genjets;
176 >      edm::Handle<BNmcparticleCollection> mcparticles;
177 >      edm::Handle<BNprimaryvertexCollection> primaryvertexs;
178 >      edm::Handle<BNbxlumiCollection> bxlumis;
179 >      edm::Handle<BNphotonCollection> photons;
180 >      edm::Handle<BNsuperclusterCollection> superclusters;
181 >      edm::Handle<BNstopCollection> stops;
182 >
183 >      edm::Handle<double> rhokt6CaloJetsHandle_;  
184 >
185 >      flagMap cumulativeFlags;
186  
187        vector<map<string, TH1D*> > oneDHists_;
188 +      vector<map<string, TH2D*> > twoDHists_;
189  
190        edm::Service<TFileService> fs_;
191  
192        CutFlow *masterCutFlow_;
193        vector<CutFlow *> cutFlows_;
194  
195 +      typedef struct {
196 +        string name;
197 +        string title;
198 +        vector<double> bins;
199 +        string inputCollection;
200 +        vector<string> inputVariables;
201 +      } histogram;
202  
203        struct cut {
204 <        string inputCollection;
205 <        string variable;
206 <        string comparativeOperator;
207 <        double cutValue;
208 <        int    numberRequired;
209 <        string eventComparativeOperator;
210 <        string name;
204 >        string inputCollection;
205 >        vector<string> functions;
206 >        vector<string> variables;
207 >        vector<string> comparativeOperators;// >, <, =, etc.
208 >        vector<double> cutValues;
209 >        vector<string> logicalOperators;//and, or
210 >        int numSubcuts;
211 >        int    numberRequired;
212 >        string eventComparativeOperator;
213 >        string name;
214        };
215  
216        struct channel {
217 <        string name;
218 <        vector<cut> cuts;
219 <        vector<string> inputCollections;
217 >        string name;
218 >        vector<string> triggers;
219 >        vector<cut> cuts;
220        };
221  
222 <      vector<string>  allNecessaryObjects;
222 >      vector<string>  objectsToGet;
223 >      vector<string>  objectsToCut;
224 >      vector<string>  objectsToPlot;
225        vector<channel> channels;
226 +      vector<histogram> histograms;
227 +
228 +      PUWeight *puWeight_;
229 +      MuonSFWeight *muonSFWeight_;
230 +      ElectronSFWeight *electronSFWeight_;
231 + #ifdef DISPLACED_SUSY
232 +      CTauWeight *cTauWeight_;
233 + #endif
234 +      double cTauScaleFactor_;
235  
236        template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string);
237 +      template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, vector<bool>, vector<bool>, string);
238 +
239 +      template <class InputCollection> void fill1DHistogram(TH1*, histogram, InputCollection, vector<bool>, double);
240 +      template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double);
241 +      template <class InputCollection> void fill2DHistogram(TH2*, histogram, InputCollection, vector<bool>, double);
242 +      template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double);
243 +      bool getPreviousCumulativeFlags(uint currentCutIndex, flagMap &individualFlags, string obj1Type, uint object1);  
244 +
245 +      template <class InputObject> int getGenMatchedParticleIndex(InputObject);
246 +      int getPdgIdBinValue(int);
247 +      int findTauMotherIndex(const BNmcparticle*);
248 +
249 +      template <class InputObject> double getGenDeltaRLowest(InputObject);
250  
251  
252 +      const BNprimaryvertex *chosenVertex ();
253 +      const BNmet *chosenMET ();
254 +      const BNelectron *chosenElectron ();
255 +      const BNmuon *chosenMuon ();
256  
257    };
258  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines