charmonium/jul10/simfast_cmp.C File Reference

Go to the source code of this file.

Functions
void	simfast_cmp (TString part="psi(2S)", TString decfile="PSI2S.DEC", int nevts=1000, double mom=15., TString outfile="")

Function Documentation

void simfast_cmp	(	TString	part = "psi(2S)",
		TString	decfile = "PSI2S.DEC",
		int	nevts = 1000,
		double	mom = 15.,
		TString	outfile = ""
	)

... generate your signal on the fly

Definition at line 6 of file charmonium/jul10/simfast_cmp.C.

References PndFastSim::AddDetector(), ctime, Double_t, PndFastSim::EnablePropagation(), fRun, gDebug, mom, outfile, primGen, printf(), rootlogon(), rtime, PndEvtGenDirect::SetStoreTree(), 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=decfile+".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);

  PndEvtGenDirect *EvtGen = new PndEvtGenDirect(part,decfile,mom);
  EvtGen->SetStoreTree();
  primGen->AddGenerator(EvtGen);
//  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);
}