Dear Paola,
Thank you very much!!! Now it is OK.
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:
> 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:51:39 CET