KVSimEvent.cpp

//Created by KVClassFactory on Wed Jun 30 15:10:22 2010
//Author: bonnet
 
#include "KVSimEvent.h"
 
ClassImp(KVSimEvent)
 
 
 
 
void KVSimEvent::Print(Option_t* t) const
{
   //Print a list of all particles in the event with some characteristics.
 
   std::cout << "SimEvent #" << GetNumber() << "  MULT=" << GetMult() << std::endl;
   fParameters.ls();
   std::cout << std::endl;
   int id = 1;
   for (auto& p : *this) {
      std::cout << id << " : ";
      p.Print();
      std::cout << std::endl << std::endl;
      ++id;
   }
}
 
 
 
 
Double_t KVSimEvent::GetTotalCoulombEnergy() const
{
   // Calculate total Coulomb potential energy of all particles in event
   // Units are MeV.
 
   Double_t Vtot = 0;
   for (KVSimEvent::Iterator it = begin(); it != end(); ++it) {
      KVSimEvent::Iterator it2(it);
      for (++it2; it2 != end(); ++it2) {
         TVector3 D12 = (*it).GetPosition() - (*it2).GetPosition();
         Vtot += (*it).GetZ() * (*it2).GetZ() * KVNucleus::e2 / D12.Mag();
      }
   }
   return Vtot;
 
}
 
 
 
void KVSimEvent::Streamer(TBuffer& R__b)
{
   //Stream an object of class KVReconstructedEvent.
   //We set the particles' angles depending on whether mean or random angles
   //are wanted (fMeanAngles = kTRUE or kFALSE)
 
   if (R__b.IsReading()) {
      UInt_t R__s, R__c;
      UInt_t R__v  = R__b.ReadVersion(&R__s, &R__c);
      if (R__v < 5) {
         // all versions written with automatic streamer
         // (version is not read correctly from TTree branches)
         Clear();
         KVBase::Streamer(R__b);
         fParticles->Streamer(R__b);
         fParameters.Streamer(R__b);
      }
      else
         R__b.ReadClassBuffer(KVSimEvent::Class(), this, R__v, R__s, R__c);
   }
   else {
      R__b.WriteClassBuffer(KVSimEvent::Class(), this);
   }
}