RE: Synchrotron radiation

Dear Alberto and all,

I want to see synchrotron radiation energy before target and i use USRBIN card for that.
But there is a red line center of my plot (see attachment).
Is it an error or what??
Can you help me..??

Thanks,
Emre

On Tue, 1 Jun 2010 01:19:43 -0700, Fasso, Alberto wrote
> Yes, you are right. I thought that I had said it sometimes in one of the ma=
> ny emails this thread has developed in.
> If I haven't, I apologize
>
> Alberto
> ________________________________________
> From: Nicholas Bolibruch [nboli_at_acanac.net]
> Sent: Monday, May 31, 2010 1:58 PM
> To: Fasso, Alberto
> Subject: Re: Synchrotron radiation
>
> Dear Alberto,
>
> In regards to the Source.f file you posted...
>
> I don't think you mentioned that the SOURCE card needs WHAT(1) set to 1
> for the OPEN statement to work for io stream 81. I spent enough time
> today trying to debug this problem, was even running it in GDB trying to
> figure out where I was getting a Floating Exception from ...
>
> --
> Nick Bolibruch
> Controls & Instrumentation Development
> Canadian Light Source Inc.
>
> Alberto Fasso' wrote:
> > Dear Emre,
> >
> > you have not specified what "some errors" are. But your BEAM card and you=
> r
> > source file are assuming that your maximum photon energy is 2 MeV (0.002 =
> GeV).
> > The spectrum in the sr.spectrum file is in MeV, which are translated to G=
> eV
> > inside the source routine:
> > DO 2 I =3D 1, 630
> > READ(81,*,END=3D999) EN(I), SPEC(I)
> > EN(I) =3D EN(I)*1.D-3 <----
> > 2 CONTINUE
> > On the other hand, the spectrum in the sr.spectrum file has a maximum of =
> 1995
> > MeV (1.995 GeV), which is very unlikely. And the sampling is done between
> > 1.D-6 (1 keV, corresponding to the first line of the spectrum) and 1.9952=
> 6D-4
> > (about 200 keV):
> > ENERGY =3D 1.D-6 + ( 1.99526D-4 - 1.D-6) * RAN01
> >
> > To conclude:
> > 1) check your spectrum file: I guess it should go from 1 keV to 200 keV, =
> and
> > not from 1 keV to 2 GeV.
> > 2) Your BEAM card is not wrong, but if the maximum is really 200 keV, set=
> ting
> > it to 2 MeV would make your cross section tabulations less precise
> > without need
> >
> > Kind regards,
> >
> > Alberto
> >
> > On Mon, 31 May 2010, Emre EROGLU wrote:
> >
> >
> >> Dear Alberto,
> >>
> >> Many thanks for your mail. But still i have some errors. I didn't find w=
> here is it. I'm
> >> sending my .inp file, sr_source.f (modified from bl1source.f) and sr.spe=
> ctrum (generated
> >> from srplot.f). And also are there any mistakes in my BEAM card?
> >> Thanks.
> >>
> >> Best Regards,
> >> Emre
> >>
> >>
> >>
> >> On Thu, 27 May 2010 06:35:42 -0700 (PDT), Alberto Fasso' wrote
> >>
> >>> I apologize. On April 28, answering a question by Emre Eroglu, I sent
> >>> an example of source routine sampling from a synchrotron radiation spec=
> trum
> >>> provided as an external file. See:
> >>> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/2879.html
> >>>
> >>> I did not realize at that time that it was an old routine, incompatible=
> with
> >>> the latest FLUKA versions. I send here in attachment a corrected versio=
> n.
> >>>
> >>> Sorry!
> >>>
> >>> Alberto
> >>>
> >>> On Wed, 26 May 2010, Emre EROGLU wrote:
> >>>
> >>>
> >>>> Dear Alberto,
> >>>>
> >>>> I try to compile bl1source.f but i have some errors. what can i do for=
> that?
> >>>>
> >>>> Regards,
> >>>> Emre
> >>>>
> >>>>
> >>>> -e g77 -o bl1source.o -c -O3 -g -mcpu=3Dpentiumpro -fexpensive-optimiz=
> ations
> >>>> -funroll-loops -fstrength-reduce -Wall -fno-automatic -fno-silent -m32
> >>>> -ffortran-bounds-check -I/home/eeroglu/fluka/flukapro bl1source.f
> >>>> `-mcpu=3D' is deprecated. Use `-mtune=3D' or '-march=3D' instead.
> >>>> source:
> >>>> bl1source.f: In subroutine `source':
> >>>> bl1source.f:11:
> >>>> INCLUDE '(AACOLL)'
> >>>> ^
> >>>> Unable to open INCLUDE file `(AACOLL)' at (^)
> >>>> bl1source.f:12:
> >>>> INCLUDE '(BEAM)'
> >>>> ^
> >>>> Unable to open INCLUDE file `(BEAM)' at (^)
> >>>> bl1source.f:13:
> >>>> INCLUDE '(EPISOR)'
> >>>> ^
> >>>> Unable to open INCLUDE file `(EPISOR)' at (^)
> >>>> bl1source.f:17:
> >>>> INCLUDE '(STACK)'
> >>>> ^
> >>>> Unable to open INCLUDE file `(STACK)' at (^)
> >>>> bl1source.f:18:
> >>>> INCLUDE '(STARS)'
> >>>> ^
> >>>> Unable to open INCLUDE file `(STARS)' at (^)
> >>>> bl1source.f:39:
> >>>> LUSSRC =3D .TRUE.
> >>>> 1 2
> >>>> Type disagreement between expressions at (1) and (2)
> >>>> bl1source.f:78:
> >>>> WT (LSTACK) =3D WEIGHT
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `wt' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:90:
> >>>> ILO (LSTACK) =3D IJHION
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `ilo' at (^) [initially =
> seen at (^)]
> >>>> bl1source.f:102:
> >>>> LO (LSTACK) =3D 1
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `lo' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:104:
> >>>> LOUSE (LSTACK) =3D 0
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `louse' at (^) [initiall=
> y seen at
> >>>> (^)]
> >>>> bl1source.f:107:
> >>>> SPAREK (ISPR,LSTACK) =3D ZERZER
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `sparek' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:111:
> >>>> ISPARK (ISPR,LSTACK) =3D 0
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `ispark' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:116:
> >>>> NUMPAR (LSTACK) =3D NPARMA
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `numpar' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:117:
> >>>> NEVENT (LSTACK) =3D 0
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `nevent' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:118:
> >>>> DFNEAR (LSTACK) =3D +ZERZER
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `dfnear' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:121:
> >>>> AGESTK (LSTACK) =3D +ZERZER
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `agestk' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:122:
> >>>> AKNSHR (LSTACK) =3D -TWOTWO
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `aknshr' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:124:
> >>>> IGROUP (LSTACK) =3D 0
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `igroup' at (^) [initial=
> ly seen at (^)]
> >>>> bl1source.f:126:
> >>>> TKE(LSTACK) =3D ENERGY
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `tke' at (^) [initially =
> seen at (^)]
> >>>> bl1source.f:128:
> >>>> PMOM (LSTACK) =3D ENERGY
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `pmom' at (^) [initially=
> seen at (^)]
> >>>> bl1source.f:132:
> >>>> TX (LSTACK) =3D TINX
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `tx' at (^) [initially s=
> een at (^)
> >>>> ]
> >>>> bl1source.f:133:
> >>>> TY (LSTACK) =3D TINY
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `ty' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:134:
> >>>> TZ (LSTACK) =3D TINZ
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `tz' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:137:
> >>>> TXPOL (LSTACK) =3D -TWOTWO
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `txpol' at (^) [initiall=
> y seen at (^)]
> >>>> bl1source.f:138:
> >>>> TYPOL (LSTACK) =3D +ZERZER
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `typol' at (^) [initiall=
> y seen at (^)]
> >>>> bl1source.f:139:
> >>>> TZPOL (LSTACK) =3D +ZERZER
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `tzpol' at (^) [initiall=
> y seen at (^)]
> >>>> bl1source.f:141:
> >>>> XA (LSTACK) =3D XINA
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `xa' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:142:
> >>>> YA (LSTACK) =3D YINA
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `ya' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:143:
> >>>> ZA (LSTACK) =3D ZINA
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `za' at (^) [initially s=
> een at (^)]
> >>>> bl1source.f:163:
> >>>> NLATTC (LSTACK) =3D MLATTC
> >>>> ^
> >>>> Invalid declaration of or reference to symbol `nlattc' at (^) [initial=
> ly seen at (^)]
> >>>> In file included from bl1source.f:11:
> >>>> In file included from bl1source.f:17:
> >>>> bl1source.f:160: warning: unused variable 'idisc'
> >>>> bl1source.f:143: warning: unused variable 'zina'
> >>>> bl1source.f:142: warning: unused variable 'yina'
> >>>> bl1source.f:141: warning: unused variable 'xina'
> >>>> bl1source.f:134: warning: unused variable 'tinz'
> >>>> bl1source.f:133: warning: unused variable 'tiny'
> >>>> bl1source.f:132: warning: unused variable 'tinx'
> >>>> bl1source.f:79: warning: unused variable 'weipri'
> >>>> bl1source.f:39: warning: unused variable 'lussrc'
> >>>> bl1source.f:21: warning: unused variable 'line'
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>
> >>>> On Wed, 28 Apr 2010 01:42:18 -0700 (PDT), Alberto Fasso' wrote
> >>>>
> >>>>> Dear Emre,
> >>>>>
> >>>>> as you know, FLUKA for the time being cannot simulate the production =
> of
> >>>>> synchrotron radiation. A source routine, therefore, must sample from =
> a
> >>>>> spectrum produced by a separate program (see
> >>>>> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/1183.html=
> )
> >>>>> I attach here an example, where the spectrum is provided as an
> >>>>> external file
> >>>>>
> >>>>> Alberto
> >>>>>
> >>>>> On Tue, 27 Apr 2010, emreeroglu_at_uludag.edu.tr wrote:
> >>>>>
> >>>>>
> >>>>>> Hi all,
> >>>>>>
> >>>>>> Does anyone have source.f which is modified for synchrotron radiatio=
> n?
> >>>>>>
> >>>>>> Thanks in advance.
> >>>>>>
> >>>>>> Regards,
> >>>>>> Emre EROGLU
> >>>>>>
> >>>> Emre EROGLU
> >>>> Uludag University
> >>>> High Energy Physcics
> >>>>
> >>>>
> >>>>
> >>>>
> >>> --
> >>> Alberto Fasso`
> >>> SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025
> >>> Phone: (1 650) 926 4762 Fax: (1 650) 926 3569
> >>> fasso_at_slac.stanford.edu
> >>>
> >> Emre EROGLU
> >> Uludag University
> >> High Energy Physcics
> >>
> >>
> >>
> >
> >=

Emre EROGLU
Uludag University
High Energy Physcics