base/TpcSimIonisationTrack.cxx

Go to the documentation of this file.

// 
// File TpcSimIonisationTrack.cxx
//
/*

 */

#include "HarpoDebug.h"
#include "TClonesArray.h"

#include "TpcSimIonisationTrack.h"


ClassImp(TpcSimIonisationTrack)

TpcSimIonisationPoint::TpcSimIonisationPoint()
{
  fNe = 0;
  fX = 0;
  fY = 0;
  fZ = 0;
  fT = 0;
  fEloss = 0;
  fEleft = 0;
}

TpcSimIonisationPoint::TpcSimIonisationPoint(const TpcSimIonisationPoint &pt) :
  TObject(pt), fNe(pt.fNe), fX(pt.fX), fY(pt.fY), fZ(pt.fZ), fT(pt.fT), fEloss(pt.fEloss), fEleft(pt.fEleft)
{
}


TpcSimIonisationPoint::TpcSimIonisationPoint(Int_t Ne, Double_t X, Double_t Y, Double_t Z, Double_t T, Double_t Eloss, Double_t Eleft)
{
  fNe = Ne;
  fX = X;
  fY = Y;
  fZ = Z;
  fT = T;
  fEloss = Eloss;
  fEleft = Eleft;
}




TpcSimIonisationTrack::TpcSimIonisationTrack()
{
  fPoints = new TClonesArray("TpcSimIonisationPoint",fkMaxNpoints);
  fNpoints = 0;
  fMCTrackId = 0;
  fEventId = 0;
  fRunId = 0;
  fParticleId = 0;
  fPx = 0;
  fPy = 0;
  fPz = 0;
  fXorig = 0;
  fYorig = 0;
  fZorig = 0;
  fTorig = 0;
  //fChamber = 0;//new TpcSimChamber();

}

TpcSimIonisationTrack::TpcSimIonisationTrack(const TpcSimIonisationTrack &tr) :
  TObject(tr),
  fNpoints(tr.fNpoints), 
  fMCTrackId(tr.fMCTrackId), 
  fEventId(tr.fEventId), 
  fRunId(tr.fRunId), 
  fParticleId(tr.fParticleId), 
  fPx(tr.fPx), 
  fPy(tr.fPy), 
  fPz(tr.fPz), 
  fXorig(tr.fXorig), 
  fYorig(tr.fYorig), 
  fZorig(tr.fZorig), 
  fTorig(tr.fTorig), 
  //fChamber(tr.fChamber),
  fPoints(tr.fPoints)
{
}

TpcSimIonisationTrack::TpcSimIonisationTrack(Double_t Xorig, Double_t Yorig, Double_t Zorig, Double_t Torig, 
                             Double_t Px, Double_t Py, Double_t Pz,
                             Int_t pID, 
                             Int_t rId, Int_t eId, Int_t tId)
{


  fNpoints = 0;
  fMCTrackId = tId;
  fEventId = eId;
  fRunId = rId;
  fParticleId = pID;
  fPx = Px;
  fPy = Py;
  fPz = Pz;
  fXorig = Xorig;
  fYorig = Yorig;
  fZorig = Zorig;
  fTorig = Torig;
  //fChamber = 0;//new TpcSimChamber();
  fPoints = new TClonesArray("TpcSimIonisationPoint",fkMaxNpoints);

}

TpcSimIonisationTrack::~TpcSimIonisationTrack()
{
  if(gHarpoDebug>1) 
    Info("~TpcSimIonisationTrack","Deleting");
  fPoints->Delete();
  delete fPoints;
}

void TpcSimIonisationTrack::AddPoint(TpcSimIonisationPoint* point)
{
  if(fNpoints>=fkMaxNpoints){
    Warning("AddPoint","Too many points in tracks (%d)",fNpoints);
    return;
  }
  new((*fPoints)[fNpoints]) TpcSimIonisationPoint(*point);
  //fPoints->AddAt(point,fNpoints);
  //fPoints[fNpoints] = point;
  fNpoints++;

}

// void TpcSimIonisationTrack::GenerateTrack(Int_t seed)
// {

//   if(!fChamber){
//     Warning("GenerateTrack","fChamber not defined");
//     return;
//   }
//   Double_t x = fXorig, y = fYorig, z = fZorig;
//   while(fChamber->GetGeom()->GetVolume(x,y,z)==1 && fNpoints<fkMaxNpoints){
//     //TpcSimIonisation* point = 
//     new((*fPoints)[fNpoints]) TpcSimIonisationPoint(1,x,y,z,0,0);
//     fNpoints++;
//     x += fPx;
//     y += fPy;
//     z += fPz;
//   }

// }