EvtDecayAmp Class Reference

#include <EvtDecayAmp.hh>

Inheritance diagram for EvtDecayAmp:

Public Member Functions
void	makeDecay (EvtParticle *p)
void	setWeight (double weight)
void	vertex (const EvtComplex &amp)
void	vertex (int i1, const EvtComplex &amp)
void	vertex (int i1, int i2, const EvtComplex &amp)
void	vertex (int i1, int i2, int i3, const EvtComplex &amp)
void	vertex (int *i1, const EvtComplex &amp)
virtual	~EvtDecayAmp ()
Public Member Functions inherited from EvtDecayBase
virtual void	getName (std::string &name)=0
virtual void	decay (EvtParticle *p)=0
virtual void	makeDecay (EvtParticle *p)=0
virtual EvtDecayBase *	clone ()=0
virtual void	init ()
virtual void	initProbMax ()
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 ()

Protected Attributes
EvtAmp	_amp2
Protected Attributes inherited from EvtDecayBase
bool	_daugsDecayedByParentModel

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 ()

Detailed Description

Definition at line 27 of file EvtDecayAmp.hh.

Constructor & Destructor Documentation

~EvtDecayAmp()

virtual EvtDecayAmp::~EvtDecayAmp ( )

inlinevirtual

Definition at line 63 of file EvtDecayAmp.hh.

{}

Member Function Documentation

makeDecay()

void EvtDecayAmp::makeDecay ( EvtParticle *  p )

virtual

Implements EvtDecayBase.

