TWiki> NA62 Web>KDecayPions (revision 7)EditAttach

Pions from in-beam Kaon Decay

Pions from in-beam Kaon Decay


This case is handled via pichdecays.F. This file is a container for charged pion decays.
The routine pichdecay_manager is called in kaon decay generators where one ore more charged pions are present in the final state (e.g. kch2pipi0.F).
Two pion decay types are implemented in pich2lnu.F and pich2lnug.F, with the possibility to choose the lepton type (electron or muon).

Software implementation

New generator files:

  • Generator/pichdecays.F
  • Generator/pich2lnu.F
  • Generator/pich2lnug.F
Generator to be modified:
  • Generator/kch2pipi0.F

Macro file to be modified:

  • macros/StandardRun.mac

C++ files to be modified:

  • src/, include/CMC.hh
  • src/, include/DatacardManager.hh
  • src/, include/DatacardMessenger.hh

How to run

To run this channel:

  • specify the charged pion decay type in the StandardRun.mac file, at the line with the tag /decay/pichDecay.
    Available charged pion decay types:
    • 0: pi is passed to Geant4;
    • 1: pi-->mu nu;
    • 2: pi-->e nu;
    • 3: pi-->mu nu g;
    • 4: pi-->e nu g;

Run the software giving the modified macro file as input parameter, e.g. on an SL6 machine:

cd /afs/
source scripts/
NA62MC macros/StandardRun.mac

In the StandardRun.mac file, one has to pay attention to the following instructions:

/decay/pichDecay 2 
as explained above. Then use

/output/fileName pluto.root
to set where the output goes.

Verbosity can be set with

/run/verbose 1
/event/verbose 1
/tracking/verbose 1

Temporary patch

  • The generation of the K+ -> pi+pi0, pi+ -> mu+ nu, pi+0 -> gg process is handled in fortran files independent from the NA62MC software. These produce a kdecayggl.txt file containing informations on three final state particles: mu+, gamma, gamma.

  • In src/ a new case has been added, in which the kdecayggl.txt file is read and the three input particles are inserted into the NA62 apparatus.
    In the macros/StandardRun.mac file this option can be selected with /beam/SetBeam PionDecayTest.
    The handling of the new parameter has been implemented in src/

How to check the output

An event analyzer is provided inside the NA62Analysis framework.

In /afs/ the "MCTruthAnalyzer" can be used to debug the MC events produced with NA62MC.

  • The code is in Analyzers/src/

After sourcing scripts/ from within UserAnalysis/

  • Compile with configFiles/MCTruthAnalyzer.conf
  • Run with ./MCTruthAnalyzer -i ../NA62MC/pluto.root -o test.root -v kUser --config configFiles/
    The kUser verbose flag has to be used to get on screen information on generated particles.


This will be documented in technical note CERN-

Edit | Attach | Watch | Print version | History: r10 < r9 < r8 < r7 < r6 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r7 - 2015-06-29 - BrunoAngelucci
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback