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.6 by lantonel, Thu Feb 7 14:52:37 2013 UTC vs.
Revision 1.22 by jbrinson, Thu Apr 11 15:52:25 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 19 | 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 38 | Line 44
44   #include "ProductArea/BNcollections/interface/BNgenjet.h"
45  
46   #include "OSUT3Analysis/AnaTools/interface/CutFlow.h"
47 + #include "OSUT3Analysis/AnaTools/interface/PUWeight.h"
48 + #include "OSUT3Analysis/AnaTools/interface/SFWeight.h"
49  
50 <
51 <
50 > #ifdef DISPLACED_SUSY
51 > #include "DisplacedSUSY/Configuration/interface/CTauWeight.h"
52 > #endif
53  
54   using namespace std;
55  
# Line 57 | Line 66 | class OSUAnalysis : public edm::EDAnalyz
66  
67        double valueLookup (const BNjet* object, string variable, string function = "");
68        double valueLookup (const BNmuon* object, string variable, string function = "");
69 +      double valueLookup (const BNmuon* object1, const BNmuon* object2, string variable, string function = "");
70        double valueLookup (const BNelectron* object, string variable, string function = "");
71 +      double valueLookup (const BNelectron* object1, const BNelectron* object2, string variable, string function = "");
72 +      double valueLookup (const BNelectron* object1, const BNmuon* object2, string variable, string function = "");
73 +      double valueLookup (const BNelectron* object1, const BNtrack* object2, string variable, string function = "");
74 +      double valueLookup (const BNmuon* object1, const BNtrack* object2, string variable, string function = "");
75 +      double valueLookup (const BNmuon* object1, const BNtau* object2, string variable, string function = "");
76 +      double valueLookup (const BNtau* object1, const BNtau* object2, string variable, string function = "");
77 +      double valueLookup (const BNtau* object1, const BNtrack* object2, string variable, string function);
78        double valueLookup (const BNevent* object, string variable, string function = "");
79        double valueLookup (const BNtau* object, string variable, string function = "");
80        double valueLookup (const BNmet* object, string variable, string function = "");
81        double valueLookup (const BNtrack* object, string variable, string function = "");
82 +      double valueLookup (const BNtrack* object1, const BNevent* object2, string variable, string function = "");
83        double valueLookup (const BNgenjet* object, string variable, string function = "");
84        double valueLookup (const BNmcparticle* object, string variable, string function = "");
85        double valueLookup (const BNprimaryvertex* object, string variable, string function = "");
86        double valueLookup (const BNbxlumi* object, string variable, string function = "");
87        double valueLookup (const BNphoton* object, string variable, string function = "");
88        double valueLookup (const BNsupercluster* object, string variable, string function = "");
89 +
90 +      int getTrkIsIso (const BNtrack* track1, const BNtrackCollection* trackColl);
91 +      double getTrkPtTrue (const BNtrack* track1, const BNmcparticleCollection* genPartColl);
92 +      double getTrkPtRes (const BNtrack* track1);
93 +      double getTrkIsIso (const BNtrack* track1);
94 +      double getTrkCaloTotRhoCorr(const BNtrack* track);
95 +      void WriteDeadEcal ();
96 +      int getTrkIsMatchedDeadEcal (const BNtrack* track1);
97 +
98        //BNskimbits
99        //BNtrigobj
100  
# Line 82 | Line 109 | class OSUAnalysis : public edm::EDAnalyz
109        //inner vector corresponds to each object in input collection
110        //bool         tells if object passes or fails cuts
111        typedef map<string, vector < vector<bool> > > flagMap;
112 <      
112 >
113        //counterMap:
114        //string holds input collection type
115        //vector corresponds to each cut
# Line 99 | Line 126 | class OSUAnalysis : public edm::EDAnalyz
126        edm::InputTag tracks_;
127        edm::InputTag genjets_;
128        edm::InputTag mcparticles_;
129 +      edm::InputTag stops_;
130        edm::InputTag primaryvertexs_;
131        edm::InputTag bxlumis_;
132        edm::InputTag photons_;
133        edm::InputTag superclusters_;
134        edm::InputTag triggers_;
135 <
136 <      TVector3 vertexPosition;
137 <
135 >      std::string puFile_;
136 >      std::string deadEcalFile_;
137 >      std::string muonSFFile_;
138 >      std::string dataPU_;
139 >      std::string electronSFID_;
140 >      std::string muonSF_;
141 >      std::string dataset_;
142 >      std::string datasetType_;
143        vector<edm::ParameterSet> channels_;
144        vector<edm::ParameterSet> histogramSets_;
145 +      bool plotAllObjectsInPassingEvents_;
146 +      bool doPileupReweighting_;
147 +      bool printEventInfo_;  
148 +      bool useTrackCaloRhoCorr_;  // to use the calo-based rho correction for the by-hand calculation of the track isolation energy
149 +      vector<double> stopCTau_;
150 +
151 +      struct DeadEcal {
152 +        double etaEcal;
153 +        double phiEcal;
154 +      };
155 +
156 +      std::vector<DeadEcal> DeadEcalVec;
157 +
158 +
159 +      //Collections
160 +      edm::Handle<BNtriggerCollection> triggers;
161 +      edm::Handle<BNjetCollection> jets;
162 +      edm::Handle<BNmuonCollection> muons;
163 +      edm::Handle<BNelectronCollection> electrons;
164 +      edm::Handle<BNeventCollection> events;
165 +      edm::Handle<BNtauCollection> taus;
166 +      edm::Handle<BNmetCollection> mets;
167 +      edm::Handle<BNtrackCollection> tracks;
168 +      edm::Handle<BNgenjetCollection> genjets;
169 +      edm::Handle<BNmcparticleCollection> mcparticles;
170 +      edm::Handle<BNprimaryvertexCollection> primaryvertexs;
171 +      edm::Handle<BNbxlumiCollection> bxlumis;
172 +      edm::Handle<BNphotonCollection> photons;
173 +      edm::Handle<BNsuperclusterCollection> superclusters;
174 +      edm::Handle<double> rhokt6CaloJetsHandle_;  
175 +
176 +      flagMap cumulativeFlags;
177  
178        vector<map<string, TH1D*> > oneDHists_;
179        vector<map<string, TH2D*> > twoDHists_;
# Line 119 | Line 184 | class OSUAnalysis : public edm::EDAnalyz
184        vector<CutFlow *> cutFlows_;
185  
186        typedef struct {
187 <        string name;
188 <        string title;
189 <        vector<double> bins;
190 <        string inputCollection;
191 <        vector<string> inputVariables;
127 <        string function;
187 >        string name;
188 >        string title;
189 >        vector<double> bins;
190 >        string inputCollection;
191 >        vector<string> inputVariables;
192        } histogram;
193  
194        struct cut {
195 <        string inputCollection;
196 <        string variable;
197 <        string comparativeOperator;
198 <        string function;
199 <        double cutValue;
200 <        int    numberRequired;
201 <        string eventComparativeOperator;
202 <        string name;
195 >        string inputCollection;
196 >        vector<string> functions;
197 >        vector<string> variables;
198 >        vector<string> comparativeOperators;// >, <, =, etc.
199 >        vector<double> cutValues;
200 >        vector<string> logicalOperators;//and, or
201 >        int numSubcuts;
202 >        int    numberRequired;
203 >        string eventComparativeOperator;
204 >        string name;
205        };
206  
207        struct channel {
208 <        string name;
209 <        vector<string> triggers;
210 <        vector<cut> cuts;
208 >        string name;
209 >        vector<string> triggers;
210 >        vector<cut> cuts;
211        };
212  
213 +      vector<string>  objectsToGet;
214 +      vector<string>  objectsToCut;
215        vector<string>  objectsToPlot;
148 –      vector<string>  allNecessaryObjects;
216        vector<channel> channels;
217        vector<histogram> histograms;
218  
219 +      PUWeight *puWeight_;
220 +      MuonSFWeight *muonSFWeight_;
221 +      ElectronSFWeight *electronSFWeight_;
222 + #ifdef DISPLACED_SUSY
223 +      CTauWeight *cTauWeight_;
224 + #endif
225 +      double cTauScaleFactor_;
226 +
227        template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string);
228 <      template <class InputCollection>  void fillHistogram(TH1*, histogram, InputCollection);
228 >      template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, vector<bool>, vector<bool>, string);
229 >
230 >      template <class InputCollection> void fill1DHistogram(TH1*, histogram, InputCollection, vector<bool>, double);
231 >      template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double);
232 >      template <class InputCollection> void fill2DHistogram(TH2*, histogram, InputCollection, vector<bool>, double);
233 >      template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double);
234 >
235 >      template <class InputObject> int getGenMatchedParticleIndex(InputObject);
236 >      int getPdgIdBinValue(int);
237 >      int findTauMotherIndex(const BNmcparticle*);
238 >
239 >      template <class InputObject> double getGenDeltaRLowest(InputObject);
240  
241  
242 +      const BNprimaryvertex *chosenVertex ();
243 +      const BNmet *chosenMET ();
244 +      const BNelectron *chosenElectron ();
245 +      const BNmuon *chosenMuon ();
246  
247    };
248  

Diff Legend

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