Message: parameterisation in kRho direction Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None parameterisation in kRho direction 

Forum: Geometry
Date: 22 Feb, 2007
From: Lingyun Chen <Lingyun Chen>

Hi, I am trying to parameterise concentric cylindrical ring, with increasing radius and constant width, in KRho direction, the problem is my run aborted with error information:

Start closing geometry.
PANIC! - G4SmartVoxelHeader::BuildNodes()
         Daughter physical volume Tally
         is entirely outside mother logical volume Phantom !!

*** G4Exception : InvalidSetup
      issued by : G4SmartVoxelHeader::BuildNodes()
Overlapping daughter with mother volume.
*** Fatal Exception *** core dump ***

*** G4Exception: Aborting execution *** Aborted

Below is my parameterisation, Anyone can point out what's wrong? Thanks a lot.

Lydia

______________________________________________________________

  G4Tubs* sPhantom = new G4Tubs("Phantom",0*cm,6*cm,
				5*cm,0.,360.);

  G4LogicalVolume* lPhantom = new G4LogicalVolume(sPhantom,phantomMaterial,"Phantom");	

  	  new G4PVPlacement(0,				//no rotation
  			   G4ThreeVector(20*cm,0*cm,0*cm),		
                           lPhantom,			//logical volume
			   "Phantom",		//name
                           lWorld,	       		//mother  volume
                           false,			//no boolean operation
                           0);				//copy number

    G4Tubs* sTally = new G4Tubs("Tally",0*mm,2*mm,1*mm,0.,360.);
    G4LogicalVolume* lTally = new G4LogicalVolume(sTally,tallyMaterial,"Tally");

  G4VPVParameterisation* tallyParam = new TallyParameterisation();
 //  G4VPhysicalVolume *tally_phys =
      new G4PVParameterised("Tally",lTally,lPhantom,kRho,20,tallyParam);

///////////////////parameterisation///////////////////////////////

void TallyParameterisation::ComputeTransformation
(const G4int,G4VPhysicalVolume *physVol) const
{
  G4ThreeVector origin;
  physVol->SetTranslation(origin);
}

void TallyParameterisation::ComputeDimensions
(G4Tubs & tallyRing, const G4int copyNo, const G4VPhysicalVolume*) const
{
  G4double tally_thick = 2.*mm;  
  G4double tally_dr = 2.*mm;
  G4double tallyTubRmin = copyNo*tally_dr; 
  tallyRing.SetInnerRadius(tallyTubRmin);
  tallyRing.SetOuterRadius(tallyTubRmin+tally_dr);
  tallyRing.SetZHalfLength(tally_thick/2);
  tallyRing.SetStartPhiAngle(0.*deg);
  tallyRing.SetDeltaPhiAngle(360.*deg);
}

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

1 None: Re: parameterisation in kRho direction   (M Tariq Siddique - 26 Sep, 2008)
 Add Message Add Message
to: "parameterisation in kRho direction"

 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 ]