EvtPhokhara_pipieta Class Reference

#include <EvtPhokhara_pipieta.hh>

Inheritance diagram for EvtPhokhara_pipieta:

Public Member Functions
	EvtPhokhara_pipieta ()
virtual	~EvtPhokhara_pipieta ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	decay (EvtParticle *p)
std::string	commandName ()
void	command (std::string cmd)
void	init ()
void	init_mode (EvtParticle *p)
void	init_evt (EvtParticle *p)
void	initProbMax ()
int	getTotalEvt ()
void	PhokharaInit (int dummy)
void	ExclusiveDecay (EvtParticle *p)
Public Member Functions inherited from EvtDecayIncoherent
void	makeDecay (EvtParticle *p)
virtual	~EvtDecayIncoherent ()
void	setDaughterSpinDensity (int daughter)
int	isDaughterSpinDensitySet (int daughter)
Public Member Functions inherited from EvtDecayBase
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 33 of file EvtPhokhara_pipieta.hh.

Constructor & Destructor Documentation

EvtPhokhara_pipieta()

EvtPhokhara_pipieta::EvtPhokhara_pipieta

(

)

Definition at line 66 of file EvtPhokhara_pipieta.cc.

{}

Referenced by clone().

~EvtPhokhara_pipieta()

EvtPhokhara_pipieta::~EvtPhokhara_pipieta ( )

virtual

Definition at line 67 of file EvtPhokhara_pipieta.cc.

                                          {
  int i;
  // the deletion of commands is really uggly!
 
  if ( nphokharadecays == 0 )
  {
    delete[] commands;
    commands = 0;
    return;
  }
 
  for ( i = 0; i < nphokharadecays; i++ )
  {
    if ( phokharadecays[i] == this )
    {
      phokharadecays[i] = phokharadecays[nphokharadecays - 1];
      nphokharadecays--;
      if ( nphokharadecays == 0 )
      {
        delete[] commands;
        commands = 0;
      }
      return;
    }
  }
 
  report( ERROR, "EvtGen" ) << "Error in destroying Phokhara model!" << endl;
}

Member Function Documentation

clone()

EvtDecayBase * EvtPhokhara_pipieta::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 101 of file EvtPhokhara_pipieta.cc.

{ return new EvtPhokhara_pipieta; }

command()

void EvtPhokhara_pipieta::command ( std::string cmd )

virtual

Reimplemented from EvtDecayBase.

Definition at line 138 of file EvtPhokhara_pipieta.cc.

                                                 {
 
  if ( ncommand == lcommand )
  {
 
    lcommand = 10 + 2 * lcommand;
 
    std::string* newcommands = new std::string[lcommand];
 
    int i;
 
    for ( i = 0; i < ncommand; i++ ) { newcommands[i] = commands[i]; }
 
    delete[] commands;
 
    commands = newcommands;
  }
 
  commands[ncommand] = cmd;
 
  ncommand++;
}

commandName()

std::string EvtPhokhara_pipieta::commandName ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 136 of file EvtPhokhara_pipieta.cc.

{ return std::string( "PhokharaPar" ); }

decay()

