EvtDecayIncoherent.cc

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 1998      Caltech, UCSB
//
// Module: EvtGen/EvtDecayIncoherent.cc
//
// Description:
//
// Modification history:
//
//    DJL/RYD     August 11, 1998         Module created
//
//------------------------------------------------------------------------
#include "EvtPatches.hh"
 
#include "EvtDecayBase.hh"
#include "EvtDecayIncoherent.hh"
#include "EvtPDL.hh"
#include "EvtParticle.hh"
#include "EvtRadCorr.hh"
#include "EvtReport.hh"
 
void EvtDecayIncoherent::makeDecay( EvtParticle* p ) {
 
  int i;
  // initialize this the hard way..
  // Lange June 26, 2000
  for ( i = 0; i < MAX_DAUG; i++ ) { spinDensitySet[i] = 0; }
  _daugsDecayedByParentModel = false;
 
  decay( p );
  p->setDecayProb( 1. );
 
  EvtSpinDensity rho;
 
  rho.SetDiag( p->getSpinStates() );
 
  p->setSpinDensityBackward( rho );
 
  if ( getPHOTOS() || EvtRadCorr::alwaysRadCorr() ) { EvtRadCorr::doRadCorr( p ); }
 
  // Now decay the daughters.
 
  if ( !daugsDecayedByParentModel() )
  {
 
    for ( i = 0; i < p->getNDaug(); i++ )
    {
      // Need to set the spin density of the daughters to be
      // diagonal.
      rho.SetDiag( p->getDaug( i )->getSpinStates() );
      // if (p->getDaug(i)->getNDaug()==0){
      // only do this if the user has not already set the
      // spin density matrix herself.
      // Lange June 26, 2000
      if ( isDaughterSpinDensitySet( i ) == 0 )
      { p->getDaug( i )->setSpinDensityForward( rho ); }
      else
      {
        // report(INFO,"EvtGen") << "spinDensitymatrix already set!!!\n";
        EvtSpinDensity temp = p->getDaug( i )->getSpinDensityForward();
        //  report(INFO,"EvtGen") <<temp<<endl;
      }
      // Now decay the daughter.  Really!
      p->getDaug( i )->decay();
    }
  }
}