HarpoSimEvent.cxx

Go to the documentation of this file.

// 
// File HarpoSimEvent.cxx
//
#include "HarpoSimEvent.h"
#include "HarpoDccMap.h"
#include "HarpoDebug.h"
#include "TpcSimIonisationTrack.h"
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <cassert>
// for swab 
#include <unistd.h> 

ClassImp(HarpoSimEvent)

HarpoSimEvent::HarpoSimEvent()
{
  if(gHarpoDebug>1)
    Info("HarpoSimEvent()","Creating Sim Event");
  fIonisationTracks = new TClonesArray("TpcSimIonisationTrack",fkMaxNtracks);
  //  fStrips = NULL;//new TClonesArray("HarpoSimStrips",fkMaxNtracks);

  fNtracks = 0; 
  fNstrips = 0;

  fMCEvent = (TpcSimMCEvent*) NULL;

  if(gHarpoDebug>1)
    Info("HarpoSimEvent()","Sim Event Created");
}

HarpoSimEvent::~HarpoSimEvent()
{
  if(gHarpoDebug>1) Info("~HarpoSimEvent","Deleting");
  fIonisationTracks->Delete();
  delete fIonisationTracks;
  //fStrips->Delete();
  if(fMCEvent) 
    fMCEvent->Delete();

}

Int_t HarpoSimEvent::AddIonisationTrack(Int_t iTr, TpcSimIonisationTrack* tr)
{

  if(fNtracks != iTr){
    return 0;
  }
  if(fNtracks>=fkMaxNtracks){
    return -1;
  }

  fMCEvent->AddTrack(tr);
  if(gHarpoDebug>1)
    Info("AddIonisationTrack","%d",fNtracks);
  new((*fIonisationTracks)[fNtracks]) TpcSimIonisationTrack(*tr);
  fNtracks++;
  //  Info("AddIonisationTrack","%d   %d",fNtracks,fIonisationTracks->GetEntries());


  return fNtracks;

}

TpcSimIonisationTrack* HarpoSimEvent::GetIonisationTrack(Int_t iTr)
{

  if(iTr < 0)
    return 0;
  if(iTr >= fNtracks)
    return 0;

  return (TpcSimIonisationTrack*)fIonisationTracks->At(iTr);

}


Int_t HarpoSimEvent::AddStrips(Int_t /* iTr */ , HarpoSimStrips* /* strips*/ )
{

  // if(fNstrips != iTr){
  //   //std::cout << "fNstrips (" << fNstrips 
  //   return 0;
  // }
  // if(fNstrips>=fkMaxNtracks){
  //   return -1;
  // }

  // new((*fStrips)[fNstrips]) HarpoSimStrips(*strips);
  fNstrips++;
  return fNstrips;

}

void          HarpoSimEvent::SetHeader(UInt_t len, UInt_t evnum)
{
  fEventHeader.eventSize=len;
  fEventHeader.eventNumb=evnum;
}


void HarpoSimEvent::print() const
{
  std::cout << "Header num : " << fEventHeader.eventNumb
            << " : size : " << fEventHeader.eventSize << std::endl;
}

ULong_t      HarpoSimEvent::GetTimeStamp(Int_t )
{
  Info("GetTimeStamp","0");
  return (ULong_t) 0;
}