macro/master/sim_complete.C File Reference

Go to the source code of this file.

Functions
int	sim_complete (Int_t nEvents=100, TString SimEngine="TGeant3", Double_t BeamMomentum=6.231552, TString prefix="evtcomplete")

Function Documentation

int sim_complete	(	Int_t	nEvents = 100,
		TString	SimEngine = "TGeant3",
		Double_t	BeamMomentum = 6.231552,
		TString	prefix = "evtcomplete"
	)

Definition at line 8 of file macro/master/sim_complete.C.

References PndMasterRunSim::AddSimTasks(), BeamMomentum, PndMasterRunSim::CreateGeometry(), PndMasterRunSim::Finish(), fRun, PndMasterRunSim::GetFilteredPrimaryGenerator(), nEvents, primGen, PndMasterRunSim::SetGenerator(), PndMasterRunSim::SetInput(), PndMasterRunSim::SetNumberOfEvents(), PndMasterRunSim::SetOptions(), PndMasterRunSim::SetParamAsciiFile(), PndMasterRunSim::Setup(), FairFilteredPrimaryGenerator::SetVerbose(), and TString.

{
  //-----User Settings:------------------------------------------------------
//  gDebug=5;
  TString parAsciiFile   = "all.par";
  
  TString  options        = "";
  // TString inputGenerator = 
  // EvtGen -> "xxxxxxxx.dec" (parses dec-file for initial particle) or "xxxxxxx.dec:initial_particle"
  // DPM    -> "dpm_xxxxx"
  // FTF    -> "ftf_xxxxx"
  // BOX    -> "box:type(pdgcode,mult):p(min,max):tht(min,max):phi(min,max)"
  // PIPI   -> "pipi:cosTheta(min,max)"
  // LEP    -> "leplep:pid(value):gegm(value):cosTheta(min,max)"

  //TString inputGenerator = "psi2s_Jpsi2pi_Jpsi_mumu.dec";
  //TString inputGenerator = "ftf";
  //TString inputGenerator = "ftf";
  //TString inputGenerator = "box:type(211,1):p(1,1):tht(10,120):phi(0,360)";
  TString inputGenerator = "box:type(211,10):p(0.5,1.5):tht(5,90):phi(0,360)";

  //-------------------------------------------------------------------------
  // -----   Create the Simulation run manager ------------------------------
  PndMasterRunSim *fRun = new PndMasterRunSim();
  fRun->SetInput(inputGenerator);
  fRun->SetName(SimEngine);
  fRun->SetParamAsciiFile(parAsciiFile);
  fRun->SetNumberOfEvents(nEvents);
  fRun->SetBeamMom(BeamMomentum);
  fRun->SetStoreTraj(kTRUE);
  fRun->SetOptions(options);
  // -----  Initialization   ------------------------------------------------
  fRun->Setup(prefix);
  // -----   Geometry   -----------------------------------------------------
  fRun->CreateGeometry();
  // -----   Event generator   ----------------------------------------------
  fRun->SetGenerator();

  // -----   Event filter setup   -------------------------------------------
  FairFilteredPrimaryGenerator *primGen = fRun->GetFilteredPrimaryGenerator();
  primGen->SetVerbose(0);
  // ---- Example configuration for the event filter ------------------------
  //FairEvtFilterOnSingleParticleCounts* chrgFilter = new FairEvtFilterOnSingleParticleCounts("chrgFilter");
  //chrgFilter->AndMinCharge(4, FairEvtFilter::kCharged);
  //primGen->AndFilter(chrgFilter);

  // -----   Add tasks   ----------------------------------------------------
  fRun->AddSimTasks();
  // -----   Intialise and run   --------------------------------------------
  fRun->Init();
  fRun->Run(nEvents);
  fRun->Finish();
  return 0;
}