#include <EvtD0Topipipi0.hh>

Inheritance diagram for EvtD0Topipipi0:

Public Member Functions
	EvtD0Topipipi0 ()

virtual	~EvtD0Topipipi0 ()

void	getName (std::string &name)

EvtDecayBase *	clone ()

void	init ()

void	initProbMax ()

void	decay (EvtParticle *p)

Public Member Functions inherited from EvtDecayProb
void	makeDecay (EvtParticle *p)

void	setProb (double prob)

double	getProb ()

void	setWeight (double weight)

virtual	~EvtDecayProb ()

Public Member Functions inherited from EvtDecayBase
virtual std::string	commandName ()

virtual void	command (std::string cmd)

double	getProbMax (double prob)

double	resetProbMax (double prob)

	EvtDecayBase ()

virtual	~EvtDecayBase ()

virtual bool	matchingDecay (const EvtDecayBase &other) const

EvtId	getParentId ()

double	getBranchingFraction ()

void	disableCheckQ ()

void	checkQ ()

int	getNDaug ()

EvtId *	getDaugs ()

EvtId	getDaug (int i)

int	getNArg ()

int	getPHOTOS ()

void	setPHOTOS ()

void	setVerbose ()

void	setSummary ()

double *	getArgs ()

std::string *	getArgsStr ()

double	getArg (int j)

std::string	getArgStr (int j)

std::string	getModelName ()

int	getDSum ()

int	summary ()

int	verbose ()

void	saveDecayInfo (EvtId ipar, int ndaug, EvtId *daug, int narg, std::vector< std::string > &args, std::string name, double brfr)

void	printSummary ()

void	setProbMax (double prbmx)

void	noProbMax ()

void	checkNArg (int a1, int a2=-1, int a3=-1, int a4=-1)

void	checkNDaug (int d1, int d2=-1)

void	checkSpinParent (EvtSpinType::spintype sp)

void	checkSpinDaughter (int d1, EvtSpinType::spintype sp)

virtual int	nRealDaughters ()

Additional Inherited Members
Static Public Member Functions inherited from EvtDecayBase
static void	findMasses (EvtParticle *p, int ndaugs, EvtId daugs[10], double masses[10])

static void	findMass (EvtParticle *p)

static double	findMaxMass (EvtParticle *p)

Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()

Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 30 of file EvtD0Topipipi0.hh.

Constructor & Destructor Documentation

◆ EvtD0Topipipi0()

EvtD0Topipipi0::EvtD0Topipipi0 ( )

inline

Definition at line 34 of file EvtD0Topipipi0.hh.

34{}

Referenced by clone().

◆ ~EvtD0Topipipi0()

EvtD0Topipipi0::~EvtD0Topipipi0 ( )

virtual

Definition at line 43 of file EvtD0Topipipi0.cc.

43{}

Member Function Documentation

◆ clone()

EvtDecayBase * EvtD0Topipipi0::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 49 of file EvtD0Topipipi0.cc.

                                   {
  return new EvtD0Topipipi0;
}

◆ decay()

void EvtD0Topipipi0::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 162 of file EvtD0Topipipi0.cc.

                                          {
/*
   double maxprob = 0.0;
   for(int ir=0;ir<=60000000;ir++){
      p->initializePhaseSpace(getNDaug(),getDaugs());
      for(int i=0; i<_nd; i++){
          _p4Lab[i]=p->getDaug(i)->getP4Lab();
          _p4CM[i]=p->getDaug(i)->getP4();
      }
      double Prob = AmplitudeSquare(charm, tagmode);
      if(Prob>maxprob) {
          maxprob=Prob;
          std::cout << "Max PDF = " << ir << " prob= " << Prob << std::endl;
      }
   }
   std::cout << "Max!!!!!!!!!!! " << maxprob<< std::endl;
*/
  p->initializePhaseSpace(getNDaug(),getDaugs());
  for(int i=0; i<_nd; i++){
          _p4Lab[i]=p->getDaug(i)->getP4Lab();
          _p4CM[i]=p->getDaug(i)->getP4();
  }
  double prob = AmplitudeSquare(charm, tagmode);
  setProb(prob);
  return;
}

◆ getName()

void EvtD0Topipipi0::getName ( std::string & name )

virtual

Implements EvtDecayBase.

Definition at line 45 of file EvtD0Topipipi0.cc.

                                                 {
  model_name="D0Topipipi0";     
}

◆ init()

