QAmacro_evtgen_2.C

Go to the documentation of this file.

// TEST 1: only check EvtGenDirect
// void qa_evtgendirect(){
int QAmacro_evtgen_2()
{
  TStopwatch timer;
  timer.Start();

  TString BaseDir =  gSystem->Getenv("VMCWORKDIR");
  TString DECAYInput = BaseDir + "/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC";
  TString PDLInput = BaseDir + "/pgenerators/EvtGen/EvtGen/Private/evt.pdl";
  TString USERInput = BaseDir + "/macro/QA/evtgen/test.dec";

  Bool_t fTest=kFALSE;
  
  FairRunSim *fRun = new FairRunSim();

  // set the MC version used
  // ------------------------

  fRun->SetName("TGeant3");
  //fRun->SetName("TGeant4");

  fRun->SetOutputFile("qa_evtgen2.root");

  // Set Material file Name
  //-----------------------
  fRun->SetMaterials("media_pnd.geo");

  // Create and add detectors
  //-------------------------
  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  fRun->AddModule(Cave); 

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

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

  // EvtGen Generator
  PndEvtGenDirect* evtGen = new PndEvtGenDirect("psi(2S)",USERInput,7,123456,DECAYInput,PDLInput);
  evtGen->SetStoreTree();
  primGen->AddGenerator(evtGen);

 
  // -----   Initialize and run   --------------------------------------------
  fRun->Init();
  fRun->Run(10);// just a few events
  // ------------------------------------------------------------------------

  fTest = kTRUE;
        

  if (fTest){
    cout << " Test Passed" << endl;
    cout << " All Ok " << endl;  
  }else{
    cout << " Test Failed" << endl;
    cout << " Not Ok " << endl;         
  }
 timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
  //return 0; 
  return 0;
}