---

RMCEvent.cpp

#include <iostream.h>
#include <math.h>
#include <string.h>
#include <cstdio>
#include "TObject.h"
#include "TObjArray.h"
#include "RMCEvent.hh"
#include "RMCSecondary.hh"
#include "ReactorEvent.hh"
#include "ReactorDetector.hh"

ClassImp(RMCEvent);

RMCEvent::RMCEvent()
{
  //Default constructor.
  fSecondaries = new TObjArray(10);
  fNSecondary = 0;

}

RMCEvent::~RMCEvent(){
  //RMCEvent destructor.

}

void RMCEvent::Build(ReactorEvent& Event, ReactorDetector& Detector, std::string EventType, int Nevent)
{
//Build event.  Information comes from ReactorEvent and ReactorDetector
// classes.

    fEventNum = Nevent;
    if (EventType == "reactor"){   
            fEventType = 0;
            int n=0;

            fKEParent = Event.GetEnu(n);

            double* rnu = new double[n];
            double* cnu = new double[n];
            double* fnu = new double[n];
            rnu[n] = Event.GetNubarVertexR(n);
            cnu[n] = Event.GetNubarVertexCosTheta(n);
            fnu[n] = Event.GetNubarVertexPhi(n);
            fInteractionVertex[0]  = rnu[n]*cos(fnu[n])*sqrt(1-cnu[n]*cnu[n]);
            fInteractionVertex[1]  = rnu[n]*sin(fnu[n])*sqrt(1-cnu[n]*cnu[n]);
            fInteractionVertex[2]  = rnu[n]*cnu[n];

            fNumSecondaries = 2;

       // clean up
            delete[] rnu;
            delete[] cnu;
            delete[] fnu;

      // Now fill the Secondary information (positron and neutron, 
      // using a TObjArray
      // Add one for each secondary.
         fNSecondary = 0;
         double energy=0.;
         //      double direction[3];
         double time=0.;
         double dircos[3];
         int ngammas = 0;
         for (int s=0; s<fNumSecondaries; s++) 
         {
                 double ETotalDeposited=0.;
                 RMCSecondary *secondary = new RMCSecondary(Detector,s);
                 if (s==0) 
                 {
                         energy = Event.GetPositronKE(n);
                         time = Event.GetPositronTime(n);
                         dircos[0]=1.0;  //These are made up directions for 
                         dircos[1]=0.0;  //now.
                         dircos[2]=0.0;
                         // Get number of `gammas' made, and total
                         // deposited energy
                         ngammas = Detector.GetNgamma("P");
                         double* eg = new double[ngammas];
                         double* tg = new double[ngammas];
                         double* rg = new double[ngammas];
                         double* cg = new double[ngammas];
                         double* fg = new double[ngammas];
                         Detector.GetGamma("P",eg,tg,rg,cg,fg);
                         for (int i=0;i < ngammas; i++) ETotalDeposited+=eg[i];
                         delete[] eg;
                         delete[] tg;
                         delete[] rg;
                         delete[] cg;
                         delete[] fg;

                 }
                 if (s==1) 
                 {
                         energy = Event.GetNeutronKE(n);
                         time = Event.GetNeutronTime(n);
                         dircos[0]=-1.0;  //These are made up directions for 
                         dircos[1]=0.0;  //now.
                         dircos[2]=0.0;
                         // Get number of gammas made, and total
                         // deposited energy
                         ngammas = Detector.GetNgamma("N");
                         double* eg = new double[ngammas];
                         double* tg = new double[ngammas];
                         double* rg = new double[ngammas];
                         double* cg = new double[ngammas];
                         double* fg = new double[ngammas];
                         Detector.GetGamma("N",eg,tg,rg,cg,fg);
                         for (int i=0;i < ngammas; i++) ETotalDeposited+=eg[i];
                         delete[] eg;
                         delete[] tg;
                         delete[] rg;
                         delete[] cg;
                         delete[] fg;
                 }
                 secondary->SetSecondaryGammas(ngammas);
                 secondary->SetSecondaryEDeposited(ETotalDeposited);
                 secondary->SetSecondaryKE(energy);
                 secondary->SetSecondaryTime(time);
                 secondary->SetSecondaryDirection(dircos);
                 fSecondaries->AddAt(secondary,s);
         }

    }
}

---