|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
|
#include "DataFormats/METReco/interface/CSCHaloData.h"
/*
[class]: CSCHaloData
[authors]: R. Remington, The University of Florida
[description]: See CSCHaloData.h
[date]: October 15, 2009
*/
using namespace reco;
CSCHaloData::CSCHaloData() {
nTriggers_PlusZ = 0;
nTriggers_MinusZ = 0;
nTracks_PlusZ = 0;
nTracks_MinusZ = 0;
HLTAccept = false;
nOutOfTimeTriggers_PlusZ = 0;
nOutOfTimeTriggers_MinusZ = 0;
nOutOfTimeHits = 0;
nTracks_Small_dT = 0;
nTracks_Small_beta = 0;
nTracks_Small_dT_Small_beta = 0;
// MLR
nFlatHaloSegments = 0;
// End MLR
nTriggers_PlusZ_TrkMuUnVeto = 0;
nTriggers_MinusZ_TrkMuUnVeto = 0;
HLTAccept = false;
segments_in_both_endcaps = false;
nFlatHaloSegments_TrkMuUnVeto = 0;
segments_in_both_endcaps_loose_TrkMuUnVeto = false;
segments_in_both_endcaps_loose_dtcut_TrkMuUnVeto = false;
segmentiscalomatched = false;
segmentisHcalomatched = false;
segmentisEBcalomatched = false;
segmentisEEcalomatched = false;
}
int CSCHaloData::NumberOfHaloTriggers(HaloData::Endcap z) const {
if (z == HaloData::plus)
return nTriggers_PlusZ;
else if (z == HaloData::minus)
return nTriggers_MinusZ;
else
return nTriggers_MinusZ + nTriggers_PlusZ;
}
int CSCHaloData::NumberOfHaloTriggers_TrkMuUnVeto(HaloData::Endcap z) const {
if (z == HaloData::plus)
return nTriggers_PlusZ_TrkMuUnVeto;
else if (z == HaloData::minus)
return nTriggers_MinusZ_TrkMuUnVeto;
else
return nTriggers_MinusZ_TrkMuUnVeto + nTriggers_PlusZ_TrkMuUnVeto;
}
short int CSCHaloData::NumberOfOutOfTimeTriggers(HaloData::Endcap z) const {
if (z == HaloData::plus)
return nOutOfTimeTriggers_PlusZ;
else if (z == HaloData::minus)
return nOutOfTimeTriggers_MinusZ;
else
return nOutOfTimeTriggers_PlusZ + nOutOfTimeTriggers_MinusZ;
}
int CSCHaloData::NumberOfHaloTracks(HaloData::Endcap z) const {
int n = 0;
for (unsigned int i = 0; i < TheTrackRefs.size(); i++) {
edm::Ref<reco::TrackCollection> iTrack(TheTrackRefs[i]);
// Does the track go through both endcaps ?
bool Traversing = (iTrack->outerPosition().z() > 0 && iTrack->innerPosition().z() < 0) ||
(iTrack->outerPosition().z() < 0 && iTrack->innerPosition().z() > 0);
// Does the track go through only +Z endcap ?
bool PlusZ = (iTrack->outerPosition().z() > 0 && iTrack->innerPosition().z() > 0);
// Does the track go through only -Z endcap ?
bool MinusZ = (iTrack->outerPosition().z() < 0 && iTrack->innerPosition().z() < 0);
if ((z == HaloData::plus) && (PlusZ || Traversing))
n++;
else if ((z == HaloData::minus) && (MinusZ || Traversing))
n++;
else if ((z == HaloData::both) && (PlusZ || MinusZ || Traversing))
n++;
}
return n;
}
|