Message: problem with primitive scorer? Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None problem with primitive scorer? 

Keywords: primitive scorer, parameterisation
Forum: Hits, Digitization and Pileup
Date: 23 Mar, 2006
From: sylvia studeny <sylvia.studeny@gsf.de>

Hi, I'm working with geant4.8.0 on a problem using phantoms. For this I took example/extended/runAndEvent/Re02. I print out an file, so that I can see the deposited energy in every voxel. For this I use the data in the primitive scorers. With the original example everything is fine and I get the expected energy distribution.Then I tried and changed the parameterisation into

#include "RE02PhantomParameterisation.hh"
#include "G4Box.hh"
#include "G4VPhysicalVolume.hh"

RE02PhantomParameterisation::RE02PhantomParameterisation(
		   const G4ThreeVector& motherSize, const G4ThreeVector& voxelSize,
		   const G4int nx, const G4int ny, const G4int nz)
//I changed this in DetectorConstruction,too!
    :G4VPVParameterisation(),fDxyzMother(motherSize),fVoxel(voxelSize),
     fNx(nx),fNy(ny),fNz(nz)
{}

RE02PhantomParameterisation::~RE02PhantomParameterisation() {}

void RE02PhantomParameterisation::ComputeTransformation(const G4int copyNo, 
				  G4VPhysicalVolume* physVol) const
{

G4int nx=0,ny=0,nz=0;
G4int nzx=fNz*fNx;
   G4int cN; 
    cN=copyNo;
    //G4cout<<"CopyNumber"<<cN<<G4endl;
    ny=cN/nzx;
    cN%=nzx;
    nx=cN/fNz;
    cN%=fNz;
    nz=cN;

  G4double Xposition= -fDxyzMother.x() + (nx+.5) * fVoxel.x(); 
  G4double Yposition= -fDxyzMother.y() + (ny+.5) * fVoxel.y();
  G4double Zposition= -fDxyzMother.z() + (nz+.5) * fVoxel.z();

  //G4double Xposition= offsetX + nx *2* fVoxel.x();
  //G4double Yposition= offsetY + ny *2* fVoxel.y();
  //G4double Zposition= offsetZ + nz *2* fVoxel.z();
  G4ThreeVector Schwerpunkt(Xposition,Yposition,Zposition);

    physVol->SetTranslation(Schwerpunkt);
    physVol->SetRotation(0);
}

which worked quite well in geant4.7.0.p1. Also when I look at the phantom with vis.mac everything looks fine. The deposited energy on the other hand looks quite weird. I double checked everything, but I don't know what I did wrong. Has anybody any idea what else I can try? Thanks a lot Sylvia

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: problem with primitive scorer?   (Makoto Asai - 28 Mar, 2006)
(_ None: Re: problem with primitive scorer?   (sylvia studeny - 29 Mar, 2006)
(_ None: Re: problem with primitive scorer?   (Sylvia Studeny - 03 Apr, 2006)
 Add Message Add Message
to: "problem with primitive scorer?"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]