#include "Parameterisation.hh"


#include "G4Box.hh"
#include "G4Sphere.hh"
#include "G4ThreeVector.hh"
#include "G4SystemOfUnits.hh"
#include "G4VPhysicalVolume.hh"

Parameterisation::Parameterisation(const G4ThreeVector& voxelsize_sphere,G4int NZ_SPHERE):
G4VPVParameterisation(),dr(voxelsize_sphere.x()),dphi(voxelsize_sphere.y()),dtheta(voxelsize_sphere.z()),nz_sphere(NZ_SPHERE)
{
pz_sphere.clear();

G4double z_sphere;

for(G4int i=0;i<nz_sphere;i++)
{
  z_sphere=i*dr;
  pz_sphere.push_back(z_sphere);

}

}

Parameterisation::~Parameterisation(){pz_sphere.clear();}

void Parameterisation::ComputeTransformation(const G4int copyNo,G4VPhysicalVolume* pv) const
{
  G4ThreeVector origin(0,0,0);
  pv->SetTranslation(origin);
  pv->SetRotation(0);
}
void Parameterisation::ComputeDimensions (G4Sphere& sphere,const G4int copyNo, const G4VPhysicalVolume*) const
{
sphere.SetInnerRadius(pz_sphere[copyNo]);
sphere.SetOuterRadius(pz_sphere[copyNo]+dr);
sphere.SetStartPhiAngle(0*deg);
sphere.SetDeltaPhiAngle(dphi);
sphere.SetStartThetaAngle(0*deg);
sphere.SetDeltaThetaAngle(dtheta);
}