G4MolecularDissociationTable.cc

Go to the documentation of this file.

//
// ********************************************************************
// * License and Disclaimer                                           *
// *                                                                  *
// * The  Geant4 software  is  copyright of the Copyright Holders  of *
// * the Geant4 Collaboration.  It is provided  under  the terms  and *
// * conditions of the Geant4 Software License,  included in the file *
// * LICENSE and available at  http://cern.ch/geant4/license .  These *
// * include a list of copyright holders.                             *
// *                                                                  *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work  make  any representation or  warranty, express or implied, *
// * regarding  this  software system or assume any liability for its *
// * use.  Please see the license in the file  LICENSE  and URL above *
// * for the full disclaimer and the limitation of liability.         *
// *                                                                  *
// * This  code  implementation is the result of  the  scientific and *
// * technical work of the GEANT4 collaboration.                      *
// * By using,  copying,  modifying or  distributing the software (or *
// * any work based  on the software)  you  agree  to acknowledge its *
// * use  in  resulting  scientific  publications,  and indicate your *
// * acceptance of all terms of the Geant4 Software license.          *
// ********************************************************************
//
//
// WARNING : This class is released as a prototype.
// It might strongly evolve or even disapear in the next releases.
//
// ----------------------------------------------------------------------
//      GEANT 4 class implementation file
//
//      History: first implementation by Alfonso Mantero 4 Mar 2009
//
// ----------------------------------------------------------------
 
#include "G4MolecularDissociationTable.hh"
#include "G4MolecularDissociationChannel.hh"
#include "G4MolecularConfiguration.hh"
 
using namespace std;
using namespace G4DNA;
 
//______________________________________________________________________________
 
G4MolecularDissociationTable::G4MolecularDissociationTable()
= default;
 
//______________________________________________________________________________
 
G4MolecularDissociationTable::~G4MolecularDissociationTable()
{
  auto it_map = fDissociationChannels.begin();
 
  for (; it_map != fDissociationChannels.end(); it_map++)
  {
    vector<const G4MolecularDissociationChannel*>& decayChannels = it_map
        ->second;
    if (!decayChannels.empty())
    {
      for (auto & decayChannel : decayChannels)
      {
        if (decayChannel != nullptr)
        {
          delete decayChannel;
          decayChannel = nullptr;
        }
      }
      decayChannels.clear();
    }
  }
  fDissociationChannels.clear();
}
 
//______________________________________________________________________________
 
G4MolecularDissociationTable::
  G4MolecularDissociationTable(const G4MolecularDissociationTable& right)
{
  *this = right;
}
 
//______________________________________________________________________________
 
G4MolecularDissociationTable&
G4MolecularDissociationTable::operator=(const G4MolecularDissociationTable& right)
{
  if(this == &right) return *this;
  fDissociationChannels = right.fDissociationChannels;
  return *this;
}
 
//______________________________________________________________________________
 
const vector<const G4MolecularDissociationChannel*>*
G4MolecularDissociationTable::
  GetDecayChannels(const G4MolecularConfiguration* conf) const
{
  auto it_exstates = fDissociationChannels.find(conf);
  if (it_exstates == fDissociationChannels.end()) return nullptr;
  return &(it_exstates->second);
}
 
//______________________________________________________________________________
 
const vector<const G4MolecularDissociationChannel*>*
G4MolecularDissociationTable::GetDecayChannels(const G4String& exState) const
{
  for(const auto & fDissociationChannel : fDissociationChannels)
  {
    if(fDissociationChannel.first->GetLabel() == exState) return &(fDissociationChannel.second);
  }
  return nullptr;
}
 
//______________________________________________________________________________
 
//void G4MolecularDissociationTable::
//  AddExcitedState(const G4String& label,
//                  const G4MolecularConfiguration* molConf);
//{
//
//}
 
//______________________________________________________________________________
 
void G4MolecularDissociationTable::
  AddChannel(const G4MolecularConfiguration* molConf,
             const G4MolecularDissociationChannel* channel)
{
  fDissociationChannels[molConf].push_back(channel);
}
 
//______________________________________________________________________________
 
void G4MolecularDissociationTable::CheckDataConsistency() const
{
  ChannelMap::const_iterator channelsIter;
 
  for(channelsIter = fDissociationChannels.begin();
      channelsIter != fDissociationChannels.end(); ++channelsIter)
  {
 
    const vector<const G4MolecularDissociationChannel*>& decayVect =
        channelsIter->second;
    G4double sum = 0;
 
    G4double max = decayVect.size();
 
    for(size_t i = 0; i < max; i++)
    {
      const G4MolecularDissociationChannel* decay = decayVect[i];
      const G4double prob = decay->GetProbability();
      sum += prob;
    }
 
    if(sum != 1)
    {
      G4ExceptionDescription errMsg;
      errMsg << "The probabilities for deecitation of molecular configuration "
             << channelsIter->first->GetName()
             << " with label :" <<  channelsIter->first->GetLabel()
             << " don't sum up to 1";
      G4Exception("G4MolecularDissociationTable::CheckDataConsistency",
                  "BRANCHING_RATIOS_CONSISTENCY",
                  FatalErrorInArgument,
                  errMsg);
    }
  }
}
 
void G4MolecularDissociationTable::Serialize(std::ostream& /*char_traits*/)
{
  // TODO
}