Dear Vahan,
about
> 1) How can I record all prticles' coordinates, energy, and moving
> directions at the point of crossing boundary in some file?
> I read fluka manual and it says it can do that with collision tape....
> but special subrutine should be written.... how program that routine?
you can read the manual in the section:
http://www.fluka.org/fluka.php?id=man_onl&sub=81 related with BXDRAW and
in the chapter 13
http://www.fluka.org/fluka.php?id=man_onl&sub=99 in the
section 13.12.13 and in particular:
* BXDRAW (boundary crossing dumping)
Argument list (all variables are input only):
ICODE : FLUKA physical compartment originating the call
= 19: call from Kaskad (hadrons and muons)
= 29: call from Emfsco (electrons, positrons and photons)
= 39: call from Kasneu (low-energy neutrons)
= 49: call from Kashea (heavy ions)
= 59: call from Kasoph (optical photons)
MREG : number of region before boundary crossing
NEWREG : number of region after boundary crossing
XSCO, YSCO, ZSCO : coordinates of crossing point
BXDRAW is called at each boundary crossing (if requested by the user with
USERDUMP, WHAT(3) < 7.0). There is no default output: any output must be
supplied by the user.
If name-based input is being used, the names corresponding to MREG and
NEWREG can be obtained via a call to routine GEOR2N:
CALL GEOR2N (NUMREG, NAMREG, IERR)
where NUMREG (input variable) is a region number, and NAMREG (returned
variable) is the corresponding region name (to be declared as
CHARACTER*8). IERR is a returned error code: if = 0, the conversion is
successful. Example:
.......................................
CHARACTER*8 MRGNAM, NRGNAM
.......................................
ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO )
CALL GEOR2N ( MREG, MRGNAM, IERR1 )
CALL GEOR2N ( NEWREG, NRGNAM, IERR2 )
IF(IERR1 .NE. 0 .OR. IERR2 .NE. 0) STOP "Error in name
conversion"
.......................................
IF(MRGNAM .EQ. "MyUpsREG" .AND. NRGNAM .EQ. "MyDwnREG") THEN
.......................................
Best Regards,
Andrea
Received on Wed Nov 13 2013 - 00:37:10 CET