interaction.cxx

Go to the documentation of this file.

#include "interaction.h"
#include <iostream>

ClassImp(interaction)

using std::endl;
using std::cout;

//____________________________________________________________________________
ostream & operator << (ostream & stream, const interaction & conf)
{
  conf.print(stream);
  return stream;
}
//____________________________________________________________________________
interaction::interaction() :
_name("default")
{

}
//____________________________________________________________________________
interaction::interaction(const char * name) :
_name(name)
{

}
//____________________________________________________________________________
interaction::interaction(const interaction * inter) 
{
  _name       = inter->_name;
  _flavor     = inter->_flavor;
  _nucleus    = inter->_nucleus;
  _ccnc       = inter->_ccnc;
  _init_state = inter->_init_state;
}
//____________________________________________________________________________
interaction::interaction(flavor_t f, nucleus_t n, ccnc_t c, init_state_t in) :
_flavor(f),_nucleus(n),_ccnc(c),_init_state(in)
{

}
//____________________________________________________________________________
interaction::~interaction()
{

}
//____________________________________________________________________________
int interaction::getProcess(void) const
{
  int icode = 0;
  int pa    = 0;
  int im    = _init_state;
  int it    = _flavor;
  int cn    = _ccnc;

  makestate_(&pa, &pa, &pa, &pa, &pa, &it, &im, &cn, &icode);

  //writestate_(&icode);  

  return icode;
}
//____________________________________________________________________________
int interaction::getProcess(final_state * final) const
{
  int proton =  final->getProton();
  int neutron = final->getNeutron();
  int piplus =  final->getPiplus();
  int piminus = final->getPiminus();
  int pizero =  final->getPizero();
  
  int icode = 0;
  int im    = _init_state;
  int it    = _flavor;
  int cn    = _ccnc;

  makestate_(&piplus, &piminus, &pizero, &proton, &neutron,
             &it, &im, &cn, &icode);

  cout << "In interaction, process (w/final)= " << icode << endl;

  return icode;
}
//____________________________________________________________________________
void interaction::print(ostream & stream) const
{
  stream << "Interaction Name:........." << _name  << endl;
  stream << "Flavor:..................." << flavor::as_string(_flavor) << endl;
  stream << "Nucleus:.................." << nucleus::as_string(_nucleus) << endl;
  stream << "CCNC:....................." << ccnc::as_string(_ccnc) << endl;
  stream << "Init_State:..............." << init_state::as_string(_init_state) << endl;  
}
//____________________________________________________________________________
        

---