PndMasterRunSim Class Reference

Class for the master simulation chain. More...

#include <PndMasterRunSim.h>

Inheritance diagram for PndMasterRunSim:

Public Member Functions
	PndMasterRunSim ()
	Default constructor. More...
virtual	~PndMasterRunSim ()
	Default destructor. More...
Bool_t	Setup (TString outprefix="")
	Initial setup. More...
void	Finish ()
	Final diagnostics. More...
void	CreateGeometry ()
	It switches between different standard geometry volumes. More...
void	CreateGeometryDefault ()
	It creates all the standard geometry volumes. More...
void	CreateGeometryPhase1 ()
	It creates the standard geometry volumes for phase1 of the experiment. More...
void	CreateGeometryDay1 ()
	It creates the standard geometry volumes for day1 phase. More...
void	AddSimTasks ()
	Add simulation tasks. More...
void	SetGenerator ()
	Set the event generator. More...
void	AddGenerator (FairGenerator *aGen)
	Add a generator to existing setup. More...
void	SetGenerator (FairGenerator *boxGen)
	Set the event generator for FairBoxGenerator. More...
void	SetGenerator (PndBoxGenerator *boxGen)
	Set the event generator for PndBoxGenerator. More...
void	SetDpmFlag (Int_t Mode)
	Set the DPM flag. More...
void	SetFtfFlag (Int_t Mode)
	Set the FTF noelastic flag. More...
void	UseDpmGenerator ()
	Use DPM as event generator. More...
void	UseFtfGenerator (TString ftfData)
	Use FTF as event generator. More...
void	UseEvtGenGenerator (TString EvtGenFile)
	Use EvtGen as event generator. More...
void	UseBoxGenerator (TString BoxConfig)
	Use BoxGen as event generator. More...
void	UsePiPiGenerator (TString pipiConfig)
	Use PiPiGenerator as event generator. More...
void	UseLepLepGenerator (TString leplepConfig)
	Use LepLepGenerator as event generator. More...
void	SetInput (TString par)
	Input of the simulation This string can be: a) the name of the dec file for EvtGen, ending w/ .dec b) "dpm" if you want to use dpm c) "ftf" if you want to use ftf d) "box:[...]" if you want to use box. More...
void	SetInputDir (TString par)
	Input directory of the simulation. More...
void	SetParamRootFile (TString par)
	Setter of the parameter root file. More...
void	SetParamAsciiFile (TString par)
	Setter of the parameter ascii file. More...
void	SetOptions (TString par)
	Setter of the simulation options This string can be: "" -> default settings full setup "day1" -> Setup for day1 experiment: no GEM, FTS1234, NO DISC, NO RICH "gem2" (added to "day1") -> Setup for day1 experiment with 2 GEM planes (first two) "gem3" (added to "day1") -> Setup for day1 experiment with 3 GEM planes "fts1256" (added to "day1") -> Setup for day1 experiment with FTS1256 insted of FTS1234 "strip" (added to "day1") -> Setup for day1 experiment with only the strip part of the MVD instead of the full MVD Example: "day1+gem+fts1256" means day1 setup + GEM planes + fst1256. More...
TString	GetOptions ()
void	SetNumberOfEvents (Int_t par)
	Setter of the number of events. More...
void	SetEventCounterRate (Int_t par)
	Setter of the event counter rate. More...
void	SetTargetMode (Short_t par)
	Setter of the target mode. More...
FairFilteredPrimaryGenerator *	GetFilteredPrimaryGenerator ()
	Getter for the primary generator, e.g. to configure the event filter. More...
PndFilteredPrimaryGenerator *	GetPndFilteredPrimaryGenerator ()

Private Member Functions
void	GetRange (TString par, double &min, double &max)

Private Attributes
TString	fInput
	Name of the input for the simulation. More...
TString	fInputDir
	Name of the input directory for the simulation. More...
TString	fOutFile
	Name of the output file. More...
TString	fParamRootFile
	Name of the parameter root file. More...
TString	fParamAsciiFile
	Name of the parameter ascii file. More...
TString	fOptions
	Options parsed to the simulation. More...
Int_t	fDpmFlag
	Flag for DPM event generator. More...
Int_t	fFtfFlag
	Flag for FTF event generator. More...
Int_t	fNEvents
	Number of events. More...
Int_t	fEventCounterRate
	After how many events the counter will print. More...
Short_t	fTargetMode
	Target mode. More...
FairRuntimeDb *	fRtdb
	Runtime DB. More...
TStopwatch	fTimer
	Timer. More...

Detailed Description

Class for the master simulation chain.

Master Simulation Class

This class is the basic for the simulation macro. It loads the environment and all the standard detectors.

Author

Stefano Spataro spata.nosp@m.ro@t.nosp@m.o.inf.nosp@m.n.it, Torino University

Version

1.0

Date

Feb 3, 2016

Definition at line 31 of file PndMasterRunSim.h.

Constructor & Destructor Documentation

PndMasterRunSim::PndMasterRunSim

(

)

Default constructor.

Definition at line 54 of file PndMasterRunSim.cxx.

References fTimer.

                                 :
  FairRunSim(), fInput(), fInputDir(""), fOutFile(), fParamRootFile(), fParamAsciiFile(), fOptions(), fDpmFlag(1), fFtfFlag(0), fNEvents(0), fEventCounterRate(100), fTargetMode(0), fRtdb(), fTimer()
{
  fTimer.Start();
}

PndMasterRunSim::~PndMasterRunSim ( )

virtual

Default destructor.

Definition at line 60 of file PndMasterRunSim.cxx.

References gGeoManager.

{
  if (gROOT->GetVersionInt() >= 60602 && gGeoManager!=NULL) {
    gGeoManager->GetListOfVolumes()->Delete();
    gGeoManager->GetListOfShapes()->Delete();
    delete gGeoManager;
  }
}

Member Function Documentation

void PndMasterRunSim::AddGenerator

(

FairGenerator *

aGen

)

Add a generator to existing setup.

Adds a generator to the existing one. BOTH generators will be active

and produce their particle content in each even simultaneously.

void PndMasterRunSim::AddSimTasks

(

)

Add simulation tasks.

Add Master simulation tasks

It adds all the standard simulation tasks

Definition at line 380 of file PndMasterRunSim.cxx.

References AddTask(), fEventCounterRate, fNEvents, fOptions, and sim().