void EvtPhokhara_pipieta::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 161 of file EvtPhokhara_pipieta.cc.

                                                {
  EvtId myvpho = EvtPDL::getId( "vpho" );
  if ( p->getId() != myvpho )
  {
    std::cout << "Parent particle is required to be vpho for Phokhara model" << std::endl;
    abort();
  }
  if ( nevtgen == 0 ) { init_mode( p ); }
  else { init_evt( p ); }
  std::cout << "PHOKHARA : pi+pi-eta mode " << std::endl;
 
  int istdheppar = EvtPDL::getStdHep( p->getId() );
  int ntrials    = 0;
  int tr_old[3];
  tr_old[0] = (int)maxima_.tr[0];
  tr_old[1] = (int)maxima_.tr[1];
  tr_old[2] = (int)maxima_.tr[2];
 
  while ( ntrials < 1000000 )
  {
    ievent++;
    RANLXDF( Ar_r, 1 );
    Ar[1] = Ar_r[0];
 
    if ( Ar[1] <=
         ( maxima_.Mmax[0] / ( maxima_.Mmax[0] + maxima_.Mmax[1] + maxima_.Mmax[2] ) ) )
    {
      maxima_.count[0] = maxima_.count[0] + 1.0;
      GEN_0PH( 2, qqmin, ctes_.Sp, cos3min, cos3max );
    }
    else if ( Ar[1] <= ( ( maxima_.Mmax[0] + maxima_.Mmax[1] ) /
                         ( maxima_.Mmax[0] + maxima_.Mmax[1] + maxima_.Mmax[2] ) ) )
    {
      maxima_.count[1] = maxima_.count[1] + 1.0;
      GEN_1PH( 2, qqmin, qqmax, cos1min, cos1max, cos3min, cos3max );
    }
    else
    {
      maxima_.count[2] = maxima_.count[2] + 1.0;
      GEN_2PH( 2, qqmin, cos1min, cos1max, cos2min, cos2max, cos3min, cos3max );
    }
 
    if ( ( (int)maxima_.tr[0] + (int)maxima_.tr[1] + (int)maxima_.tr[2] ) >
         ( tr_old[0] + tr_old[1] + tr_old[2] ) ) // event accepted after cuts
    { goto storedEvents; }
    ntrials++;
  }
  std::cout << "FATAL: Could not satisfy cuts after " << ntrials << "trials. Terminate."
            << std::endl;
  //----
storedEvents:
  int         more      = 0;
  int         numstable = 0;
  int         numparton = 0;
  EvtId       evtnumstable[100]; //
  EvtVector4R p4[20];
 
  // except ISR photos, products depending on channel
  if ( flags_.pion == 0 )
  { // mu+ mu-
    // mu+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -13 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // mu -
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 13 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
  }
  if ( flags_.pion == 1 )
  { // pi+ pi-
    // pi+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // pi -
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
  }
  if ( flags_.pion == 2 )
  { // pi+ pi-2pi0
    // pi0
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 111 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // pi0
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 111 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
    // pi-
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][7], ctes_.momenta[1][7], ctes_.momenta[2][7],
                       ctes_.momenta[3][7] );
    numstable++;
    // pi +
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][8], ctes_.momenta[1][8], ctes_.momenta[2][8],
                       ctes_.momenta[3][8] );
    numstable++;
  }
  if ( flags_.pion == 3 )
  { // 2(pi+ pi-)
    // pi+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // pi-
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
    // pi+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][7], ctes_.momenta[1][7], ctes_.momenta[2][7],
                       ctes_.momenta[3][7] );
    numstable++;
    // pi -
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][8], ctes_.momenta[1][8], ctes_.momenta[2][8],
                       ctes_.momenta[3][8] );
    numstable++;
  }
  if ( flags_.pion == 4 )
  { // ppbar
    // pbar
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -2212 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // p
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 2212 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
  }
  if ( flags_.pion == 5 )
  { // nnbar
    // pbar
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -2112 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // p
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 2112 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
  }
  if ( flags_.pion == 6 )
  { // K+ K-
    // K+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 321 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // K -
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -321 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
  }
  if ( flags_.pion == 7 )
  { // K0K0bar
    // Kbar
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 311 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // K0
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -311 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
  }
  if ( flags_.pion == 8 )
  { // pi+ pi-pi0
    // pi+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // pi-
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
    // pi0
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 111 );
    p4[numstable].set( ctes_.momenta[0][7], ctes_.momenta[1][7], ctes_.momenta[2][7],
                       ctes_.momenta[3][7] );
    numstable++;
  }
  if ( flags_.pion == 9 )
  { // Lambda Lambdabar-> pi+ pi- ppbar
    // pi+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][7], ctes_.momenta[1][7], ctes_.momenta[2][7],
                       ctes_.momenta[3][7] );
    numstable++;
    // pbar
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -2212 );
    p4[numstable].set( ctes_.momenta[0][8], ctes_.momenta[1][8], ctes_.momenta[2][8],
                       ctes_.momenta[3][8] );
    numstable++;
    // pi-
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][9], ctes_.momenta[1][9], ctes_.momenta[2][9],
                       ctes_.momenta[3][9] );
    numstable++;
    // p
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 2212 );
    p4[numstable].set( ctes_.momenta[0][10], ctes_.momenta[1][10], ctes_.momenta[2][10],
                       ctes_.momenta[3][10] );
    numstable++;
  }
  if ( flags_.pion == 10 )
  { // pi+ pi-eta
    // pi+
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 211 );
    p4[numstable].set( ctes_.momenta[0][5], ctes_.momenta[1][5], ctes_.momenta[2][5],
                       ctes_.momenta[3][5] );
    numstable++;
    // pi-
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( -211 );
    p4[numstable].set( ctes_.momenta[0][6], ctes_.momenta[1][6], ctes_.momenta[2][6],
                       ctes_.momenta[3][6] );
    numstable++;
    // eta
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 221 );
    p4[numstable].set( ctes_.momenta[0][7], ctes_.momenta[1][7], ctes_.momenta[2][7],
                       ctes_.momenta[3][7] );
    numstable++;
  }
 
  // ISR gamma
  evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 22 );
  p4[numstable].set( ctes_.momenta[0][2], ctes_.momenta[1][2], ctes_.momenta[2][2],
                     ctes_.momenta[3][2] );
  numstable++;
  if ( ctes_.momenta[0][3] != 0 ) // second photon exists
  {
    evtnumstable[numstable] = EvtPDL::evtIdFromStdHep( 22 );
    p4[numstable].set( ctes_.momenta[0][3], ctes_.momenta[1][3], ctes_.momenta[2][3],
                       ctes_.momenta[3][3] );
    numstable++;
  }
 
  int channel = EvtDecayTable::inChannelList( p->getId(), numstable, evtnumstable );
  more        = ( channel != -1 );
  if ( more )
  {
    std::cout << "Existence of mode " << channel
              << " in exclusive decay list has the same final state as this one" << std::endl;
    abort();
  }
 
  p->makeDaughters( numstable, evtnumstable );
  // double weight = p->initializePhaseSpace(getNDaug(),getDaugs());
 
  int         ndaugFound = 0;
  EvtVector4R SUMP4( 0, 0, 0, 0 );
  for ( int i = 0; i < numstable; i++ )
  {
    p->getDaug( i )->init( evtnumstable[i], p4[i] );
    ndaugFound++;
  }
  if ( ndaugFound == 0 )
  {
    report( ERROR, "EvtGen" ) << "Phokhara has failed to do a decay ";
    report( ERROR, "EvtGen" ) << EvtPDL::name( p->getId() ).c_str() << " " << p->mass()
                              << endl;
    assert( 0 );
  }
 
  nevtgen++;
  return;
}

