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

Idea Re: Ice 

Forum: Geometry
Re: Question Ice (Zack Snyder)
Re: Feedback Re: Ice (Joel Long)
Date: 14 Jun, 2013
From: Zack Snyder <Zack Snyder>

Hello Joel,

Thank you for your reply. From what I can tell, the same things appear to be occurring in both water and ice, which leads me to believe that everything is now working fine. I was misunderstanding the output slightly, so I am partially to blame. One of my concerns is about the molecular structure of ice versus water. Ice has a crystal structure and I do not know if this is taken into account by Geant4 when building ice via the code I quoted in a previous message or if Geant4 must be told to do so in some fashion. In addition, I do not know whether this will have a noticeable effect on neutron interaction or not.

On Thu, 13 Jun 2013 22:19:44 GMT, Joel Long wrote:

> Hello Zach,
> 
> I do not immediately see anything obviously wrong with your ice
> definition, but there is an easy way to check: run your simulation with
> NIST water, and see if the results there are what you expect. If
> unexpected things are still happening, then your material is probably
> not your issue (though I would leave it as water until you were getting
> the results you expected, _then_ change it back to ice; debugging is all
> about isolating possible issues to examine individually).
> 
> It is difficult to make any suggestions beyond that without knowing A.
> what neutron physics you are trying to use and B. what type of anomalies
> you are seeing. If you haven't already, I suggest posting a new thread
> in the Physics List or Hadronic Processes forums with this information.
> 
> Hope this helps, Joel Long
> 
> On Tue, 16 Apr 2013 21:55:28 GMT, Zack Snyder wrote:
> 
> > 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 detectorconstruction.cc:
> > 
> > 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,         
> >                         shape1_mat,          
> >                         "Shape1");           
> > 
> >   new G4PVPlacement(0,                      
> >                     pos1,                    
> >                     logicShape1,             
> >                     "Shape1",                
> >                     logicEnv,                
> >                     false,                   
> >                     0,                       
> >                     checkOverlaps); 
> > 
> > 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 More: Re: Ice   (Joel Long - 14 Jun, 2013)
(_ None: Re: Ice   (Michael H. Kelsey - 16 Jun, 2013)
 Add Message Add Message
to: "Re: 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 ]