Referenced by prod_sim(), sim_box(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_mvd(), sim_radlength_complete(), and tut_sim().

{
  // -----   Event Counter   --------------------------------
  AddTask(new PndEventCounterTask("Event Counter", fNEvents, fEventCounterRate));

  PndMasterSimTask *sim = new PndMasterSimTask(fOptions);
  AddTask(sim);

  return;
}

void PndMasterRunSim::CreateGeometry

(

)

It switches between different standard geometry volumes.

Master Geometry List

According to fOptions, it creates the standard geometry volumes for the full setup, or dedicated geometries (such as day1)

Definition at line 126 of file PndMasterRunSim.cxx.

References CreateGeometryDay1(), CreateGeometryDefault(), CreateGeometryPhase1(), and fOptions.

Referenced by prod_sim(), sim_box(), sim_complete(), sim_day1(), and tut_sim().

{
  if (fOptions.Contains("phase1")) CreateGeometryPhase1();
  else if (fOptions.Contains("day1")) CreateGeometryDay1();
  else CreateGeometryDefault();
}

void PndMasterRunSim::CreateGeometryDay1

(

)

It creates the standard geometry volumes for day1 phase.

Master Geometry List

It creates all the standard geometry volumes which have to be used in simulation, with the setup for day1 experiments. All the MCPoint will be stored, except for EMC

Definition at line 284 of file PndMasterRunSim.cxx.

References AddModule(), Cave, Dipole, Drc, Emc, fOptions, FTof, Fts, Gem, Muo, Mvd, Pipe, PndMdt::SetBarrel(), PndMdt::SetEndcap(), PndMdt::SetForward(), PndEmc::SetGeometryVersion(), PndMdt::SetMdtCoil(), PndMdt::SetMdtMagnet(), PndMdt::SetMdtMFIron(), PndMdt::SetMuonFilter(), PndDrc::SetRunCherenkov(), PndEmc::SetStorageOfData(), and Stt.

Referenced by CreateGeometry().

{
  //-------------------------  CAVE      -----------------
  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  AddModule(Cave);
  //-------------------------  Magnet   -----------------
  // This part is commented because the MDT geometry contains the magnet now
  //FairModule *Magnet= new PndMagnet("MAGNET");
  //Magnet->SetGeometryFileName("FullSolenoid_V842.root");
  //Magnet->SetGeometryFileName("FullSuperconductingSolenoid_v831.root");
  //AddModule(Magnet);
  FairModule *Dipole= new PndMagnet("MAGNET");
  Dipole->SetGeometryFileName("dipole.geo");
  AddModule(Dipole);
  //-------------------------  Pipe     -----------------
  FairModule *Pipe= new PndPipe("PIPE");
  Pipe->SetGeometryFileName("beampipe_201309.root");
  AddModule(Pipe);
  //-------------------------  STT       -----------------
  FairDetector *Stt= new PndStt("STT", kTRUE);
  Stt->SetGeometryFileName("straws_skewed_blocks_35cm_pipe.geo");
  AddModule(Stt);
  //-------------------------  MVD       -----------------
  if (fOptions.Contains("strip")||fOptions.Contains("nopixels")){
          FairDetector *Mvd = new PndMvdDetector("MVD", kTRUE);
          Mvd->SetGeometryFileName("Mvd-2.1-Strips.root");
          AddModule(Mvd);
  } else {
          FairDetector *Mvd = new PndMvdDetector("MVD", kTRUE);
          Mvd->SetGeometryFileName("Mvd-2.1_FullVersion.root");
          AddModule(Mvd);
  }
  //-------------------------  EMC       -----------------
  PndEmc *Emc = new PndEmc("EMC",kTRUE);
  Emc->SetGeometryVersion(1);
  Emc->SetStorageOfData(kFALSE);
  AddModule(Emc);
  //-------------------------  SCITIL    -----------------
  FairDetector *SciT = new PndSciT("SCIT",kTRUE);
  SciT->SetGeometryFileName("SciTil_201601.root");
  AddModule(SciT);
  //-------------------------  DRC       -----------------
  PndDrc *Drc = new PndDrc("DIRC", kTRUE);
  Drc->SetGeometryFileName("dirc_e3_b3_l6_m40.root");
  Drc->SetRunCherenkov(kFALSE);
  AddModule(Drc);
  //-------------------------  MDT       -----------------
  PndMdt *Muo = new PndMdt("MDT",kTRUE);
  Muo->SetBarrel("fast");
  Muo->SetEndcap("fast");
  Muo->SetMuonFilter("fast");
  Muo->SetForward("fast");
  Muo->SetMdtMagnet(kTRUE);
  Muo->SetMdtCoil(kTRUE);
  Muo->SetMdtMFIron(kTRUE);
  AddModule(Muo);
  //-------------------------  FTOF      -----------------
  FairDetector *FTof = new PndFtof("FTOF",kTRUE);
  FTof->SetGeometryFileName("ftofwall.root");
  AddModule(FTof);

  if (fOptions.Contains("nogem")||fOptions.Contains("gem0")) {
    // do nothing
  }
  else if (fOptions.Contains("gem3")) // GEM 3 Stations
  {
          FairDetector *Gem = new PndGemDetector("GEM", kTRUE);
          Gem->SetGeometryFileName("gem_3Stations_realistic_v2.root");
          AddModule(Gem);
  }
  else   //GEM 2 Stations
  {
    //-------------------------  GEM       -----------------
    FairDetector *Gem = new PndGemDetector("GEM", kTRUE);
    Gem->SetGeometryFileName("gem_2Stations_realistic_v2.root");
    AddModule(Gem);
  }

  if (fOptions.Contains("fts1256"))
    {
      //-------------------------  FTS       -----------------
      FairDetector *Fts= new PndFts("FTS", kTRUE);
      Fts->SetGeometryFileName("fts_1256.geo");
      AddModule(Fts);
    }
  else
    {
      //-------------------------  FTS       -----------------
      FairDetector *Fts= new PndFts("FTS", kTRUE);
      Fts->SetGeometryFileName("fts_reduced.geo");
      AddModule(Fts);
    }
}

void PndMasterRunSim::CreateGeometryDefault

(

)

It creates all the standard geometry volumes.

Master Geometry List

It creates all the standard geometry volumes which have to be used in simulation. All the MCPoint will be stored, expect for EMC.

Definition at line 134 of file PndMasterRunSim.cxx.

References AddModule(), Cave, Dipole, Drc, Dsk, Emc, FTof, Fts, Gem, Muo, Mvd, Pipe, PndMdt::SetBarrel(), PndMdt::SetEndcap(), PndMdt::SetForward(), PndEmc::SetGeometryVersion(), PndMdt::SetMdtCoil(), PndMdt::SetMdtMagnet(), PndMdt::SetMdtMFIron(), PndMdt::SetMuonFilter(), PndDrc::SetRunCherenkov(), PndEmc::SetStorageOfData(), PndDsk::SetStoreCerenkovs(), PndDsk::SetStoreTrackPoints(), and Stt.

Referenced by CreateGeometry().

{
  //-------------------------  CAVE      -----------------
  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  AddModule(Cave);
  //-------------------------  Magnet   -----------------
  // This part is commented because the MDT geometry contains the magnet now
  //FairModule *Magnet= new PndMagnet("MAGNET");
  //Magnet->SetGeometryFileName("FullSolenoid_V842.root");
  //Magnet->SetGeometryFileName("FullSuperconductingSolenoid_v831.root");
  //AddModule(Magnet);
  FairModule *Dipole= new PndMagnet("MAGNET");
  Dipole->SetGeometryFileName("dipole.geo");
  AddModule(Dipole);
  //-------------------------  Pipe     -----------------
  FairModule *Pipe= new PndPipe("PIPE");
  Pipe->SetGeometryFileName("beampipe_201309.root");
  AddModule(Pipe);
  //-------------------------  STT       -----------------
  FairDetector *Stt= new PndStt("STT", kTRUE);
  Stt->SetGeometryFileName("straws_skewed_blocks_35cm_pipe.geo");
  AddModule(Stt);
  //-------------------------  MVD       -----------------
  FairDetector *Mvd = new PndMvdDetector("MVD", kTRUE);
  Mvd->SetGeometryFileName("Mvd-2.1_FullVersion.root");
  AddModule(Mvd);
  //-------------------------  GEM       -----------------
  FairDetector *Gem = new PndGemDetector("GEM", kTRUE);
  Gem->SetGeometryFileName("gem_3Stations_realistic_v2.root");
  AddModule(Gem);
  //-------------------------  EMC       -----------------
  PndEmc *Emc = new PndEmc("EMC",kTRUE);
  Emc->SetGeometryVersion(1);
  Emc->SetStorageOfData(kFALSE);
  AddModule(Emc);
  //-------------------------  SCITIL    -----------------
  FairDetector *SciT = new PndSciT("SCIT",kTRUE);
  SciT->SetGeometryFileName("SciTil_201601.root");
  AddModule(SciT);
  //-------------------------  DRC       -----------------
  PndDrc *Drc = new PndDrc("DIRC", kTRUE);
  Drc->SetGeometryFileName("dirc_e3_b3_l6_m40.root");
  Drc->SetRunCherenkov(kFALSE);
  AddModule(Drc);
  //-------------------------  DISC      -----------------
  PndDsk* Dsk = new PndDsk("DSK", kTRUE);
  Dsk->SetStoreCerenkovs(kFALSE);
  Dsk->SetStoreTrackPoints(kFALSE);
  AddModule(Dsk);
  //-------------------------  MDT       -----------------
  PndMdt *Muo = new PndMdt("MDT",kTRUE);
  Muo->SetBarrel("fast");
  Muo->SetEndcap("fast");
  Muo->SetMuonFilter("fast");
  Muo->SetForward("fast");
  Muo->SetMdtMagnet(kTRUE);
  Muo->SetMdtCoil(kTRUE);
  Muo->SetMdtMFIron(kTRUE);
  AddModule(Muo);
  //-------------------------  FTS       -----------------
  FairDetector *Fts= new PndFts("FTS", kTRUE);
  Fts->SetGeometryFileName("fts.geo");
  AddModule(Fts);
  //-------------------------  FTOF      -----------------
  FairDetector *FTof = new PndFtof("FTOF",kTRUE);
  FTof->SetGeometryFileName("ftofwall.root");
  AddModule(FTof);
  //-------------------------  RICH       ----------------
  PndRich *Rich= new PndRich("RICH",kTRUE);
  Rich->SetGeometryFileName("rich_v313.root");
  AddModule(Rich);
}

void PndMasterRunSim::CreateGeometryPhase1

(

)

It creates the standard geometry volumes for phase1 of the experiment.

Master Geometry List

It creates all the standard geometry volumes which have to be used in simulation, with the setup for phase1 experiments. All the MCPoint will be stored, except for EMC

Definition at line 209 of file PndMasterRunSim.cxx.

References AddModule(), Cave, Dipole, Drc, Emc, FTof, Fts, Gem, Muo, Mvd, Pipe, PndMdt::SetBarrel(), PndMdt::SetEndcap(), PndMdt::SetForward(), PndEmc::SetGeometryVersion(), PndMdt::SetMdtCoil(), PndMdt::SetMdtMagnet(), PndMdt::SetMdtMFIron(), PndMdt::SetMuonFilter(), PndDrc::SetRunCherenkov(), PndEmc::SetStorageOfData(), and Stt.

Referenced by CreateGeometry().

{
  //-------------------------  CAVE      -----------------
  FairModule *Cave= new PndCave("CAVE");
  Cave->SetGeometryFileName("pndcave.geo");
  AddModule(Cave);
  //-------------------------  Magnet   -----------------
  // This part is commented because the MDT geometry contains the magnet now
  //FairModule *Magnet= new PndMagnet("MAGNET");
  //Magnet->SetGeometryFileName("FullSolenoid_V842.root");
  //Magnet->SetGeometryFileName("FullSuperconductingSolenoid_v831.root");
  //AddModule(Magnet);
  FairModule *Dipole= new PndMagnet("MAGNET");
  Dipole->SetGeometryFileName("dipole.geo");
  AddModule(Dipole);
  //-------------------------  Pipe     -----------------
  FairModule *Pipe= new PndPipe("PIPE");
  Pipe->SetGeometryFileName("beampipe_201309.root");
  AddModule(Pipe);
  //-------------------------  STT       -----------------
  FairDetector *Stt= new PndStt("STT", kTRUE);
  Stt->SetGeometryFileName("straws_skewed_blocks_35cm_pipe.geo");
  AddModule(Stt);
  //-------------------------  MVD       -----------------
  FairDetector *Mvd = new PndMvdDetector("MVD", kTRUE);
  Mvd->SetGeometryFileName("Mvd-2.1_FullVersion.root");
  AddModule(Mvd);
  //-------------------------  GEM       -----------------
  FairDetector *Gem = new PndGemDetector("GEM", kTRUE);
  Gem->SetGeometryFileName("gem_3Stations_realistic_v2.root");
  AddModule(Gem);
  //-------------------------  EMC       -----------------
  PndEmc *Emc = new PndEmc("EMC",kTRUE);
  Emc->SetGeometryVersion(1);
  Emc->SetStorageOfData(kFALSE);
  AddModule(Emc);
  //-------------------------  SCITIL    -----------------
  FairDetector *SciT = new PndSciT("SCIT",kTRUE);
  SciT->SetGeometryFileName("SciTil_201601.root");
  AddModule(SciT);
  //-------------------------  DRC       -----------------
  PndDrc *Drc = new PndDrc("DIRC", kTRUE);
  Drc->SetGeometryFileName("dirc_e3_b3_l6_m40.root");
  Drc->SetRunCherenkov(kFALSE);
  AddModule(Drc);
  //-------------------------  DISC      -----------------
  //PndDsk* Dsk = new PndDsk("DSK", kTRUE);
  //Dsk->SetStoreCerenkovs(kFALSE);
  //Dsk->SetStoreTrackPoints(kFALSE);
  //AddModule(Dsk);
  //-------------------------  MDT       -----------------
  PndMdt *Muo = new PndMdt("MDT",kTRUE);
  Muo->SetBarrel("fast");
  Muo->SetEndcap("fast");
  Muo->SetMuonFilter("fast");
  Muo->SetForward("fast");
  Muo->SetMdtMagnet(kTRUE);
  Muo->SetMdtCoil(kTRUE);
  Muo->SetMdtMFIron(kTRUE);
  AddModule(Muo);
  //-------------------------  FTS       -----------------
  FairDetector *Fts= new PndFts("FTS", kTRUE);
  Fts->SetGeometryFileName("fts.geo");
  AddModule(Fts);
  //-------------------------  FTOF      -----------------
  FairDetector *FTof = new PndFtof("FTOF",kTRUE);
  FTof->SetGeometryFileName("ftofwall.root");
  AddModule(FTof);
  //-------------------------  RICH       ----------------
  //PndRich *Rich= new PndRich("RICH",kTRUE);
  //Rich->SetGeometryFileName("rich_v313.root");
  //AddModule(Rich);
}

void PndMasterRunSim::Finish

(

)

Final diagnostics.

Master Final diagnostics

It prints CPU time, memory usage, used parameters, and eventually send the information to CDash

Definition at line 728 of file PndMasterRunSim.cxx.

References ctime, Double_t, fOptions, fOutFile, fParamAsciiFile, fParamRootFile, fRtdb, fTimer, gGeoManager, outfile, and rtime.

Referenced by prod_sim(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_mvd(), sim_radlength_complete(), and tut_sim().

{
  fRtdb->saveOutput();

  cout<<"PndMasterRunAna::Finish(): Tasks that ran just now:"<<endl;
  TFile* outfile=fRootManager->GetOutFile();
  bool wasopen=outfile->IsOpen ();
  if (!wasopen)
  {
    cout<<"file is "<< ((wasopen) ? "" : "not " ) <<"open" <<endl;
    outfile=new TFile(outfile->GetName(),"UPDATE");
  }
  outfile->cd();

  // write the summary of event filter to output root file
  if (!strcmp(fGen->ClassName(),"PndFilteredPrimaryGenerator"))
        ((PndFilteredPrimaryGenerator*)fGen)->WriteEvtFilterStatsToRootFile(outfile);
  else
        ((FairFilteredPrimaryGenerator*)fGen)->WriteEvtFilterStatsToRootFile(outfile);

  FairFileHeader* outheader=(FairFileHeader*)outfile->Get("FileHeader");
  cout<<"Task tha ran just now:"<<endl;
  for(const auto&& os : *(outheader->GetListOfTasks()) ) cout<<" - "<<((TObjString*)os)->GetString().Data()<<endl;

  TObjString outoptions(fOptions);
  outoptions.Write("PndOptions",kOverwrite);

  outfile->Write();
  if(!wasopen) outfile->Close();

  cout << endl;
  if (gROOT->GetVersionInt() >= 60602) {
    gGeoManager->GetListOfVolumes()->Delete();
    gGeoManager->GetListOfShapes()->Delete();
    delete gGeoManager;
  }

  // Extract the maximal used memory an add is as Dart measurement
  // This line is filtered by CTest and the value send to CDash
  FairSystemInfo sysInfo;
  Float_t maxMemory=sysInfo.GetMaxMemory();
  cout << "<DartMeasurement name=\"MaxMemory\" type=\"numeric/double\">";
  cout << maxMemory;
  cout << "</DartMeasurement>" << endl;

  fTimer.Stop();
  Double_t rtime = fTimer.RealTime();
  Double_t ctime = fTimer.CpuTime();

  Float_t cpuUsage=ctime/rtime;
  cout << "<DartMeasurement name=\"CpuLoad\" type=\"numeric/double\">";
  cout << cpuUsage;
  cout << "</DartMeasurement>" << endl;

  cout << endl;
  cout << "Output file is\t\t"    << fOutFile << endl;
  cout << "Parameter ROOT file is\t" << fParamRootFile << endl;
  cout << "Parameter ASCII file is\t" << fParamAsciiFile << endl;
  cout << "Real time " << rtime << " s, CPU time " << ctime
            << "s" << endl;
  cout << "CPU usage " << cpuUsage*100. << "%" << endl;
  cout << "Max Memory " << maxMemory << " MB" << endl;

  cout << "Macro finished successfully." << endl;

}

FairFilteredPrimaryGenerator* PndMasterRunSim::GetFilteredPrimaryGenerator ( )

inline

Getter for the primary generator, e.g. to configure the event filter.

Definition at line 273 of file PndMasterRunSim.h.

Referenced by sim_complete(), sim_complete_newSTT(), sim_day1(), sim_radlength_complete(), and tut_sim().

{return (FairFilteredPrimaryGenerator*)fGen;}

TString PndMasterRunSim::GetOptions ( )

inline

Definition at line 248 of file PndMasterRunSim.h.

References fOptions.

{ return fOptions;}

PndFilteredPrimaryGenerator* PndMasterRunSim::GetPndFilteredPrimaryGenerator ( )

inline

Definition at line 274 of file PndMasterRunSim.h.

Referenced by prod_sim().

{return (PndFilteredPrimaryGenerator*)fGen;}

void PndMasterRunSim::GetRange ( TString  par, double &  min, double &  max  )

private

Definition at line 798 of file PndMasterRunSim.cxx.

References par, and TString.

Referenced by UseBoxGenerator(), UseLepLepGenerator(), and UsePiPiGenerator().

{
        par.ReplaceAll(" ","");
        par = par(par.Index("(")+1, par.Length()-par.Index("(")-2);

        TString smin=par, smax=par;

        if (par.Contains(","))
        {
                smin = par(0,par.Index(","));
                smax = par(par.Index(",")+1,1000);
        }

        min = smin.Atof();
        max = smax.Atof();
}

void PndMasterRunSim::SetDpmFlag ( Int_t  Mode )

inline

Set the DPM flag.

Parameters

Mode	= 0. - DPM - No elastic scattering, only inelastic
Mode	= 1. - DPM - Elastic and inelastic interactions (default)
Mode	= 2. - DPM - Only elastic scattering, no inelastic one

Definition at line 145 of file PndMasterRunSim.h.

References fDpmFlag.

Referenced by prod_sim().

{ fDpmFlag = Mode; };

void PndMasterRunSim::SetEventCounterRate ( Int_t  par )

inline

Setter of the event counter rate.

Definition at line 258 of file PndMasterRunSim.h.

References fEventCounterRate, and par.

{ fEventCounterRate = par;}

void PndMasterRunSim::SetFtfFlag ( Int_t  Mode )

inline

Set the FTF noelastic flag.

Parameters

Mode	= 0. - FTF - Elastic and inelastic interactions (default)
Mode	= 1. - FTF - No elastic scattering, only inelastic

Definition at line 152 of file PndMasterRunSim.h.

References fFtfFlag.

Referenced by prod_sim().

{ fFtfFlag = Mode; };

void PndMasterRunSim::SetGenerator

(

)

Set the event generator.

Master event generator

This call set the event generator according to the input name. If the input name contains "dpm" it uses dpm, if it contains "ftf" then ftf, if ".dec" it runs evtgen using the input name as namefile of the .dec file. If the input name contains "box" the macro breaks, since in that case the SetGenerator(PndBoxGenerator *boxGen) function must be used.

Definition at line 392 of file PndMasterRunSim.cxx.

References fInput, fOptions, fTargetMode, TString, UseBoxGenerator(), UseDpmGenerator(), UseEvtGenGenerator(), UseFtfGenerator(), and UsePiPiGenerator().

Referenced by prod_sim(), sim_box(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_mvd(), sim_radlength_complete(), and tut_sim().

{
  if (fOptions.Contains("pndfiltprim"))
        fGen = new PndFilteredPrimaryGenerator();
  else
        fGen = new FairFilteredPrimaryGenerator();

  switch (fTargetMode)
    {
    case 0:
      LOG(INFO) << "Using no Vertex smearing" << FairLogger::endl;
      break;
    case 1:
      LOG(INFO) << "Using Cluster Jet Target" << FairLogger::endl;
      // a cluster-jet beam at the interaction zone with a horizontal width
      // of e.g. dx = 1 mm and a length in accelerator beam direction of dz = 10 mm.
      // Target TDR, page 44
      fGen->SetTarget(0., 1./2.355); // From FWHM to sigma
      fGen->SmearVertexZ(kTRUE);
      fGen->SmearGausVertexZ(kTRUE);
      fGen->SetBeam(0., 0., 0.1, 0.1);
      fGen->SmearVertexXY(kTRUE);
      break;
    case 2:
      LOG(INFO) << "Using Pellet Target" << FairLogger::endl;
      // At PANDA a beam diameter around 3mm is needed and one may want even
      // smaller size when PTR is not possible.
      // Target TDR, page 61
      fGen->SetTarget(0., 0.3);
      fGen->SmearVertexZ(kTRUE);
      fGen->SmearGausVertexZ(kTRUE);
      fGen->SetBeam(0., 0., 0.3, 0.3);
      fGen->SmearVertexXY(kTRUE);
      break;
    case 3:
      LOG(INFO) << "Using Pellet Tracking Target" << FairLogger::endl;
      // A position resolution σ(x, y, z) < 0.2 mm in the in- teraction position
      // is desirable for event reconstruc- tion.
      // Target TDR, page 61
      fGen->SetTarget(0., 0.02);
      fGen->SmearVertexZ(kTRUE);
      fGen->SmearGausVertexZ(kTRUE);
      fGen->SetBeam(0., 0., 0.02, 0.02);
      fGen->SmearVertexXY(kTRUE);
      break;
    default:
      LOG(INFO) << "Unkwown target mode - Using no vertex smearing" << FairLogger::endl;
    }

  TString input = fInput;
  input.ToLower();

  if (input.EndsWith(".dec") || input.Contains(".dec:"))
    {
      UseEvtGenGenerator(fInput);
    }
  else if (input.BeginsWith("dpm"))
    {
      UseDpmGenerator();
    }
  else if (input.BeginsWith("ftf"))
    {
      UseFtfGenerator(fInput);
    }
  else if (input.BeginsWith("box"))
    {
      UseBoxGenerator(fInput);
    }
  else if (input.BeginsWith("pipi"))
  {
          UsePiPiGenerator(fInput);
  }
  else
    {
      LOG(WARNING)<< "Generator could not be identified from input '"<<fInput.Data()<<"'!!" <<  FairLogger::endl;
    }

}

void PndMasterRunSim::SetGenerator

(

FairGenerator *

boxGen

)

Set the event generator for FairBoxGenerator.

Master event generator for FairBoxGenerator

This call set the FairBoxGenerator as event generator. The user should create a FairBoxGenerator object with all the settings, and pass it as argument to the function.

Definition at line 631 of file PndMasterRunSim.cxx.

References fOptions.

 {
  LOG(INFO) << "Using FairBoxGenerator generator" << FairLogger::endl;
  if (fOptions.Contains("PndFiltPrim"))
        fGen = new PndFilteredPrimaryGenerator();
  else
        fGen = new FairFilteredPrimaryGenerator();
  fGen->AddGenerator(boxGen);
}

void PndMasterRunSim::SetGenerator

(

PndBoxGenerator *

boxGen

)

Set the event generator for PndBoxGenerator.

Master event generator for PndBoxGenerator

This call set the PndBoxGenerator as event generator. The user should create a PndBoxGenerator object with all the settings, and pass it as argument to the function.

Definition at line 620 of file PndMasterRunSim.cxx.

References fOptions.

{
  LOG(INFO) << "Using PndBoxGenerator generator" << FairLogger::endl;
  if (fOptions.Contains("PndFiltPrim"))
        fGen = new PndFilteredPrimaryGenerator();
  else
        fGen = new FairFilteredPrimaryGenerator();
  fGen->AddGenerator(boxGen);
}

void PndMasterRunSim::SetInput ( TString  par )

inline

Input of the simulation This string can be: a) the name of the dec file for EvtGen, ending w/ .dec b) "dpm" if you want to use dpm c) "ftf" if you want to use ftf d) "box:[...]" if you want to use box.

