TpcSimEventGenerator.cxx

Go to the documentation of this file.

#include "TClonesArray.h"
#include "TRandom3.h"

#include "TpcSimMCEvent.h"
#include "TpcSimEventGenerator.h"

ClassImp(TpcSimEventGenerator)



TpcSimEventGenerator::TpcSimEventGenerator()
{
}


TpcSimEventGenerator::~TpcSimEventGenerator()
{
}


TpcSimMCEvent* TpcSimEventGenerator::GenerateEvent(Int_t /* seed */)
{

  TpcSimMCEvent* event = new TpcSimMCEvent();

  Int_t ntracks = 1;
  Double_t x0, y0, z0, t0;
  gRandom->Sphere(x0,y0,z0,50);
  z0 -= 100;
  t0 = 2;
  for(Int_t iTr = 0; iTr<ntracks; iTr++){
    if(event->GetNtracks()>=fkMaxNtracks) break;
    Double_t px0, py0, pz0;
    Int_t pid = 11;
    // x0 = gRandom->Gaus(0,1);
    // y0 = gRandom->Gaus(0,1);
    // z0 = 10 + gRandom->Gaus(0,1);
    px0 = 0.1*gRandom->Gaus(0,1);
    py0 = 0.1*gRandom->Gaus(0,1);
    pz0 = 1;
    if(gHarpoDebug>0) 
      Info("GenerateEvent","%f %f %f %f | %f %f %f | %d",
           x0,y0,z0,t0,
           px0,py0,pz0,
           pid);
    TpcSimMCTrack* tr = new TpcSimMCTrack(x0,y0,z0,t0,
                                          px0,py0,pz0,
                                          pid); // electron
    
    event->AddTrack(tr);
  }
  
  return event;
}