QAmacro_fastsim_1.C File Reference

Go to the source code of this file.

Functions
int	QAmacro_fastsim_1 (Int_t nEvents=1000)

Function Documentation

int QAmacro_fastsim_1

(

Int_t

nEvents = 1000

)

Macro for running fast simulation it will only fill the PndData without any transport(M. Al-Turany)

Definition at line 4 of file QAmacro_fastsim_1.C.

References PndFastSim::AddDetector(), Bool_t, boxGen, Cave, ctime, Double_t, PndFastSim::EnablePropagation(), evtGen, fField, fRun, gDebug, kParameterMerged, nEvents, output, primGen, printf(), rtdb, rtime, PndFastSim::SetVerbosity(), timer, and TString.

{
  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;
}