flyopunta/src/src_dev/devclrt.cpp

/***************************************************************************
                          devclrt.cpp  -  description
                             -------------------
    begin                : Sat May 26 2001
    copyright            : (C) 2001 by Giuseppe Pierazzini
    email                : peppe@unipi.it
  ***************************************************************************
  *                                                                         *
  *   NA48  simulation program.                                             *
  *                                                                         *
  ***************************************************************************/
#include "parm.h"
#include "devclrt.h"
 
// using namespace std;

void DevClRt::Prgeom()
{
 Gout<<"==================  R i v e l a t o r e =================";
 Gout<<"\n --> < "<<nome <<" >  [ "<<fun
     <<" ]  Cilindrico con buco di x "<<Lin.x<<" e y " <<Lin.y<<" rettangolare ";
  Device::Prgeom();
}
//======================================================
//----------------- -----P o s i z i o n e ---
int DevClRt::Posizione()
{
 if( DevCln::Posizione()>0) return 1; //  il punto e' esterno
 // se sono interno al rivelatore verifico se siamo nel buco
 P_Cin=X_dev-Centrin;
 Dind=P_Cin-Lin;
 Dins=-(Lin+P_Cin);
 if(Dins[2]>0.0 || Dind[2]>0.0 || Dins[0]>0.0 ||
    Dind[0]>0.0 || Dins[1]>0.0 || Dind[1]>0.0   ) return 0;
// ritorno -1 se il punto e' nel buco altrimento ritorno 0
   return -1;
}
//----------------------C a m E s t e r --------------------
// usa quella di Cilindro per eredita'
//===========================================================
//----------------------C a m I n t e r --------------------

double DevClRt::CamInter()
{
//uso quella di Cilindro  per calcolare il cammino
// interno senza il buco...
  double tratto=DevCln::CamInter();
//poi verifico se colpisce le pareti del buco .....
  double nv,t;
  camm=-1.;
// calcolo il cammino positivo verso le sei superfici del buco. Per dist>0
// scelgo il cammino maggiore

  if((nv=V_dev[0])>0.0 && Dins[0]>0.0)  // piani x
      {t=Dins[0]/nv;  if(t>camm)camm=t;}
  else if(nv<0.0 && Dind[0]>0.0)
      {t=-Dind[0]/nv; if(t>camm)camm=t;}

  if((nv=V_dev[1])>0.0 && Dins[1]>0.0)  // piani y
      {t=Dins[1]/nv;  if(t>camm)camm=t;}
  else if(nv<0.0 && Dind[1]>0.0)
      {t=-Dind[1]/nv; if(t>camm)camm=t;}

  if((nv=V_dev[2])>0.0 && Dins[2]>0.0)  // piani z
      {t=Dins[2]/nv;  if(t>camm)camm=t;}
  else if(nv<0.0 && Dind[2]>0.0)
      {t=-Dind[2]/nv; if(t>camm)camm=t;}

  // quindi confronto camm con tratto e accetto il cammino piu' breve......
  if(camm<0.0) {camm=tratto; return camm;}
  if(camm>0.0 && tratto<camm) camm=tratto;
  return camm;
}
//===============================================================
//--------------- C a m B u c o ---------------------------------
double DevClRt::CamBuco()
{
  double nv,t;
  camm=1.0e+10;
// calcolo il cammino verso le sei superfici
// le distanze sono per definizione negative...
// e scelgo il cammino positivo  minore

  if((nv=V_dev[0])>0.0)  // piani x
      {t=-Dind[0]/nv; if(t<camm)camm=t; }
  else if(nv<0.0)
      {t=Dins[0]/nv;  if(t<camm)camm=t; }

  if((nv=V_dev[1])>0.0)  // piani y
      {t=-Dind[1]/nv; if(t<camm)camm=t; }
  else if(nv<0.0)
      {t=Dins[1]/nv;  if(t<camm)camm=t; }

  if((nv=V_dev[2])>0.0)  // piani z
      {t=-Dind[2]/nv; if(t<camm)camm=t; }
  else if(nv<0.0)
      {t=Dins[2]/nv;  if(t<camm)camm=t; }
return camm;
}