[fluka-discuss]: RE: Primary protons absorbed in a nuclear interaction vs elastic large angle scattered

From: Anton Lechner <Anton.Lechner_at_cern.ch>
Date: Thu, 8 Oct 2015 19:43:50 +0000

Dear Ana,

You could determine the number of inelastic collisions of primary protons in the relevant regions. One way to do that is to score with USRBIN (per region scoring) the number of "stars" of primary protons (using BEAMPART).

On the other hand, one could also count manually the collisions in the USDRAW entry of the MGDRAW routine, which is called after each interaction. This is maybe more bothersome than the above approach, but allows to access additional information (like the position, secondaries produced etc.). With following code snippet you could filter out the events which are inelastic collisions of primaries (using the generation number LTRACK, which is 1 for primaries and making use of the ICODE variable, which is 101 for inelastic nuclear collisions):

* Inelastic collision of a primary particle
      IF ( LTRACK .EQ. 1 .AND. ICODE .EQ. 101 ) THEN

* Writing position of inelastic collisions
            WRITE( 98,* ) "Inelastic coll at", XSCO, YSCO, ZSCO
      END IF

Cheers, Anton

From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Ferreira De Almeida Lourenco, Ana [am.lourenco_at_ucl.ac.uk]
Sent: 07 October 2015 17:57
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: Primary protons absorbed in a nuclear interaction vs elastic large angle scattered

Dear Fluka experts,

I am studying the number of primary protons that are removed from a proton beam. For that, I would like to distinguish between the protons that are absorbed in a nuclear interaction from the ones that are elastic large angle scattered. How can I get this information from FLUKA?

Many thanks,


You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Oct 08 2015 - 23:26:21 CEST

This archive was generated by hypermail 2.3.0 : Thu Oct 08 2015 - 23:26:26 CEST