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/
> >
>
>
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:28 EST