void EvtD0Topipipi0::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 53 of file EvtD0Topipipi0.cc.

                         {
  checkNArg(2);
  checkNDaug(3);
  charm   = getArg(0);
  tagmode = getArg(1);
  //std::cout << "Initializing EvtD0Topipipi0: charm= "<<charm<<" tagmode= "<<tagmode<<std::endl;
 
 
    g_uv.clear();
    for(int i=0; i<4; i++){
        for(int j=0; j<4; j++){
            if(i!=j){
                g_uv.push_back(0.0);
            }else if(i<3){
                g_uv.push_back(-1.0);
            }else if(i==3){
                g_uv.push_back(1.0);
            }
        }
    }
 
    epsilon_uvmn.clear();
    for(int i=0; i<4; i++){
        for(int j=0; j<4; j++){
            for(int k=0; k<4; k++){
                for(int l=0; l<4; l++){
                    if(i==j || i==k || i==l || j==k || j==l || k==l){
                        epsilon_uvmn.push_back(0.0);
                    }else{
                        if(i==0 && j==1 && k==2 && l==3) epsilon_uvmn.push_back(1.0);
                        if(i==0 && j==1 && k==3 && l==2) epsilon_uvmn.push_back(-1.0);
                        if(i==0 && j==2 && k==1 && l==3) epsilon_uvmn.push_back(-1.0);
                        if(i==0 && j==2 && k==3 && l==1) epsilon_uvmn.push_back(1.0);
                        if(i==0 && j==3 && k==1 && l==2) epsilon_uvmn.push_back(1.0);
                        if(i==0 && j==3 && k==2 && l==1) epsilon_uvmn.push_back(-1.0);
 
                        if(i==1 && j==0 && k==2 && l==3) epsilon_uvmn.push_back(-1.0);
                        if(i==1 && j==0 && k==3 && l==2) epsilon_uvmn.push_back(1.0);
                        if(i==1 && j==2 && k==0 && l==3) epsilon_uvmn.push_back(1.0);
                        if(i==1 && j==2 && k==3 && l==0) epsilon_uvmn.push_back(-1.0);
                        if(i==1 && j==3 && k==0 && l==2) epsilon_uvmn.push_back(-1.0);
                        if(i==1 && j==3 && k==2 && l==0) epsilon_uvmn.push_back(1.0);
 
                        if(i==2 && j==0 && k==1 && l==3) epsilon_uvmn.push_back(1.0);
                        if(i==2 && j==0 && k==3 && l==1) epsilon_uvmn.push_back(-1.0);
                        if(i==2 && j==1 && k==0 && l==3) epsilon_uvmn.push_back(-1.0);
                        if(i==2 && j==1 && k==3 && l==0) epsilon_uvmn.push_back(1.0);
                        if(i==2 && j==3 && k==0 && l==1) epsilon_uvmn.push_back(1.0);
                        if(i==2 && j==3 && k==1 && l==0) epsilon_uvmn.push_back(-1.0);
 
                        if(i==3 && j==0 && k==1 && l==2) epsilon_uvmn.push_back(-1.0);
                        if(i==3 && j==0 && k==2 && l==1) epsilon_uvmn.push_back(1.0);
                        if(i==3 && j==1 && k==0 && l==2) epsilon_uvmn.push_back(1.0);
                        if(i==3 && j==1 && k==2 && l==0) epsilon_uvmn.push_back(-1.0);
                        if(i==3 && j==2 && k==0 && l==1) epsilon_uvmn.push_back(-1.0);
                        if(i==3 && j==2 && k==1 && l==0) epsilon_uvmn.push_back(1.0);
 
                    }
                }
            }
        }
    }
 
  _nd = 3;
    math_pi = 3.1415926;
    mass_Pion = 0.13957;
 
    rRes = 3.0*0.197321;
    rD   = 5.0*0.197321;
    m_Pi  = mass_Pion;
    m2_Pi = m_Pi*m_Pi;
    m_Pi0  = 0.134977;
    m2_Pi0 = m_Pi0*m_Pi0;
 
    m0_rho7700 = 0.77526;
    w0_rho7700 = 0.1478;
 
    m0_rho770p = 0.77511;
    w0_rho770p = 0.1491;
 
    m0_f21270 = 1.2755;
    w0_f21270 = 0.1867;
 
    s0_prod = -5.0;
 
    fitpara.clear();
    fitpara.push_back(complex<double>(100,0));
    fitpara.push_back(complex<double>(69.8939*cos(3.14983),69.8939*sin(3.14983)));
    fitpara.push_back(complex<double>(58.521*cos(-2.90685),58.521*sin(-2.90685)));
    fitpara.push_back(complex<double>(483.035*cos(-0.679009),483.035*sin(-0.679009)));
    fitpara.push_back(complex<double>(441.921*cos(-0.879847),441.921*sin(-0.879847)));
    fitpara.push_back(complex<double>(1356.95*cos(-0.206653),1356.95*sin(-0.206653)));
    fitpara.push_back(complex<double>(559.218*cos(0.501728),559.218*sin(0.501728)));
    fitpara.push_back(complex<double>(3165.25*cos(3.39939),3165.25*sin(3.39939)));
    fitpara.push_back(complex<double>(1422.93*cos(3.05347),1422.93*sin(3.05347)));
    fitpara.push_back(complex<double>(2399.8*cos(2.24983),2399.8*sin(2.24983)));
    fitpara.push_back(complex<double>(4601.19*cos(2.74388),4601.19*sin(2.74388)));
    fitpara.push_back(complex<double>(1684.1*cos(1.99894),1684.1*sin(1.99894)));
    fitpara.push_back(complex<double>(678.674*cos(-2.510691),678.674*sin(-2.510691)));
    fitpara.push_back(complex<double>(2.19068*cos(0.991805),2.19068*sin(0.991805)));
 
 
  return;
}

◆ initProbMax()

void EvtD0Topipipi0::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 158 of file EvtD0Topipipi0.cc.

                                {
  setProbMax(60000000);
}

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EvtD0Topipipi0()

◆ ~EvtD0Topipipi0()

Member Function Documentation

◆ clone()

◆ decay()

◆ getName()

◆ init()

◆ initProbMax()