flyopunta/src/src_par/part3pi0.cpp

/***************************************************************************
                          part3pi0.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;

// constructor
    Part3pi0::Part3pi0(double pant)
    { // si usa questa procedura se il decadimento e'mediato dalla
      // risonanza virtuale pi02 che va esplicitamente scritta nella catenata
      // della reazione in Epc.
      mpadre=0.;
      slope=pant;    // definisce la slope
      Gout<<"\n Pi02 done --- Slope ="<<slope<<std::endl;
    };
//--------------------------------------------------------------*/
//       D e c P a r m   ==>   vedi PartFase::DecParm()         */
//--------------------------------------------------------------*/
//------------------------------
 void Part3pi0::Get_PhaseParm()
{
//-------------------------------------------
// define the particle links.....
   PartFase::Get_PhaseParm();
   if(p2->Get_Nome()!=p3->Get_Nome())
    { Gout<<"\n Particle error in 3 pions decay; pi02 -> "
          << p2->Get_Nome()<<" "<<p3->Get_Nome()
          <<"  not equal ...malformed.... exit!"<<std::endl;
      exit(0);
    }
 }
//---------------
 void  Part3pi0::Get_MatrixParm()
{
// matrix element max weight  guess
    if(slope >0.0) { wfase_max=1.+slope*2.63;} else wfase_max=1.;
    Gout<<"\n Matrix prob max "<< wfase_max<<" "<<SelRea->Get_Titolo()<<std::endl;
}

//---------------
 int  Part3pi0::Get_MatrixValue()
{
// attenzioine questa matrix non è simmetrica sulla permutazione????? 
double su=(mresq-s0)/MPIOQ;
double sv=(2.*mpadre*(e2-e1))/MPIOQ;
double sr=sqrt(su*su+sv*sv);
 prb= 1. +slope*sr;

 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.0001; // aggiorno ma informo.....non dovrebbe mai accadere!
  }
  
 if((wmtr_max*Pran()) <prb) return 1;   
  else return 0;
}