Hi,<br><br>You are wright, the problem was already mentioned and is caused by the new random engine.<br>In fact this new engine takes different seeds than the last one but the gjs program was not updated to produce these new seeds ...<br>
<br>To fix this out you have to change the definition of the random engine at the beginning of the main program to set the old one instead of the new one.<br>I send you with this mail the file Gate.cc where I did the modification. Just replace the original Gate.cc file by this one and recompile, then it will use the old random engine with the old seed files generated by the gjs program and it will work fine.<br>
<br>A new patch is coming soon with which you can manage all the random engines and seeds aspects with macro commands.<br><br>Sorry for your CPU time lost !<br>Cheers,<br>Simon<br><br><div class="gmail_quote">On Wed, Dec 10, 2008 at 4:01 PM, KB <span dir="ltr"><<a href="mailto:kd8had@gmail.com">kd8had@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hello All,<br><br>GATE version name: gate_v4.0.0 (20-August-2008). I am running SGE. <br><br>
Using gjs, I split the main macro into 25 parts and submit them as a array job using qsub:<br><br>"qsub -t 1-25 -V /home/data/.Gate/Hits/Hits.sh"<br>
<br>This produced 25 root files which I am supposed to merge with gjm. I did so and noticed that the resultant root file has some strange behavior when some parameters were histogrammed. There was a periodicity of 25. So I randomly picked 2 of the 25 root files and using root printed out the events. They are unfortunately identical. This makes the culster version of GATE almost useless (at least for me). Does anyone else notice the same thing? <br>
<br>Here is some output from running root:<br>I open the first file:<br>root [0] TFile f2("CylinderHitsSingleGamma3.root")<br>root [1] Hits->Show(1)<br>======> EVENT:1<br> PDGEncoding = 22<br> trackID = 4<br>
parentID = 1<br> time = 0.400002<br> edep = 0.0018285<br> stepLength = 0.00627867<br> posX = 455.453<br> posY = 75.4412<br> posZ = -30.2955<br> localPosX = -1.67338<br>
localPosY = -0.606797<br> localPosZ = 0.267013<br> gantryID = 0<br> blockID = 0<br> crystalID = 339<br> unused3ID = -1<br> unused4ID = -1<br> unused5ID = -1<br> photonID = 1<br>
nPhantomCompton = 6<br> nCrystalCompton = 0<br> nPhantomRayleigh = 0<br> nCrystalRayleigh = 0<br> primaryID = 1<br> sourcePosX = 17.0792<br> sourcePosY = -100.416<br> sourcePosZ = -49.8924<br> sourceID = 0<br>
eventID = 3<br> runID = 0<br> axialPos = 0<br> rotationAngle = 0<br> volumeID = 0, <br> 0, 0, 339, -1, -1, -1, -1, -1, -1<br> processName = LowEnPhotoElec<br> comptVolName = wholecylinder_P<br>
RayleighVolName = NULL<br><br><br>Here is the 1st event from another root file:<br>root [0] TFile f1("CylinderHitsSingleGamma10.root")<br>root [1] Hits->Show(1)<br>======> EVENT:1<br> PDGEncoding = 22<br>
trackID = 4<br> parentID = 1<br> time = 1.8<br> edep = 0.0018285<br> stepLength = 0.00627867<br> posX = 455.453<br> posY = 75.4412<br> posZ = -30.2955<br>
localPosX = -1.67338<br> localPosY = -0.606797<br> localPosZ = 0.267013<br> gantryID = 0<br> blockID = 0<br> crystalID = 339<br> unused3ID = -1<br> unused4ID = -1<br> unused5ID = -1<br>
photonID = 1<br> nPhantomCompton = 6<br> nCrystalCompton = 0<br> nPhantomRayleigh = 0<br> nCrystalRayleigh = 0<br> primaryID = 1<br> sourcePosX = 17.0792<br> sourcePosY = -100.416<br> sourcePosZ = -49.8924<br>
sourceID = 0<br> eventID = 3<br> runID = 1<br> axialPos = 0<br> rotationAngle = 0<br> volumeID = 0,<br> 0, 0, 339, -1, -1, -1, -1, -1, -1<br> processName = LowEnPhotoElec<br>
comptVolName = wholecylinder_P<br> RayleighVolName = NULL<br><br>The only thing that changes is the time. Everything else is identical. Printint out the other events shows the same similarities. <br><br>Karthik.<br>
<br>_______________________________________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.healthgrid.org">Gate-users@lists.healthgrid.org</a><br>
<a href="http://lists.healthgrid.org/mailman/listinfo/gate-users" target="_blank">http://lists.healthgrid.org/mailman/listinfo/gate-users</a><br>
<br></blockquote></div><br>