BranchProbability.cpp

Go to the documentation of this file.

// BranchProbability.cpp: implementation of the BranchProbability class.
//
 
#include "BranchProbability.h"
#include "../Scheduler/CommandControlCenter.h"
#include "../Scheduler/AnalysisInfoContainer.h"
#include "../Crew/OperatorDM.h"
#include <boost/serialization/export.hpp>
#include "../Global/SimulationContext.h"
#include <boost/archive/binary_oarchive.hpp>
#include <boost/archive/binary_iarchive.hpp>
BOOST_CLASS_EXPORT_IMPLEMENT(BranchProbability)
using namespace std;
 
// Construction/Destruction
 
BranchProbability::BranchProbability() {
  ID_of_this_branch = -999;
  ID_of_the_associated_branch_cluster = -999;
  ID_of_the_child_branch_cluster = -999;
  unconditional_branch_probability = -1.0;
  type_of_child = sim::None;//Branch Or ES
}
 
BranchProbability::~BranchProbability() {
  CommunicationPlace::clearCommunicationPlaceAll();
  child_IDs.clear();
}
void BranchProbability::calculateProbability() {
  double b_probability = 1.0;
  for (auto &e : goal) { b_probability *= e->getProbability(); }
  for (auto &e : strategy) { b_probability *= e->getProbability(); }
  for (auto &e : subgoal) { b_probability *= e->getProbability(); }
  for (auto &e : diagnosis) { b_probability *= e->getProbability(); }
  for (auto &e : action_objective) { b_probability *= e->getProbability(); }
  for (auto &e : action_pack) { b_probability *= e->getProbability(); }
  for (auto &e : action) { b_probability *= e->getProbability(); }
  for (auto &e : alarm_state) { b_probability *= e->getProbability(); }
  for (auto &e : component_state) { b_probability *= e->getProbability(); }
  for (auto &e : support_system_state) { b_probability *= e->getProbability(); }
  for (auto &e : parameter_trend) { b_probability *= e->getProbability(); }
  for (auto &e : parameter_value) { b_probability *= e->getProbability(); }
  for (auto &e : parameter_difference) { b_probability *= e->getProbability(); }
  for (auto &e : procedure_step) { b_probability *= e->getProbability(); }
  for (auto &e : instinctive_response_procedure) { b_probability *= e->getProbability(); }
  for (auto &e : info_gather_mode) { b_probability *= e->getProbability(); }
  for (auto &e : HW_reliability) { b_probability *= e->getProbability(); }
  for (auto &e : SS_reliability) { b_probability *= e->getProbability(); }
  for (auto &e : mental_belief) { b_probability *= e->getProbability(); }
  for (auto &e : scanned_parameter) { b_probability *= e->getProbability(); }
  Probability::setProbability(b_probability);
}
 
