static const float ppm=2107.84, ovr=OVR*OVR*OVER;
static const float crf[4]={ ppm*4.889*0.894, // em shower
			    ppm*4.076*0.860, // hdr m/GeV
			    ppm*1.1720,      // bare
			    ppm*0.0324 };    // muon

photon p;
void output();
deque<string> flnz;
unsigned int flnb=0, flne=0;

int hcmp(const void *a, const void *b){
  return ((hit *) a)->n - ((hit *) b)->n;
}

void print(){
  if(flnb<flne) qsort(d.hits, e.hidx, sizeof(hit), hcmp);

  for(int i=0; i<e.hidx; i++){
    hit & h = d.hits[i];

    for(; flnb<h.n; flnb++){
      cout<<flnz.front()<<endl;
      flnz.pop_front();
    }

    name n=d.names[h.i];
    printf("HIT %d %d %f %f\n", (int) n.str, (int) n.dom, h.t, h.z);
  }

  for(; flnb<flne; flnb++){
    cout<<flnz.front()<<endl;
    flnz.pop_front();
  }
  p.q=0; flnb=0; flne=0;
}

unsigned long long f2ki, f2kn;

void cascade(float rx, float ry, float rz, float t, float E, int type){
  p.t=t; p.r[0]=rx; p.r[1]=ry; p.r[2]=rz;
  unsigned long long num=llroundf(crf[type]*E);
  f2kn+=num;
  for(; f2ki<f2kn; f2ki+=ovr){
    pz[pn++]=p; if(pn==pmax/OVER) output();
  }
}

void muon(float rx, float ry, float rz, float t, float E, float dr){
  unsigned long long num=llroundf(max(0.0f, dr*(crf[2]+crf[3]*logf(E))));
  f2kn+=num;
  for(unsigned long long i=0; f2ki<f2kn; f2ki+=ovr, i+=ovr){
    p.t=t; p.r[0]=rx; p.r[1]=ry; p.r[2]=rz;
    float a=(dr*i)/num; p.t+=a*d.ocv;
    for(int i=0; i<3; i++) p.r[i]+=p.n[i]*a;
    pz[pn++]=p; if(pn==pmax/OVER) output();
  }
}

void f2k(){
  ini(0);

  string in;
  while(getline(cin, in)){
    flnz.push_back(in); flne++;
    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(), "%31s", name)==1){ if(strcmp(name, "EM")==0){ f2ki=0; f2kn=0; } }
    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*=fpi/180; ph*=fpi/180;
      float costh=cosf(th), sinth=sinf(th), cosph=cosf(ph), sinph=sinf(ph);
      p.q=flne; p.n[0]=sinth*cosph; p.n[1]=sinth*sinph; p.n[2]=costh;
      if(0==strcmp(name, "amu+") || 0==strcmp(name, "amu-") || 0==strcmp(name, "amu")) muon(x, y, z, 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, t, E, 0);
      else if(0==strcmp(name, "munu") || 0==strcmp(name, "hadr")) cascade(x, y, z, t, E, 1);
    }
  }
  output();

  fin();
}