Message: Geant4 BUG???wrong material with G4PhantomParameterisation to get dose Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Geant4 BUG???wrong material with G4PhantomParameterisation to get dose 

Forum: Event and Track Management
Date: 22 Jul, 2013
From: Geng <Geng>

Dear experts: I am meeting a very confused problem now. When I use G4PhantomParamerisation to construct my geometry and G4PSDosedeposit to score dose.Then, I want to check the right material with input ID. So, I modified G4psdosedeposit with:

///////////////////////////////////////////////////////////////////////////////
  G4TouchableHandle theTouchable = aStep->GetTrack()->GetStep()->GetPreStepPoint()->GetTouchableHandle();
  G4int index=theTouchable->GetCopyNumber();
      G4String materialname=aStep->GetTrack()->GetStep()->GetPreStepPoint()->GetMaterial()->GetName();
       const DicomDetectorConstruction* detector =
     (const DicomDetectorConstruction*)
     (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
if(detector->fMateIDs[index]==0 )
  {if(materialname!="1a") G4cout<<"error"<<materialname<<G4endl;}
/////////////////////////////////////////////////////////////////////////////

"1a" is the material name of ID=0. In my opinion, there should not occur any output of "error". But actually, it does.

So, I have no idea about the problem of this case.I think this will be a great mistake to get dose in simulation, because the mass of the volume might be wrong.

And I have no where to go to ask this problem except here,so I really need your help.

Any advice or solution will be great appreciate,thanks.

Geng

 Add Message Add Message
to: "Geant4 BUG???wrong material with G4PhantomParameterisation to get dose"

 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 ]