Message: Ice Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Ice 

Forum: Geometry
Date: 16 Apr, 2013
From: Zack Snyder <Zack Snyder>

Good evening,

First, I apologize if this is a simple/dumb question.

I am trying to define ice in geant4 and then shoot neutrons into the ice and get the track length. I have defined what I believe to be ice (and compared against water imported from nist material builder, with not much of a difference), however the professor that I am working under is questioning what I show him to be happening.

Here is what I have in my

density = 920*mg/cm3;
 temperature = 216.15*kelvin;
 G4Material* shape1_mat = new G4Material(name="shape1_mat", density, ncomponents=2, kStateSolid,temperature);

shape1_mat->AddElement(elH, natoms=2);
shape1_mat->AddElement(elO, natoms=1);
G4ThreeVector pos1 = G4ThreeVector(0*m, 0*m, 0*m);
G4double shape1_dx = 100*m, shape1_dy = 100*m, shape1_dz = 100*m;
  G4Box* solidShape1 =
    new G4Box("Shape1",                    
         1*shape1_dx, 1*shape1_dy, 1*shape1_dz)
  G4LogicalVolume* logicShape1 =
    new G4LogicalVolume(solidShape1,         

  new G4PVPlacement(0,                      

No errors arise when I run make, I just believe there to be an error somewhere in here that is causing an improper definition of ice, and therefore improper interactions with neutrons.

Also, perhaps it is because I am forgetting to include a library necessary for neutrons in one of the other .cc or .hh files?

Again, sorry if this question is dumb.

Thanks, Zack

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

1 Feedback: Re: Ice   (Joel Long - 13 Jun, 2013)
(_ Idea: Re: Ice   (Zack Snyder - 14 Jun, 2013)
(_ More: Re: Ice   (Joel Long - 14 Jun, 2013)
(_ None: Re: Ice   (Michael H. Kelsey - 16 Jun, 2013)
 Add Message Add Message
to: "Ice"

 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 ]