---

RMCSecondary.cpp

#include "RMCSecondary.hh"
#include "RMCPMT.hh"
#include "ReactorDetector.hh"
#include "TObject.h"

ClassImp(RMCSecondary);

RMCSecondary::RMCSecondary() : TObject()
{
  //Default constructor.
   fPMTs = new TObjArray(2000);  //assume there are not more than 2000
                                 //PMT hits for each secondary
   fNPMT = 0;

}
RMCSecondary::RMCSecondary(ReactorDetector& Detector, int s) : TObject()
{
        fPMTs = new TObjArray(2000); //assume there are not more than 2000
                                     //PMT hits for each secondary
        fSecondaryType = s;
        std::string parent;
        if (fSecondaryType==0) {                      //positron
                fNumPMTs = Detector.GetHitPMTsPos();
                parent = "P";
        }
        if (fSecondaryType==1) {                     //neutron 
                fNumPMTs = Detector.GetHitPMTsNeu();
                parent = "N";
        }
        fNPMT = 0;
        int npmt=Detector.GetNpmt();

// Now loop over all PMT hits, adding a new RCPMT object each time
// and filling it with information from the PMT class.
        for(int pmtindex=0;pmtindex<npmt;pmtindex++) 
        {
                int npe=0;
                int hits=0;
                double Q=0.;
                double Time=0.;
                double PH=0;
                double XYZpmt[3];
                for (int i=0;i<3;i++){
                        XYZpmt[i]=0;
                }
                Detector.GetPMTdata(parent,pmtindex,XYZpmt,&hits,
                                &PH,&Q,&Time,&npe);
                if (hits > 0){
                  RMCPMT *pmt = new RMCPMT();
                  fPMTs->AddAt(pmt,pmtindex);
                  pmt->SetPMTid(pmtindex);
                  pmt->SetPMTQ(Q);
                  pmt->SetPMTtime(Time);
                  pmt->SetPMTPH(PH);
                  pmt->SetPMTnpe(npe);
                  pmt->SetPMThits(hits);
                  pmt->SetPMTxyz(XYZpmt);
                }
        }
}       



RMCSecondary::~RMCSecondary(){
  //RMCSecondaries destructor.

}

---