ElectronEfficiencyHelpers Namespace Reference

Functions
CP::CorrectionCode	getEleTruthCharge (const xAOD::Electron &ele, int &truthcharge)
	Get the charge of the original electron.
CP::CorrectionCode	isGoodEle (const xAOD::Electron &ele, bool &goodEle)
	Return true if it's good ele for charge flip measurements.

Function Documentation

getEleTruthCharge()

CP::CorrectionCode ElectronEfficiencyHelpers::getEleTruthCharge	(	const xAOD::Electron &	ele,
		int &	truthcharge )

Get the charge of the original electron.

Definition at line 14 of file ElectronEfficiencyHelpers.cxx.

{
 
  // Define an AsgMessaging instance
  static const asg::AsgMessaging msg("ElectronEfficiencyHelpers");
 
  static const SG::ConstAccessor<int> firstEgMotherPdgIdAcc("firstEgMotherPdgId");
  if ( !(firstEgMotherPdgIdAcc.isAvailable(ele)) ) {
    msg.msg(MSG::ERROR) << "Link not available for firstEgMotherPdgId...BAD!!!" << endmsg;
    msg.msg(MSG::ERROR) << "Need to have present: ( !(ele.isAvailable<int>('firstEgMotherPdgId')) )" << endmsg;
    return CP::CorrectionCode::OutOfValidityRange;
  }
 
  truthcharge = (-1)*firstEgMotherPdgIdAcc(ele);
  // Make truthcharge -1, 0, +1
  truthcharge = (0 < truthcharge) - (truthcharge < 0);
 
  return CP::CorrectionCode::Ok;
 
}

isGoodEle()

CP::CorrectionCode ElectronEfficiencyHelpers::isGoodEle	(	const xAOD::Electron &	ele,
		bool &	goodEle )

Return true if it's good ele for charge flip measurements.

Definition at line 36 of file ElectronEfficiencyHelpers.cxx.

{
 
  // good ele => (firstEgMotherPdgId) == 11 
  // this is the case for type 2,3 and 4 electrons, so is not a real truth match!!!
  // but we believe that it doesn't hurt to apply scale factors for all these cases
  goodEle = false;
  int firstEgPdgId = -9999;
 
  // Define an AsgMessaging instance
  static const asg::AsgMessaging msg("ElectronEfficiencyHelpers");
 
  static const SG::ConstAccessor<int> firstEgMotherPdgIdAcc("firstEgMotherPdgId");
  if ( !(firstEgMotherPdgIdAcc.isAvailable(ele)) ) {
    msg.msg(MSG::ERROR) << "firstEgMotherPdgId IS NOT AVAILABLE!!" << endmsg;
    return CP::CorrectionCode::OutOfValidityRange;
  }
  else {
 
    firstEgPdgId = firstEgMotherPdgIdAcc(ele);
 
    if ( std::abs(firstEgPdgId) != 11) {
 
      goodEle = false;
      return CP::CorrectionCode::Ok;
 
    }
    else {
 
      goodEle = true;
      return CP::CorrectionCode::Ok;
 
    }
  }
 
  return CP::CorrectionCode::OutOfValidityRange;
}