void BranchProbability::setIDOfThisBranch(int n) { ID_of_this_branch = n; }
void BranchProbability::setIDOfTheAssociatedBranchCluster(int n) { ID_of_the_associated_branch_cluster = n; }
//void BranchProbability::setIThBranchInTheAssociatedBranchCluster(int n){i_th_branch_in_the_associated_branch_cluster = n;}
void BranchProbability::setUnconditionalBranchProbability(double x) { unconditional_branch_probability = x; }
void BranchProbability::calUnconditionalBranchProbability(double x) {
  unconditional_branch_probability = x * Probability::getProbability();
}
int BranchProbability::getIDOfThisBranch() { return ID_of_this_branch; }
int BranchProbability::getIDOfTheAssociatedBranchCluster() { return ID_of_the_associated_branch_cluster; }
//int BranchProbability::getIThBranchInTheAssociatedBranchCluster(){return i_th_branch_in_the_associated_branch_cluster;}
double BranchProbability::getUnconditionalBranchProbability() { return unconditional_branch_probability; }
void BranchProbability::sendBranchInfoToCommandCenter(int _associate_branch_ID) {
  // sets branch timing and adds branch to Command_handler
  double ODM_action_time_multiplier, OAT_action_time_multiplier, OCT_action_time_multiplier, action_time_multiplier;
  double time;
  int sender, receiver;
  string str;
  //iterators: range-based for loops used below
//
 
// TODO: disabled for now
//  if (Decision_Maker.m_RM
//      != NULL) { Decision_Maker.m_RM->m_event_schema_base->UpdateTheBranchingHistory((void *) this); }
//
//  if (Action_Taker.m_RM != NULL) { Action_Taker.m_RM->m_event_schema_base->UpdateTheBranchingHistory((void *) this); }
  ODM_action_time_multiplier = g_sim.SScheduler.getODMActionTimeMultiplier();
  OAT_action_time_multiplier = g_sim.SScheduler.getOATActionTimeMultiplier();
  OCT_action_time_multiplier = g_sim.SScheduler.getOCTActionTimeMultiplier();
  for (auto &e : goal) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addGoal(e, _associate_branch_ID);
  }
  for (auto &e : strategy) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addStrategy(e, _associate_branch_ID);
  }
  for (auto &e : subgoal) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addSubGoal(e, _associate_branch_ID);
  }
  for (auto &e : diagnosis) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addDiagnosis(e, _associate_branch_ID);
  }
  for (auto &e : action_objective) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addActionObjective(e, _associate_branch_ID);
  }
  for (auto &e : action_pack) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addActionPack(e, _associate_branch_ID);
  }
 
  for (auto &e : scanned_parameter) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addScannedParameter(e, _associate_branch_ID);
  }
 
  for (auto &e : action) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    //
    if (receiver == sim::Source::Plant) {
      time = e->getActionTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    } else if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else {//
      //
      // sender is the receiver - used for "Waiting" action steps to simulate time delay
      e->setBothTotalTimeRequiredAndRemainingTime(
          (e->getTotalTimeRequired()) * action_time_multiplier);
    }
    //
    g_sim.Command_handler.addAction(e, _associate_branch_ID);
  }
  for (auto &e : alarm_state) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == receiver) {
      e->setBothTotalTimeRequiredAndRemainingTime(0.0);
    } else if (sender == sim::Operator::ODM && receiver == sim::Source::Plant)
    {}// the timer is already set when the branch is created
    else {
      time = e->calculateMeanRequiredTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    }
    g_sim.Command_handler.addExpectedAlarmState(e, _associate_branch_ID);
  }
  for (auto &e : component_state) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == receiver) {
      e->setBothTotalTimeRequiredAndRemainingTime(0.0);
    } else if (sender == sim::Operator::ODM && receiver == sim::Source::Plant)
    {}// the timer is already set when the branch is created
    else {
      time = e->calculateMeanRequiredTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    }
    g_sim.Command_handler.addExpectedComponentState(e, _associate_branch_ID);
  }
  for (auto &e : support_system_state) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == receiver) {
      e->setBothTotalTimeRequiredAndRemainingTime(0.0);
    } else if (sender == sim::Operator::ODM && receiver == sim::Source::Plant)
    {}// the timer is already set when the branch is created
    else {
      time = e->calculateMeanRequiredTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    }
    g_sim.Command_handler.addExpectedSupportSystemState(e, _associate_branch_ID);
  }
  for (auto &e : parameter_trend) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == receiver) {
      e->setBothTotalTimeRequiredAndRemainingTime(0.0);
    } else {
      time = e->calculateMeanRequiredTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    }
    g_sim.Command_handler.addExpectedParameterTrend(e, _associate_branch_ID);
  }
  for (auto &e : parameter_value) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == receiver) {
      e->setBothTotalTimeRequiredAndRemainingTime(0.0);
    } else if (sender == sim::Operator::ODM && receiver == sim::Source::Plant)
    {}// the timer is already set when the branch is created
    else {
      time = e->calculateMeanRequiredTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    }
    g_sim.Command_handler.addExpectedParameterValue(e, _associate_branch_ID);
  }
  for (auto &e : parameter_difference) {
    sender = e->getSender();
    receiver = e->getIntendedReceiver();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    if (sender == sim::Operator::ODM && receiver == sim::Operator::OAT) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Decision_Maker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == sim::Operator::OAT && receiver == sim::Operator::ODM) {
      e->setBothTotalTimeRequiredAndRemainingTime(
          g_sim.Action_Taker.getNominalCommunicationTime() * action_time_multiplier);
    } else if (sender == receiver) {
      e->setBothTotalTimeRequiredAndRemainingTime(0.0);
    } else {
      time = e->calculateMeanRequiredTime();
      e->setBothTotalTimeRequiredAndRemainingTime(time * action_time_multiplier);
    }
    g_sim.Command_handler.addExpectedParameterDifference(e, _associate_branch_ID);
  }
  for (auto &e : procedure_step) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addProcedureStep(e, _associate_branch_ID);
  }
  for (auto &e : instinctive_response_procedure) {
    sender = e->getSender();
    if (sender == sim::Operator::ODM) { action_time_multiplier = ODM_action_time_multiplier; }
    else if (sender == sim::Operator::OAT) { action_time_multiplier = OAT_action_time_multiplier; }
    else if (sender == sim::Operator::OCT) { action_time_multiplier = OCT_action_time_multiplier; }
    else { action_time_multiplier = 1.0; }
    e->setBothTotalTimeRequiredAndRemainingTime(e->getRemainingTime() * action_time_multiplier);
    g_sim.Command_handler.addInstinctiveResponseProcedure(e, _associate_branch_ID);
  }
  for (auto &e : info_gather_mode) {
    e->setBothTotalTimeRequiredAndRemainingTime(-1.0);
    g_sim.Command_handler.addInfoGatherMode(e, _associate_branch_ID);
  }
  for (auto &e : mental_belief) {
    g_sim.Command_handler.addMentalBelief(e, _associate_branch_ID);
  }
  for (auto &e : HW_reliability) {
    g_sim.Command_handler.addHWReliability(e, _associate_branch_ID);
  }
  for (auto &e : SS_reliability) {
    g_sim.Command_handler.addSSReliability(e, _associate_branch_ID);
  }
  for (auto &e : new_alarm) {
    g_sim.Command_handler.addNewAlarm(e, _associate_branch_ID);
  }
}
void BranchProbability::setTypeOfChild(int n) { type_of_child = n; }
int BranchProbability::getTypeOfChild() { return type_of_child; }
void BranchProbability::saveInfoToFile(ostream &fout) {
  CommunicationPlace::saveInfoToFileAsBranchForADSV1(fout);
}
void BranchProbability::setIDOfTheChildBranchCluster(int n) { ID_of_the_child_branch_cluster = n; }
int BranchProbability::getIDOfTheChildBranchCluster() { return ID_of_the_child_branch_cluster; }
void BranchProbability::saveBranchProbabilityToFile(ofstream &fout) {
  fout << ID_of_this_branch << ' ' << ID_of_the_associated_branch_cluster << ' ' << ID_of_the_child_branch_cluster
       << ' '
       << unconditional_branch_probability << ' ' << type_of_child << '\n';
  this->saveProbabilityToFile(fout);
  this->saveCommunicationPlaceToFile(fout);
  this->saveStateOfActivationToFile(fout);
}
void BranchProbability::retrieveBranchProbabilityFromFile(ifstream &fin) {
  fin >> ID_of_this_branch >> ID_of_the_associated_branch_cluster >> ID_of_the_child_branch_cluster
      >> unconditional_branch_probability >> type_of_child;
  this->retrieveProbabilityFromFile(fin);
  this->retrieveCommunicationPlaceFromFile(fin);
  this->retrieveStateOfActivationFromFile(fin);
}
void BranchProbability::copyBranchProbability(BranchProbability &_bp) {
  ID_of_this_branch = _bp.getIDOfThisBranch();
  ID_of_the_associated_branch_cluster = _bp.getIDOfTheAssociatedBranchCluster();
  ID_of_the_child_branch_cluster = _bp.getIDOfTheChildBranchCluster();
  unconditional_branch_probability = _bp.getUnconditionalBranchProbability();
  type_of_child = _bp.getTypeOfChild();//Branch Or ES
  this->copyProbability(_bp);
  this->copyStateOfActivation(_bp);
  this->copyCommunicationPlace(_bp);
 
}