| 3024 |
|
} |
| 3025 |
|
value = trkMuonDeltaRMin; |
| 3026 |
|
} |
| 3027 |
< |
|
| 3028 |
< |
else if(variable == "isPassMuonLooseIDVeto") { |
| 3029 |
< |
|
| 3030 |
< |
// calculate minimum deltaR between track and any other loose-Id |
| 3031 |
< |
double trkMuonDeltaRMin = 99.; |
| 3032 |
< |
if (!muons.product()) clog << "ERROR: cannot find deltaRMinMuonLooseId because muons collection is not initialized." << endl; |
| 3033 |
< |
for (uint imuon = 0; imuon<muons->size(); imuon++) { |
| 3034 |
< |
string empty = ""; |
| 3035 |
< |
double isLooseIdMuon = valueLookup(&muons->at(imuon), "looseID", "", empty); |
| 3036 |
< |
if (!isLooseIdMuon) continue; // only consider muons that pass the looseID |
| 3037 |
< |
double muonEta = valueLookup(&muons->at(imuon), "eta", "", empty); |
| 3038 |
< |
double muonPhi = valueLookup(&muons->at(imuon), "phi", "", empty); |
| 3039 |
< |
double trkMuonDeltaR = deltaR(object->eta, object->phi, muonEta, muonPhi); |
| 3040 |
< |
if (trkMuonDeltaR < trkMuonDeltaRMin) trkMuonDeltaRMin = trkMuonDeltaR; |
| 3041 |
< |
} |
| 3042 |
< |
if (trkMuonDeltaRMin < 0.15) value = 1; |
| 3043 |
< |
if (trkMuonDeltaRMin > 0.15) value = 0; |
| 3027 |
> |
else if(variable == "isPassMuonLooseID") { |
| 3028 |
> |
string empty = ""; |
| 3029 |
> |
double trkMuonDeltaRMin = valueLookup(object, "deltaRMinMuonLooseId", "", empty); |
| 3030 |
> |
if (trkMuonDeltaRMin < 0.15) value = 1; |
| 3031 |
> |
else value = 0; |
| 3032 |
|
} |
| 3045 |
– |
|
| 3046 |
– |
|
| 3047 |
– |
|
| 3033 |
|
else if(variable == "deltaRMinSecMuonLooseIdGlobal") { |
| 3034 |
|
// calculate minimum deltaR between track and any other loose-Id muon |
| 3035 |
|
double trkMuonDeltaRMin = 99.; |
