G4RootHitContainer.cc

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4RootHitContainer.cc

#include "G4RootHitContainer.h"

#include <g4main/PHG4Hit.h>      // for PHG4Hit
#include <g4main/PHG4HitEval.h>

#include <TClonesArray.h>

#include <cmath>                // for NAN
#include <ostream>               // for basic_ostream::operator<<, operator<<

using namespace std;

static const int NMAX = 100000;

G4RootHitContainer::G4RootHitContainer()
  : etotal(NAN)
  , eion(NAN)
  , leakage(NAN)
  , event(0)
{
  SnglHits = new TClonesArray("PHG4HitEval", NMAX);
}

G4RootHitContainer::~G4RootHitContainer()
{
  SnglHits->Clear();
  delete SnglHits;
}

void G4RootHitContainer::Reset()
{
  etotal = NAN;
  leakage = NAN;
  event = 0;
  SnglHits->Clear();
  if (SnglHits->GetSize() > NMAX)
  {
    SnglHits->Expand(NMAX);
  }
  return;
}

PHG4Hit *
G4RootHitContainer::AddHit(const PHG4Hit *g4hit)
{
  TClonesArray &cl = *SnglHits;
  int nextindex = SnglHits->GetLast() + 1;
  if (nextindex == SnglHits->GetSize())
  {
    SnglHits->Expand(SnglHits->GetSize() + 10000);
  }
  new (cl[nextindex]) PHG4HitEval(g4hit);
  return (static_cast<PHG4Hit *>(cl[nextindex]));
}

void G4RootHitContainer::identify(ostream &os) const
{
  os << "Number of Hits: " << SnglHits->GetLast() << endl;
  return;
}