Re: monte carlo details

From: Timothy Smith (tim_smith@MIT.EDU)
Date: Thu Sep 05 2002 - 15:34:14 EDT


Hello Aaron,

        I just read through the fort.99 . It is very well laid out and
explained - thank you . I thin it is the last word on blastmc - until
the next revolution in blastmc.

        Thanks!
                                Tim Smith

On Thu, 5 Sep 2002, Aaron Joseph Maschinot wrote:

>
> Hello, Adrian and anyone else interested in the Monte Carlo:
>
> About a month ago, I checked in a fort file (called fort.99) to the
> Blast_Params directory in CVS. In this file, you will find a
> thoroughly-detailed exaplanation of EVERY fort file command valid for our
> Monte Carlo. It also explains in detail all of the options for all of the
> commands (the file is roughly 900 lines long). I think that, if you look
> at the file, all fort file questions will be answered -- including the
> vertex questions you have been asking. (The triangularly-distributed
> vertex is the default vertex distribution style but ONLY under the NEW
> Monte Carlo files (I also checked in a bunch of "corrections" to
> blast_init.f and blast_proc.f); the old "default" (which you are probably
> using), assuming you didn't specify a vertex distribution style, was the
> point origin.).
>
> Aaron
>
>
>
>
> On Thu, 5 Sep 2002, Timothy Smith wrote:
>
> > Hello Adrian,
> >
> > The size and distribution of you target depends upon the input
> > parameters in your fort.3* file. for example,
> >
> >
> > VERT 2 0. 0. -20.0 0. 0. 40.0
> >
> > means that I am using vertex generator 2
> > which then means it is uniform 0. < x < 0. + 0.
> > 0. < y < 0. + 0.
> > -20. < z < -20. + 40. = +20.
> >
> > All the different vertex generators are in the subroutine KIVERT located
> > in the file "blast_proc.f" . There is a good explanation in the top of
> > that subroutine which is below.
> >
> > Tim
> >
> > ____________________________________________________________________
> > Timothy Paul Smith Assistant Research Professor
> > still using my MIT email for now -> tim_smith@mit.edu
> > Dartmouth College
> > Department of Physics and Astronomy
> > 6127 Wilder Lab.
> > Hanover, New Hampshire 03755
> >
> > =============================================
> >
> >
> > C ******************************************************************
> > C * SUBROUTINE: KIVERT *
> > C * FUNCTION: Generates vertex coordinates according to the data *
> > C * card VERT. The format of the data card is: *
> > C * VERT IFFVER FFVERT(1) FFVERT(2) FFVERT(3) ... *
> > C * The first parameter (IFFVER) gives the type of *
> > C * vertex distribution: *
> > C * IFFVER = 0: (default vertex style) Distributed *
> > C * with a triangular z distribution *
> > C * 40.0cm long and a flat distribution *
> > C * in r and theta out to 0.5mm (1.0mm *
> > C * in diameter for the beam X target *
> > C * profile). *
> > C * IFFVER = 1: Normal distribution centered at *
> > C * (X, Y, Z) = (FFVERT(1), FFVERT(2), *
> > C * FFVERT(3)) and having (SIGMAX, *
> > C * SIGMAY, SIGMAZ) = (FFVERT(4), *
> > C * FFVERT(5), FFVERT(6)) *
> > C * IFFVER = 2: Uniform distribution such that: *
> > C * FFVERT(1) <= X <= FFVERT(4) *
> > C * FFVERT(2) <= Y <= FFVERT(5) *
> > C * FFVERT(3) <= Z <= FFVERT(6) *
> > C * IFFVER = 3: Uniform distribution within a tube *
> > C * of radius FFVERT(1) and of length *
> > C * FFVERT(2) in z direction, centered *
> > C * at FFVERT(3,4,5) in Cartesian *
> > C * coordinates *
> > C * IFFVER = else: The point origin is the origin *
> > C * CALLED BY: GUKINE *
> > C * AUTHOR: M. Guckes (July 1986) *
> > C * M. Guckes (October 1987) *
> > C * A. Maschinot (July 2002) *
> > C * commented, corrected, and updated for the BLAST *
> > C * 2002 experiments *
> > C ******************************************************************
> > SUBROUTINE KIVERT(NVTX, UBUF, NUBUF)
> >
> >
> >
> >
> >
> >
> > On Thu, 5 Sep 2002, Adrian T Sindile wrote:
> >
> > > Hi, Jason, everybody!
> > > Before leaving on vacation, I left Tancredi a couple of CODA files (Monte
> > > Carlo generated).
> > > He told me that it was possible that I did not have the extended target.
> > > To the best of my knowledge, that has always been hardwired (with the
> > > corresponding triangular distribution). I will also look at my code and
> > > doublecheck, but Tancredi suggested I contact you.
> > > If anybody knows otherwise about the extended target, please let me know!
> > >
> > > Adrian
> > >
> > > -------------------------------
> > > Adrian Sindile
> > > Research Assistant
> > > Nuclear Physics Group
> > > University of New Hampshire
> > > phone: (603)862-1691
> > > FAX: (603)862-2998
> > > email: asindile@alberti.unh.edu
> > > http://einstein.unh.edu/~adrian/
> > >
> >
> >
>

-- 
____________________________________________________________________
 Timothy Paul Smith                     Assistant Research Professor
      still using my MIT email for now ->          tim_smith@mit.edu
      no voice mail yet                ->       tel:  (603) 646-9346
 Dartmouth College
 Department of Physics and Astronomy
 6127 Wilder Lab. 
 Hanover, New Hampshire 03755



This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:28 EST