#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"

  float xrnd(){
    float rnd;
    do rnd=rand(); while(rnd==0);
    const float a=1./(1ll+RAND_MAX);
    return a*rnd;
  }

  union name{
    unsigned short k;
    struct{
      unsigned char dom;
      unsigned char str;
      unsigned short one;
    } n;
    float f;

    name(){
      f=1.0;
    }
  };

  struct DOM{
    name t;
    float r[3];

    DOM(){
      for(int i=0; i<3; i++) r[i]=0;
    }
  };

  vector<DOM> oms;

#ifdef ACCL2
  static const int cx=10, cy=10, cz=10;
  int cn[3]={cx, cy, cz};
  vector<int> cells[cx][cy][cz];
  float cl[3], ch[3], crst[3];
#endif

  void load(){
    ifstream inFile;
    inFile.open("geo-f2k", ifstream::in);

    if(!inFile.fail()){
      DOM dom;
      string mbid;
      int om, str;
      unsigned long long omid;
      while(inFile>>mbid>>hex>>omid>>dec>>dom.r[0]>>dom.r[1]>>dom.r[2]>>str>>om) if(str>0 && om>=1 && om<=60){
	dom.r[2]+=zoff;
	dom.t.n.dom=om;
	dom.t.n.str=str;
	oms.push_back(dom);
      }
      inFile.close();
    }

#ifdef ACCL2
    for(int i=0; i<cx; i++) for(int j=0; j<cy; j++) for(int k=0; k<cz; k++)
      cells[i][j][k].erase(cells[i][j][k].begin(), cells[i][j][k].end());

    for(int n=0; n<oms.size(); n++) for(int m=0; m<3; m++){
      if(n==0 || oms[n].r[m]<cl[m]) cl[m]=oms[n].r[m];
      if(n==0 || oms[n].r[m]>ch[m]) ch[m]=oms[n].r[m];
    }
    for(int m=0; m<3; m++){
      double diff=ch[m]-cl[m];
      if(diff<R){ cerr<<"Error in cell initialization"<<endl; exit(1); }
      crst[m]=(cn[m]-1)/diff;
    }

    for(int n=0; n<oms.size(); n++){
      int i[3];
      for(int m=0; m<3; m++){
	i[m]=lroundf((oms[n].r[m]-cl[m])*crst[m]);
	if(i[m]<0 && i[m]>=cn[m]){ cerr<<"Error in cell initialization"<<endl; exit(1); }
      }
      cells[i[0]][i[1]][i[2]].push_back(n);
    }
#endif

    cerr<<oms.size()<<" DOMs loaded"<<endl;
  }

  struct SRC{
    float t, r[3], n[3];
    int type;

    SRC(){
      load();
      for(int i=0; i<3; i++){ r[i]=0; n[i]=i==0?1:0; }
    }

    void cascade(){
      type=1;
      m.next(xrnd());
    }

    void advance(float a){
      t+=a*m.ocv;
      for(int i=0; i<3; i++) r[i]+=n[i]*a;
    }

    void init_flasher(int str, int om){
      type=0;
      t=0;
      for(int n=0; n<oms.size(); n++) if(oms[n].t.n.str==str && oms[n].t.n.dom==om){
	for(int m=0; m<3; m++) r[m]=oms[n].r[m]; break;
      }
    }

    void flasher(){
      static const float rms=9.2*sqrtf(3)*M_PI/180;
      float xi=xrnd();
      xi=rms*(2*xi-1);
      float nz=sinf(xi), np=cosf(xi);

      xi=xrnd();
      xi*=2*M_PI;
      n[0]=np*cosf(xi);
      n[1]=np*sinf(xi);
      n[2]=nz;
    }
  } x;

  struct photon{
    float t, r[3], n[3];

    photon(){
      t=x.t;
      for(int i=0; i<3; i++){
	r[i]=x.r[i]; n[i]=x.n[i];
      }

      if(x.type==1) rotate(m.w->coschr, m.w->sinchr);
    }

    void rotate(float cs, float si){
      float u[3];
      {
	float p1[3], p2[3];
	{
	  int i0=0;
	  float n0=n[0];
	  for(int i=1; i<3; i++) if(abs(n[i])>abs(n0)){ i0=i; n0=n[i]; }
	  int i1=(i0+1)%3, i2=(i0+2)%3;
	  float n1=n[i1], n2=n[i2], nq=n0*n0;
	  float r1=1/sqrtf(nq+n1*n1), r2=1/sqrtf(nq+n2*n2);

	  p1[i0]=-n1*r1; p1[i1]=n0*r1; p1[i2]=0;
	  p2[i0]=-n2*r2; p2[i1]=0; p2[i2]=n0*r2;
	}
	{
	  float q1[3], q2[3], r1=0, r2=0;
	  for(int i=0; i<3; i++){
	    float a1=p1[i]-p2[i]; q1[i]=a1; r1+=a1*a1;
	    float a2=p1[i]+p2[i]; q2[i]=a2; r2+=a2*a2;
	  }

	  r1=1/sqrtf(r1); r2=1/sqrtf(r2);

	  for(int i=0; i<3; i++){
	    p1[i]=q1[i]*r1;
	    p2[i]=q2[i]*r2;
	  }
	}

	{
	  float zi=xrnd(); zi*=2*M_PI;
	  float px=cosf(zi), py=sinf(zi);
	  for(int i=0; i<3; i++) u[i]=px*p1[i]+py*p2[i];
	}
      }

      {
	float r=0;
	for(int i=0; i<3; i++){
	  float a=n[i]*cs+u[i]*si;
	  n[i]=a; r+=a*a;
	}
	{
	  r=1/sqrtf(r);
	  for(int i=0; i<3; i++) n[i]*=r;
	}
      }
    }

    void rotate(){
      float xi=xrnd();
      xi=2*xi-1;
      float g=m.g;
      if(g!=0){
	float g2=g*g;
	float ga=(1-g2)/(1+g*xi);
	xi=(1+g2-ga*ga)/(2*g);
      }
      if(xi>1) xi=1; else if(xi<-1) xi=-1;

      float si=sqrtf(1-xi*xi);
      rotate(xi, si);
    }

    void advance(float a){
      t+=a*m.w->ocm;
      for(int i=0; i<3; i++) r[i]+=n[i]*a;
    }

    void propagate(){
      float TOT=-logf(xrnd());
      while(true){
	float SCA=-logf(xrnd());
	float sca=m.sX(SCA, r[2], n[2]);
	float tot=m.aX(TOT, r[2], n[2]);
	if(tot<sca) sca=tot;

	name om=sphere(sca);
	TOT-=m.xA(sca, r[2], n[2]);
	advance(sca);

	if(om.k!=0){
	  cout<<"HIT "<<(int)om.n.str<<" "<<(int)om.n.dom<<" "<<t<<" "<<n[2]<<endl;
	  break;
	}

	if(TOT<xx) break;
	rotate();
      }
    }

    name sphere(float &sca){
      name imin;

#ifdef ACCL0
      float dx[3], dn[3];
      for(int m=0; m<3; m++){
	float ini=r[m];
	float fin=ini+sca*n[m];
	if(n[m]>0){
	  dn[m]=ini;
	  dx[m]=fin;
	}
	else{
	  dn[m]=fin;
	  dx[m]=ini;
	}
	dn[m]-=R;
	dx[m]+=R;
      }
#endif

#ifdef ACCL2
      int xl[3], xh[3];
      for(int m=0; m<3; m++){
	xl[m]=min(max((int) lroundf((dn[m]-cl[m])*crst[m]), 0), cn[m]);
	xh[m]=max(min((int) lroundf((dx[m]-cl[m])*crst[m]), cn[m]-1), -1);
      }
      for(int i=xl[0]; i<=xh[0]; i++) for(int j=xl[1]; j<=xh[1]; j++) for(int k=xl[2]; k<=xh[2]; k++) for(int l=0; l<cells[i][j][k].size(); l++){
	const DOM & om=oms[cells[i][j][k][l]];
#else
      for(int l=0; l<oms.size(); l++){
	const DOM & om=oms[l];
#endif
#ifdef ACCL1
	bool pass=false;
	for(int i=0; i<3; i++) if(om.r[i]<dn[i] || om.r[i]>dx[i]){ pass=true; break; }
	if(pass) continue;
#endif

	{
	  float b=0, c=0;
	  for(int i=0; i<3; i++){
	    float dr=om.r[i]-r[i];
	    b+=n[i]*dr; c+=dr*dr;
	  }
	  c-=R*R;
	  float D=b*b-c, h=-1;
	  if(D>0) h=b-sqrtf(D);

	  if(h>0 && h<sca){ imin=om.t; sca=h; }
	}
      }

      return imin;
    }
  };

  void flasher(int str, int dom, unsigned long long num){
    x.init_flasher(str, dom);
    for(unsigned long long i=0; i<num; i++){
      x.flasher();
      photon p;
      p.propagate();
    }
  }

  void cascade(float rx, float ry, float rz, float nx, float ny, float nz, float t, float E, int type){
    x.t=t;
    x.r[0]=rx; x.r[1]=ry; x.r[2]=rz;
    x.n[0]=nx; x.n[1]=ny; x.n[2]=nz;
    unsigned long long num=(unsigned long long) roundl(crf[type]*E);
    for(unsigned long long i=0; i<num; i++){
      x.cascade();
      photon p;
      p.propagate();
    }
  }

  void muon(float rx, float ry, float rz, float nx, float ny, float nz, float t, float E, float dr){
    x.n[0]=nx; x.n[1]=ny; x.n[2]=nz;
    unsigned long long num=(unsigned long long) roundl(max((float) 0, dr*(crf[2]+crf[3]*logf(E))));
    for(unsigned long long i=0; i<num; i++){
      x.t=t;
      x.r[0]=rx; x.r[1]=ry; x.r[2]=rz;
      x.advance(dr*i/num);
      x.cascade();
      photon p;
      p.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]);
    srand(1+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]);

    srand(1+seed);
    flasher(str, dom, num);
  }
}