ExclusiveDecay()

void EvtPhokhara_pipieta::ExclusiveDecay

(

EvtParticle *

p

)

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 96 of file EvtPhokhara_pipieta.cc.

                                                         {
 
  model_name = "PHOKHARA_pipieta";
}

getTotalEvt()

int EvtPhokhara_pipieta::getTotalEvt ( )

inline

Definition at line 51 of file EvtPhokhara_pipieta.hh.

{ return nevt; }

init()

void EvtPhokhara_pipieta::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 114 of file EvtPhokhara_pipieta.cc.

                               {
  checkNArg( 0 );
 
  std::string locvp = getenv( "BESEVTGENROOT" );
  system( "cat $BESEVTGENROOT/share/phokhara_10.0.param>phokhara_10.0.param" );
  system( "cat $BESEVTGENROOT/share/phokhara_10.0.fferr>phokhara_10.0.fferr" );
  system( "cat $BESEVTGENROOT/share/phokhara_10.0.ffwarn>phokhara_10.0.ffwarn" );
 
  if ( getParentId().isAlias() )
  {
 
    report( ERROR, "EvtGen" ) << "EvtPhokhara finds that you are decaying the" << endl
                              << " aliased particle " << EvtPDL::name( getParentId() ).c_str()
                              << " with the Phokhara model" << endl
                              << " this does not work, please modify decay table." << endl;
    report( ERROR, "EvtGen" ) << "Will terminate execution!" << endl;
    ::abort();
  }
 
  store( this );
}

init_evt()

void EvtPhokhara_pipieta::init_evt

(

EvtParticle *

p

)

======== list parameters to be initialized

Definition at line 474 of file EvtPhokhara_pipieta.cc.

                                                   {
  m_pion = 10;
  // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),
  // 2pi+2pi-(3),ppbar(4),nnbar(5),
  // K+K-(6),K0K0bar(7),pi+pi-pi0(8),
  // Lamb Lambbar->pi-pi+ppbar(9)
#include "Phokhara_init_evt.txt"
}

Referenced by decay().

init_mode()

void EvtPhokhara_pipieta::init_mode

(

EvtParticle *

p

)

======== list parameters to be initialized

Definition at line 105 of file EvtPhokhara_pipieta.cc.

                                                    {
  m_pion = 10;
  // mu+mu-(0),pi+pi-(1),2pi0pi+pi-(2),
  // 2pi+2pi-(3),ppbar(4),nnbar(5),
  // K+K-(6),K0K0bar(7),pi+pi-pi0(8),
  // Lamb Lambbar->pi-pi+ppbar(9)
#include "Phokhara_init_mode.txt"
}

Referenced by decay().

initProbMax()

void EvtPhokhara_pipieta::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 103 of file EvtPhokhara_pipieta.cc.

{ noProbMax(); }

PhokharaInit()

void EvtPhokhara_pipieta::PhokharaInit

(

int

dummy

)

Definition at line 463 of file EvtPhokhara_pipieta.cc.

                                                  {
  static int first = 1;
  if ( first )
  {
 
    first = 0;
    // for(int i=0;i<ncommand;i++)
    //  lugive0_(commands[i].c_str(),strlen(commands[i].c_str()));
  }
}

---

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

• EvtPhokhara_pipieta.hh
• EvtPhokhara_pipieta.cc