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.23 by wulsin, Wed Apr 17 15:39:43 2013 UTC vs.
Revision 1.31 by ahart, Tue May 7 02:41:34 2013 UTC

# Line 9 | Line 9
9   #include <fstream>
10   #include <iostream>
11   #include <algorithm>
12 + #include <regex>
13   #include "TH1.h"
14   #include "TH2.h"
15   #include "TH1D.h"
# Line 43 | Line 44
44   #include "ProductArea/BNcollections/interface/BNtau.h"
45   #include "ProductArea/BNcollections/interface/BNgenjet.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
51 < #include "DisplacedSUSY/Configuration/interface/CTauWeight.h"
52 < #endif
52 > #include "OSUT3Analysis/AnaTools/interface/StopCTauWeight.h"
53  
54   using namespace std;
55  
# Line 61 | Line 61 | class OSUAnalysis : public edm::EDAnalyz
61  
62        void analyze (const edm::Event &, const edm::EventSetup &);
63        bool evaluateComparison (double, string, double);
64 +      bool evaluateComparison (string, string, string);
65        bool evaluateTriggers   (vector<string>,const BNtriggerCollection*);
66        double applyFunction(string, double);
67  
68 <      double valueLookup (const BNjet* object, string variable, string function = "");
69 <      double valueLookup (const BNmuon* object, string variable, string function = "");
70 <      double valueLookup (const BNmuon* object1, const BNmuon* object2, string variable, string function = "");
71 <      double valueLookup (const BNelectron* object, string variable, string function = "");
72 <      double valueLookup (const BNelectron* object1, const BNelectron* object2, string variable, string function = "");
73 <      double valueLookup (const BNelectron* object1, const BNmuon* object2, string variable, string function = "");
74 <      double valueLookup (const BNelectron* object1, const BNtrack* object2, string variable, string function = "");
75 <      double valueLookup (const BNmuon* object1, const BNtrack* object2, string variable, string function = "");
76 <      double valueLookup (const BNmuon* object1, const BNtau* object2, string variable, string function = "");
77 <      double valueLookup (const BNtau* object1, const BNtau* object2, string variable, string function = "");
78 <      double valueLookup (const BNtau* object1, const BNtrack* object2, string variable, string function);
79 <      double valueLookup (const BNevent* object, string variable, string function = "");
80 <      double valueLookup (const BNtau* object, string variable, string function = "");
81 <      double valueLookup (const BNmet* object, string variable, string function = "");
82 <      double valueLookup (const BNtrack* object, string variable, string function = "");
83 <      double valueLookup (const BNtrack* object1, const BNevent* object2, string variable, string function = "");
84 <      double valueLookup (const BNgenjet* object, string variable, string function = "");
85 <      double valueLookup (const BNmcparticle* object, string variable, string function = "");
86 <      double valueLookup (const BNprimaryvertex* object, string variable, string function = "");
87 <      double valueLookup (const BNbxlumi* object, string variable, string function = "");
88 <      double valueLookup (const BNphoton* object, string variable, string function = "");
89 <      double valueLookup (const BNsupercluster* object, string variable, string function = "");
68 >      double valueLookup (const BNjet* object, string variable, string function, string &stringValue);
69 >      double valueLookup (const BNmuon* object, string variable, string function, string &stringValue);
70 >      double valueLookup (const BNmuon* object1, const BNmuon* object2, string variable, string function, string &stringValue);
71 >      double valueLookup (const BNelectron* object, string variable, string function, string &stringValue);
72 >      double valueLookup (const BNelectron* object1, const BNelectron* object2, string variable, string function, string &stringValue);
73 >      double valueLookup (const BNelectron* object1, const BNmuon* object2, string variable, string function, string &stringValue);
74 >      double valueLookup (const BNelectron* object1, const BNjet* object2, string variable, string function, string &stringValue);
75 >      double valueLookup (const BNelectron* object1, const BNtrack* object2, string variable, string function, string &stringValue);
76 >      double valueLookup (const BNmuon* object1, const BNjet* object2, string variable, string function, string &stringValue);
77 >      double valueLookup (const BNmuon* object1, const BNtrack* object2, string variable, string function, string &stringValue);
78 >      double valueLookup (const BNmuon* object1, const BNtau* object2, string variable, string function, string &stringValue);
79 >      double valueLookup (const BNtau* object1, const BNtau* object2, string variable, string function, string &stringValue);
80 >      double valueLookup (const BNtau* object1, const BNtrack* object2, string variable, string function, string &stringValue);
81 >      double valueLookup (const BNjet* object1, const BNjet* object2, string variable, string function, string &stringValue);
82 >      double valueLookup (const BNevent* object, string variable, string function, string &stringValue);
83 >      double valueLookup (const BNtau* object, string variable, string function, string &stringValue);
84 >      double valueLookup (const BNmet* object, string variable, string function, string &stringValue);
85 >      double valueLookup (const BNtrack* object, string variable, string function, string &stringValue);
86 >      double valueLookup (const BNtrack* object1, const BNevent* object2, string variable, string function, string &stringValue);
87 >      double valueLookup (const BNgenjet* object, string variable, string function, string &stringValue);
88 >      double valueLookup (const BNmcparticle* object, string variable, string function, string &stringValue);
89 >      double valueLookup (const BNprimaryvertex* object, string variable, string function, string &stringValue);
90 >      double valueLookup (const BNbxlumi* object, string variable, string function, string &stringValue);
91 >      double valueLookup (const BNphoton* object, string variable, string function, string &stringValue);
92 >      double valueLookup (const BNsupercluster* object, string variable, string function, string &stringValue);
93 >      double valueLookup (const BNtrigobj* object, string variable, string function, string &stringValue);
94 >      double valueLookup (const BNelectron* object1, const BNtrigobj* object2, string variable, string function, string &stringValue);
95 >      double valueLookup (const BNmuon* object1, const BNtrigobj* object2, string variable, string function, string &stringValue);
96 >      double valueLookup (const BNstop* object, string variable, string function, string &stringValue);
97  
98        int getTrkIsIso (const BNtrack* track1, const BNtrackCollection* trackColl);
99        double getTrkPtTrue (const BNtrack* track1, const BNmcparticleCollection* genPartColl);
# Line 99 | Line 107 | class OSUAnalysis : public edm::EDAnalyz
107        //BNtrigobj
108  
109        vector<string> splitString (string);
110 <      void getTwoObjs(string tempInputCollection, string& obj1, string& obj2);
111 <      string getObjToGet(string obj);  
110 >      void getTwoObjs(string tempInputCollection, string& obj1, string& obj2);
111 >      string getObjToGet(string obj);
112  
113      private:
114  
# Line 133 | Line 141 | class OSUAnalysis : public edm::EDAnalyz
141        edm::InputTag photons_;
142        edm::InputTag superclusters_;
143        edm::InputTag triggers_;
144 <      std::string puFile_;
145 <      std::string deadEcalFile_;
146 <      std::string muonSFFile_;
147 <      std::string dataPU_;
148 <      std::string electronSFID_;
149 <      std::string muonSF_;
150 <      std::string dataset_;
151 <      std::string datasetType_;
144 >      edm::InputTag trigobjs_;
145 >      string puFile_;
146 >      string deadEcalFile_;
147 >      string muonSFFile_;
148 >      string dataPU_;
149 >      string electronSFID_;
150 >      string muonSF_;
151 >      string dataset_;
152 >      string datasetType_;
153        vector<edm::ParameterSet> channels_;
154        vector<edm::ParameterSet> histogramSets_;
155        bool plotAllObjectsInPassingEvents_;
156        bool doPileupReweighting_;
157 <      bool printEventInfo_;  
158 <      bool useTrackCaloRhoCorr_;  // to use the calo-based rho correction for the by-hand calculation of the track isolation energy
157 >      bool applyLeptonSF_;
158 >      bool printEventInfo_;
159 >      bool useTrackCaloRhoCorr_;  // to use the calo-based rho correction for the by-hand calculation of the track isolation energy
160        vector<double> stopCTau_;
161  
162        struct DeadEcal {
# Line 154 | Line 164 | class OSUAnalysis : public edm::EDAnalyz
164          double phiEcal;
165        };
166  
167 <      std::vector<DeadEcal> DeadEcalVec;
167 >      vector<DeadEcal> DeadEcalVec;
168  
169  
170        //Collections
171        edm::Handle<BNtriggerCollection> triggers;
172 +      edm::Handle<BNtrigobjCollection> trigobjs;
173        edm::Handle<BNjetCollection> jets;
174        edm::Handle<BNmuonCollection> muons;
175        edm::Handle<BNelectronCollection> electrons;
# Line 172 | Line 183 | class OSUAnalysis : public edm::EDAnalyz
183        edm::Handle<BNbxlumiCollection> bxlumis;
184        edm::Handle<BNphotonCollection> photons;
185        edm::Handle<BNsuperclusterCollection> superclusters;
186 <      edm::Handle<double> rhokt6CaloJetsHandle_;  
186 >      edm::Handle<BNstopCollection> stops;
187 >
188 >      edm::Handle<double> rhokt6CaloJetsHandle_;
189  
190        flagMap cumulativeFlags;
191  
# Line 188 | Line 201 | class OSUAnalysis : public edm::EDAnalyz
201          string name;
202          string title;
203          vector<double> bins;
204 +        vector<double> variableBinsX;
205 +        vector<double> variableBinsY;
206          string inputCollection;
207          vector<string> inputVariables;
208        } histogram;
# Line 198 | Line 213 | class OSUAnalysis : public edm::EDAnalyz
213          vector<string> variables;
214          vector<string> comparativeOperators;// >, <, =, etc.
215          vector<double> cutValues;
216 +        vector<string> cutStringValues;
217          vector<string> logicalOperators;//and, or
218          int numSubcuts;
219          int    numberRequired;
# Line 220 | Line 236 | class OSUAnalysis : public edm::EDAnalyz
236        PUWeight *puWeight_;
237        MuonSFWeight *muonSFWeight_;
238        ElectronSFWeight *electronSFWeight_;
239 < #ifdef DISPLACED_SUSY
240 <      CTauWeight *cTauWeight_;
241 < #endif
226 <      double cTauScaleFactor_;
239 >
240 >      StopCTauWeight *stopCTauWeight_;
241 >      double stopCTauScaleFactor_;
242  
243        template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string);
244        template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, vector<bool>, vector<bool>, string);
# Line 232 | Line 247 | class OSUAnalysis : public edm::EDAnalyz
247        template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double);
248        template <class InputCollection> void fill2DHistogram(TH2*, histogram, InputCollection, vector<bool>, double);
249        template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double);
250 <      bool getPreviousCumulativeFlags(uint currentCutIndex, flagMap &individualFlags, string obj1Type, uint object1);  
250 >      bool getPreviousCumulativeFlags(uint currentCutIndex, flagMap &individualFlags, string obj1Type, uint object1);
251  
252        template <class InputObject> int getGenMatchedParticleIndex(InputObject);
253        int getPdgIdBinValue(int);

Diff Legend

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