TWiki> NA62 Web>KDecayPions (revision 6)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

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 | r8 < r7 < r6 < r5 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r6 - 2015-04-20 - DanProtopopescu
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2021 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback