ggrc

//*****************************************************************************
int GetGenieReactionCode(const std::string& reactionCode){
//*****************************************************************************

  int code = 0;

  //if(!p_silent) cout << "3: " << reactionCode << endl;
  if(reactionCode.find("charm") == string::npos &&
    reactionCode.find("Weak[CC],QES") != string::npos){
    code = 1;
  }
  else if(reactionCode.find("Weak[CC],RES") != string::npos){
    code = 11;
  }
  else if(reactionCode.find("charm") == string::npos &&
    reactionCode.find("Weak[CC],DIS") != string::npos){
    code = 21;
  }
  else if(reactionCode.find("Weak[CC],COH") != string::npos){
    code = 16;
  }
  else if(reactionCode.find("Weak[NC],RES") != string::npos){
    code = 31;
  }
  else if(reactionCode.find("Weak[NC],QES") != string::npos){
    code = 51;
  }
  else if(reactionCode.find("Weak[NC],COH") != string::npos){
    code = 36;
  }
  else if(reactionCode.find("Weak[NC],DIS") != string::npos){
    code = 41;
  }
  return code;
}