Definition at line 35 of file EvtDecayAmp.cc.

                                         {
 
  int ntimes=10000;
 
  int i,more;
 
  EvtSpinDensity rho;
  double prob,prob_max;
 
  _amp2.init(p->getId(),getNDaug(),getDaugs());
//  report(INFO,"EvtGen") << "Decaying " << EvtPDL::name(p->getId()) << endl;
  do{
 
    _daugsDecayedByParentModel=false;
    _weight = 1.0;
    decay(p);
 
    rho=_amp2.getSpinDensity();
    prob=p->getSpinDensityForward().NormalizedProb(rho);
    if (prob<0.0) {
 
      report(ERROR,"EvtGen")<<"Negative prob:"<<p->getId().getId()
                <<" "<<p->getChannel()<<endl;
 
      report(ERROR,"EvtGen") << "rho_forward:"<<endl;
      report(ERROR,"EvtGen") << p->getSpinDensityForward();
      report(ERROR,"EvtGen") << "rho decay:"<<endl;
      report(ERROR,"EvtGen") << rho <<endl;
    }
 
    if (prob!=prob) {
 
      report(DEBUG,"EvtGen") << "Forward density matrix:"<<endl;
      report(ERROR,"EvtGen") << p->getSpinDensityForward();
 
      report(DEBUG,"EvtGen") << "Decay density matrix:"<<endl;
      report(ERROR,"EvtGen") << rho;
 
      report(DEBUG,"EvtGen") << "prob:"<<prob<<endl;
      
      report(DEBUG,"EvtGen") << "Particle:"
                 <<EvtPDL::name(p->getId()).c_str()<<endl;
      report(DEBUG,"EvtGen") << "channel        :"<<p->getChannel()<<endl;
      report(DEBUG,"EvtGen") << "Momentum:" << p->getP4() << " " << p->mass() << endl;
      if( p->getParent()!=0){
    report(DEBUG,"EvtGen") << "parent:"
                   <<EvtPDL::name(
                p->getParent()->getId()).c_str()<<endl;
    report(DEBUG,"EvtGen") << "parent channel        :"
                   <<p->getParent()->getChannel()<<endl;
 
        int i;
    report(DEBUG,"EvtGen") << "parent daughters  :";
        for (i=0;i<p->getParent()->getNDaug();i++){
      report(DEBUG,"") << EvtPDL::name(
                p->getParent()->getDaug(i)->getId()).c_str()
                 << " ";
        }
    report(DEBUG,"") << endl;
 
    report(DEBUG,"EvtGen") << "daughters  :";
        for (i=0;i<p->getNDaug();i++){
      report(DEBUG,"") << EvtPDL::name(
                p->getDaug(i)->getId()).c_str()
                 << " ";
        }
    report(DEBUG,"") << endl;
 
    report(DEBUG,"EvtGen") << "daughter momenta  :" << endl;;
        for (i=0;i<p->getNDaug();i++){
      report(DEBUG,"") << p->getDaug(i)->getP4() << " " << p->getDaug(i)->mass();
      report(DEBUG,"") << endl;
        }
 
      }
    }
 
 
    prob/=_weight;
 
 
    prob_max = getProbMax(prob);
    p->setDecayProb(prob/prob_max);
 
    //report(INFO,"EvtGen") << "Prob,prob_max,weight:"<<prob<<" "<<prob_max<<" "<<_weight<<endl;
 
    more=prob<EvtRandom::Flat(prob_max);
 
    ntimes--;
 
  }while(ntimes&&more);
  //report(INFO,"EvtGen") << "Done\n";
 
  if (ntimes==0){
    report(DEBUG,"EvtGen") << "Tried accept/reject:10000"
               <<" times, and rejected all the times!"<<endl;
    report(DEBUG,"")<<p->getSpinDensityForward()<<endl;
    report(DEBUG,"EvtGen") << "Is therefore accepting the last event!"<<endl;
    report(DEBUG,"EvtGen") << "Decay of particle:"<<
      EvtPDL::name(p->getId()).c_str()<<"(channel:"<<
      p->getChannel()<<") with mass "<<p->mass()<<endl;
    
    int ii;
    for(ii=0;ii<p->getNDaug();ii++){
      report(DEBUG,"EvtGen") <<"Daughter "<<ii<<":"<<
    EvtPDL::name(p->getDaug(ii)->getId()).c_str()<<" with mass "<<
    p->getDaug(ii)->mass()<<endl;
    }                  
  }
 
 
  EvtSpinDensity rho_list[10];
 
  rho_list[0]=p->getSpinDensityForward();
   
 
  EvtAmp ampcont;
 
  if (_amp2._pstates!=1){
    ampcont=_amp2.contract(0,p->getSpinDensityForward());   // J2BB2 bugging here
 
  }
  else{
    ampcont=_amp2;
  }
 
          
 
  // it may be that the parent decay model has already
  // done the decay - this should be rare and the
  // model better know what it is doing..
  
  if ( !daugsDecayedByParentModel() ){
    
 
 //   report(INFO,"EvtGen") << "Found " << p->getNDaug() << " daughters\n";
    for(i=0;i<p->getNDaug();i++){
    
      rho.SetDim(_amp2.dstates[i]);
 
      if (_amp2.dstates[i]==1) {
    rho.Set(0,0,EvtComplex(1.0,0.0));
      }
      else{
    rho=ampcont.contract(_amp2._dnontrivial[i],_amp2);
      }
      
      if (!rho.Check()) {
    
    report(ERROR,"EvtGen") << "-------start error-------"<<endl;
    report(ERROR,"EvtGen")<<"forward rho failed Check:"<<
      EvtPDL::name(p->getId()).c_str()<<" "<<p->getChannel()<<" "<<i<<endl;
    
    report(ERROR,"EvtGen")<<"Parent:"<<EvtPDL::name(p->getParent()->getId()).c_str()<<endl;
    report(ERROR,"EvtGen")<<"GrandParent:"<<EvtPDL::name(p->getParent()->getParent()->getId()).c_str()<<endl;
    report(ERROR,"EvtGen")<<"GrandGrandParent:"<<EvtPDL::name(p->getParent()->getParent()->getParent()->getId()).c_str()<<endl;
    
    report(ERROR,"EvtGen") << rho;
    int ii; 
    _amp2.dump();
        
    for(ii=0;ii<i+1;ii++){
      report(ERROR,"EvtGen") << rho_list[ii];
    }
    report(ERROR,"EvtGen") << "-------Done with error-------"<<endl;  
      }
 
 
      
      p->getDaug(i)->setSpinDensityForward(rho);
 
      //debugging
      // p->printTree();
     p->getDaug(i)->decay();
 
     rho_list[i+1]= p->getDaug(i)->getSpinDensityBackward();
 
      if (_amp2.dstates[i]!=1){
    ampcont=ampcont.contract(_amp2._dnontrivial[i],rho_list[i+1]);
      }
  
      
    }
 
    p->setSpinDensityBackward(_amp2.getBackwardSpinDensity(rho_list));
   
    
    if (!p->getSpinDensityBackward().Check()) {
      
      report(ERROR,"EvtGen")<<"rho_backward failed Check"<<
    p->getId().getId()<<" "<<p->getChannel()<<endl;
      
      report(ERROR,"EvtGen") << p->getSpinDensityBackward();
      
    }
  }
 
  if (getPHOTOS() || EvtRadCorr::alwaysRadCorr()) {
    EvtRadCorr::doRadCorr(p);
 
  }
 
}

setWeight()

void EvtDecayAmp::setWeight ( double  weight )

inline

Definition at line 32 of file EvtDecayAmp.hh.

{_weight=weight;}

Referenced by EvtbTosllAli::decay(), and EvtbTosllBall::decay().

vertex() [1/5]

void EvtDecayAmp::vertex ( const EvtComplex &  amp )

inline

sets the amplitudes calculated in the decay objects

Definition at line 37 of file EvtDecayAmp.hh.

{_amp2.vertex(amp);}

