KVSimEvent.cpp

Go to the documentation of this file.

//Created by KVClassFactory on Wed Jun 30 15:10:22 2010
//Author: bonnet
 
#include "KVSimEvent.h"
 
#include <KVSimNucleus.h>
using namespace std;
 
ClassImp(KVSimEvent)
 
//___________________________
 
 
KVSimEvent::KVSimEvent(Int_t mult, const char* classname): KVEvent(mult, classname)
{
   // Default constructor
}
 
 
 
 
KVSimEvent::~KVSimEvent()
{
   //destructeur
}
 
 
 
 
void KVSimEvent::Print(Option_t* t) const
{
   //Print a list of all particles in the event with some characteristics.
   //Optional argument t can be used to select particles (="ok", "groupname", ...)
 
   cout << "\nKVSimEvent with " << ((KVSimEvent*) this)->GetMult(t) << " particles :" << endl;
   cout << "------------------------------------" << endl;
   fParameters.Print();
   KVSimNucleus* frag = 0;
   const_cast<KVSimEvent*>(this)->ResetGetNextParticle();
   while ((frag = (KVSimNucleus*)((KVSimEvent*) this)->GetNextParticle(t))) {
      frag->Print();
      cout << "   Position: (" << frag->GetPosition().x() << "," << frag->GetPosition().y() << "," << frag->GetPosition().z() << ")" << endl;
   }
 
}
 
 
 
 
Double_t KVSimEvent::GetTotalCoulombEnergy() const
{
   // Calculate total Coulomb potential energy of all particles in event
   // Units are MeV.
 
   Double_t Vtot = 0;
   for (KVEvent::Iterator it = begin(); it != end(); ++it) {
      KVEvent::Iterator it2(it);
      for (++it2; it2 != end(); ++it2) {
         TVector3 D12 = it.get_reference<KVSimNucleus>().GetPosition() -
                        it2.get_reference<KVSimNucleus>().GetPosition();
         Vtot += (*it).GetZ() * (*it2).GetZ() * KVNucleus::e2 / D12.Mag();
      }
   }
   return Vtot;
 
}