Ejemplo n.º 1
0
def get_channel(operator):
    PSQ = operator.getMomentumSquared()
    irrep = operator.getLGIrrep()
    isospin = operator.getIsospin()
    strangeness = operator.getStrangeness()
    channel = defs.Channel(PSQ, irrep, isospin, strangeness)
    return channel
Ejemplo n.º 2
0
def get_channel(correlator):
    op = correlator.getSink().getBasicLapH()
    P = (op.getXMomentum(), op.getYMomentum(), op.getZMomentum())
    irrep = op.getLGIrrep()
    irrep_row = op.getLGIrrepRow()
    isospin = op.getIsospin()
    strangeness = op.getStrangeness()
    channel = defs.Channel(P, irrep, irrep_row, isospin, strangeness)
    return channel
Ejemplo n.º 3
0
def get_binfile(correlator, ensemble_name, base_output_dir):
    operator = correlator.getSink().getBasicLapH()
    P = (operator.getXMomentum(), operator.getYMomentum(),
         operator.getZMomentum())
    irrep = operator.getLGIrrep()
    irrep_row = operator.getLGIrrepRow()
    isospin = operator.getIsospin()
    strangeness = operator.getStrangeness()
    channel = defs.Channel(P, irrep, irrep_row, isospin, strangeness)

    output_dir = os.path.join(base_output_dir, ensemble_name)
    os.makedirs(output_dir, exist_ok=True)

    bin_file = os.path.join(output_dir, f"{channel.iso_strange_str()}.bin")
    return bin_file
Ejemplo n.º 4
0
def get_channel(correlator):
  op = correlator.getSink().getBasicLapH()
  P = (op.getXMomentum(), op.getYMomentum(), op.getZMomentum())
  irrep = op.getLGIrrep()
  irrep_row = op.getLGIrrepRow()
  isospin = op.getIsospin()
  if isospin=="pion":
    isospin="triplet"

  if isospin!="triplet":
    print("shit")
    exit()

  strangeness = op.getStrangeness()
  channel = defs.Channel(P, irrep, irrep_row, isospin, strangeness)
  return channel