ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/dhidas/OSUAnalysis/Tools/test/TestElectron.h
Revision: 1.1.1.1 (vendor branch)
Committed: Thu Dec 1 16:28:47 2011 UTC (13 years, 5 months ago) by dhidas
Content type: text/plain
Branch: dhidas, MAIN
CVS Tags: START, HEAD
Changes since 1.1: +0 -0 lines
Log Message: 
osu copy modified

File Contents

# User Rev Content
1 dhidas 1.1 #include "cute/cute.h"
2     #include "cute/cute_suite.h"
3    
4     #include "../interface/RecoObjects/Electron.h"
5     #include "TestObjectFactory.h"
6    
7     using namespace BAT;
8    
9     struct TestElectron {
10     private:
11     ElectronPointer isolatedElectron;
12     ElectronPointer goodElectron;
13     ElectronPointer badEtElectron;
14     ElectronPointer badEtaElectron;
15     ElectronPointer badInCrackElectron;
16     ElectronPointer badD0Electron;
17     ElectronPointer badElectronWithMissingInnerLayerHit;
18     ElectronPointer badElectronWithPartnerTrack;
19     ElectronPointer looseElectron;
20     ElectronPointer badLooseElectronNoID;
21     ElectronPointer badElectronNoID;
22    
23     public:
24     TestElectron() :
25     isolatedElectron(TestObjectFactory::goodIsolatedElectron()),
26     goodElectron(TestObjectFactory::goodCaloElectron()),
27     badEtElectron(TestObjectFactory::badEtElectron()),
28     badEtaElectron(TestObjectFactory::badEtaElectron()),
29     badInCrackElectron(TestObjectFactory::badInCrackElectron()),
30     badD0Electron(TestObjectFactory::badD0Electron()),
31     badElectronWithMissingInnerLayerHit(TestObjectFactory::electronWithMissingInnerLayerHit()),
32     badElectronWithPartnerTrack(TestObjectFactory::electronWithMatchingPartnerTrack()),
33     looseElectron(TestObjectFactory::goodLooseElectron()),
34     badLooseElectronNoID(TestObjectFactory::badLooseElectronNoID()),
35     badElectronNoID(TestObjectFactory::badElectronNoID()) {
36     }
37    
38     public:
39     void testEcalIsolation() {
40     ASSERT_EQUAL_DELTA(0.3, isolatedElectron->ecalIsolation(), 0.01);
41     }
42    
43     void testHcalIsolation() {
44     ASSERT_EQUAL_DELTA(0.5, isolatedElectron->hcalIsolation(), 0.01);
45     }
46    
47     void testTrackerIsolation() {
48     ASSERT_EQUAL_DELTA(0.4, isolatedElectron->trackerIsolation(), 0.01);
49     }
50    
51     void testBadEtElectron() {
52     ASSERT(badEtElectron->isGood() == false);
53     }
54    
55     void testBadEtaElectron() {
56     ASSERT(badEtaElectron->isGood() == false);
57     }
58    
59     void testBadInCrackElectron() {
60     ASSERT(badInCrackElectron->isGood() == false);
61     }
62    
63     void testElectronFromConversion() {
64     ASSERT(badElectronWithMissingInnerLayerHit->isFromConversion());
65     }
66    
67     void testElectronFromConversionWithPartnerTrack() {
68     ASSERT(badElectronWithPartnerTrack->isTaggedAsConversion(0.2,0.2));
69     }
70    
71     void testBadD0Electron() {
72     ASSERT(badD0Electron->isGood() == false);
73     }
74    
75     void testLooseElectron() {
76     ASSERT_EQUAL(true, badD0Electron->isLoose());
77     ASSERT_EQUAL(true, looseElectron->isLoose());
78     }
79    
80     void testRelativeIsolation() {
81     float isoSum = isolatedElectron->trackerIsolation() + isolatedElectron->ecalIsolation()
82     + isolatedElectron->hcalIsolation();
83     ASSERT_EQUAL_DELTA(isoSum / isolatedElectron->et(), isolatedElectron->relativeIsolation(), 0.01);
84     ASSERT(isolatedElectron->isIsolated());
85     }
86    
87     void testGoodElectron() {
88     ASSERT(goodElectron->isGood());
89     }
90    
91     void testBadLooseElectronNoID() {
92     ASSERT(badLooseElectronNoID->isLoose() == false);
93     }
94    
95     void testBadElectronNoID() {
96     ASSERT(badElectronNoID->isGood() == false);
97     }
98    
99     void testVBTFW70Barrel() {
100     Electron electron = Electron(50, 40, 30, 0);
101     electron.setSigmaIEtaIEta(0.01 - 0.01
102     * 0.01);
103     electron.setDPhiIn(0.03 - 0.01 * 0.03);
104     electron.setDEtaIn(0.004 - 0.01 * 0.004);
105     electron.setHadOverEm(0.025 - 0.01 * 0.025);
106     electron.setSuperClusterEta(0);
107     assert(electron.isInBarrelRegion());
108     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
109     }
110    
111     void testVBTFW70BarrelNegativeValues() {
112     Electron electron = Electron(50, 40, 30, 0);
113     electron.setSigmaIEtaIEta(-0.01 + 0.01
114     * 0.01);
115     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
116     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
117     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
118     electron.setSuperClusterEta(0);
119     assert(electron.isInBarrelRegion());
120     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
121     }
122    
123     void testVBTFW70BarrelFailSigmaIEtaIEta() {
124     Electron electron = Electron(50, 3, 0, 40);
125     electron.setSigmaIEtaIEta(0.01 + 0.01
126     * 0.01);
127     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
128     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
129     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
130     electron.setSuperClusterEta(0);
131     assert(electron.isInBarrelRegion());
132     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
133     }
134    
135     void testVBTFW70BarrelNegativeSigmaIEtaIEta() {
136     Electron electron = Electron(50, 3, 0, 40);
137     electron.setSigmaIEtaIEta(-0.01 - 0.01
138     * 0.01);
139     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
140     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
141     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
142     electron.setSuperClusterEta(0);
143     assert(electron.isInBarrelRegion());
144     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
145     }
146    
147     void testVBTFW70BarrelFailDPhiIn() {
148     Electron electron = Electron(50, 3, 0, 40);
149     electron.setSigmaIEtaIEta(-0.01 + 0.01
150     * 0.01);
151     electron.setDPhiIn(0.03 + 0.01 * 0.03);
152     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
153     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
154     electron.setSuperClusterEta(0);
155     assert(electron.isInBarrelRegion());
156     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
157     }
158    
159     void testVBTFW70BarrelFailNegativeDPhiIn() {
160     Electron electron = Electron(50, 3, 0, 40);
161     electron.setSigmaIEtaIEta(-0.01 + 0.01
162     * 0.01);
163     electron.setDPhiIn(-0.03 - 0.01 * 0.03);
164     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
165     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
166     electron.setSuperClusterEta(0);
167     assert(electron.isInBarrelRegion());
168     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
169     }
170    
171     void testVBTFW70BarrelFailDEtaIn() {
172     Electron electron = Electron(50, 3, 0, 40);
173     electron.setSigmaIEtaIEta(-0.01 + 0.01
174     * 0.01);
175     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
176     electron.setDEtaIn(0.004 + 0.01 * 0.004);
177     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
178     electron.setSuperClusterEta(0);
179     assert(electron.isInBarrelRegion());
180     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
181     }
182    
183     void testVBTFW70BarrelFailNegativeDEtaIn() {
184     Electron electron = Electron(50, 3, 0, 40);
185     electron.setSigmaIEtaIEta(-0.01 + 0.01
186     * 0.01);
187     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
188     electron.setDEtaIn(-0.004 - 0.01 * 0.004);
189     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
190     electron.setSuperClusterEta(0);
191     assert(electron.isInBarrelRegion());
192     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
193     }
194    
195     void testVBTFW70BarrelFailHadOverEm() {
196     Electron electron = Electron(50, 3, 0, 40);
197     electron.setSigmaIEtaIEta(-0.01 + 0.01
198     * 0.01);
199     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
200     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
201     electron.setHadOverEm(0.025 + 0.01 * 0.025);
202     electron.setSuperClusterEta(0);
203     assert(electron.isInBarrelRegion());
204     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
205     }
206    
207     void testVBTFW70BarrelNegativeHadOverEm() {
208     Electron electron = Electron(50, 3, 0, 40);
209     electron.setSigmaIEtaIEta(-0.01 + 0.01
210     * 0.01);
211     electron.setDPhiIn(-0.03 + 0.01 * 0.03);
212     electron.setDEtaIn(-0.004 + 0.01 * 0.004);
213     electron.setHadOverEm(-0.025 - 0.01 * 0.025);
214     electron.setSuperClusterEta(0);
215     assert(electron.isInBarrelRegion());
216     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
217     }
218    
219     void testVBTFW70Endcap() {
220     Electron electron = Electron(50, 3, 0, 40);
221     electron.setSigmaIEtaIEta(0.03 - 0.01
222     * 0.03);
223     electron.setDPhiIn(0.02 - 0.01 * 0.02);
224     electron.setDEtaIn(0.005 - 0.01 * 0.005);
225     electron.setHadOverEm(0.025 - 0.01 * 0.025);
226     electron.setSuperClusterEta(2.6);
227     assert(electron.isInEndCapRegion());
228     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
229     }
230    
231     void testVBTFW70EndcapNegativeValues() {
232     Electron electron = Electron(50, 3, 0, 40);
233     electron.setSigmaIEtaIEta(-0.03 + 0.01
234     * 0.03);
235     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
236     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
237     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
238     electron.setSuperClusterEta(2.6);
239     assert(electron.isInEndCapRegion());
240     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
241     }
242    
243     void testVBTFW70EndcapFailSigmaIEtaIEta() {
244     Electron electron = Electron(50, 3, 0, 40);
245     electron.setSigmaIEtaIEta(0.03 + 0.01
246     * 0.03);
247     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
248     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
249     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
250     electron.setSuperClusterEta(2.6);
251     assert(electron.isInEndCapRegion());
252     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
253     }
254    
255     void testVBTFW70EndcapNegativeSigmaIEtaIEta() {
256     Electron electron = Electron(50, 3, 0, 40);
257     electron.setSigmaIEtaIEta(-0.03 - 0.01
258     * 0.03);
259     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
260     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
261     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
262     electron.setSuperClusterEta(2.6);
263     assert(electron.isInEndCapRegion());
264     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
265     }
266    
267     void testVBTFW70EndcapFailDPhiIn() {
268     Electron electron = Electron(50, 3, 0, 40);
269     electron.setSigmaIEtaIEta(-0.03 + 0.01
270     * 0.03);
271     electron.setDPhiIn(0.02 + 0.01 * 0.02);
272     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
273     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
274     electron.setSuperClusterEta(2.6);
275     assert(electron.isInEndCapRegion());
276     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
277     }
278    
279     void testVBTFW70EndcapFailNegativeDPhiIn() {
280     Electron electron = Electron(50, 3, 0, 40);
281     electron.setSigmaIEtaIEta(-0.03 + 0.01
282     * 0.03);
283     electron.setDPhiIn(-0.02 - 0.01 * 0.02);
284     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
285     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
286     electron.setSuperClusterEta(2.6);
287     assert(electron.isInEndCapRegion());
288     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
289     }
290    
291     void testVBTFW70EndcapDEtaIn() {
292     Electron electron = Electron(50, 3, 0, 40);
293     electron.setSigmaIEtaIEta(-0.03 + 0.01
294     * 0.03);
295     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
296     electron.setDEtaIn(0.005 - 0.01 * 0.005);
297     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
298     electron.setSuperClusterEta(2.6);
299     assert(electron.isInEndCapRegion());
300     //expect true, dEtaIn cut is not used in the endcap at the moment
301     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
302     }
303    
304     void testVBTFW70EndcapNegativeDEtaIn() {
305     Electron electron = Electron(50, 3, 0, 40);
306     electron.setSigmaIEtaIEta(-0.03 + 0.01
307     * 0.03);
308     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
309     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
310     electron.setHadOverEm(-0.025 + 0.01 * 0.025);
311     electron.setSuperClusterEta(2.6);
312     assert(electron.isInEndCapRegion());
313     //expect true, dEtaIn cut is not used in the endcap at the moment
314     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
315     }
316    
317     void testVBTFW70EndcapFailHadOverEm() {
318     Electron electron = Electron(50, 3, 0, 40);
319     electron.setSigmaIEtaIEta(-0.03 + 0.01
320     * 0.03);
321     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
322     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
323     electron.setHadOverEm(0.025 + 0.01 * 0.025);
324     electron.setSuperClusterEta(2.6);
325     assert(electron.isInEndCapRegion());
326     ASSERT_EQUAL(false, electron.VBTF_W70_ElectronID());
327     }
328    
329     void testVBTFW70EndcapNegativeHadOverEm() {
330     Electron electron = Electron(50, 3, 0, 40);
331     electron.setSigmaIEtaIEta(-0.03 + 0.01
332     * 0.03);
333     electron.setDPhiIn(-0.02 + 0.01 * 0.02);
334     electron.setDEtaIn(-0.005 + 0.01 * 0.005);
335     electron.setHadOverEm(-0.025 - 0.01 * 0.025);
336     electron.setSuperClusterEta(2.6);
337     assert(electron.isInEndCapRegion());
338     ASSERT_EQUAL(true, electron.VBTF_W70_ElectronID());
339     }
340    
341     void testVBTFW95Barrel() {
342     Electron electron = Electron(50, 40, 30, 0);
343     electron.setSigmaIEtaIEta(0.01 - 0.01
344     * 0.01);
345     electron.setDPhiIn(0.8 - 0.01 * 0.8);
346     electron.setDEtaIn(0.007 - 0.01 * 0.007);
347     electron.setHadOverEm(0.15 - 0.01 * 0.15);
348     electron.setSuperClusterEta(0);
349     assert(electron.isInBarrelRegion());
350     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
351     }
352    
353     void testVBTFW95BarrelNegativeValues() {
354     Electron electron = Electron(50, 40, 30, 0);
355     electron.setSigmaIEtaIEta(-0.01 + 0.01
356     * 0.01);
357     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
358     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
359     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
360     electron.setSuperClusterEta(0);
361     assert(electron.isInBarrelRegion());
362     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
363     }
364    
365     void testVBTFW95BarrelFailSigmaIEtaIEta() {
366     Electron electron = Electron(50, 3, 0, 40);
367     electron.setSigmaIEtaIEta(0.01 + 0.01
368     * 0.01);
369     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
370     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
371     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
372     electron.setSuperClusterEta(0);
373     assert(electron.isInBarrelRegion());
374     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
375     }
376    
377     void testVBTFW95BarrelNegativeSigmaIEtaIEta() {
378     Electron electron = Electron(50, 3, 0, 40);
379     electron.setSigmaIEtaIEta(-0.01 - 0.01
380     * 0.01);
381     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
382     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
383     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
384     electron.setSuperClusterEta(0);
385     assert(electron.isInBarrelRegion());
386     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
387     }
388    
389     void testVBTFW95BarrelFailDPhiIn() {
390     Electron electron = Electron(50, 3, 0, 40);
391     electron.setSigmaIEtaIEta(-0.01 + 0.01
392     * 0.01);
393     electron.setDPhiIn(0.8 + 0.01 * 0.8);
394     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
395     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
396     electron.setSuperClusterEta(0);
397     assert(electron.isInBarrelRegion());
398     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
399     }
400    
401     void testVBTFW95BarrelFailNegativeDPhiIn() {
402     Electron electron = Electron(50, 3, 0, 40);
403     electron.setSigmaIEtaIEta(-0.01 + 0.01
404     * 0.01);
405     electron.setDPhiIn(-0.8 - 0.01 * 0.8);
406     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
407     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
408     electron.setSuperClusterEta(0);
409     assert(electron.isInBarrelRegion());
410     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
411     }
412    
413     void testVBTFW95BarrelFailDEtaIn() {
414     Electron electron = Electron(50, 3, 0, 40);
415     electron.setSigmaIEtaIEta(-0.01 + 0.01
416     * 0.01);
417     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
418     electron.setDEtaIn(0.007 + 0.01 * 0.007);
419     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
420     electron.setSuperClusterEta(0);
421     assert(electron.isInBarrelRegion());
422     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
423     }
424    
425     void testVBTFW95BarrelFailNegativeDEtaIn() {
426     Electron electron = Electron(50, 3, 0, 40);
427     electron.setSigmaIEtaIEta(-0.01 + 0.01
428     * 0.01);
429     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
430     electron.setDEtaIn(-0.007 - 0.01 * 0.007);
431     electron.setHadOverEm(-0.15 + 0.01 * 0.15);
432     electron.setSuperClusterEta(0);
433     assert(electron.isInBarrelRegion());
434     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
435     }
436    
437     void testVBTFW95BarrelFailHadOverEm() {
438     Electron electron = Electron(50, 3, 0, 40);
439     electron.setSigmaIEtaIEta(-0.01 + 0.01
440     * 0.01);
441     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
442     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
443     electron.setHadOverEm(0.15 + 0.01 * 0.15);
444     electron.setSuperClusterEta(0);
445     assert(electron.isInBarrelRegion());
446     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
447     }
448    
449     void testVBTFW95BarrelNegativeHadOverEm() {
450     Electron electron = Electron(50, 3, 0, 40);
451     electron.setSigmaIEtaIEta(-0.01 + 0.01
452     * 0.01);
453     electron.setDPhiIn(-0.8 + 0.01 * 0.8);
454     electron.setDEtaIn(-0.007 + 0.01 * 0.007);
455     electron.setHadOverEm(-0.15 - 0.01 * 0.15);
456     electron.setSuperClusterEta(0);
457     assert(electron.isInBarrelRegion());
458     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
459     }
460    
461     void testVBTFW95Endcap() {
462     Electron electron = Electron(50, 3, 0, 49);
463     electron.setSigmaIEtaIEta(0.03 - 0.01
464     * 0.03);
465     electron.setDPhiIn(0.7 - 0.01 * 0.7);
466     electron.setDEtaIn(0.01 - 0.01 * 0.01);
467     electron.setHadOverEm(0.07 - 0.01 * 0.07);
468     electron.setSuperClusterEta(2.4);
469    
470     assert(electron.isInEndCapRegion());
471     assert(electron.sigmaIEtaIEta() < 0.03);
472     assert(fabs(electron.dPhiIn()) < 0.7);
473     assert(fabs(electron.dEtaIn()) < 0.01);
474     assert(electron.HadOverEm() < 0.07);
475     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
476     }
477    
478     void testVBTFW95EndcapNegativeValues() {
479     Electron electron = Electron(50, 3, 0, 40);
480     electron.setSigmaIEtaIEta(-0.03 + 0.01
481     * 0.03);
482     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
483     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
484     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
485     electron.setSuperClusterEta(2.4);
486    
487     assert(electron.sigmaIEtaIEta() < 0.03);
488     assert(fabs(electron.dPhiIn()) < 0.7);
489     assert(fabs(electron.dEtaIn()) < 0.01);
490     assert(electron.HadOverEm() < 0.07);
491     assert(electron.isInEndCapRegion());
492     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
493     }
494    
495     void testVBTFW95EndcapFailSigmaIEtaIEta() {
496     Electron electron = Electron(50, 3, 0, 40);
497     electron.setSigmaIEtaIEta(0.03 + 0.01
498     * 0.03);
499     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
500     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
501     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
502     electron.setSuperClusterEta(2.4);
503     assert(electron.isInEndCapRegion());
504     assert(electron.sigmaIEtaIEta() > 0.03);
505     assert(fabs(electron.dPhiIn()) < 0.7);
506     assert(fabs(electron.dEtaIn()) < 0.01);
507     assert(electron.HadOverEm() < 0.07);
508     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
509     }
510    
511     void testVBTFW95EndcapNegativeSigmaIEtaIEta() {
512     Electron electron = Electron(50, 3, 0, 40);
513     electron.setSigmaIEtaIEta(-0.03 - 0.01
514     * 0.03);
515     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
516     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
517     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
518     electron.setSuperClusterEta(2.6);
519     assert(electron.isInEndCapRegion());
520     assert(electron.sigmaIEtaIEta() < 0.03);
521     assert(fabs(electron.dPhiIn()) < 0.7);
522     assert(fabs(electron.dEtaIn()) < 0.01);
523     assert(electron.HadOverEm() < 0.07);
524     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
525     }
526    
527     void testVBTFW95EndcapFailDPhiIn() {
528     Electron electron = Electron(50, 3, 0, 40);
529     electron.setSigmaIEtaIEta(-0.03 + 0.01
530     * 0.03);
531     electron.setDPhiIn(0.7 + 0.01 * 0.7);
532     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
533     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
534     electron.setSuperClusterEta(2.6);
535     assert(electron.sigmaIEtaIEta() < 0.03);
536     assert(fabs(electron.dPhiIn()) > 0.7);
537     assert(fabs(electron.dEtaIn()) < 0.01);
538     assert(electron.HadOverEm() < 0.07);
539     assert(electron.isInEndCapRegion());
540     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
541     }
542    
543     void testVBTFW95EndcapFailNegativeDPhiIn() {
544     Electron electron = Electron(50, 3, 0, 40);
545     electron.setSigmaIEtaIEta(-0.03 + 0.01
546     * 0.03);
547     electron.setDPhiIn(-0.7 - 0.01 * 0.7);
548     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
549     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
550     electron.setSuperClusterEta(2.6);
551     assert(electron.isInEndCapRegion());
552     assert(electron.sigmaIEtaIEta() < 0.03);
553     assert(fabs(electron.dPhiIn()) > 0.7);
554     assert(fabs(electron.dEtaIn()) < 0.01);
555     assert(electron.HadOverEm() < 0.07);
556     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
557     }
558    
559     void testVBTFW95EndcapDEtaIn() {
560     Electron electron = Electron(50, 3, 0, 40);
561     electron.setSigmaIEtaIEta(-0.03 + 0.01
562     * 0.03);
563     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
564     electron.setDEtaIn(0.01 - 0.01 * 0.01);
565     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
566     electron.setSuperClusterEta(2.6);
567     assert(electron.isInEndCapRegion());
568     assert(electron.sigmaIEtaIEta() < 0.03);
569     assert(fabs(electron.dPhiIn()) < 0.7);
570     assert(fabs(electron.dEtaIn()) < 0.01);
571     assert(electron.HadOverEm() < 0.07);
572     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
573     }
574    
575     void testVBTFW95EndcapNegativeDEtaIn() {
576     Electron electron = Electron(50, 3, 0, 40);
577     electron.setSigmaIEtaIEta(-0.03 + 0.01
578     * 0.03);
579     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
580     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
581     electron.setHadOverEm(-0.07 + 0.01 * 0.07);
582     electron.setSuperClusterEta(2.6);
583     assert(electron.isInEndCapRegion());
584     assert(electron.sigmaIEtaIEta() < 0.03);
585     assert(fabs(electron.dPhiIn()) < 0.7);
586     assert(fabs(electron.dEtaIn()) < 0.01);
587     assert(electron.HadOverEm() < 0.07);
588     //expect true, dEtaIn cut is not used in the endcap at the moment
589     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
590     }
591    
592     void testVBTFW95EndcapFailHadOverEm() {
593     Electron electron = Electron(50, 3, 0, 40);
594     electron.setSigmaIEtaIEta(-0.03 + 0.01
595     * 0.03);
596     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
597     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
598     electron.setHadOverEm(0.07 + 0.01 * 0.07);
599     electron.setSuperClusterEta(2.6);
600     assert(electron.isInEndCapRegion());
601     assert(electron.sigmaIEtaIEta() < 0.03);
602     assert(fabs(electron.dPhiIn()) < 0.7);
603     assert(fabs(electron.dEtaIn()) < 0.01);
604     assert(electron.HadOverEm() > 0.07);
605     ASSERT_EQUAL(false, electron.VBTF_W95_ElectronID());
606     }
607    
608     void testVBTFW95EndcapNegativeHadOverEm() {
609     Electron electron = Electron(50, 3, 0, 40);
610     electron.setSigmaIEtaIEta(-0.03 + 0.01
611     * 0.03);
612     electron.setDPhiIn(-0.7 + 0.01 * 0.7);
613     electron.setDEtaIn(-0.01 + 0.01 * 0.01);
614     electron.setHadOverEm(-0.07 - 0.01 * 0.07);
615     electron.setSuperClusterEta(2.6);
616     assert(electron.isInEndCapRegion());
617     assert(electron.sigmaIEtaIEta() < 0.03);
618     assert(fabs(electron.dPhiIn()) < 0.7);
619     assert(fabs(electron.dEtaIn()) < 0.01);
620     assert(electron.HadOverEm() < 0.07);
621     ASSERT_EQUAL(true, electron.VBTF_W95_ElectronID());
622     }
623    
624     void testElectronCharge() {
625     Electron electron = Electron(50, 3, 0, 40);
626     electron.setCharge(2);
627     ASSERT_EQUAL(2, electron.charge());
628     }
629    
630     void testElectronChargeInCopy() {
631     Electron electron = Electron(50, 3, 0, 40);
632     electron.setCharge(2);
633     Electron copyElectron(electron);
634     ASSERT_EQUAL(2, copyElectron.charge());
635     }
636    
637     void testElectronInCollection() {
638     ElectronCollection coll;
639     ElectronPointer electron(new Electron(50, 3, 0, 40));
640     ElectronPointer electron2(new Electron(50, 3, 0, 40));
641     electron->setCharge(2);
642     electron2->setCharge(2);
643     coll.push_back(electron);
644     coll.push_back(electron2);
645     ASSERT_EQUAL(2, coll.front()->charge());
646     }
647    
648     void testElectronInSTDCollection() {
649     std::vector<Electron> coll;
650     Electron electron(50, 3, 0, 40);
651     Electron electron2(50, 3, 0, 40);
652     electron.setCharge(2);
653     electron2.setCharge(2);
654     coll.push_back(electron);
655     coll.push_back(electron2);
656     ASSERT_EQUAL(2, coll.front().charge());
657     }
658    
659     void testElectronSetMass() {
660     Electron electron = Electron(50, 3, 0, 40);
661     electron.setMass(500);
662     ASSERT_EQUAL(500, electron.mass());
663     }
664    
665     void testGSFTrack() {
666     TrackPointer track = TrackPointer(new Track(20, 20, 0, 0));
667     goodElectron->setGSFTrack(track);
668     ASSERT_EQUAL(track, goodElectron->GSFTrack());
669     }
670    
671     };
672     extern cute::suite make_suite_TestElectron() {
673     cute::suite s;
674     s.push_back(CUTE_SMEMFUN(TestElectron, testRelativeIsolation));
675     s.push_back(CUTE_SMEMFUN(TestElectron, testGoodElectron));
676     s.push_back(CUTE_SMEMFUN(TestElectron, testBadEtElectron));
677     s.push_back(CUTE_SMEMFUN(TestElectron, testBadEtaElectron));
678     s.push_back(CUTE_SMEMFUN(TestElectron, testBadD0Electron));
679     s.push_back(CUTE_SMEMFUN(TestElectron, testBadInCrackElectron));
680     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronFromConversion));
681     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronFromConversionWithPartnerTrack));
682     s.push_back(CUTE_SMEMFUN(TestElectron, testEcalIsolation));
683     s.push_back(CUTE_SMEMFUN(TestElectron, testHcalIsolation));
684     s.push_back(CUTE_SMEMFUN(TestElectron, testTrackerIsolation));
685     s.push_back(CUTE_SMEMFUN(TestElectron, testLooseElectron));
686     s.push_back(CUTE_SMEMFUN(TestElectron, testBadLooseElectronNoID));
687     s.push_back(CUTE_SMEMFUN(TestElectron, testBadElectronNoID));
688    
689     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70Barrel));
690     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelNegativeValues));
691     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelFailSigmaIEtaIEta));
692     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelNegativeSigmaIEtaIEta));
693     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelFailDPhiIn));
694     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelFailNegativeDPhiIn));
695     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelFailDEtaIn));
696     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelFailNegativeDEtaIn));
697     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelFailHadOverEm));
698     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70BarrelNegativeHadOverEm));
699    
700     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70Endcap));
701     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapNegativeValues));
702     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapFailSigmaIEtaIEta));
703     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapNegativeSigmaIEtaIEta));
704     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapFailDPhiIn));
705     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapFailNegativeDPhiIn));
706     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapDEtaIn));
707     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapNegativeDEtaIn));
708     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapFailHadOverEm));
709     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW70EndcapNegativeHadOverEm));
710    
711     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95Barrel));
712     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelNegativeValues));
713     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelFailSigmaIEtaIEta));
714     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelNegativeSigmaIEtaIEta));
715     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelFailDPhiIn));
716     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelFailNegativeDPhiIn));
717     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelFailDEtaIn));
718     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelFailNegativeDEtaIn));
719     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelFailHadOverEm));
720     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95BarrelNegativeHadOverEm));
721    
722     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95Endcap));
723     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapNegativeValues));
724     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapFailSigmaIEtaIEta));
725     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapNegativeSigmaIEtaIEta));
726     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapFailDPhiIn));
727     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapFailNegativeDPhiIn));
728     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapDEtaIn));
729     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapNegativeDEtaIn));
730     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapFailHadOverEm));
731     s.push_back(CUTE_SMEMFUN(TestElectron, testVBTFW95EndcapNegativeHadOverEm));
732    
733     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronCharge));
734     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronChargeInCopy));
735     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronInCollection));
736     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronSetMass));
737     s.push_back(CUTE_SMEMFUN(TestElectron, testElectronInSTDCollection));
738     // s.push_back(CUTE_SMEMFUN(TestElectron, testSwissCrossBarrel));
739     // s.push_back(CUTE_SMEMFUN(TestElectron, testSwissCrossEndcap));
740     // s.push_back(CUTE_SMEMFUN(TestElectron, testSwissCrossNotEcalDriven));
741     s.push_back(CUTE_SMEMFUN(TestElectron, testGSFTrack));
742     return s;
743     }