LmvmSimParam.h

Go to the documentation of this file.

/* Copyright (C) 2021 Justus-Liebig-Universitaet Giessen, Giessen
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Semen Lebedev [committer] */
 
#ifndef LMVM_SIM_PARAM_H
#define LMVM_SIM_PARAM_H
 
#include "Logger.h"
#include "TObject.h"
 
#include <string>
 
class LmvmSimParam {
 public:
  static Double_t GetWeight(const std::string& energy, const std::string& particle)
  {
    // Weight = Multiplicity * Branching Ratio
    if (energy == "8gev") {
      if (particle == "omegaepem") return 2.28721 * 7.36e-5;
      if (particle == "omegadalitz") return 2.28721 * 7.7e-4;
      if (particle == "phi") return 0.311619 * 2.97e-4;
      if (particle == "inmed" || particle == "rho0") return 0.0304706;
      if (particle == "qgp" || particle == "qgp_epem") return 4.52941e-4;
    }
    else if (energy == "12gev") {
      if (particle == "omegaepem") return 5.82093 * 7.36e-5;
      if (particle == "omegadalitz") return 5.82093 * 7.7e-4;
      if (particle == "phi") return 1.01808 * 2.97e-4;
      if (particle == "inmed" || particle == "rho0") return 0.041625;
      if (particle == "qgp" || particle == "qgp_epem") return 6.1875e-4;
    }
    else if (energy == "25gev") {
      if (particle == "rho0") return 23 * 4.7e-5;
      if (particle == "omegaepem") return 38 * 7.28e-5;
      if (particle == "omegadalitz") return 38 * 7.7e-4;
      if (particle == "phi") return 1.28 * 2.97e-4;
      if (particle == "inmed") return 4.45e-2;
      if (particle == "qgp") return 1.15e-2;
    }
    else if (energy == "3.5gev") {
      if (particle == "rho0") return 1.0 * 4.7e-5;
      if (particle == "omegaepem") return 1.2 * 7.28e-5;
      if (particle == "omegadalitz") return 1.2 * 7.7e-5;
      if (particle == "phi") return 0.1 * 2.97e-4;
    }
    else if (energy == "4.5gev") {
      if (particle == "omegadalitz") return 1.2 * 7.7e-5;
      if (particle == "omegaepem") return 1.2 * 7.28e-6;
      if (particle == "phi") return 1.2 * 2.97e-6;
      if (particle == "inmed") return 2.4 * 10e-3;
    }
    else {
      LOG(fatal) << "LmvmSimParam::SetEnergyAndParticle(): energy or particle is not correct, energy:" << energy
                 << " particle:" << particle;
      return 0.;
    }
    return 0.;
  }
};
 
 
#endif