float ReconstructionMCMatching::DrBHad(Interpretation& i){
float dr=999;
if(b_had.Pt()>0.01 && i.BHad().Pt()>0.01){
dr=b_had.DeltaR(i.BHad());
}
return dr;
}
bool ReconstructionMCMatching::MatchBHad(Interpretation& i){
if(setupMC==0){
cerr << "need to setup mc matching first!" << endl;
return false;
}
bool matchBhad=b_had.E()>1.&&i.BHad().E()>1.&&b_had.DeltaR(i.BHad())<maxdr;
// cout << "dr bhad " << b_had.DeltaR(i.BHad()) << endl;
return matchBhad;
}
float ReconstructionMCMatching::SumDrTTH(Interpretation& i){
if(setupMC<2){
cerr << "need to setup mc matching (incl Higgs) first!" << endl;
return false;
}
float sumDr=0;
sumDr+=b_had.DeltaR(i.BHad());
sumDr+=q1.DeltaR(i.Q1());
sumDr+=q2.DeltaR(i.Q2());
sumDr+=b_lep.DeltaR(i.BLep());
sumDr+=lep.DeltaR(i.Lep());
sumDr+=nu.DeltaR(i.Nu());
sumDr+=b1.DeltaR(i.B1());
sumDr+=b2.DeltaR(i.B1());
return sumDr;
}