1485 |
|
else if(variable == "puJetId_loose_simple") value = object->puJetId_loose_simple; |
1486 |
|
else if(variable == "puJetId_loose_cutbased") value = object->puJetId_loose_cutbased; |
1487 |
|
|
1488 |
+ |
//user defined variable |
1489 |
+ |
else if(variable == "disappTrkLeadingJetID") { |
1490 |
+ |
value = object->pt > 110 |
1491 |
+ |
&& fabs(object->eta) < 2.4 |
1492 |
+ |
&& object->chargedHadronEnergyFraction > 0.2 |
1493 |
+ |
&& object->neutralHadronEnergyFraction < 0.7 |
1494 |
+ |
&& object->chargedEmEnergyFraction < 0.5 |
1495 |
+ |
&& object->neutralEmEnergyFraction < 0.7; |
1496 |
+ |
} |
1497 |
+ |
|
1498 |
+ |
else if(variable == "disappTrkSubLeadingJetID") { |
1499 |
+ |
value = object->pt > 30 |
1500 |
+ |
&& fabs(object->eta) < 4.5 |
1501 |
+ |
&& object->neutralHadronEnergyFraction < 0.7 |
1502 |
+ |
&& object->chargedEmEnergyFraction < 0.5; |
1503 |
+ |
} |
1504 |
+ |
|
1505 |
|
|
1506 |
|
else{cout << "WARNING: invalid variable '" << variable << "'\n"; value = -999;} |
1507 |
|
|
1668 |
|
else if(variable == "time_ndof") value = object->time_ndof; |
1669 |
|
|
1670 |
|
//user-defined variables |
1671 |
+ |
else if(variable == "looseID") { |
1672 |
+ |
value = object->pt > 10 |
1673 |
+ |
&& object->isGlobalMuon > 0 |
1674 |
+ |
&& object->isTrackerMuon > 0; |
1675 |
+ |
} |
1676 |
|
else if(variable == "correctedD0VertexErr") value = hypot (object->tkD0err, hypot (chosenVertex ()->xError, chosenVertex ()->yError)); |
1677 |
|
else if(variable == "correctedD0VertexSig") value = object->correctedD0Vertex / hypot (object->tkD0err, hypot (chosenVertex ()->xError, chosenVertex ()->yError)); |
1678 |
|
else if(variable == "detIso") value = (object->trackIsoDR03) / object->pt; |
2143 |
|
} |
2144 |
|
} |
2145 |
|
|
2146 |
+ |
else if(variable == "looseID"){ |
2147 |
+ |
value = object->pt > 10 |
2148 |
+ |
&& object->mvaNonTrigV0 > 0; |
2149 |
+ |
} |
2150 |
|
else if(variable == "correctedD0VertexInEBPositiveCharge"){ |
2151 |
|
if(fabs(object->eta) < 0.8 && object->charge > 0) value = object->correctedD0Vertex; |
2152 |
|
else value = -999; |
2417 |
|
else if(variable == "HPSdecayModeFinding") value = object->HPSdecayModeFinding; |
2418 |
|
else if(variable == "leadingTrackValid") value = object->leadingTrackValid; |
2419 |
|
|
2420 |
< |
|
2420 |
> |
else if (variable == "looseHadronicID") { |
2421 |
> |
value = object->pt > 10 |
2422 |
> |
&& object->eta < 2.3 |
2423 |
> |
&& object->HPSbyLooseCombinedIsolationDeltaBetaCorr > 0 |
2424 |
> |
&& object->HPSdecayModeFinding > 0 |
2425 |
> |
&& object->HPSagainstElectronLoose > 0 |
2426 |
> |
&& object->HPSagainstMuonTight > 0; |
2427 |
> |
} |
2428 |
|
|
2429 |
|
else if(variable == "genDeltaRLowest") value = getGenDeltaRLowest(object); |
2430 |
|
|
2622 |
|
pt = object->pt; |
2623 |
|
value = vz - (vx * px + vy * py)/pt * (pz/pt); |
2624 |
|
} |
2625 |
+ |
|
2626 |
+ |
|
2627 |
|
else if(variable == "genDeltaRLowest") value = getGenDeltaRLowest(object); |
2628 |
|
|
2629 |
|
else if(variable == "genMatchedPdgId"){ |