[Halld-offline] Random number seeds in mcsmear

Mark M. Ito marki at jlab.org
Mon Jan 9 13:37:17 EST 2012


David,

Don't we want to save the hdgeant random seeds so they can be used for 
another run of hdgeant later? Whey does mcsmear use the random seeds 
used by hdgeant? I would think that mcsmear needs to save its own seeds 
in case mcsmear needs to be re-run again.

Shouldn't the default remain the same as it was before: use a new random 
sequence? In other words the default should be to ignore any stored 
seeds when generating random numbers, I would think.

My 2 cents.

   -- Mark

On 01/09/2012 11:24 AM, David Lawrence wrote:
> Hi All,
>
>     I have committed changes to the repository this morning to implement
> the storage and retrieval of the random number seeds used by mcsmear.
> This required several changes:
>
> 1.) The random number generator used in mcsmear was changed from
> TRandom3 to TRandom2. This is because the former uses 624 random number
> seeds and the latter only 3. TRandom2 does have a period of 10^26 so
> this should still be sufficient.
>
> 2.) The HDDM data model was modified to include slots for the 3 seeds
> that need to be stored every event. This means old HDDM files will cause
> warnings when being read into executables linked with the new data
> model. They should still work except for the new random seed features
> may not function properly.
>
> 3.) An option "-i" was added to mcsmear to allow one to specify that any
> seeds found in the input HDDM file should be ignored. Thus, allowing one
> to re-run mcsmear using the same input file and get a different result.
> This was previously the default behavior with no option to change it.
> Note that if you are using a file with no mcsmear seeds in it (the most
> common case), specifying "-i" will have no effect and the sequence will
> be different each time you run mcsmear with that input file.
>
> Please let me know if you run into any problems with hdgeant or mcsmear
> related to these changes or if you have any questions.
>
> Regards,
> -David
>
>
>
> _______________________________________________
> Halld-offline mailing list
> Halld-offline at jlab.org
> https://mailman.jlab.org/mailman/listinfo/halld-offline


-- 
Mark M. Ito
Jefferson Lab (www.jlab.org)
(757)269-5295




More information about the Halld-offline mailing list