QAmacro_fastsim_1.C

Go to the documentation of this file.

int QAmacro_fastsim_1(Int_t nEvents = 1000 )
{
  TString BaseDir =  gSystem->Getenv("VMCWORKDIR");

 //-----User Settings:-----------------------------------------------
  gDebug                  = 0;

  // choose your event generator
  Bool_t UseEvtGen            = kTRUE;
  Bool_t UseDpm               = kFALSE;
  Bool_t UseBoxGenerator  = kFALSE;

  TString EvtInput = BaseDir + "/input/psi2s_jpsi2pi_1k.evt"; //  Input EvtGen

  Double_t MomDpm  = 7.24; // pbar momentum for DPM generator; matches psi(2S) energy

  Double_t MomMin  = 0.5;  // minimum momentum for box generator
  Double_t MomMax  = 2.0;  // maximum   "       "


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

  // Create the Simulation run manager--------------------------------
  FairRunSim *fRun = new FairRunSim();
  fRun->SetOutputFile("sim_fast.root");
  fRun->SetBeamMom(MomDpm);
  fRun->SetGenerateRunInfo(kFALSE);
  fRun->SetMaterials("media_pnd.geo");
  FairRuntimeDb *rtdb=fRun->GetRuntimeDb();
  Bool_t kParameterMerged=kTRUE;
  FairParRootFileIo* output=new FairParRootFileIo(kParameterMerged);
  output->open("dummypar.root");
  rtdb->setOutput(output);

  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  fRun->AddModule(Cave);

// Create and Set Event Generator
  //-------------------------------
  FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
  fRun->SetGenerator(primGen);
  fRun->SetName("TGeant3");

  if(UseBoxGenerator){  // Box Generator
     FairBoxGenerator* boxGen = new FairBoxGenerator(211, 5); // 211 = pion; 1 = multipl.
     boxGen->SetPtRange(MomMin,MomMax); // GeV/c
     boxGen->SetPhiRange(0., 360.); // Azimuth angle range [degree]
     boxGen->SetThetaRange(0., 90.); // Polar angle in lab system range [degree]
     boxGen->SetXYZ(0., 0., 0.); // mm o cm ??
     primGen->AddGenerator(boxGen);
  }
  if(UseDpm){
          PndDpmDirect *Dpm= new PndDpmDirect(MomDpm,0);
          primGen->AddGenerator(Dpm);
  }
  if(UseEvtGen){
          FairEvtGenGenerator* evtGen = new FairEvtGenGenerator(EvtInput.Data());
          primGen->AddGenerator(evtGen);
  }

  // ------------- switch off the transport of particles
  primGen->DoTracking(kFALSE);



  PndMultiField *fField= new PndMultiField("FULL");
  fRun->SetField(fField);


  //-------- Setup the Fast Simulation Task  --------------
  //-------------------------------------------------------
  PndFastSim* fastSim = new PndFastSim();
  fastSim->SetVerbosity(0);
  fastSim->AddDetector("CmpDet");
  fastSim->EnablePropagation();

  fRun->AddTask(fastSim);


  //-------------------------  Initialize the RUN  -----------------
  fRun->Init();
  //-------------------------  Run the Simulation  -----------------
  fRun->Run(nEvents);
 //-------------------------  Save the parameters -----------------
  rtdb->saveOutput();
  //------------------------Print some info and exit----------------
  timer.Stop();

  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);

  cout << " Test passed" << endl;
  cout << " All ok " << endl;
  return 0;
}