Definition at line 219 of file PndMasterRunSim.h.

References fInput, and par.

Referenced by prod_sim(), sim_box(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_radlength_complete(), and tut_sim().

{ fInput          = par;}

void PndMasterRunSim::SetInputDir ( TString  par )

inline

Input directory of the simulation.

Definition at line 224 of file PndMasterRunSim.h.

References fInputDir, and par.

Referenced by sim_complete(), and sim_day1().

{ fInputDir          = par;}

void PndMasterRunSim::SetNumberOfEvents ( Int_t  par )

inline

Setter of the number of events.

Definition at line 253 of file PndMasterRunSim.h.

References fNEvents, and par.

Referenced by prod_sim(), sim_box(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_mvd(), sim_radlength_complete(), and tut_sim().

{ fNEvents = par;}

void PndMasterRunSim::SetOptions ( TString  par )

inline

Setter of the simulation options This string can be: "" -> default settings full setup "day1" -> Setup for day1 experiment: no GEM, FTS1234, NO DISC, NO RICH "gem2" (added to "day1") -> Setup for day1 experiment with 2 GEM planes (first two) "gem3" (added to "day1") -> Setup for day1 experiment with 3 GEM planes "fts1256" (added to "day1") -> Setup for day1 experiment with FTS1256 insted of FTS1234 "strip" (added to "day1") -> Setup for day1 experiment with only the strip part of the MVD instead of the full MVD Example: "day1+gem+fts1256" means day1 setup + GEM planes + fst1256.

Definition at line 247 of file PndMasterRunSim.h.

References fOptions, and par.

Referenced by prod_sim(), sim_complete(), and sim_day1().

{ fOptions = par; fOptions.ToLower();}

void PndMasterRunSim::SetParamAsciiFile ( TString  par )

inline

Setter of the parameter ascii file.

Definition at line 234 of file PndMasterRunSim.h.

References fParamAsciiFile, and par.

Referenced by prod_sim(), sim_box(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_mvd(), sim_radlength_complete(), and tut_sim().

{ fParamAsciiFile = par;}

void PndMasterRunSim::SetParamRootFile ( TString  par )

inline

Setter of the parameter root file.

Definition at line 229 of file PndMasterRunSim.h.

References fParamRootFile, and par.

Referenced by Setup().

{ fParamRootFile  = par;}

void PndMasterRunSim::SetTargetMode ( Short_t  par )

inline

Setter of the target mode.

Target mode 0 - No IP smearing (default) 1 - Cluster Jet 2 - Pellet target 3 - Pellet Tracking target

Definition at line 268 of file PndMasterRunSim.h.

References fTargetMode, and par.

{ fTargetMode = par;}

Bool_t PndMasterRunSim::Setup

(

TString

outprefix = ""

)

Initial setup.

Master Inital setup

This command set the source files, load the proper parameters, and set the relevant flags. If something fails, it returns a kFALSE value.

Definition at line 70 of file PndMasterRunSim.cxx.

References allDigiFile, Bool_t, creator, Double_t, fIni, fInput, fOptions, fOutFile, fParamAsciiFile, fParamRootFile, fRtdb, fTargetMode, kParameterMerged, mom, parIo1, parOutput, SetField(), SetMaterials(), SetOutputFile(), SetParamRootFile(), SetUseFairLinks(), sqrt(), and TString.

Referenced by prod_sim(), sim_box(), sim_complete(), sim_complete_newSTT(), sim_day1(), sim_mvd(), sim_radlength_complete(), and tut_sim().

{
  TString inputName = outprefix;

  // If no prefix is given, we create one from fInput and force lower-case
  if (inputName=="")
  {
    inputName = fInput;
    inputName.ToLower();
  }

  if (inputName.EndsWith(".dec")) inputName.Remove(inputName.Length()-4,4);
  inputName.ReplaceAll(":","_");

  PndFileNameCreator creator(inputName.Data());
  SetOutputFile(creator.GetSimFileName().data());
  fOutFile = creator.GetSimFileName().data();
  SetParamRootFile(creator.GetParFileName().data());
  SetMaterials("media_pnd.geo");
  SetGenerateRunInfo(kFALSE);
  SetUseFairLinks(kTRUE);

  // -----  Parameter database   --------------------------------------------
  TString allDigiFile = gSystem->Getenv("VMCWORKDIR");
  allDigiFile += "/macro/params/";
  allDigiFile += fParamAsciiFile;

  fRtdb = this->GetRuntimeDb();
  Bool_t kParameterMerged=kFALSE; // No use until now
  FairParRootFileIo* parOutput = new FairParRootFileIo(kParameterMerged);
  parOutput->open(fParamRootFile,"RECREATE");

  FairParAsciiFileIo* parIo1 = new FairParAsciiFileIo();
  parIo1->open(allDigiFile.Data(),"in");

  fRtdb->setFirstInput(parIo1);
  fRtdb->setOutput(parOutput);

  // -----  Create and Set the Field(s) ------------------------------------
  PndMultiField *field= new PndMultiField("AUTO");
  SetField(field);

  // ---- Defining PANDA particles -----------------------------------------
  Double_t mom = GetBeamMom();
  TLorentzVector fIni(0, 0, mom, sqrt(mom*mom+9.3827203e-01*9.3827203e-01)+9.3827203e-01);
  TDatabasePDG::Instance()->AddParticle("ppSystem" ,"ppSystem",        fIni.M(), kFALSE, 0.1, 6, "", 98888);
  TDatabasePDG::Instance()->AddParticle("pbarpSystem" ,"pbarpSystem",  fIni.M(), kFALSE, 0.1, 0, "", 88888);
  TDatabasePDG::Instance()->AddParticle("pbarpSystem0","pbarpSystem0", fIni.M(), kFALSE, 0.1, 0, "", 88880);
  TDatabasePDG::Instance()->AddParticle("pbarpSystem1","pbarpSystem1", fIni.M(), kFALSE, 0.1, 0, "", 88881);
  TDatabasePDG::Instance()->AddParticle("pbarpSystem2","pbarpSystem2", fIni.M(), kFALSE, 0.1, 0, "", 88882);

  if (fOptions.Contains("day1")) fTargetMode = 1;
  return kTRUE;
}

void PndMasterRunSim::UseBoxGenerator

(

TString

BoxConfig

)

Use BoxGen as event generator.

Box event generator

This call sets BoxGenerator as event generator. The format of the config string is for isotrop events in theta: 'BOX:type(pdg,mult):p(min,max):phi(min,max):tht(min,max)' for isotrop events in cos(theta): 'BOX:type(pdg,mult):p(min,max):phi(min,max):ctht(min,max)' Instead of range 'var(min,max)' also a fixed value can be set with 'var(value)' All variables left out are set to defaults.

Parameters

BoxConfig

configuration string of the BOX generator

Definition at line 471 of file PndMasterRunSim.cxx.

References Bool_t, boxGen, Double_t, GetRange(), mult, PndBoxGenerator::SetCosTheta(), PndBoxGenerator::SetDebug(), PndBoxGenerator::SetPhiRange(), PndBoxGenerator::SetPRange(), PndBoxGenerator::SetPtRange(), PndBoxGenerator::SetThetaRange(), PndBoxGenerator::SetXYZ(), and TString.

Referenced by SetGenerator().

{
  // use BOX generator; defaults

  Double_t BoxMomMin  = 0.05;   // minimum momentum for box generator
  Double_t BoxMomMax  = 10.;    // maximum   "       "
  Double_t BoxThtMin  = 0. ;    // minimum theta for box generator
  Double_t BoxThtMax  = 180.;   // maximum   "       "
  Double_t BoxPhiMin  = 0. ;    // minimum phi for box generator
  Double_t BoxPhiMax  = 360.;   // maximum   "       "
  Bool_t   BoxCosTht  = false;  // isotropic in cos(theta) instead theta
  Bool_t   BoxPt          = false;  // is pt given instead of p

  Int_t    BoxType    = 13;     // default particle muon
  Int_t    BoxMult    = 1;      // default particle multiplicity
  Double_t type=0,mult=0;       // ref. parameters for range function

  BoxConfig.ToLower();

  if (BoxConfig!="box")
  {
    BoxConfig.ReplaceAll("box","");
    BoxConfig.ReplaceAll(" ","");
    BoxConfig += ":";

    while (BoxConfig.Contains(":"))
    {
      TString curpar = BoxConfig(0,BoxConfig.Index(":"));
      BoxConfig = BoxConfig(BoxConfig.Index(":")+1,1000);
      curpar.ReplaceAll("[","(");
      curpar.ReplaceAll("]",")");

      if (curpar.BeginsWith("type(")) {GetRange(curpar,type,mult); BoxType = (Int_t)type; BoxMult = (Int_t)mult; }
      if (curpar.BeginsWith("p("))    GetRange(curpar,BoxMomMin,BoxMomMax);
      if (curpar.BeginsWith("pt("))    {GetRange(curpar,BoxMomMin,BoxMomMax); BoxPt = true;}
      if (curpar.BeginsWith("tht("))   GetRange(curpar,BoxThtMin,BoxThtMax);
      if (curpar.BeginsWith("ctht(")) {GetRange(curpar,BoxThtMin,BoxThtMax); BoxCosTht=true;}
      if (curpar.BeginsWith("phi("))   GetRange(curpar,BoxPhiMin,BoxPhiMax);
    }
  }

  PndBoxGenerator* boxGen = new PndBoxGenerator(BoxType, BoxMult);
  boxGen->SetDebug(0);

  if (BoxPt == true){
          boxGen->SetPtRange(BoxMomMin, BoxMomMax);
  } else {
          boxGen->SetPRange(BoxMomMin,BoxMomMax);      // GeV/c
  }
  boxGen->SetPhiRange(BoxPhiMin, BoxPhiMax);   // Azimuth angle range [degree]
  boxGen->SetThetaRange(BoxThtMin, BoxThtMax); // Polar angle in lab system range [degree]

  if (BoxCosTht) boxGen->SetCosTheta();

  boxGen->SetXYZ(0., 0., 0.); //cm

  LOG(INFO) << "Using PndBoxGenerator(" << GetBeamMom() <<", pdg="<<BoxType<<" mult="<<BoxMult
            <<" ) generator with range p["<<BoxMomMin<<","<<BoxMomMax<<"]  tht["<<BoxThtMin<<","<<BoxThtMax<<"]"<<(BoxCosTht?"*":"")<<"  phi["<<BoxPhiMin<<","<<BoxPhiMax<<"]" << FairLogger::endl;

  //  cout <<"BOX generator range: p["<<BoxMomMin<<","<<BoxMomMax<<"]  tht["<<BoxThtMin<<","<<BoxThtMax<<"]"<<(BoxCosTht?"*":"")<<"  phi["<<BoxPhiMin<<","<<BoxPhiMax<<"]"<<endl;

  fGen->AddGenerator(boxGen);
}

void PndMasterRunSim::UseDpmGenerator

(

)

Use DPM as event generator.

Definition at line 642 of file PndMasterRunSim.cxx.

References fDpmFlag.

Referenced by SetGenerator().

{
  LOG(INFO) << "Using PndDpmDirect(" << GetBeamMom() << ", " << fDpmFlag << ") generator" << FairLogger::endl;
  PndDpmDirect *Dpm= new PndDpmDirect(GetBeamMom(), fDpmFlag);
  fGen->AddGenerator(Dpm);
}

void PndMasterRunSim::UseEvtGenGenerator

(

TString

EvtGenFile

)

Use EvtGen as event generator.

EvtGen event generator

This call set EvtGen as event generator. The user should set the .dec file, and the function will retrieve automaticall beam momentum and initial state.

Parameters

EvtGenFile

Filename of the .dec file

Definition at line 665 of file PndMasterRunSim.cxx.

References fInputDir, s, PndEvtGenDirect::SetStoreTree(), and TString.

Referenced by SetGenerator().

{

  TString IniRes="";

  if (EvtGenFile.Contains(":")) // is the initial resonance provide as <decfile>.dec:iniRes ?
  {
    IniRes = EvtGenFile(EvtGenFile.Index(":")+1,1000);
    EvtGenFile = EvtGenFile(0,EvtGenFile.Index(":"));
  }

  if (IniRes=="") // we need to search the decay file
  {
    TString fnamepath=fInputDir+EvtGenFile;
    std::ifstream fs(fnamepath.Data());
    char line[250];

    while (fs)
    {
      fs.getline(line,249);
      TString s(line);
      s.ReplaceAll("\r","");
      if (IniRes=="" && s.Contains("Decay "))
      {
        if (s.Contains("#")) s=s(0,s.Index("#"));
        s.ReplaceAll("Decay ","");
        s.ReplaceAll(" ","");
        IniRes = s;
      }
    }
    fs.close();
  }
  /*
  // Looping over the dec file trying to find the first string "Decay", in order to find the initai
  // state as the following string
  FILE *dec = fopen(fInputDir+EvtGenFile,"r");
  if (dec==NULL) LOG(FATAL) << "The EvtGen dec file does not exist!! " << EvtGenFile << FairLogger::endl;

  char temp[6], particle[20];
  Bool_t found = kFALSE;
  while(fgets(temp, 6, dec) !=NULL)
    {
      if((strstr(temp, "Decay")) != NULL)
        {
          fscanf(dec, "%s",particle);
          LOG(INFO) << "It was found a " << particle << " as initial state." << FairLogger::endl;
          found = kTRUE;
          break;
        }
    }
  */
  if (IniRes=="") LOG(FATAL) << "The input file is not a proper .dec!! " << FairLogger::endl;

    //   TString  EvtInput =gSystem->Getenv("VMCWORKDIR");
  //   EvtInput+="/macro/run/psi2s_Jpsi2pi_Jpsi_mumu.dec";
  LOG(INFO) << "Using PndEvtGenDirect(" <<IniRes << ", " << (fInputDir+EvtGenFile).Data() << ", " << GetBeamMom() << ") generator" << FairLogger::endl;
  PndEvtGenDirect *EvtGen = new PndEvtGenDirect(IniRes, (fInputDir+EvtGenFile).Data(), GetBeamMom());
  EvtGen->SetStoreTree(kTRUE);
  fGen->AddGenerator(EvtGen);

}

void PndMasterRunSim::UseFtfGenerator

(

TString

ftfData

)

Use FTF as event generator.

FTF event generator

This call set FTF as event generator. If no input file is given FtfDirect is used. With input file the events are taken from ftfData.

Definition at line 650 of file PndMasterRunSim.cxx.

References fFtfFlag.

Referenced by SetGenerator().

{
  //if ( strncmp(fName,"TGeant4",7 ) == 0 ) LOG(FATAL) << "FTF does not run with Geant4 !!!"  << FairLogger::endl;
  if (ftfData.Contains(".root")){
          LOG(INFO) << "Using PndFtfGenerator with input file " << ftfData << FairLogger::endl;
          PndFtfGenerator* Ftf = new PndFtfGenerator(ftfData);
          fGen->AddGenerator(Ftf);
  } else {
          LOG(INFO) << "Using PndFtfDirect(anti_proton, G4_H, 1, ftfp, " << GetBeamMom() << ", " << gRandom->GetSeed() <<", "<<fFtfFlag<< ") generator" << FairLogger::endl;
          PndFtfDirect *Ftf = new PndFtfDirect("anti_proton", "G4_H", 1, "ftfp", GetBeamMom(), gRandom->GetSeed(), fFtfFlag);
          fGen->AddGenerator(Ftf);
  }
}

void PndMasterRunSim::UseLepLepGenerator

(

TString

leplepConfig

)

Use LepLepGenerator as event generator.

LepLep event generator

This call sets the PndLepLepGenerator as event generator. The format of the config string is: 'LEPLEP:pid(value):gegm(value):cosTheta(min, max)'

Parameters

leplepConfig

configuration string of the LepLep generator

Definition at line 574 of file PndMasterRunSim.cxx.

References Double_t, GetRange(), pid(), PndLepLepGenerator::SetBeamMom(), PndLepLepGenerator::SetCosThetaMax(), PndLepLepGenerator::SetCosThetaMin(), and TString.

{
  // use BOX generator; defaults

  Double_t cosThetaMin = 0.0;
  Double_t cosThetaMax = 1.0;
  Double_t pid = -1.0;
  Double_t GeGmRatio = -1.0;
  Double_t dummy = 1.0;

  leplepConfig.ToLower();

  if (leplepConfig!="")
  {
    leplepConfig.ReplaceAll("leplep","");
    leplepConfig.ReplaceAll(" ","");
    leplepConfig += ":";

    while (leplepConfig.Contains(":"))
    {
      TString curpar = leplepConfig(0,leplepConfig.Index(":"));
      leplepConfig = leplepConfig(leplepConfig.Index(":")+1,1000);
      curpar.ReplaceAll("[","(");
      curpar.ReplaceAll("]",")");

      if (curpar.BeginsWith("pid(")) GetRange(curpar, pid, dummy);
      if (curpar.BeginsWith("gegm(")) GetRange(curpar, GeGmRatio, dummy);
      if (curpar.BeginsWith("cosTheta(")) GetRange(curpar, cosThetaMin, cosThetaMax);
    }
  }

  PndLepLepGenerator* leplepGen = new PndLepLepGenerator();
  leplepGen->SetBeamMom(GetBeamMom());

  leplepGen->SetCosThetaMin(cosThetaMin);
  leplepGen->SetCosThetaMax(cosThetaMax);


  LOG(INFO) << "Using PndleplepGenerator(BeamMom = " << GetBeamMom() << " GeV/c, cosThetaRange = "<< cosThetaMin << " : " << cosThetaMax << ")" << FairLogger::endl;

  //  cout <<"BOX generator range: p["<<BoxMomMin<<","<<BoxMomMax<<"]  tht["<<BoxThtMin<<","<<BoxThtMax<<"]"<<(BoxCosTht?"*":"")<<"  phi["<<BoxPhiMin<<","<<BoxPhiMax<<"]"<<endl;

  fGen->AddGenerator(leplepGen);
}

void PndMasterRunSim::UsePiPiGenerator

(

TString

pipiConfig

)

Use PiPiGenerator as event generator.

PiPi event generator

This call sets the PndPiPiGenerator as event generator. The format of the config string is: 'PIPI:cosTheta(min, max)'

Parameters

pipiConfig

configuration string of the PiPi generator

Definition at line 535 of file PndMasterRunSim.cxx.

References Double_t, GetRange(), PndPiPiGenerator::SetBeamMom(), PndPiPiGenerator::SetCosThetaMax(), PndPiPiGenerator::SetCosThetaMin(), and TString.

Referenced by SetGenerator().

{
  // use BOX generator; defaults

  Double_t cosThetaMin = 0.0;
  Double_t cosThetaMax = 1.0;

  pipiConfig.ToLower();

  if (pipiConfig!="pipi")
  {
    pipiConfig.ReplaceAll("pipi","");
    pipiConfig.ReplaceAll(" ","");
    pipiConfig += ":";

    while (pipiConfig.Contains(":"))
    {
      TString curpar = pipiConfig(0,pipiConfig.Index(":"));
      pipiConfig = pipiConfig(pipiConfig.Index(":")+1,1000);
      curpar.ReplaceAll("[","(");
      curpar.ReplaceAll("]",")");

      if (curpar.BeginsWith("cosTheta(")) GetRange(curpar, cosThetaMin, cosThetaMax);
    }
  }

  PndPiPiGenerator* pipiGen = new PndPiPiGenerator();
  pipiGen->SetBeamMom(GetBeamMom());
  pipiGen->SetCosThetaMin(cosThetaMin);
  pipiGen->SetCosThetaMax(cosThetaMax);


  LOG(INFO) << "Using PndPiPiGenerator(BeamMom = " << GetBeamMom() << " GeV/c, cosThetaRange = "<< cosThetaMin << " : " << cosThetaMax << ")" << FairLogger::endl;

  //  cout <<"BOX generator range: p["<<BoxMomMin<<","<<BoxMomMax<<"]  tht["<<BoxThtMin<<","<<BoxThtMax<<"]"<<(BoxCosTht?"*":"")<<"  phi["<<BoxPhiMin<<","<<BoxPhiMax<<"]"<<endl;

  fGen->AddGenerator(pipiGen);
}

Member Data Documentation

Int_t PndMasterRunSim::fDpmFlag

private

Flag for DPM event generator.

Definition at line 287 of file PndMasterRunSim.h.

Referenced by SetDpmFlag(), and UseDpmGenerator().

Int_t PndMasterRunSim::fEventCounterRate

private

After how many events the counter will print.

Definition at line 290 of file PndMasterRunSim.h.

Referenced by AddSimTasks(), and SetEventCounterRate().

Int_t PndMasterRunSim::fFtfFlag

private

Flag for FTF event generator.

Definition at line 288 of file PndMasterRunSim.h.

Referenced by SetFtfFlag(), and UseFtfGenerator().

TString PndMasterRunSim::fInput

private

Name of the input for the simulation.

Definition at line 280 of file PndMasterRunSim.h.

Referenced by SetGenerator(), SetInput(), and Setup().

TString PndMasterRunSim::fInputDir

private

Name of the input directory for the simulation.

Definition at line 281 of file PndMasterRunSim.h.

Referenced by SetInputDir(), and UseEvtGenGenerator().

Int_t PndMasterRunSim::fNEvents

private

Number of events.

Definition at line 289 of file PndMasterRunSim.h.

Referenced by AddSimTasks(), and SetNumberOfEvents().

TString PndMasterRunSim::fOptions

private

Options parsed to the simulation.

Definition at line 285 of file PndMasterRunSim.h.

Referenced by AddSimTasks(), CreateGeometry(), CreateGeometryDay1(), Finish(), GetOptions(), SetGenerator(), SetOptions(), and Setup().

TString PndMasterRunSim::fOutFile

private

Name of the output file.

Definition at line 282 of file PndMasterRunSim.h.

Referenced by Finish(), and Setup().

TString PndMasterRunSim::fParamAsciiFile

private

Name of the parameter ascii file.

Definition at line 284 of file PndMasterRunSim.h.

Referenced by Finish(), SetParamAsciiFile(), and Setup().

TString PndMasterRunSim::fParamRootFile

private

Name of the parameter root file.

Definition at line 283 of file PndMasterRunSim.h.

Referenced by Finish(), SetParamRootFile(), and Setup().

FairRuntimeDb* PndMasterRunSim::fRtdb

private

Runtime DB.

Definition at line 293 of file PndMasterRunSim.h.

Referenced by Finish(), and Setup().

Short_t PndMasterRunSim::fTargetMode

private

Target mode.

Definition at line 291 of file PndMasterRunSim.h.

Referenced by SetGenerator(), SetTargetMode(), and Setup().

TStopwatch PndMasterRunSim::fTimer

private

Timer.

Definition at line 294 of file PndMasterRunSim.h.

Referenced by Finish(), and PndMasterRunSim().

---

The documentation for this class was generated from the following files:

• PndMasterRunSim.h
• PndMasterRunSim.cxx