From: paolo maestro (paolo.maestro@pi.infn.it)
Date: Wed May 12 2004 - 12:40:18 CEST
Hi All,
I wrote an input card describing a simple detector consisting of 6
pixels obtained by cloning
a basic pixel unit by means of the LATTICE card and the routine
lattice.f .
I made two versions of the card, both in fixed and free format.
The first one is working , while using the second the program
crashes when it reads the LATTICE card reporting the following error :
Abort called from LATINP reason IERR>0 Run stopped !
STOP IERR>0
I put in attachment both the free format input card and the output
file.
Can you explain which is the meaning of the error message and why
it happens only with the free format ?
I have some more general questions:
- My complete problem is to simulate the response of a 2500 pixels
silicon matrix.
My original idea was to use the LATTICE capability as I did in my
trial card.
But I wonder if I could get the same results in a cheaper way
simulating a complete not segmented silicon plane and then
reconstructing the
pixel response by using the USRBIN card .
In other words, is this the philosophy of USRBIN and EVENTBIN
cards?
- Moreover, I would like to know if
it is possible to have a modular geometry in which the basic region
in its turn
consists of many identical elemental regions (nested modular
geometries).
For istance, consider a detector made of thousand strips, each one
consisting
of ten glued scintillating fibers. Can I use the LATTICE card for
such a problem ?
Do you suggest an alternative approach ?
If the answer to the first question is positive, I have a further
doubt.
At the tracking level, in the mgdraw.f subroutine it is possible to
know the
current region and lattice cell number traversed by a particle
track
from the MREG and LT1TRK variables .
Then, considering the strips/fibers example,
I could get the fiber number but how could I know the number
of the strip the fiber belongs to ?
Thanks for your help.
Paolo Maestro
TITLE
Test of simple calorimeter
GLOBAL 3000.0 0.0 0.0 2.0 1.0 0.
FREE
DEFAULTS 0. 0. 0. 0. 0. 0. CALORIME
*
*
*BEAM 20. 0. 0. 0. 0. 1.0 RAY
BEAM 20.0 0. 0. 0. 0. 1.0 ELECTRON
*BEAM 20. 0. 0. 0. 0. 1.0 PROTON
*BEAM 20. 0. 0. 0. 0. 1.0 HEAVYION
* Heavy ion case: BEAM WHAT(1) is kin energy (momentum) per nucleon NOT total
*HI-PROPE 5.0 10.0 0.0
*EVENTYPE 0. 0. 2. 0. 0. 0. DPMJET
*
* Particles must injected in vacuum or material not in black hole
*
BEAMPOS 0. 0. 0. 0. 0. 0.
SOURCE 0. 0. 30.0 0. 0. -1.
*SOURCE -50.0 0.0 9.975 1. 0. 0.
*SOURCE 2.0 0.0 8.0 0. 0. 1.
***********************************************
*
* Geometry description
*
***********************************************
GEOBEGIN 0. 0. 5.0 11.0 0. 0. COMBINAT
0 0 Cream
SPH World 0.0 0.0 0.0 +5000.0
RPP Exph -500.0 +500.0 -500.0 +500.0 -500.0 +500.0
RPP Calo -25.0 +25.0 -25.0 +25.0 -7.0 0.0
RPP Pixel0 -0.75 +0.75 -0.75 +0.75 9.965 10.0
RPP Pixel1 +0.75 +2.25 -0.75 +0.75 9.965 10.0
RPP Pixel2 +2.25 +3.75 -0.75 +0.75 9.965 10.0
RPP Pixel3 +3.75 +5.25 -0.75 +0.75 9.965 10.0
RPP Pixel4 +5.25 +6.75 -0.75 +0.75 9.965 10.0
RPP Pixel5 +6.75 +8.25 -0.75 +0.75 9.965 10.0
END
*
BKHL 5 +World -Exph
VACUUM 5 +Exph -Calo -Pixel0 -Pixel1 -Pixel2 -Pixel3 -Pixel4 -Pixel5
CALO 1 +Calo
SIPIX0 1 +Pixel0
SIPIX1 1 +Pixel1
SIPIX2 1 +Pixel2
SIPIX3 1 +Pixel3
SIPIX4 1 +Pixel4
SIPIX5 1 +Pixel5
END
*
LATTICE 5.0 9.0 1. 101. 105. 1.
GEOEND 0. 0. 0. 0. 0. 0.
***********************************************
*
* Material definition
*
MATERIAL 1.0 1.0079 0.0000899 3.0 0. 1.0 HYDROGEN
MATERIAL 6.0 12.011 2.265 6.0 0. 0. CARBON
MATERIAL 8.0 15.999 0.001429 8.0 0. 0. OXYGEN
MATERIAL 12.0 24.305 1.738 9.0 0. 0. MAGNESIU
MATERIAL 19.0 39.102 0.031165 11.0 0. 0. POTASSIU
MATERIAL 20.0 40.08 1.54 12.0 0. 0. CALCIUM
MATERIAL 11.0 22.990 0.038349 13.0 0. 0. SODIUM
MATERIAL 17.0 35.4529 0.0029947 14.0 0. 0. CHLORINE
MATERIAL 16.0 32.066 2.070 15.0 0. 0. SULFUR
MATERIAL 74.0 183.84 19.3 23.0 0. 0. TUNGSTEN
MATERIAL 14.0 28.085 2.33 26.0 0. 0. SILICON
*
EMFCUT 0.0001 0.0001 0. 3. 3. 1.0
EMFCUT 0.00001 0.00001 0. 4. 4. 1.0
*
* Here assigns materials to regions
ASSIGNMA 1.0 1.0 1.0 1.0 0. 0.
ASSIGNMA 2.0 2.0 2.0 1.0 0. 0.
ASSIGNMA 23.0 3.0 3.0 1.0 0. 0.
ASSIGNMA 26.0 4.0 4.0 1.0 0. 0.
*
* Setting SCORE and Ouptut levels.
SCORE 208.0 211.0 0. 0. 0. 0.
OUTLEVEL 1.0 7.0 0. 0. 0. 0.
EVENTDAT -23. 0.0 0. 0. 0. 0. EVT.SCOR
USERDUMP 101.0 49.0 0. 1.0 0. 0. CALOR
*
* Activates Random Number initialization
RANDOMIZE 1. 54217137 0. 0. 0. 0.
*
* Activates user dependent initialization
*
USRICALL 0. 0. 0. 0. 0. 0.
*
* Starts Job asking for a given number of events
*
START 5. 99999999. 80 0.0 0.0 100.
*
* Activates user dependent output
*
USROCALL 0. 0. 0. 0. 0. 0.
STOP 0. 0. 0. 0. 0. 0.
1 *====================================================================*
* *
* *
* *
* *
* FFFFF L U U K K AAA 222 000 000 333 *
* F L U U K K A A 2 2 0 0 0 0 3 3 *
* F L U U K K A A 2 0 0 0 0 3 *
* FFFF L U U KK AAAAA == 2 0 0 0 0 33 *
* F L U U K K A A 2 0 0 0 0 3 *
* F L U U K K A A 2 0 0 0 0 3 3 *
* F LLLLL UUU K K A A 22222 000 000 333 *
* *
* *
* *
* Version of Fluka2003 for GNU/Linux operating system *
* *
* By : *
* *
* Alfredo Ferrari & Paola Sala INFN Milan *
* *
* ALFREDO.FERRARI@MI.INFN.IT PAOLA.SALA@MI.INFN.IT *
* *
* *
* *
* This version includes all the features of the Fluka2003 *
* *
* package together with many further options, which, *
* *
* are not yet included in the Fluka93-2003 releases *
* *
* *
* *
* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *
* *
* !!! THIS VERSION IS THE MAIN DEVELOPMENT VERSION OF !!! *
* *
* !!! FLUKA AND IT IS FOLLOWED BY THE MILAN GROUP, IT !!! *
* *
* !!! IS A STRICTLY PRIVATE ONE. IT SHOULD CONTAIN THE !!! *
* *
* !!! BEST PHYSICS ... AND A CORRESPONDING LARGE NUM- !!! *
* *
* !!! BER OF BUGS ...... YOU ARE WARNED !!!! !!! *
* *
* !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *
* *
* *
* *
*====================================================================*
1
Draft of the FLUKA User license:
By downloading and/or using the FLUKA distribution files, you agree with the
following conditions:
- Title to the downloaded materials and all copies thereof remain with the
Italian National Institute for Nuclear Physics (INFN for brief),
the European Organization for Nuclear Research (CERN for brief)
and the FLUKA authors (A.Fasso', A.Ferrari, J.Ranft and P.R.Sala)
- INFN, CERN and the FLUKA authors are the only official source of
distribution of the code, bug fixes and documentation, as well as the unique
recipients of the code intellectual property.
- The FLUKA software can be dowloaded after explicit acceptance of these
conditions only from the dedicated INFN WEB server (http://www.fluka.org).
- The downloaded software cannot be redistributed in toto or in part.
- The downloaded software can be used free of charge by scientific and
non profit organisation. For commercial use ad hoc agreements have to be
negotiated.
- The FLUKA code or parts of it can be inserted into other codes only after
an official permission is granted by INFN, CERN, and the authors through
the INFN/CERN FLUKA steering committee.
The combined software, after the permission is obtained, can be used free
of charge by scientific and non profit organisations. For commercial use
ad hoc agreements have to be negotiated.
- The authors welcome user contribution to the code, that might be included in
the code, acknowledging and properly referencing the originators. Possible
contributions to the code must be sent separately and documented to the
FLUKA authors, who will be the sole responsible for accepting them for
inclusion in the distributed code.
- The use of the Fluka code must be acknowledged explicitly by quoting
the following set of references:
- A.Fasso', A.Ferrari, P.R.Sala, "Electron-photon transport in
FLUKA: status", Proceedings of the MonteCarlo 2000 Conference,
Lisbon, October 23-26 2000, A.Kling, F.Barao, M.Nakagawa,
L.Tavora, P.Vaz - eds., Springer-Verlag Berlin, p.159-164 (2001).
- A.Fasso', A.Ferrari, J.Ranft, P.R.Sala, "FLUKA: Status and
Prospective for Hadronic Applications", Proceedings of
the MonteCarlo 2000 Conference, Lisbon, October 23-26 2000,
A.Kling, F.Barao, M.Nakagawa, L.Tavora, P.Vaz - eds. ,
Springer-Verlag Berlin, p.955-960 (2001).
For questions/clarifications/problems about the license please
contact the chairman of the joint INFN/CERN FLUKA steering committee:
Prof. Giuseppe Battistoni
INFN - Via Celoria 16 I-20133 Milano (Italy)
Giuseppe.Battistoni@mi.infn.it
1
DISCLAIMER
THIS SOFTWARE IS PROVIDED BY INFN, CERN, THE FLUKA AUTHORS, AND CONTRIBUTORS
"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL INFN, CERN, THE FLUKA AUTHORS, OR
CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;
OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF
ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
1 ****************************************************************************************************
FLUKA2003 Version 1.0b Mar-04 by A. Ferrari DATE: 5/12/ 4 TIME: 10:27: 8
****************************************************************************************************
*****NEXT CONTROL CARD ***** TITLE 0.000 0.000 0.000 0.000 0.000 0.000
Test of simple calorimeter
*****NEXT CONTROL CARD ***** GLOBAL 3000. 0.000 0.000 2.000 1.000 0.000
*-----------------------------------------
--- Free-Format Input -------------------------------------------------*
*****NEXT CONTROL CARD ***** DEFAULTS 0.000 0.000 0.000 0.000 0.000 0.000 CALORIME
*---------------- ----------------*
*---------------- ----------------*
*---------------- BEAM 20. 0. 0. 0. 0. 1.0 RAY ----------------*
*****NEXT CONTROL CARD ***** BEAM 20.00 0.000 0.000 0.000 0.000 1.000 ELECTRON
*---------------- BEAM 20. 0. 0. 0. 0. 1.0 PROTON ----------------*
*---------------- BEAM 20. 0. 0. 0. 0. 1.0 HEAVYION ----------------*
*---------------- Heavy ion case: BEAM WHAT(1) is kin energy (momentum) per nucleon NOT total ----------------*
*---------------- HI-PROPE 5.0 10.0 0.0 ----------------*
*---------------- EVENTYPE 0. 0. 2. 0. 0. 0. DPMJET ----------------*
*---------------- ----------------*
*---------------- Particles must injected in vacuum or material not in black hole ----------------*
*---------------- ----------------*
*****NEXT CONTROL CARD ***** BEAMPOS 0.000 0.000 0.000 0.000 0.000 0.000
*****NEXT CONTROL CARD ***** SOURCE 0.000 0.000 30.00 0.000 0.000 -1.000
*---------------- SOURCE -50.0 0.0 9.975 1. 0. 0. ----------------*
*---------------- SOURCE 2.0 0.0 8.0 0. 0. 1. ----------------*
*---------------- ********************************************** ----------------*
*---------------- ----------------*
*---------------- Geometry description ----------------*
*---------------- ----------------*
*---------------- ********************************************** ----------------*
*****NEXT CONTROL CARD ***** GEOBEGIN 0.000 0.000 5.000 11.00 0.000 0.000 COMBINAT
1
Cream
IVOPT = 0 IDBG = 0
BODY DATA
SPH World 1 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.50000000E+04 0.00000000E+00 0.00000000E+00 5
RPP Exph 2 -0.50000000E+03 0.50000000E+03 -0.50000000E+03 0.50000000E+03 -0.50000000E+03 0.50000000E+03 15
RPP Calo 3 -0.25000000E+02 0.25000000E+02 -0.25000000E+02 0.25000000E+02 -0.70000000E+01 0.00000000E+00 25
RPP Pixel0 4 -0.75000000E+00 0.75000000E+00 -0.75000000E+00 0.75000000E+00 0.99650000E+01 0.10000000E+02 35
RPP Pixel1 5 0.75000000E+00 0.22500000E+01 -0.75000000E+00 0.75000000E+00 0.99650000E+01 0.10000000E+02 45
RPP Pixel2 6 0.22500000E+01 0.37500000E+01 -0.75000000E+00 0.75000000E+00 0.99650000E+01 0.10000000E+02 55
RPP Pixel3 7 0.37500000E+01 0.52500000E+01 -0.75000000E+00 0.75000000E+00 0.99650000E+01 0.10000000E+02 65
RPP Pixel4 8 0.52500000E+01 0.67500000E+01 -0.75000000E+00 0.75000000E+00 0.99650000E+01 0.10000000E+02 75
RPP Pixel5 9 0.67500000E+01 0.82500000E+01 -0.75000000E+00 0.75000000E+00 0.99650000E+01 0.10000000E+02 85
END 10 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 95
NUMBER OF BODIES 9
LENGTH OF FPD-ARRAY 100
INPUT REGION DATA
*
*--------------- ---------------*
*
BKHL 5 + World - Exph
VACUUM 5 + Exph - Calo - Pixel0 -
Pixel1 - Pixel2 - Pixel3 - Pixel4 - Pixel5
CALO 1 + Calo
SIPIX0 1 + Pixel0
SIPIX1 1 + Pixel1
SIPIX2 1 + Pixel2
SIPIX3 1 + Pixel3
SIPIX4 1 + Pixel4
SIPIX5 1 + Pixel5
END
NUMBER OF INPUT REGIONS 9
NUMBER OF CODE ZONES 9
LENGTH OF INTEGER ARRAY 175
CODE ZONE INPUT REGION ZONE DATA LOC. NO. OF BODIES REGION NO.
1 1 64 2 1
2 2 73 8 2
3 3 106 1 3
4 4 111 1 4
5 5 116 1 5
6 6 121 1 6
7 7 126 1 7
8 8 131 1 8
9 9 136 1 9
I IR1(I) IR2(I)
1 1 1
2 2 2
3 3 3
4 4 4
5 5 5
6 6 6
7 7 7
8 8 8
9 9 9
Interpreted region echo
Region n. 1
1 -2
Region n. 2
2 -3 -4 -5 -6 -7 -8
-9
Region n. 3
3
Region n. 4
4
Region n. 5
5
Region n. 6
6
Region n. 7
7
Region n. 8
8
Region n. 9
9
1 OPTION 0 WAS USED IN CALCULATING VOLUMES, FOR 9 REGIONS
3: INPUT VOLUMES, ANYTHING ELSE: VOLUMES = 1.0
VOLUMES (CM**3)
1 REG 1 2 3 4 5 6 7 8 9
VOLUME 1.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00 1.000E+00
*---------------- ----------------*
*** Geometry lattice capabilities activated ***
*** Geometry lattice capabilities activated ***
Abort called from LATINP reason IERR>0 Run stopped!
STOP IERR>0
This archive was generated by hypermail 2.1.6 : Wed May 12 2004 - 14:43:42 CEST