Beispiel #1
0
 void particle() {
   String outpar;
   outpar =
       "x=" + Output.f(p.p.x) + " " + "y=" + Output.f(p.p.y) + " " + "z=" + Output.f(p.p.z) + "\n";
   outpar += "t=" + Output.f(p.p.t) + " r=" + Output.f(p.p.r) + "\n";
   outpar += "theta=" + Output.f(p.p.theta) + " phi=" + Output.f(p.p.phi) + "\n";
   outpar += "energy=" + Output.f(p.p.e) + " momentum=" + Output.f(p.p.p);
   par.append(outpar + "\n");
 }
Beispiel #2
0
  public void actionPerformed(ActionEvent event) {
    double e = 10, x = 10, result;
    String line;

    double vcut, ecut;
    String med, type;
    boolean lpm, cont, scat, time, intr;
    int bs, ph, sh, bb;
    boolean flag, iflag;

    if (event.getSource() == panel2.button) {
      type = panel3.type.getSelectedItem();
      if (type.equals("mpl")) type = "monopole";
      if (type.equals("all")) {
        Output.DEBUG = false;
        panel5.debug.setState(false);
        almc lmc =
            new almc(
                panel4.text1.getText()
                    + " -tdir="
                    + getCodeBase()
                    + " -user -raw "
                    + panel4.text2.getText());
        reset = true;
        return;
      } else stdReset();

      e = (new Double(panel4.text1.getText())).doubleValue();
      x = (new Double(panel4.text2.getText())).doubleValue();

      out.append("x=" + Output.f(x) + " e=" + Output.f(e) + " ");
      par.append("\nplease wait ...\n");

      lpm = panel5.lpm.getState();
      cont = panel5.cont.getState();
      scat = panel5.scat.getState();
      time = panel5.time.getState();
      debug = panel5.debug.getState();
      intr = panel5.inter.getState();
      switch (panel5.cuts.getSelectedIndex()) {
        case 0:
          vcut = 1;
          ecut = -1;
          break;
        case 1:
          vcut = 0.05;
          ecut = -1;
          break;
        case 2:
          vcut = 0.01;
          ecut = -1;
          break;
        case 3:
          vcut = 1.e-3;
          ecut = -1;
          break;
        default:
        case 4:
          vcut = -1;
          ecut = 500;
          break;
      }
      switch (panel5.brem.getSelectedIndex()) {
        case 0:
          bs = 1;
          break;
        case 1:
          bs = 2;
          break;
        case 2:
          bs = 3;
          break;
        default:
        case 3:
          bs = 4;
          break;
      }
      switch (panel5.phnu.getSelectedIndex()) {
        case 0:
          ph = 1;
          sh = 1;
          bb = 3;
          break;
        case 1:
          ph = 1;
          sh = 1;
          bb = 4;
          break;
        case 2:
          ph = 2;
          sh = 1;
          bb = 3;
          break;
        case 3:
          ph = 3;
          sh = 1;
          bb = 1;
          break;
        case 4:
          ph = 3;
          sh = 1;
          bb = 2;
          break;
        default:
        case 5:
          ph = 4;
          sh = 2;
          bb = 1;
          break;
      }
      med = panel3.medi.getSelectedItem();
      switch (panel3.romb.getSelectedIndex()) {
        case 0:
          romb = 3;
          break;
        case 1:
          romb = 4;
          break;
        default:
        case 2:
          romb = 5;
          break;
      }

      flag = ini;
      ini = true;
      iflag = this.intr;
      if (this.vcut != vcut) {
        this.vcut = vcut;
        flag = false;
      }
      if (this.ecut != ecut) {
        this.ecut = ecut;
        flag = false;
      }
      if (!this.med.equals(med)) {
        this.med = med;
        flag = false;
      }
      if (!this.type.equals(type)) {
        this.type = type;
        flag = false;
      }

      if (!flag) p = new Propagate(med, ecut, vcut, type);
      Output.DEBUG = debug;
      if (!flag) iflag = false;

      if (this.lpm != lpm) {
        this.lpm = lpm;
        flag = false;
      }
      if (this.cont != cont) {
        this.cont = cont;
        flag = false;
      }
      if (this.scat != scat) {
        this.scat = scat;
        flag = false;
      }
      if (this.time != time) {
        this.time = time;
        flag = false;
      }
      if (this.intr != intr) {
        this.intr = intr;
        flag = false;
      }
      if (this.bs != bs) {
        this.bs = bs;
        flag = false;
      }
      if (this.ph != ph) {
        this.ph = ph;
        flag = false;
      }
      if (this.sh != sh) {
        this.sh = sh;
        flag = false;
      }
      if (this.bb != bb) {
        this.bb = bb;
        flag = false;
      }

      p.sdec = true;
      p.recc = true;
      p.contiCorr = cont;
      p.exactTime = time;
      p.s.lpm = lpm;
      p.s.b.form = bs;
      p.s.n.form = ph;
      p.s.n.shadow = sh;
      p.s.n.bb = bb;
      p.molieScat = scat;

      Propagate.g = 5;
      if (!flag) {
        if (intr) {
          err.append("\n");
          p.interpolate("all");
        } else if (iflag) {
          err.append("\nSwitching off parameterizations:\n");
          p.interpolate("none");
        }
      }
      Propagate.g = romb;

      result = p.propagateTo(x * 1.e2, e * 1.e3);
      line = "ef=" + Output.f(result > 0 ? result * 1.e-3 : result * 1.e-2);
      if (time) line += " tf=" + Output.f(p.p.t);
      out.append(line + "\n");
      particle();
    }
  }