Re: Activity value 0.0000E+00???

Hello again Angelo,

I think it is a problem of physics cards
Quoting the releases notes which comes with FLUKA :

***********************************************
-- IMPORTANT WARNINGS FOR THE USERS --

- Whenever residual nuclei (and residual dose rates) scoring is of
  importance, the heavy residue emission ("fragmentation") and the
  coalescence emission of fast complex particles should be switched on,
  through the following data cards:

PHYSICS 3.0
EVAPORAT
PHYSICS 1.0
COALESCE

  and (as a consequence of coalescence) it would be wise to link with
  rQMD-2.4 (and DPMJET) and activate ion transport and interactions. These
  suggestions are mandatory for residual nuclei benchmarking and validation.

  Those options are not on by default because the heavy evaporation
carries a
  big CPU penalty which would be a waste for most problems when
residuals are
  not a issue.
*************************************************************************

Retry your input file with those options on, and maybe a longer
irradiation time, 1 second is maybe a little short.

Hope it helps.
Regards.

Sebastien.

Angelo Infantino wrote :

> Sorry again, I try to use your advice but i don't have results.
> Try to look this .inp file. Is a simulation on a real thing.
> Resnuclei is empty and the value on the USRBIN is 0.
> If you read the sum.list file of the Resnuclei you can see that the
> beam intensity is 0.0166 pr/s...Why?
> My value of primary is 1000000 and my particle current is in the order
> of 10^15 p/s.
> Regards
>
> Angelo
>
> ------------------------------------------------------------------------
> Date: Wed, 25 Nov 2009 13:08:04 +0100
> From: wurth_at_ipno.in2p3.fr
> To: angeloinfa_at_hotmail.com
> CC: fluka-discuss_at_fluka.org
> Subject: Re: Activity value 0.0000E+00???
>
> Hello Angelo,
>
> I began my answer before I saw your second mail, that's the reason why
> I started the explanation from the beginning.
>
> To understand properly you should have some real example : I have an
> electron beam of 6.25e13 e- per second, for say a week (convert in
> seconds 86400*7), here are the two values for "irrprofile". Then you
> want to calculate activity in a given region (resnuclei), one minute
> after the end of irradiation, here is your cooling time : 60.0.
> You should have :
> * a raddecay card to activate the radioactive decay
> * an irrprofi card to set your irradiation profile
> * a decay time card
> * a resnuclei card
> * a dcyscore card to link the last two
>
> And of course, you have to create some nuclear reaction in that given
> region !
> Make sure your physical properties are set in the proper way.
>
> When DCYSCORE is applied to a detector, results are given per unit
> time (for instance Bq for a resnuclei card).
> The number of primaries put in the START card has nothing to do with
> the number of particles involved in the physical process whereas the
> one given in the IRRPROFIle card has.
>
> I think I never tried witout defining a profile, I don't know if it
> could give some result.
>
> Hope it helps,
> Regards.
>
> Sebastien.
>
>
> Angelo Infantino wrote :
>
> Sorry, but probably i don't have understand very good.
> I tryed to do a simplest simulation to understand better the use
> of this card but
> i have the same results.
> Activity 0.0000E+00 and Resnuclei empty. What is wrong?
> Another question. What mean have the number of primary when i put
> a particle current in the IRRPROFI card?
> For example if i put primary=100 and p/s=100000, what is the
> number of particle that i'm simulating??? 100 or 100000?
> Thank you
>
> Angelo Infantino
>
>
>
> > Date: Tue, 24 Nov 2009 12:50:42 +0100
> > From: wurth_at_ipno.in2p3.fr <mailto:wurth_at_ipno.in2p3.fr>
> > To: angeloinfa_at_hotmail.com <mailto:angeloinfa_at_hotmail.com>
> > CC: fluka-discuss_at_fluka.org <mailto:fluka-discuss_at_fluka.org>
> > Subject: Re: Activity value 0.0000E+00???
> >
> > Hello Angelo,
> >
> > You should define some cooling time (card DCYTIMES) and some link
> > between those cooling times (card DCYSCORE) and your detectors
> > (RESNUCLEi or USRBIN) and some irradiation profile (card
> IRRPROFIle)
> > (see manual for those options). The RADDECAY card only is not
> enough to
> > get radioactive activites.
> >
> > Hope it helps.
> >
> > Regards.
> > Sebastien.
> >
> >
> > Angelo Infantino wrote :
> >
> > > Good morning to everyone,
> > > i'm trying to calculate the activity in the several regions of my
> > > target. I know that i produce some radionuclides during the
> > > irradiation by protons of my target but the value of the
> activity in
> > > the several regions is always 0.0000E+00. How it's possible?
> I'm also
> > > tryed to active the RADDECAY card but i obtain the same
> result. What
> > > is wrong in my input file? I attach the .inp file.
> > > Thank you very much for yours help.
> > > Regards
> > >
> > > Angelo Infantino
> > > <http://www.windowslive.it/moviemaker.aspx>
> >
>
> ------------------------------------------------------------------------
> Troppi account di posta? Riuniscili con Windows Live Mail
> <http://www.windowslive.it/mail.aspx>
>
>
>
> ------------------------------------------------------------------------
> Personalizza il tuo spazio online: ritagliati il tuo Spaces
> <http://www.windowslive.it/spaces.aspx>
> ------------------------------------------------------------------------
>
>------------------------------------------------------------------------
>
>TITLE
>Prova di irraggiamento
>DEFAULTS NEW-DEFA
>BEAM -.0165 1.0PROTON
>BEAMPOS 0.0 0.0 0.0 0.0 0.0
>GEOBEGIN COMBNAME
> 0 0
>* Black body
>SPH BLK 0.0 0.0 0.0 10000000.0
>RPP VOI -500. 500. -500. 500. -500. 500.
>ZCC TARG 0.0 0.0 1.
>XYP P1 0.0
>XYP P2 0.1
>END
>* Black hole
>BLKBODY 5 +BLK -VOI
>* Void around
>VOID 5 +VOI -TARG | +VOI +P1 | +VOI -P2
>* Target
>TARGET 5 +TARG +P2 -P1
>END
>GEOEND
>* AIR
>* Composition (%):
>* 21 O, 79 N. Density: 0.0013 g/cm3
>MATERIAL .0013 AIR
>COMPOUND -.21 OXYGEN -.79 NITROGEN AIR
>* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
>ASSIGNMA BLCKHOLE BLKBODY
>ASSIGNMA AIR VOID
>ASSIGNMA NICKEL TARGET
>RADDECAY 2.
>USRBIN 2. ACTIVITY 50. TARGET Act0
>USRBIN TARGET 1. &
>DCYSCORE 1. Act0 USRBIN
>USRBIN 2. ACTIVITY 51. TARGET Act600
>USRBIN TARGET 1. &
>DCYSCORE 3. Act600 USRBIN
>RESNUCLE 3. -52. TARGET Res0
>DCYSCORE 2. Res0 RESNUCLE
>IRRPROFI 1.1000000000
>DCYTIMES 0.0 60. 600.
>RANDOMIZ 1.0
>START 1000000. 0.0
>STOP
>
>
Received on Wed Nov 25 2009 - 15:40:43 CET