ReaKld Class Reference

#include <reakld.h>

Inheritance diagram for ReaKld:

List of all members.

Public Member Functions
virtual double	Produci ()
Protected Attributes
double	ak [5]

---

Detailed Description

decay Kl beam production

Author:

GmP

Definition at line 17 of file reakld.h.

---

Member Function Documentation

double ReaKld::Produci

(

)

[virtual]

Produzione da proton di Kl

Reimplemented from Reazione.

Definition at line 180 of file reaprod.cpp.

{

  gvet vtet= Vbeam&Vinc_ppp; vtet.Norma();
  double  tetq= asin ( vtet.norma );tetq*=tetq;



  //  old Niels code...end
  //  double A=0.16, B=8.5, C=3.0;
  /*
  double ppbs=pow(ppb,2.);
  double A=0.16, B=10.5, C=3.0;
  double c1= ( B/inc_pinc ) *exp ( -B*ppb/inc_pinc );
  double c2= ( C*ppbs/PiGreco ) *exp ( -C*ppb*ppb*tetq );
  double Kpfluss= A*c1*c2 ;

  // A=0.1, B=13., C=3.5;
  A=0.1, B=13., C=3.5;Inizio

  // spettro*=1. + 8.8*exp ( -.05*ppb );
  // spettro*=exp ( -5.*(ppb-60.)/inc_pinc  );

  c1= ( B/inc_pinc ) *exp ( -B*ppb/inc_pinc );
  c2= ( C*ppbs/PiGreco ) *exp ( -C*ppb*ppb*tetq );
  double Kmfluss= A*c1*c2 ;
  double spettro=(Kpfluss+3.*Kmfluss)*.25;

  */


  //  double C=3.5;


  // italo..

  //double slope =0.032840722 == 1./30.45;   org power ===  1.671;
  double slope =1./30.45;
  double power=1.671;          //org power ===  1.671;
  double spettro =  pow ( ppb,power ) *exp ( -slope*ppb );

  //correzione per cammino....
  double Inizio,Fine;
  Inizio=  Beg_fid - Ctarg.z;
  Fine=    End_fid - Ctarg.z;
  double dzdec=End_fid-Beg_fid;
  double camvol = ctau_avo_su_mass * ppb ;
  double  pathcor= exp ( -Inizio/camvol ) * ( 1.-exp ( -dzdec/camvol ) ) ;
  spettro*=pathcor;

  //  spettro*= ( 0.6325-0.2324e-2*ppb+0.2244e-4*ppb*ppb );  //standard...
  //  spettro*=(0.6325-0.2324e-2*ppb+0.2144e-4*ppb*ppb);


  // spettro*=(0.6325-0.1900e-2*ppb+0.2450e-4*ppb*ppb);  // quasi perfetto
  spettro*= ( 0.6325-0.1900e-2*ppb+0.27e-4*ppb*ppb );  // forse più perfetto
  return spettro;
}

---

The documentation for this class was generated from the following files:

• flyopunta/src/src_reaz/reakld.h
• flyopunta/src/src_reaz/reaprod.cpp