flyopunta/src/src_dev/devrtrt.cpp

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


DevRtRt::~DevRtRt(){
}
//=================DevRtRt==============================
//----------------- --------
void DevRtRt::Prgeom()
 {
  Gout<<"==================  R i v e l a t o r e =================";
  Gout<<"\n --> < "<<nome <<" >  [ "<<fun
  <<" ]  Rettang. con buco rettang.";
  Device::Prgeom();
  Gout<<std::endl;
}
//======================================================
//----------------- -----P o s i z i o n e ---
int DevRtRt::Posizione()
{
  if( DevRt::Posizione()>0) return 1; //  il punto e' esterno

  // se sono interno al rivelatore verifico se siamo nel buco

  P_Cin=X_dev-Centrin;   // mi riferisco al centro del buco interno
  Dind=P_Cin-Lin;
  Dins=-(Lin+P_Cin);
  
//  if(idev==22){
//   X_dev.print( evento_.Event,"XDev"); }


//   P_Cin.print(evento_.Event,"PCIN");
//   Dins.print(evento_.Event,nome);
//   Dind.print(evento_.Event,nome) ;    }
   
  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 I n t e r --------------------

double DevRtRt::CamInter()
{
//uso quella di DevRt per calcolare il cammino
// interno 
  tratto=DevRt::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 DevRtRt::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;
}