EvtOpenCharm Class Reference

#include <EvtOpenCharm.hh>

Inheritance diagram for EvtOpenCharm:

Public Member Functions
	EvtOpenCharm ()
virtual	~EvtOpenCharm ()
void	getName (std::string &name)
EvtDecayBase *	clone ()
void	decay (EvtParticle *p)
std::string	commandName ()
void	command (std::string cmd)
void	init ()
void	initProbMax ()
int	getTotalEvt ()
bool	isbelong (EvtId myid)
int	which_mode (EvtId myid)
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 ()

Static Public Member Functions
static void	OpencrmInit (int f)
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)

Static Public Attributes
static int	myiter = 1

Additional Inherited Members
Protected Member Functions inherited from EvtDecayBase
bool	daugsDecayedByParentModel ()
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

Detailed Description

Definition at line 34 of file EvtOpenCharm.hh.

Constructor & Destructor Documentation

EvtOpenCharm()

EvtOpenCharm::EvtOpenCharm

(

)

Definition at line 62 of file EvtOpenCharm.cc.

{}

Referenced by clone().

~EvtOpenCharm()

EvtOpenCharm::~EvtOpenCharm ( )

virtual

Definition at line 64 of file EvtOpenCharm.cc.

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

Member Function Documentation

clone()

EvtDecayBase * EvtOpenCharm::clone ( )

virtual

Implements EvtDecayBase.

Definition at line 97 of file EvtOpenCharm.cc.

{ return new EvtOpenCharm; }

command()

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

virtual

Reimplemented from EvtDecayBase.

Definition at line 121 of file EvtOpenCharm.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 EvtOpenCharm::commandName ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 119 of file EvtOpenCharm.cc.

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

decay()

void EvtOpenCharm::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 144 of file EvtOpenCharm.cc.

                                         {
 
  int istdheppar = EvtPDL::getStdHep( p->getId() );
 
  if ( istdheppar != 9000443 && istdheppar != 9010443 && istdheppar != 9030443 &&
       istdheppar != 9020443 )
  {
    std::cout << "EvtGen: EvtOpenCharm cann't not decay the particle pid= " << istdheppar
              << endl;
    ::abort();
  }
 
  double mp    = p->mass();
  float  xmp   = mp;
  double totEn = 0;
 
  // debugging
  //   std::cout<<"parent "<<EvtPDL::name(p->getId())<<"float xmp="<<xmp<<"
  //   "<<p->getP4Lab()<<std::endl;
 
  EvtVector4R p4[20];
 
  int i, more;
  int ip = EvtPDL::getStdHep( p->getId() );
  int ndaugjs;
 
  static int    myflag;
  EvtPsi3Sdecay theIni;
  EvtId         pid = p->getId();
 
  static int themode;
  if ( getNArg() == 1 )
  {
    themode = getArg( 0 );
    theIni.setMode( themode );
  }
 
  int count = 0;
  do {
 
    theIni.PHSPDecay( p );
    std::vector<EvtVector4R> v_p4 = theIni.getDaugP4();
    std::vector<EvtId>       Vid  = theIni.getDaugId();
    ndaugjs                       = Vid.size();
 
    EvtId myId[3];
 
    for ( int i = 0; i < ndaugjs; i++ ) { myId[i] = Vid[i]; }
 
    if ( p->getNDaug() != 0 ) p->resetNDaug();
    p->makeDaughters( ndaugjs, myId );
 
    for ( int i = 0; i < ndaugjs; i++ )
    {
      // std::cout<<"Vid "<<EvtPDL::name(Vid[i])<<" p4: "<<v_p4[i]<<std::endl;
      p->getDaug( i )->init( Vid[i], v_p4[i] );
    }
 
    theMode = theIni.getMode();
    p->setGeneratorFlag( theMode );
 
    totEn = 0;
    for ( i = 0; i < ndaugjs; i++ )
    {
      totEn += p->getDaug( i )->getP4().get( 0 );
      if ( p->getDaug( i )->getId() == EvtId( -1, -1 ) )
      {
        report( ERROR, "EvtGen" )
            << "OpenCharm returned particle:" << EvtPDL::name( p->getDaug( i )->getId() )
            << endl;
        report( ERROR, "EvtGen" ) << "This can not be translated to evt number" << endl;
        report( ERROR, "EvtGen" ) << "and the decay will be rejected!" << endl;
        report( ERROR, "EvtGen" ) << "The decay was of particle:" << ip << endl;
      }
    }
    int channel = EvtDecayTable::inChannelList( p->getId(), ndaugjs, myId );
 
    // std::cout<<"channel= "<<channel<<std::endl;
    more = ( channel != -1 );
    // if(more){std::cout<<"count= "<<count<<" inchannel= "<<channel<<std::endl;}
 
    count++;
 
  } while ( more && ( count < 10000 ) );
 
  if ( fabs( xmp - totEn ) > 0.01 )
  {
    std::cout << "Warning:OPENCHARM generate incomplet final state, " << mp << " " << totEn
              << endl;
    ::abort();
  }
 
  if ( count > 9999 )
  {
    report( INFO, "EvtGen" ) << "Too many loops in EvtOpenCharm!!!" << endl;
    report( INFO, "EvtGen" ) << "Parent:" << EvtPDL::name( getParentId() ).c_str() << endl;
  }
 
  fixPolarizations( p );
 
  return;
}

getName()

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

virtual

Implements EvtDecayBase.

Definition at line 95 of file EvtOpenCharm.cc.

{ model_name = "OPENCHARM"; }

getTotalEvt()

int EvtOpenCharm::getTotalEvt ( )

inline

Definition at line 50 of file EvtOpenCharm.hh.

{ return nevt; }

init()

void EvtOpenCharm::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 101 of file EvtOpenCharm.cc.

                        {
 
  //  checkNArg(1);
 
  if ( getParentId().isAlias() )
  {
 
    report( ERROR, "EvtGen" ) << "EvtOpenCharm finds that you are decaying the" << endl
                              << " aliased particle " << EvtPDL::name( getParentId() ).c_str()
                              << " with the OpenCharm model" << endl
                              << " this does not work, please modify decay table." << endl;
    report( ERROR, "EvtGen" ) << "Will terminate execution!" << endl;
    ::abort();
  }
 
  store( this );
}

initProbMax()

void EvtOpenCharm::initProbMax ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 99 of file EvtOpenCharm.cc.

{ noProbMax(); }

isbelong()

bool EvtOpenCharm::isbelong

(

EvtId

myid

)

Definition at line 314 of file EvtOpenCharm.cc.

                                        {
  for ( int i = 0; i < mypar.size(); i++ )
  {
    if ( myid == mypar[i] )
    {
      _index = i;
      return true;
    }
  }
  return false;
}

OpencrmInit()

void EvtOpenCharm::OpencrmInit ( int f )

static

Definition at line 312 of file EvtOpenCharm.cc.

{}

which_mode()

int EvtOpenCharm::which_mode

(

EvtId

myid

)

Definition at line 326 of file EvtOpenCharm.cc.

                                         {
  for ( int i = 0; i < mypar.size(); i++ )
  {
    if ( myid == mypar[i] )
    {
      _index = i;
      return i;
    }
  }
  std::cout << "EvtOpenCharm::which_mode() fails to find element" << std::endl;
  abort();
}

Member Data Documentation

myiter

int EvtOpenCharm::myiter = 1

static

Definition at line 53 of file EvtOpenCharm.hh.

---

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

• EvtOpenCharm.hh
• EvtOpenCharm.cc