---

decay.cxx

Go to the documentation of this file.

/*
Test particle decays in Gismo

  $Id: decay.cxx,v 1.2 2001/08/09 17:18:55 burnett Exp $
*/
#include "facilities/error.h"
#include "gismo/GParticle.h"
#include "gismo/PDT.h"
#include "gismo/Material.h"
#include <iostream>

#if 0 //defined(_DEBUG) && defined (_WIN32)
#  include <CRTDBG.H>
#  pragma comment(lib, "msvcrtd.lib")
#endif
//#ifdef _MSC_DEV
using namespace std;
//#endif
const char* particleName= "pi0";
int   number = 1;

void WARNING(const char * text){std::cerr << "WARNING: " << text << '\n';}
void FATAL(const char * text){
    std::cerr << "FATAL error: " <<text << '\n';
    std::cerr << "type help for instructions\n";
    exit(-1);}

static const char* instructions =
"Format for the gismo test program decay:\n"
"\n"
" decay [[particle-name ] count ]\n\n"
"Where particle-name is the name of a particle in the Particle data table\n"
"                    or 'list' to get a list of particles\n"
"                    or 'help' for these instructions\n"
"      count is the number of particles to decay\n";

int main(int argc, char* argv[])
{
    if( argc>1 ) particleName = argv[1];
    if( argc>2 ) number = atoi(argv[2]);

    if( strcmp(particleName,"list")==0 )
    {
        GParticle::thePDT()->printOn(cout);
        return 0;
    }

    if( strcmp(particleName, "help")==0) {
        cout << instructions;
        return 0;
    }

    std::cout << " Decay test  \n";

    PData* pd = GParticle::thePDT()->lookUp(particleName);
        if( pd==0) {
                std::cout 
                        << "Particle "<<particleName 
                        << " not found in the PDT. Use the 'list' option for a list of particles" 
                        << std::endl;
                return -1;
        }
    std::cout << "                                      -------- decays ---------\n";
    std::cout << "   id  name     mass     width        fraction  type   products";

    pd->printOn(cout);


    for(int n=0; n<number; n++)
    {
        GParticle p(pd);
        p.decayOnce();
        cout  << p << '\n';
        cout.flush();
    }
#if 0 //defined(_DEBUG) && defined (_WIN32)
    // this to verify that memory has been released
    _CrtSetDbgFlag( _CrtSetDbgFlag( _CRTDBG_REPORT_FLAG ) | _CRTDBG_LEAK_CHECK_DF);

    delete GParticle::thePDT();
    Material::deleteAll();

#endif
    return 0;
}

---