Referenced by EvtAV2GV::decay(), EvtBHadronic::decay(), EvtBto2piCPiso::decay(), EvtBTo3piCP::decay(), EvtBTo4piCP::decay(), EvtBtoKD3P::decay(), EvtBtoKpiCPiso::decay(), EvtBToKpipiCP::decay(), EvtCBTo3piMPP::decay(), EvtCBTo3piP00::decay(), EvtChi0BB1::decay(), EvtChi0BB2::decay(), EvtChi1BB1::decay(), EvtChi1BB2::decay(), EvtChi2BB1::decay(), EvtChi2BB2::decay(), EvtD0mixDalitz::decay(), EvtDDalitz::decay(), EvtDeBD::decay(), EvtEDM::decay(), EvtEtaDalitz::decay(), EvtHelPPJ::decay(), EvtHypWK::decay(), EvtIntervalDecayAmp< T >::decay(), EvtJ2BB1::decay(), EvtJ2BB2::decay(), EvtJ2BB3::decay(), EvtJPE::decay(), EvtJpipi::decay(), EvtJTO3P::decay(), EvtKstarnunu::decay(), EvtKstarstargamma::decay(), EvtLambdaP_BarGamma::decay(), EvtLNuGamma::decay(), EvtMultibody::decay(), EvtOmegaDalitz::decay(), EvtP2GC0::decay(), EvtP2GC1::decay(), EvtP2GC2::decay(), EvtPBB1::decay(), EvtPBB2::decay(), EvtPhiDalitz::decay(), EvtS2GV::decay(), EvtSll::decay(), EvtSLN::decay(), EvtSPL::decay(), EvtSSDCP::decay(), EvtSSSCP::decay(), EvtSSSCPpng::decay(), EvtSSSCPT::decay(), EvtSTS::decay(), EvtSTSCP::decay(), EvtSVPCP::decay(), EvtSVPHelAmp::decay(), EvtSVS::decay(), EvtSVSCP::decay(), EvtSVSCPiso::decay(), EvtSVSCPLH::decay(), EvtSVSNONCPEIGEN::decay(), EvtT2GV::decay(), EvtTauGamMu::decay(), EvtTauHadnu::decay(), EvtTaulnunu::decay(), EvtTauScalarnu::decay(), EvtTauVectornu::decay(), EvtTSS::decay(), EvtTVSPwave::decay(), EvtVll::decay(), EvtVPHOtoVISR::decay(), EvtVPHOtoVISRHi::decay(), EvtVSPPwave::decay(), EvtVSS::decay(), EvtVSSBMixCPT::decay(), EvtVSSMix::decay(), EvtVVP::decay(), EvtVVpipi::decay(), EvtVVPIPI_WEIGHTED::decay(), EvtVVSPwave::decay(), and EvtTauGamMu::init().

vertex() [2/5]

void EvtDecayAmp::vertex ( int *  i1, const EvtComplex &  amp  )

inline

sets the amplitudes calculated in the decay objects

Definition at line 60 of file EvtDecayAmp.hh.

  { _amp2.vertex(i1,amp);}

vertex() [3/5]

void EvtDecayAmp::vertex ( int  i1, const EvtComplex &  amp  )

inline

sets the amplitudes calculated in the decay objects

Definition at line 42 of file EvtDecayAmp.hh.

{_amp2.vertex(i1,amp);}

vertex() [4/5]

void EvtDecayAmp::vertex ( int  i1, int  i2, const EvtComplex &  amp  )

inline

sets the amplitudes calculated in the decay objects

Definition at line 47 of file EvtDecayAmp.hh.

  {_amp2.vertex(i1,i2,amp);}

vertex() [5/5]

void EvtDecayAmp::vertex ( int  i1, int  i2, int  i3, const EvtComplex &  amp  )

inline

sets the amplitudes calculated in the decay objects

Definition at line 54 of file EvtDecayAmp.hh.

  {_amp2.vertex(i1,i2,i3,amp);}

Member Data Documentation

_amp2

EvtAmp EvtDecayAmp::_amp2

protected

Definition at line 66 of file EvtDecayAmp.hh.

Referenced by EvtbTosllAli::decay(), EvtbTosllBall::decay(), EvtHelAmp::decay(), EvtHQET::decay(), EvtHQET2::decay(), EvtHypNonLepton::decay(), EvtISGW::decay(), EvtISGW2::decay(), EvtKKLambdaC::decay(), EvtMelikhov::decay(), EvtPartWave::decay(), EvtSLBKPole::decay(), EvtSLPole::decay(), EvtSVVCP::decay(), EvtSVVCPLH::decay(), EvtSVVHelAmp::decay(), EvtSVVNONCPEIGEN::decay(), makeDecay(), and vertex().

---

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

• EvtDecayAmp.hh
• EvtDecayAmp.cc