Hi
my fault, I did not look carefully to the cards, and the manual is not
clear:
When using SPECSOUR, one should input the TOTAL momentum of the ion, not
the momentum/nucleon.
We'll correct the manual accordingly.
Yes, fragments are in the same stack, with negative particle ID coding
Z*100000+A*100 +internal index
Paola
On 11/15/2017 10:07 AM, Mauro Gallio wrote:
> Dear Paola,
> as you can see from the output in attachment, I got the message that
> RQMD and DPMJET is loaded.
> From the output usrein.out I can see (like you wrote) that sometimes a
> lot of particles is produced in the collision, but:
> - I do not see any neutron with a momentum value of 2760 GeV/c or
> something like that (the so called spectator neutrons)
> I tried 100 events and all the forward neutrons have around 13 GeV/c
> which it seems to me a bit strange.
> - Another question: are the nuclear fragments produced in the collision
> written in the same stack?
> Thanks again
> Mauro
>
>
>
> Mauro Gallio
> Dip. di Fisica
> Via P. Giuria,1
> 10125 Torino (Italy)
> phone : +39 011 6707361
>
>
> On Wed, 15 Nov 2017, paola sala wrote:
>
>> Hello
>>
>> could you, please, check again? look in the .out file, you should
>> find messages like
>>
>> "requested particle production model RQMD"
>>
>> and
>>
>> "requested particle production model DPMJET"
>>
>> otherwise the executable does not include the heavy ion generators.
>>
>> (**** always look in the output file*******)
>>
>> Note that the executable produced by ldpmqmd is called flukadpm3, not
>> flukahp,
>>
>> one has to specify it in the rfluka command line
>>
>> rfluka -e flukadpm3
>>
>>
>> running your input file without modifications, I get, for instance
>>
>> *** Event no. 1 fluka particles in stack = 322
>>
>> or
>>
>> *** Event no. 3 fluka particles in stack = 16
>> x,y,z collision vertex = 0.00000 0.00000 0.00000
>> id = -8220807 E,ux,uy,uz,T 2613.94169 0.00000 0.00000
>> -1.00000 2436.93212
>> id = 8 E,ux,uy,uz,T 15.64084 -0.01480
>> -0.01393 0.99979 14.72947
>> id = 8 E,ux,uy,uz,T 13.56707 -0.00653
>> -0.00071 0.99998 12.66000
>> id = 8 E,ux,uy,uz,T 12.60993 -0.00053
>> 0.00417 0.99999 11.70532
>>
>> Paola
>>
>> On 11/14/2017 10:19 AM, Mauro Gallio wrote:
>>> Dear Paola,
>>> yes, I created the executable with the following command:
>>> $FLUPRO/flutil/ldpmqmd -o flukamy usrein.o
>>> and I have seen that there is no difference betwwen the two scripts:
>>> ldpmqmd and ldpm3qmd
>>> (diff ldpmqmd ldpm3qmd ==> no difference)
>>> Are there errors in the input data cards (may be concerning units)?
>>> Even if e.m. dissociation only was activated, I would expect
>>> neutrons with greater energies !
>>> Thanks
>>> Mauro
>>>
>>>
>>> Mauro Gallio
>>> Dip. di Fisica
>>> Via P. Giuria,1
>>> 10125 Torino (Italy)
>>> phone : +39 011 6707361
>>>
>>>
>>> On Tue, 14 Nov 2017, Paola Sala wrote:
>>>
>>>> Dera Mauro,
>>>> are you sure that you linked the heavy ions event generators
>>>> (ldpm3qmd)?
>>>> Your output events are too simple to be leed-lead collisions at 7TeV!
>>>> looks like they only perform EM dissociation, no inelastic
>>>> interactions.
>>>> Paola
>>>>>
>>>>> Dear Fluka experts,
>>>>> following an old discussions concerning the production of spectator
>>>>> neutrons in heavy ions collisions, I wrote a simple user routine
>>>>> usrein.f
>>>>> (in attachment together with the input file) which should write the
>>>>> kinematics properties of the particles generated in the collision.
>>>>>
>>>>> Looking to the output from usrein.f, I see that there are no
>>>>> neutrons emitted in the forward direction with a momentum value close
>>>>> to that of the nucleons of the colliding nuclei).
>>>>> Experimentally we see them with a Zero Degree Calorimeter.
>>>>>> From the output I can see the two residual nuclei and also some
>>>>>> neutrons
>>>>> but with a momentum value around 13 GeV/c, no neutrons with momentum
>>>>> around 2760 GeV/c, neither nuclear fragments in the case of hadronic
>>>>> interaction (see in the last attachment a sample of the output from
>>>>> usrein.f).
>>>>> What could be the reason?
>>>>> Thanks in advance
>>>>> Mauro
>>>>>
>>>>> Mauro Gallio
>>>>> Dip. di Fisica
>>>>> Via P. Giuria,1
>>>>> 10125 Torino (Italy)
>>>>> phone : +39 011 6707361
>>>>>
>>>>>
>>>>> On Tue, 7 Nov 2017, Francesco Cerutti wrote:
>>>>>
>>>>>>
>>>>>> Hello Mauro
>>>>>>
>>>>>> continuation cards must still start with the card name (SPECSOUR
>>>>>> in this
>>>>>> case)
>>>>>>
>>>>>> Best
>>>>>>
>>>>>> Francesco
>>>>>>
>>>>>> **************************************************
>>>>>> Francesco Cerutti
>>>>>> CERN-EN/STI
>>>>>> CH-1211 Geneva 23
>>>>>> Switzerland
>>>>>> tel. +41 22 7678962
>>>>>> fax +41 22 7668854
>>>>>>
>>>>>> On Tue, 7 Nov 2017, Mauro Gallio wrote:
>>>>>>
>>>>>>>
>>>>>>> I forgot the attachment!!
>>>>>>> Sorry
>>>>>>> Mauro
>>>>>>>
>>>>>>>
>>>>>>> Mauro Gallio
>>>>>>> Dip. di Fisica
>>>>>>> Via P. Giuria,1
>>>>>>> 10125 Torino (Italy)
>>>>>>> phone : +39 011 6707361
>>>>>>>
>>>>>>>
>>>>>>> ---------- Forwarded message ----------
>>>>>>> Date: Tue, 7 Nov 2017 12:47:26 +0100 (CET)
>>>>>>> From: Mauro Gallio <gallio_at_to.infn.it>
>>>>>>> To: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
>>>>>>> Cc: Alfredo Ferrari <alfredo.ferrari_at_cern.ch>,
>>>>>>> Trinh Ngoc-Duy <Ngoc-Duy.Trinh_at_ganil.fr>,
>>>>>>> "fluka-discuss_at_fluka.org" <fluka-discuss_at_fluka.org>
>>>>>>> Subject: RE: [fluka-discuss]: Error OBSOLETE CODEWORD trying to
>>>>>>> simulate
>>>>>>> Pb-Pb
>>>>>>> at LHC
>>>>>>>
>>>>>>>
>>>>>>> Dear Fluka expert,
>>>>>>> Following the suggestion of Alfredo Ferrari,I tried to insert a
>>>>>>> very
>>>>>>> simple
>>>>>>> dummy geometry following the example exfixed.inp (which runs
>>>>>>> without
>>>>>>> errors). The input file is in attachment.
>>>>>>> I'm interested essentially in the collision physics, for the
>>>>>>> moment.
>>>>>>>
>>>>>>> The program still crashes with the same error OBSOLETE CODEWORD.
>>>>>>> The command which runs the progam and the corresponding output
>>>>>>> follows:
>>>>>>> [gallio_at_alz117xl fluka]$ ./rfluka -e flukadpm3 -M1 gene0
>>>>>>> $TARGET_MACHINE = Linux
>>>>>>> $FLUPRO = /usr/local/fluka
>>>>>>>
>>>>>>> Initial seed copied from /usr/local/fluka
>>>>>>> Running fluka in /home/gallio/fluka/fluka_24860
>>>>>>>
>>>>>>> ======================= Running FLUKA for cycle # 1
>>>>>>> =======================
>>>>>>> ./rfluka: line 359: 24887 Aborted (core dumped)
>>>>>>> "${EXE}" <
>>>>>>> "$INPN" 2> "$LOGF" > "$LOGF"
>>>>>>>
>>>>>>> Any hint about the error?
>>>>>>> Thanks
>>>>>>> Mauro
>>>>>>>
>>>>>>>
>>>>>>> Mauro Gallio
>>>>>>> Dip. di Fisica
>>>>>>> Via P. Giuria,1
>>>>>>> 10125 Torino (Italy)
>>>>>>> phone : +39 011 6707361
>>>>>>>
>>>>>>>
>>>>>>> On Fri, 3 Nov 2017, Francesco Cerutti wrote:
>>>>>>>
>>>>>>>>
>>>>>>>> Note in addition that, despite the dedicated PHYSICS card, you
>>>>>>>> did not
>>>>>>>> activate electromagnetic-dissociation, to that purpose you need
>>>>>>>> WHAT(1)=2.0
>>>>>>>>
>>>>>>>> Best
>>>>>>>>
>>>>>>>> Francesco
>>>>>>>>
>>>>>>>> **************************************************
>>>>>>>> Francesco Cerutti
>>>>>>>> CERN-EN/STI
>>>>>>>> CH-1211 Geneva 23
>>>>>>>> Switzerland
>>>>>>>> tel. +41 22 7678962
>>>>>>>> fax +41 22 7668854
>>>>>>>>
>>>>>>>> On Fri, 3 Nov 2017, Alfredo Ferrari wrote:
>>>>>>>>
>>>>>>>>> ... you must have a valid geometry in your input, even a dummy
>>>>>>>>> one.
>>>>>>>>> If not the code revert to the 30 year old cylindrical geometry
>>>>>>>>> which
>>>>>>>>> do no longer exists and prints that message.
>>>>>>>>>
>>>>>>>>> Alfredo
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> +----------------------------------------------------------------------+
>>>>>>>>> | Alfredo Ferrari || Tel.: +41.22.76.76119
>>>>>>>>> |
>>>>>>>>> | CERN-EN/STI || Fax.: +41.22.76.69474
>>>>>>>>> |
>>>>>>>>> | 1211 Geneva 23 || e-mail:
>>>>>>>>> Alfredo.Ferrari_at_cern.ch
>>>>>>>>> |
>>>>>>>>> | Switzerland ||
>>>>>>>>> |
>>>>>>>>> +----------------------------------------------------------------------+
>>>>>>>>>
>>>>>>>>> On Fri, 3 Nov 2017, Mauro Gallio wrote:
>>>>>>>>>
>>>>>>>>>> Dear Trinh,
>>>>>>>>>> unfortunately I downloaded exactly the last version from the
>>>>>>>>>> official
>>>>>>>>>> site:
>>>>>>>>>>
>>>>>>>>>> fluka2011.2c-linuxAA.tar.gz
>>>>>>>>>> FLUKA Package Version 2011.2c-6 5th of June, 2017
>>>>>>>>>> 147 MB
>>>>>>>>>>
>>>>>>>>>> Is there an other repository for the last version?
>>>>>>>>>> Best regards
>>>>>>>>>> Mauro
>>>>>>>>>>
>>>>>>>>>> Mauro Gallio
>>>>>>>>>> Dip. di Fisica
>>>>>>>>>> Via P. Giuria,1
>>>>>>>>>> 10125 Torino (Italy)
>>>>>>>>>> phone : +39 011 6707361
>>>>>>>>>>
>>>>>>>>>>
>>>>>>>>>> On Fri, 3 Nov 2017, Trinh Ngoc-Duy wrote:
>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>> Dear Mauro,
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> My first guess is that your Fluka version is out of date.
>>>>>>>>>>> Please
>>>>>>>>>>> try
>>>>>>>>>>> the
>>>>>>>>>>> latest version Fluka2011.2c.6 (June 2017).
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Your faithfully.
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Duy, GANIL.
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> -----Message d'origine-----
>>>>>>>>>>> De : owner-fluka-discuss_at_mi.infn.it
>>>>>>>>>>> [mailto:owner-fluka-discuss_at_mi.infn.it]
>>>>>>>>>>> De la part de Mauro Gallio
>>>>>>>>>>> Envoyé : Friday, November 03, 2017 2:57 PM
>>>>>>>>>>> � : fluka-discuss_at_fluka.org
>>>>>>>>>>> Objet : [fluka-discuss]: Error OBSOLETE CODEWORD trying to
>>>>>>>>>>> simulate
>>>>>>>>>>> Pb-Pb at
>>>>>>>>>>> LHC
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Dear Fluka experts,
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> I am interested in relativistic heavy ions collisions at
>>>>>>>>>>> LHC. In
>>>>>>>>>>> particular
>>>>>>>>>>> I would like to simulate the formation of nuclear fragments
>>>>>>>>>>> and the
>>>>>>>>>>> number
>>>>>>>>>>> of free spectator neutrons and protons.
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> I downloaded the fluka version
>>>>>>>>>>>
>>>>>>>>>>> -Â Fluka2011.2c (October 2014) released on: October 15th
>>>>>>>>>>> 2014 and
>>>>>>>>>>> I
>>>>>>>>>>> installed fluka on a Linux Desktop with SLC4 (gcc and g77
>>>>>>>>>>> version
>>>>>>>>>>>
>>>>>>>>>>> 3.4.6)
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> I tried to run the executable flukadpm3 with a simple .inp data
>>>>>>>>>>> cards
>>>>>>>>>>> in
>>>>>>>>>>> attachment, but I got the following message:
>>>>>>>>>>>
>>>>>>>>>>> Â Â ******** !!!!!!! OBSOLETE CODEWORDÂ Â Â !!!!!!! ********
>>>>>>>>>>>
>>>>>>>>>>> (see the output file in the second attachment).
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> I did not implement any geometry, since I am interested for the
>>>>>>>>>>> moment
>>>>>>>>>>> only
>>>>>>>>>>> in the collision tape.
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Thanks in advance for your help
>>>>>>>>>>>
>>>>>>>>>>> Best regards
>>>>>>>>>>>
>>>>>>>>>>> Â Â Â Â Â Â Â Â Â Â Â Mauro
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Â
>>>>>>>>>>>
>>>>>>>>>>> Mauro Gallio
>>>>>>>>>>>
>>>>>>>>>>> Dip. di Fisica
>>>>>>>>>>>
>>>>>>>>>>> Via P. Giuria,1
>>>>>>>>>>>
>>>>>>>>>>> 10125 Torino (Italy)
>>>>>>>>>>>
>>>>>>>>>>> phone : +39 011 6707361
>>>>>>>>>>>
>>>>>>>>>>> ____________________________________________________________________________
>>>>>>>>>>> Préservons notre environnement, nâ??imprimez ce mail que si
>>>>>>>>>>> nécessaire.
>>>>>>>>>>> Preserve our environment, print this email only if necessary.
>>>>>>>>>>>
>>>>>>>>>>>
>>>>>>>>>>
>>>>>>
>>>>
>>>>
>>>> Paola Sala
>>>> INFN Milano
>>>> tel. Milano +39-0250317374
>>>> tel. CERN +41-227679148
>>>>
>>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Wed Nov 15 2017 - 12:11:37 CET