---

BaseApp.cxx

Go to the documentation of this file.

// $Id: BaseApp.cxx,v 1.2 2001/10/06 14:18:53 burnett Exp $
//
// file BaseApp.cc, a part of the Gismo project
//
//

#include "BaseApp.h"

#include "geometry/Box.h"
#include "gismo/MCParticle.h"
#include "gismo/PDT.h"
#include "gismo/Field.h"
#include <signal.h>
#include <stdio.h>
#ifdef _MSC_VER
  using namespace std;
#endif

BaseApp::BaseApp()
: eventNumber(0)
,  _world(0)
, _generator(0)
{}

BaseApp& 
BaseApp::setBeam(const char* name, const Vector& p, const Point& r)
{  
  return setGenerator( new Generator(name,p,r));
}

BaseApp&
BaseApp::setWorld(World* w)
{ 
  if( _world )
    delete _world;
  _world=w;  return *this;
}

BaseApp&
BaseApp::setGenerator(Generator * newgen)
{
  if( _generator )
         delete _generator;
  _generator = newgen;
  return *this;
} 
BaseApp&
BaseApp::addInteractor(Interactor* inter)
{
 MCParticle::addInteractor(inter); return *this;
}
BaseApp::~BaseApp()
{
  if(_world)delete _world;
  if(_generator)delete _generator;
}
World& BaseApp::world()
{  if(!_world) FATAL("BaseApp: No world has been set!");
   return *_world;
}
Generator& BaseApp::generator()
{
   if(!_generator) FATAL("BaseApp: no generator has been defined!");
   return *_generator;
}
//------------------------------------------------------------
//  Some stuff for signals, that should be eventually done with
//  C++ 
#include <signal.h>
#include <setjmp.h>

static jmp_buf env;             /* holds environment to get back after trap */
static int valid = 0;
static int gotControlC=0;


typedef void (*fptr)(int); // really need the (int)?
void catcher( int sig, int /* code*/, struct sigcontext * /*scp*/ )
{
  if( sig==SIGINT )
  {  // if single control-C, just set flag to stop run
     if( !gotControlC )
     {  gotControlC = 1;
        return;
     }
  }
  if(valid)
     longjmp(env, sig);

  FATAL("Invalid interrupt");
}

void setup_signals()
{           // took off (ftpr) cast for catcher ??
        if( (int)signal(SIGINT, (fptr)catcher) == -1 )
                 FATAL("failure to setup signal SIGINT");
        if( (int)signal(SIGFPE, (fptr)catcher) == -1 )
                 FATAL("failure to setup signal SIGFPE");
}

static const char * error;

const char*
BaseApp::run(unsigned long nevent)
{
  setup_signals();
  eventNumber=0;
  gotControlC = 0;
  while( nevent-- > 0)
  {  valid=1;
     int val;
     if( (val = setjmp(env)) != 0 )
     {  /* got here because of longjump with value val */
        if( val == SIGINT )
                 return "User abort";
        else
                 return "Run terminated: hardware trap ";
     }
     if( (error=trigger())!=0 ) return error;
     if( gotControlC ){
         WARNING("user pause of run");
         break;
     }
  }
  valid = 0;
  return 0;

}
void
BaseApp::stop()
{
  gotControlC=1;
}

const char *
BaseApp::trigger()
{
  if( !_generator ) return "no event source!";
  eventNumber++;
  clearEvent();
  clearResponse();
  MCParticle* thisEvent = _generator->getEvent();
  if( thisEvent==0 )
          return "end of data";
  beginEvent();
  const Medium* whereToStart =  _world->inside(_generator->position());
  if( ! whereToStart )
     FATAL("Inital position is outside known world\n");
  _generator->propagate(whereToStart);
  _world->generateResponse();
  endEvent();
  return 0;
}
void
BaseApp::printMaterials(ostream& os)
{
    os << "\n Print materials not yet implemented \n";
};

void
BaseApp::clearEvent()
{
}
void
BaseApp::printPDT(ostream& os)
{
  os << "============================ PDT ================================\n";
  GParticle::thePDT()->printOn(os);
}

void
BaseApp::printEvent(ostream& os)
{
  os << "=========================== event "
         << _generator->eventNumber() << "=====================\n";
  if( _generator)
                os << *_generator;
  os << "\n===============================================================\n";
}

void 
BaseApp::printResponse(ostream& os)
{
  if(_world)
     _world->printResponse(os);
}

void 
BaseApp::printDetector(ostream& os)
{                    
  os <<  "======================= Detector ==============================\n";
  if(_world)
     _world->printOn(os);
  os << "\n===============================================================\n";

}

---