699 |
|
if( ZCandidates[z1index].first == ZCandidates[z2index].second ) continue; |
700 |
|
if( ZCandidates[z1index].second == ZCandidates[z2index].first ) continue; |
701 |
|
if( ZCandidates[z1index].second == ZCandidates[z2index].second ) continue; |
702 |
+ |
//ghost removal |
703 |
+ |
if(mithep::MathUtils::DeltaR(ZCandidatesLeptons[z1index].first.vec.Phi(),ZCandidatesLeptons[z1index].first.vec.Eta(),ZCandidatesLeptons[z2index].first.vec.Phi(),ZCandidatesLeptons[z2index].first.vec.Eta()) < 0.02) continue; |
704 |
+ |
if(mithep::MathUtils::DeltaR(ZCandidatesLeptons[z1index].first.vec.Phi(),ZCandidatesLeptons[z1index].first.vec.Eta(),ZCandidatesLeptons[z2index].second.vec.Phi(),ZCandidatesLeptons[z2index].second.vec.Eta()) < 0.02) continue; |
705 |
+ |
if(mithep::MathUtils::DeltaR(ZCandidatesLeptons[z1index].second.vec.Phi(),ZCandidatesLeptons[z1index].second.vec.Eta(),ZCandidatesLeptons[z2index].first.vec.Phi(),ZCandidatesLeptons[z2index].first.vec.Eta()) < 0.02) continue; |
706 |
+ |
if(mithep::MathUtils::DeltaR(ZCandidatesLeptons[z1index].second.vec.Phi(),ZCandidatesLeptons[z1index].second.vec.Eta(),ZCandidatesLeptons[z2index].second.vec.Phi(),ZCandidatesLeptons[z2index].second.vec.Eta()) < 0.02) continue; |
707 |
+ |
|
708 |
|
ZZCandidates.push_back(std::pair<int,int> (z1index,z2index)); |
709 |
|
// Z2type = abs(lepvec[ZCandidates[z2index].first].type); |
710 |
|
Z2type = abs(ZCandidatesLeptons[z2index].first.type); |