Added custom DIME generator

alisw · Jul 21, 2017 · a4314cf · a4314cf
1 parent 7f4ec9a
commit a4314cf
Showing 1 changed file with 28 additions and 0 deletions.
diff --git a/MC/CustomGenerators/PWGUD/Dime.C b/MC/CustomGenerators/PWGUD/Dime.C
@@ -0,0 +1,28 @@
+AliGenerator* GeneratorCustom() {
+  printf("Starting DIME generator\n");
+
+  TString proc;
+  if      (processConfig.Contains("Pipm")) proc = "pipm";
+  else if (processConfig.Contains("Kpkm")) proc = "kpkm";
+  else printf("ERROR: uknown process\n");
+
+  TString ff;
+  if      (processConfig.Contains("Orexp")) ff = "orexp";
+  else if (processConfig.Contains("Exp"))   ff = "exp";
+  else if (processConfig.Contains("Power")) ff = "power";
+  else printf("ERROR: uknown process\n");
+
+  gSystem->Load("libEVGEN");
+  gSystem->Load("libTDime");
+  gSystem->Load("libdime");
+  AliGenDime* dime = new AliGenDime(1000);
+  AliDimeRndm::GetDimeRandom()->SetSeed(seedConfig);
+  printf("Setting new random seed: %i\n",AliDimeRndm::GetDimeRandom()->GetSeed());
+  dime->GetTDime()->SetEnergyCMS(energyConfig);       // sqrt(s)
+  dime->GetTDime()->SetProcess(proc.Data());          // Process {"pipm", "pi0", "kpkm", "ks", "rho"}
+  dime->GetTDime()->SetFormf(ff.Data());              // Meson-Pomeron form factor {"orexp", "exp", "power"}
+  dime->GetTDime()->SetFsi("true");                   // Exclusive supression {"true", "false"}
+  dime->GetTDime()->SetYRange(yminConfig,ymaxConfig); // Set rapidity range of mesons
+  dime->GetTDime()->SetMinPt(0.);                     // Minimum pT of mesons
+  return dime;
+}