flyopunta/src/src_par/partmunu.cpp

/***************************************************************************
                          partmunu.cpp  -  description
                             -------------------
    begin                : Sun May 9 2004
    copyright            : (C) 2004 by gmp
    email                : peppe@newpeppe
  ***************************************************************************
  *                                                                         *
  *   NA48  simulation program.                                             *
  *                                                                         *
  ***************************************************************************/
//--------------------------------
//            P a r t M u n u    for    K->3l decay
//----------------------------------------------

#include "flyoh.h"

using namespace GmpConst;
void PartMunu::
 DecParm()
 {
// Attenzione questa routine vale solo per
// il decadimento semileptonico dei K->l3
// va scritta tutta.....

   P[3] = P.m=massa =  0.126+ Pran()*0.239;   //da definire


// calcolo la direzione di decadimento dei figli.
   double cost = 1.0 -  2.0 * Pran();

//mettere la distribuzione in cost giusta.. p(costet)= (1.+costet)/2.

   double   fi = DuePiGreco * Pran();
   dw->cstar=cost;     dw->fistar=fi;
   dw->rg->cstar=-cost; dw->rg->fistar=fi+PiGreco;
 }