13 |
|
BaseMod(name,title), |
14 |
|
fMuonBranchName(Names::gkMuonBrn), |
15 |
|
fCleanMuonsName(ModNames::gkCleanMuonsName), |
16 |
+ |
fVertexName("PrimaryVertexesBeamSpot"), |
17 |
|
fMuonIDType("Loose"), |
18 |
|
fMuonIsoType("TrackCaloSliding"), |
19 |
|
fMuonClassType("Global"), |
21 |
|
fCaloIsolationCut(3.0), |
22 |
|
fCombIsolationCut(5.0), |
23 |
|
fMuonPtMin(10), |
24 |
< |
fMuons(0) |
24 |
> |
fD0Cut(0.025), |
25 |
> |
fMuIDType(kIdUndef), |
26 |
> |
fMuIsoType(kIsoUndef), |
27 |
> |
fMuClassType(kClassUndef), |
28 |
> |
fMuons(0), |
29 |
> |
fVertices(0), |
30 |
> |
fMuonTools(0) |
31 |
|
{ |
32 |
|
// Constructor. |
33 |
|
} |
38 |
|
// Process entries of the tree. |
39 |
|
|
40 |
|
LoadBranch(fMuonBranchName); |
41 |
+ |
LoadBranch(fVertexName); |
42 |
|
|
43 |
|
MuonOArr *CleanMuons = new MuonOArr; |
44 |
|
CleanMuons->SetName(fCleanMuonsName); |
127 |
|
if (!isopass) |
128 |
|
continue; |
129 |
|
|
130 |
+ |
// d0 cut |
131 |
+ |
double d0_real = 99999; |
132 |
+ |
for(uint i0 = 0; i0 < fVertices->GetEntries(); i0++) { |
133 |
+ |
double pD0 = mu->GlobalTrk()->D0Corrected(*fVertices->At(i0)); |
134 |
+ |
if(TMath::Abs(pD0) < TMath::Abs(d0_real)) d0_real = TMath::Abs(pD0); |
135 |
+ |
} |
136 |
+ |
if(d0_real >= fD0Cut) continue; |
137 |
+ |
|
138 |
|
// add good muon |
139 |
|
CleanMuons->Add(mu); |
140 |
|
} |
153 |
|
// we just request the muon collection branch. |
154 |
|
|
155 |
|
ReqBranch(fMuonBranchName, fMuons); |
156 |
+ |
ReqBranch(fVertexName, fVertices); |
157 |
|
|
158 |
|
fMuonTools = new MuonTools; |
159 |
|
|