flyopunta/src/src_par/partpvv.cpp

/***************************************************************************
 *   Copyright (C) 2005 by Giuseppe Pierazzini          Pisa 20.02.05      *
 *   pierazzini@unipi.it                                                   *
 *                                     F l y o                             *
 *   Epsi/NA48                                                             *
 ***************************************************************************/


#include "flyoh.h"

using namespace std;
//-------------------------------------------------------
// constructor
PartPvv::PartPvv ( )
{
  // si usa questa procedura se il decadimento e'mediato dalla
  // risonanza virtuale pvv che va esplicitamente scritta nella catenata
  // della reazione in Epc.
  mas_pip=ptr_Pip->Get_Massa();
  mas_kp= ptr_Kapp->Get_Massa();
  masq_kp=mas_kp*mas_kp;
  masq_pip=mas_pip*mas_pip;
  Gout<<"\n PartPvv done ";


}

//--------------------------------------------------------------*/
//       D e c P a r m   ==>   vedi PartFase::DecParm()         */
//--------------------------------------------------------------*/
//------------------------------------------
void PartPvv::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  PartPvv::Get_MatrixParm()
{
  // calcolare qui il massimo della matrice
  wmtr_max=0.013;
  Gout<<"\nPvv  done:   prob_max "<<setprecision ( 3 ) <<setw ( 8 ) << wmtr_max;
}

//

int  PartPvv::Get_MatrixValue()
{
  // Compute weight
  double Eprime,a,q2,v;
  Eprime = ( masq_kp +masq_pip ) / ( 2.*mas_kp )-e0;
  a      = mas_kp* ( 2.*e1*e2-mas_kp*Eprime );
  q2     = masq_kp +masq_pip-2.*mas_kp*e0;
  v      = 1.+0.0286*q2/ ( masq_pip );
  prb = a*pow ( v,2 );

  // Check against wmtr_max
  if ( prb>wmtr_max )
  {
    Gout<<"\n Wmtr off: Ev "
    <<setprecision ( 4 ) <<fixed
    <<setw ( 6 ) <<evento_.Gen.Event
    <<" 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;
  }
  if ( ( wmtr_max*Pran() ) <prb ) return 1;
  else return 0;
}