TEvtGen/EvtGenBase/EvtParticleDecay.cxx

   1 //--------------------------------------------------------------------------
   2 //
   3 // Environment:
   4 //      This software is part of the EvtGen package developed jointly
   5 //      for the BaBar and CLEO collaborations.  If you use all or part
   6 //      of it, please give an appropriate acknowledgement.
   7 //
   8 // Copyright Information: See EvtGen/COPYRIGHT
   9 //      Copyright (C) 1998      Caltech, UCSB
  10 //
  11 // Module: EvtDecayParm.cc
  12 //
  13 // Description: Store decay parameters for one decay.
  14 //
  15 // Modification history:
  16 //
  17 //    RYD     April 5, 1997         Module created
  18 //
  19 //------------------------------------------------------------------------
  20 //
  21 #include "EvtGenBase/EvtPatches.hh"
  22 #include <iostream>
  23 #include <fstream>
  24 #include <stdlib.h>
  25 #include <ctype.h>
  26 #include "EvtGenBase/EvtParticleDecay.hh"
  27 #include "EvtGenBase/EvtParticle.hh"
  28 #include "EvtGenBase/EvtRandom.hh"
  29 #include "EvtGenBase/EvtReport.hh"
  30 #include "EvtGenBase/EvtPDL.hh"
  31 #include "EvtGenBase/EvtId.hh"
  32 #include <string>
  33 #include <vector>
  34 using std::fstream;
  35 void EvtParticleDecay::printSummary(){
  36
  37   if (_decay!=0) _decay->printSummary();
  38
  39 }
  40
  41 void EvtParticleDecay::chargeConj(EvtParticleDecay *decay){
  42
  43   _brfrsum=decay->_brfrsum;
  44   _massmin=decay->_massmin;
  45
  46   _decay=decay->_decay->clone();
  47
  48   int ndaug=decay->_decay->getNDaug();
  49   int narg=decay->_decay->getNArg();
  50   double brfr=decay->_decay->getBranchingFraction();
  51   std::string name=decay->_decay->getName();
  52   EvtId ipar=EvtPDL::chargeConj(decay->_decay->getParentId());
  53   int i;
  54   EvtId* daug=new EvtId[ndaug];
  55   for(i=0;i<ndaug;i++){
  56     daug[i]=EvtPDL::chargeConj(decay->_decay->getDaug(i));
  57   }
  58   //Had to add 1 to make sure the vector is not empty!
  59   std::vector<std::string> args;
  60   for(i=0;i<narg;i++){
  61     args.push_back(decay->_decay->getArgStr(i));
  62   }
  63
  64   _decay->saveDecayInfo(ipar,ndaug,daug,narg,args,name,brfr);
  65
  66   if (decay->_decay->getPHOTOS()) _decay->setPHOTOS();
  67
  68   delete [] daug;
  69
  70 }
  71
  72
  73
  74
  75
  76
  77