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 36 of file EvtOpenCharm.hh.

Constructor & Destructor Documentation

EvtOpenCharm()

EvtOpenCharm::EvtOpenCharm

(

)

Definition at line 65 of file EvtOpenCharm.cc.

{}

Referenced by clone().

~EvtOpenCharm()

EvtOpenCharm::~EvtOpenCharm ( )

virtual

Definition at line 67 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 104 of file EvtOpenCharm.cc.

                                 {
 
  return new EvtOpenCharm;
 
}

command()

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

virtual

Reimplemented from EvtDecayBase.

Definition at line 147 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 141 of file EvtOpenCharm.cc.

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

decay()

void EvtOpenCharm::decay ( EvtParticle * p )

virtual

Implements EvtDecayBase.

Definition at line 176 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 98 of file EvtOpenCharm.cc.

                                               {
 
  model_name="OPENCHARM";     
 
}

getTotalEvt()

int EvtOpenCharm::getTotalEvt ( )

inline

Definition at line 53 of file EvtOpenCharm.hh.

{return nevt;}

init()

void EvtOpenCharm::init ( )

virtual

Reimplemented from EvtDecayBase.

Definition at line 118 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 111 of file EvtOpenCharm.cc.

                              {
 
  noProbMax();
 
}

isbelong()

bool EvtOpenCharm::isbelong

(

EvtId

myid

)

Definition at line 343 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 339 of file EvtOpenCharm.cc.

                                       {
}

which_mode()

int EvtOpenCharm::which_mode

(

EvtId

myid

)

Definition at line 350 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 56 of file EvtOpenCharm.hh.

---

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

• EvtOpenCharm.hh
• EvtOpenCharm.cc