BDSParticle.cc

/* 
Beam Delivery Simulation (BDSIM) Copyright (C) Royal Holloway, 
University of London 2001 - 2018.

This file is part of BDSIM.

BDSIM is free software: you can redistribute it and/or modify 
it under the terms of the GNU General Public License as published 
by the Free Software Foundation version 3 of the License.

BDSIM is distributed in the hope that it will be useful, but 
WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with BDSIM.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "BDSParticle.hh"

#include <ostream>

BDSParticle::BDSParticle():
  position(0,0,0),
  momentum(0,0,0),
  totalEnergy(0.0),
  T(0.0),
  weight(1.0),
  trackID(-1),
  parentID(-1)
{;}

BDSParticle::BDSParticle(G4double x,
                         G4double y,
                         G4double z,
                         G4double xp,
                         G4double yp,
                         G4double zp,
                         G4double totalEnergyIn,
                         G4double globalTimeIn,
                         G4double weightIn,
                         G4int trackIDIn,
                         G4int parentIDIn):
  position(x,y,z),
  momentum(xp,yp,zp),
  totalEnergy(totalEnergyIn),
  T(globalTimeIn),
  weight(weightIn),
  trackID(trackIDIn),
  parentID(parentIDIn)
{;}

BDSParticle::BDSParticle(G4ThreeVector pos,
                         G4ThreeVector mom,
                         G4double      totalEnergyIn,
                         G4double      globalTimeIn,
                         G4double      weightIn,
                         G4int         trackIDIn,
                         G4int         parentIDIn):
  position(pos),
  momentum(mom),
  totalEnergy(totalEnergyIn),
  T(globalTimeIn),
  weight(weightIn),
  trackID(trackIDIn),
  parentID(parentIDIn)
{;}

std::ostream& operator<< (std::ostream& out, BDSParticle const& p)
{
  out << "Total E: " << p.GetTotalEnergy() << " MeV" << G4endl;
  out << "Position: (" << p.GetX() << ", " << p.GetY() << ", " << p.GetZ() << ")" << G4endl;
  out << "Momentum: (" << p.GetXp() << ", " << p.GetYp() << ", " << p.GetZp() << ")" << G4endl;
  return out;
}