analysis/simfast_cmp.C File Reference

Go to the source code of this file.

Functions
void	simfast_cmp (TString infile="dsdsj20k.evt", int nevts=0, TString outfile="")

Function Documentation

void simfast_cmp	(	TString	infile = "dsdsj20k.evt",
		int	nevts = 0,
		TString	outfile = ""
	)

Definition at line 6 of file analysis/simfast_cmp.C.

References PndFastSim::AddDetector(), ctime, Double_t, PndFastSim::EnablePropagation(), evtGen, fRun, gDebug, infile(), outfile, primGen, printf(), rootlogon(), rtime, PndFastSim::SetVerbosity(), and timer.

{
  TStopwatch timer;
  timer.Start();
  gDebug=0;

  // Load basic libraries
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/rootlogon.C");
  rootlogon();

  // Load this example libraries
  gSystem->Load("libRho");
  gSystem->Load("libfastsim");
 
  if (outfile=="") outfile=infile+".fast.root";

  FairRunSim *fRun = new FairRunSim();
  fRun->SetOutputFile(outfile.Data());
  fRun->SetName("TGeant3");

  // Create and Set Event Generator
  //-------------------------------

  FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
  fRun->SetGenerator(primGen);

  FairEvtGenGenerator* evtGen = new FairEvtGenGenerator(infile.Data());
  primGen->AddGenerator(evtGen);

  primGen->DoTracking(kFALSE);

  
  //Setup the Fast Simulation Task
  //-----------------------------
  PndFastSim* fastSim = new PndFastSim();
  fastSim->SetVerbosity(0);
  fastSim->AddDetector("CmpDet");
  fastSim->EnablePropagation();
  
  fRun->AddTask(fastSim);
  
  fRun->Init();
  if (nevts==0) nevts=10;
  fRun->Run(nevts);             
  
  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
}