#include <cmath>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <fstream>
#include <cstring>

#define ACCL1
#define ACCL2

#ifdef ACCL1
#define ACCL0
#endif
#ifdef ACCL2
#define ACCL0
#endif

namespace ppc{
  using namespace std;
  struct setprc{
    short store;
    setprc(){
      short byte=127, end;
      __asm__("fstcw %0\nfldcw %1\nfstcw %2"::"m" (store), "m" (byte), "m" (end));
      // cerr<<store<<" "<<byte<<" "<<end<<endl;
    }
    ~setprc(){
      short byte, end;
      __asm__("fstcw %0\nfldcw %1\nfstcw %2"::"m" (byte), "m" (store), "m" (end));
      // cerr<<byte<<" "<<store<<" "<<end<<endl;
    }
  } f;

#include "ice.cxx"

  void flasher(int str, int dom, unsigned long long num){
    d.type=0;
    d.xt=0;
    for(int n=0; n<oms.size(); n++) if(oms[n].t.str==str && oms[n].t.dom==dom){
      for(int m=0; m<3; m++) d.xr[m]=oms[n].r[m]; break;
    }

    for(unsigned long long i=0; i<num; i++) d.propagate();
  }

  void cascade(float rx, float ry, float rz, float nx, float ny, float nz, float t, float E, int type){
    d.type=1;
    d.xt=t;
    d.xr[0]=rx; d.xr[1]=ry; d.xr[2]=rz;
    d.xn[0]=nx; d.xn[1]=ny; d.xn[2]=nz;
    unsigned long long num=llroundf(crf[type]*E);
    for(unsigned long long i=0; i<num; i++){
      m.next(d.xrnd());
      d.propagate();
    }
  }

  void muon(float rx, float ry, float rz, float nx, float ny, float nz, float t, float E, float dr){
    d.type=1;
    d.xn[0]=nx; d.xn[1]=ny; d.xn[2]=nz;
    unsigned long long num=llroundf(max((float) 0, dr*(crf[2]+crf[3]*logf(E))));
    for(unsigned long long i=0; i<num; i++){
      d.xt=t;
      d.xr[0]=rx; d.xr[1]=ry; d.xr[2]=rz;

      float a=dr*i/num;
      d.xt+=a*d.ocv;
      for(int i=0; i<3; i++) d.xr[i]+=d.xn[i]*a;

      m.next(d.xrnd());
      d.propagate();
    }
  }
}

using namespace ppc;

void f2k(){
  string in;
  while(getline(cin, in)){
    cout<<in<<endl;
    char name[32];
    int gens, igen;
    float x, y, z, th, ph, l, E, t;
    char * str = "TR %d %d %31s %f %f %f %f %f %f %f %f";
    if(sscanf(in.c_str(), str, &gens, &igen, &name, &x, &y, &z, &th, &ph, &l, &E, &t)==11){
      th=180-th; ph=ph<180?ph+180:ph-180;
      th*=M_PI/180; ph*=M_PI/180;
      float costh=cosf(th), sinth=sinf(th), cosph=cosf(ph), sinph=sinf(ph);
      float nx=sinth*cosph, ny=sinth*sinph, nz=costh;
      if(0==strcmp(name, "amu+") || 0==strcmp(name, "amu-") || 0==strcmp(name, "amu")) muon(x, y, z, nx, ny, nz, t, E, l);
      else if(0==strcmp(name, "delta") || 0==strcmp(name, "brems") || 0==strcmp(name, "epair")
	      || 0==strcmp(name, "e+") || 0==strcmp(name, "e-") || 0==strcmp(name, "e")) cascade(x, y, z, nx, ny, nz, t, E, 0);
      else if(0==strcmp(name, "munu") || 0==strcmp(name, "hadr")) cascade(x, y, z, nx, ny, nz, t, E, 1);
    }
  }
}

main(int arg_c, char *arg_a[]){
  if(arg_c<=1) cerr<<"Use: "<<arg_a[0]<<" [str] [om] [num] [seed]"<<endl;
  else if(arg_c<=2){
    unsigned int seed=atoi(arg_a[1]);
    d.setr(seed);
    f2k();
  }
  else{
    int str=0, dom=0;
    unsigned int seed=0;
    unsigned long long num=1000000ULL;

    if(arg_c>1) str=atoi(arg_a[1]);
    if(arg_c>2) dom=atoi(arg_a[2]);
    if(arg_c>3) num=(unsigned long long) atof(arg_a[3]);
    if(arg_c>4) seed=atoi(arg_a[4]);

    d.setr(seed);
    flasher(str, dom, num);
  }
}