flyopunta/src/src_dev/devot.cpp

/***************************************************************************
                          devot.cpp  -  description
                             -------------------
    begin                : Sat Dec 6 2003
    copyright            : (C) 2003 by gmp
    email                : peppe@newpeppe
  ***************************************************************************
  *                                                                         *
  *   NA48  simulation program.                                             *
  *                                                                         *
  ***************************************************************************/
#include "parm.h"
#include "devot.h"
#include "particella.h"

DevOt::DevOt()
{
  devtype=TypDevOt;
  devclass="Ottagon";
  Firot=.707107;
};
DevOt::~DevOt()
{
}


void DevOt::Prgeom()
{

  Gout<<"==================  R i v e l a t o r e =================";
  Gout<<"\n --> < "<<nome <<" >  [ "<<fun
  <<" ]  ottagonale  senza buco "<<std::endl;
  Device::Prgeom();
}
//----------------- -----P o s i z i o n e ---
int DevOt::Posizione()
{
//     Posizione della particella rispetto questo dev
// Le distanze di un punto interno dai piani del
// rivelatore sono tutte negative o nulle data
// la definizione delle normali ai piani ( verso l'esterno)
// Le distanze di un punto esterno sono sia positive che negative:
// ma almeno  una distanza e' positiva.
// X_dev   e' la posizione  della particella rispeto al sistema del dev
// Ruoto di -45 gradi il punto e lo paragono con i lati inclinati....
// ripeto la stessa procedura di prima per simmetria (Lout non cambia);
// double Firot=0.707107= 1/sqrt(2.);
// X_rot   e' la posizione rotata di 45 gradi

  X_rot.setvn ( Firot* ( X_dev[0]+X_dev[1] ),Firot* ( X_dev[1]-X_dev[0] ), X_dev[2] );
  Doutrd=X_rot-Lout;
  Doutrs=- ( Lout+X_rot );

//Test now
  if ( DevRt::Posizione() >0 ) return 1;
  if ( Doutrs[0]>0.0 ||
          Doutrd[0]>0.0 || Doutrs[1]>0.0 || Doutrd[1]>0.0 ) return 1;

  return 0;  // 0 se interno, 1 se esterno. Il buco qui non c'e'.

}
//----------------------------------------------------------
//----------------------C a m E s t e r --------------------
// calcola il cammino esterno per raggiungere il rivelatore...
// vedi la routine nella classe base rettangolo.
// V_dev e' il versore della traiettoria nel dev
// V_rot e' il versore ruotato di 45 gradi per i lati inclinati
//
double DevOt::CamEster()
{
  double nv,t;
//   test sulla classe base
  camm=DevRt::CamEster();

  /*
     Doutd.print("doutd");
     Douts.print("douts");
     Doutrd.print("doutrd");
     Doutrs.print("doutrs");
  */

// roto di 45 gradi per controllare il moto rispetto ai piani inclinati...
  V_rot.setvn ( Firot* ( V_dev[0]+V_dev[1] ),Firot* ( V_dev[1]-V_dev[0] ), V_dev[2] );

  if ( ( nv=V_rot[0] ) >0.0 && Doutrs[0]>0.0 )  // piani x inclinati
    {t=Doutrs[0]/nv; if ( t>camm ) camm=t;}
  else if ( nv<0.0 && Doutrd[0]>0.0 )
    {t=-Doutrd[0]/nv;if ( t>camm ) camm=t;}



  if ( ( nv=V_rot[1] ) >0.0 && Doutrs[1]>0.0 )  // piani y inclinati
  {
    t=Doutrs[1]/nv;if ( t>camm ) camm=t;
  }
  else if ( nv<0.0 && Doutrd[1]>0.0 )
    {t=-Doutrd[1]/nv;if ( t>camm ) camm=t;}

  return camm;
}

//----------------------C a m I n t e r --------------------
// calcola il cammino interno per uscire dal rivelatore...
// vedi la routine nella classe base rettangolo.
// V_dev e' il versore della traiettoria nel dev
// V_rot e' il versore ruotato di 45 gradi per i lati inclinati
//
double DevOt::CamInter()
{
  double nv,t;
  // calcolo il cammino verso le otto +2 superfici
  // le distanze sono per definizione negative...
  // e scelgo il cammino positivo  maggiore
  // camm in uscita e' sempre definito
  // prima la classe base....
  camm=DevRt::CamInter();

// roto di 45 gradi per controllare il moto rispetto ai piani inclinati...
  V_rot.setvn ( Firot* ( V_dev[0]+V_dev[1] ),Firot* ( V_dev[1]-V_dev[0] ), V_dev[2] );

  if ( ( nv=V_rot[0] ) >0.0 )  // piani x    rotato
    {t=-Doutrd[0]/nv; if ( t<camm ) camm=t; }
  else if ( nv<0.0 )
    {t=Doutrs[0]/nv;  if ( t<camm ) camm=t; }
  if ( ( nv=V_rot[1] ) >0.0 )  // piani y    rotato
    {t=-Doutrd[1]/nv; if ( t<camm ) camm=t; }
  else if ( nv<0.0 )
    {t=Doutrs[1]/nv;  if ( t<camm ) camm=t; }

  return camm;
}