Re: LET spectrum in a spacecraft Si component

From: Ercan Pilicer (ercanpilicer@gmail.com)
Date: Sat Nov 10 2007 - 13:48:02 CET

  • Next message: Niels Bassler: "USRTRACK fluence spectrum scoring in GeV/nucleon instead of GeV?"

    dear all

    i have seen that there is a previous post about LET spectra..well..

    but the problem is still the distribution of particles..most of the
    source particles escape from reg2 to reg1..but very very few of them
    enter to reg3...it it normal?

    source definition is like;

    * Kinetic energy of the particle (GeV)
          TKEFLK (NPFLKA) = ENERGY

    * Particle momentum
          PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
         & + TWOTWO * AM (ILOFLK(NPFLKA)) ) )

          rsph = whasou(1)

    * Isotropic and uniform distribution
          CALL SFLOOD ( XXX,YYY,ZZZ,UXXX, VYYY,WZZZ )

    * Cosines (tx,ty,tz)
          TXFLK (NPFLKA) = UXXX
          TYFLK (NPFLKA) = VYYY
          TZFLK (NPFLKA) = WZZZ

    * Polarization cosines:
          TXPOL (NPFLKA) = -TWOTWO
          TYPOL (NPFLKA) = +ZERZER
          TZPOL (NPFLKA) = +ZERZER

    * Particle coordinates
          XFLK (NPFLKA) = XXX * rsph
          YFLK (NPFLKA) = YYY * rsph
          ZFLK (NPFLKA) = ZZZ * rsph

    when i dump them within mgdraw.f

          write(91,*) mreg,newreg,icode,jtrack,etrack,am(jtrack)

    i see most of them (%99) goes to the reg1, but less to reg3.

    and in the input i wrote;

    ...
    BEAM -0.4 PROTON
    SOURCE 100.0
    *
    GEOBEGIN COMBNAME
        0 0 experimental hall geometry
    SPH B1 0.0 0.0 0.0 10000.0
    SPH B2 0.0 0.0 0.0 100.0
    SPH B3 0.0 0.0 0.0 5.1
    SPH B4 0.0 0.0 0.0 5.0
    ...

    something i have done wrong? should i define radius (whasou(1)) less than 100?

    any help would be apperciated

    ciao
    e.

    On Nov 7, 2007 6:27 PM, Ercan Pilicer <ercanpilicer@gmail.com> wrote:
    >
    > dear All,
    >
    > the reason of that work is to simulate a LET (MeV cm2 g-1) spectrum
    > vs. differential flux (m-2 sr-1 s-1 (MeV/n)-1 ) as can be done by
    > spenvis. but i dont know if i can :)
    >
    > anyone managed to do that with fluka?
    >
    > what i have is;
    > a spectrum of energy (MeV) vs. trapped proton differential flux (cm-2
    > MeV-1 s-1) from spenvis. i modified source.f from the discussions by
    > Hamide, Ferrari and Mario.
    >
    > what i encounter is;
    > primary source particles go out of the vacuum void to external
    > blackhole (from region2 to region1) which is not the case i assumed.
    > i am missing a point i assume. what could that be?
    >
    > thanks in advance.
    > e.
    >
    >
    > ps.
    > attached files are source.f, mgdraw.f, pbbsi.inp
    >
    >
    >
    > ==============================
    > ()".'."() Ercan Pilicer
    > ( (T) ) Uludag University, Turkey
    > (ö) (ö) High Energy Physics Department
    > ==============================
    > 

    -- 
==============================
()".'."()    Ercan Pilicer
(  (T)  )   Uludag University, Turkey
(ö)   (ö)   High Energy Physics Department
==============================
  • Next message: Niels Bassler: "USRTRACK fluence spectrum scoring in GeV/nucleon instead of GeV?"

    This archive was generated by hypermail 2.1.6 : Sun Nov 11 2007 - 18:29:55 CET