|Message: Re: Ice||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
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:||Outline Depth:||Add message:|