flyopunta/src/src_par/part3pic.cpp

/***************************************************************************
                          part3pic.cpp  -  description
                             -------------------
    begin                : Sun May 9 2004
    copyright            : (C) 2004 by gmp
    email                : peppe@newpeppe
  ***************************************************************************
  *                                                                         *
  *   NA48  simulation program.                                             *
  *                                                                         *
  ***************************************************************************/

#define MPIC   0.13956
#define MPICQ  0.019477
#define MPIO   0.13497
#define MPIOQ  0.0182169
 
#include "flyoh.h"

using namespace std;
//-------------------------------------------------------
//        K -> pic pic2 -> 3 pic  D e c a y
//-------------------------------------------------------

// constructor
    Part3pic::Part3pic(double pant)
    { 
      // si usa questa procedura se il decadimento e'mediato dalla
      // risonanza virtuale pic2 che va esplicitamente scritta nella catenata
      // della reazione in Epc.
      mpadre=0.;
      slope=pant;
    }

//--------------------------------------------------------------*/
//       D e c P a r m   ==>   vedi PartFase::DecParm()         */
//--------------------------------------------------------------*/
//------------------------------------------
 void Part3pic::Get_PhaseParm()
{
//-------------------------------------------
// define the particle links.....
   PartFase::Get_PhaseParm();
   if(p2->Get_Nome()!=p3->Get_Nome())
    { Gout<<"\n Particle error in 3 pions decay; pic2 -> "
          << p2->Get_Nome()<<" "<<p3->Get_Nome()
          <<"   not equal ..  malformed.... exit!"<<std::endl;
      exit(0);
    }
 }
//---------------
 void  Part3pic::Get_MatrixParm()
{
  // matrix element max weight  guess
  // for K+ and K- ...but you may define different parameters
  // if(strncmp(dw->Get_Nome(),"Kp,2)==0) {....} else Km {}..
  gslope= -0.2154;
  hcurv=   0.012;
  kcurv= -0.0101;
  Gout<<"\nPic2 done:  Slope= "<<setprecision(3)<<setw(8)<< slope
      <<" g,h,k= "<<setw(8)<<gslope
      <<" "<<setw(8)<<hcurv
      <<" "<<setw(8)<<kcurv;
  wmtr_max=1.+gslope*1.2+ hcurv*1.44; 
}

//---------------
 int  Part3pic::Get_MatrixValue()
{
 double su=(mresq-s0)/MPICQ;
 double sv=(2.*mpadre*(e2-e1))/MPICQ;  //ricorda le energie sono e0,e1,e2
 prb=1.+gslope*su+hcurv*su*su + kcurv*sv*sv;
//  Gout<<"\nPart3 matrix ok prb "<<setprecision(4)<<prb;
 if(prb>wmtr_max)
    {
       Gout<<"\n Wmtr off: Ev "
        <<setprecision(4)<<fixed
        <<setw(6)<<evento_.Gen.Event<<" iter="
        <<setw(4)<<itermtr<<" pb="
        <<setw(7)<<prb<<" pbx="
        <<setw(7)<<wmtr_max <<" e0="
        <<setprecision(2)<<fixed
        <<setw(6)<<e0<<" e1="
        <<setw(6)<<e1<<" "
        <<setw(25)<<SelRea->Get_Titolo()<<std::endl;  
      wmtr_max=prb*1.001; // aggiorno ma informo.... dovrebbe  accadere poche volte!
    }

 if( (wmtr_max*Pran()) <prb) return 1;  else return 0;
}