flyopunta/src/src_dev/devsfera.cpp

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

//using namespace std;

DevSfera::~DevSfera()
{
}
//==================  Sfera=======Erede di RettEdro
DevSfera::DevSfera()
{
  devtype=TypDevSfera;
  devclass="Sfera";
  rout=Lout.z;
  routq=rout*rout;
  Lout.Norma();
  Rqsize=Lout.normaq;
  rin=Lin.z;
  rinq=rin*rin;
  rbuco=Lin.x=Lin.y;  // attenzione a definire le dimensione del buco di ingresso del beam...
  rbucoq=rbuco*rbuco;
  zbucoout= -sqrt ( routq-rbucoq ); // superficie interna del buco.
  zbucoin = -sqrt ( rinq-rbucoq ); // superficie interna del buco.

}

//----------------- --------
void DevSfera::Prgeom()
{
  Gout<<"==================  R i v e l a t o r e =================";
  Gout<<"\n --> < "<<nome <<" >  [ "<<fun
  <<" ]  Sfera di rout "<< rout<<" cm   con un buco cilin. di "<<rbuco<<" cm";
  Device::Prgeom();

}

//======================================================
//----------------- -----P o s i z i o n e ---
int DevSfera::Posizione()
{

  if ( X_dev.normaq>routq ) return 1;//e' fuori

  if ( X_dev.normaq>rinq ) // è interno tra i due gusci sferici
    { // è forse anche nel buco all'ingresso?
      if ( X_dev.z>0.0 ) return 0; // è interno nel guscio della semisfera dowstream
      if ( ( X_dev.x*X_dev.x+ X_dev.y*X_dev.y ) >rbucoq ) return 0;// è interno comunque anche  nella semisfera upstream
      if ( X_dev.z<zbucoout ) return 1; // è  fuori dall'ingresso upstream
    }
  return -1;// e se arivva qui è nel buco !
}
//----------------------------------------------------------
//----------------------C a m E s t e r --------------------
double DevSfera::CamEster()
{
// calcolo il cammino positivo verso la  superfice della sfera
// e scelgo il cammino maggiore
// calcolo i coefficienti della equazione intersezione di una retta con la sfera.
  camm=-1.;

// se mi allontano  ritorno -1.

//con la sfera

  double a=X_dev%V_dev;
  if ( a>0.0 ) return -1.;
  double b=X_dev.normaq-routq;
  double d=a*a-b;
  if ( d<0.0 ) return -1.;
  camm= -a - sqrt ( d );
  //
// distanza dal piano del buco di ingresso del beam
  if ( abs ( V_dev.z ) >0.00001 )
    {
      double t= ( zbucoout-X_dev.z ) /V_dev.z;
      if ( t<0.0 ) return -1;
      if ( t>camm ) camm=t;
    }

  return camm;
}
//===========================================================
//----------------------C a m I n t e r --------------------
double DevSfera::CamInter()
{
  camm=0.01;
// sono tra le due sfere
  double a=X_dev%V_dev;
  double b=X_dev.normaq-rinq;  // sempre positivo
  double d=a*a-b;
  if ( a<=0.0&&d>=0.0)
    { // va verso la sfera interna
      camm= -a - sqrt ( d );
    }
  else if ( a>0.0 )
    {// dovrebbe colpire la sfera esterna
      b=X_dev.normaq-routq;  // sempre negativo
      camm= -a + sqrt ( a*a-b ); // è positivo
    }
   return camm; 
}
//===========================================================
//----------------------C a m I n t e r --------------------
double DevSfera::CamBuco()
{
  camm=1.0e+10;
// calcolo il cammino verso la  superfice sferica esterna
// scelgo il cammino positivo minore...
// calcolo il cammino positivo verso la superfice del Sfera interno.

//con la sfera

  double a=X_dev%V_dev;
  double b=X_dev.normaq-rinq;  // sempre negativo
  camm= -a + sqrt ( a*a-b ); // è positivo
// distanza dal piano di ingresso
  if ( abs ( V_dev.z ) >0.00001 )
    {
      double t= ( zbucoout-X_dev.z ) /V_dev.z;
      if ( t<camm ) camm=t;
    }
  return camm;
}