flyopunta/src/src_dev/devfluka.cpp

/***************************************************************************
 *   Copyright (C) 2005 by Giuseppe Pierazzini          Pisa 20.02.05      *
 *   pierazzini@unipi.it                                                   *
 *                                     F l y o                             *
 *   Epsi/NA48                                                             *
 ***************************************************************************/

#include "flyoh.h"
#define ESTERNO  1
#define INTERNO  0
#define INBUCO   -1

using namespace std;

DevFluka::DevFluka()
{
  //tipico device per lo studio delle interazioni co il gasss
  devtype=TypDevGas;
  devclass="Dec_Tubo";
  fun =Gasint;
  Lato_z=2.*Lout.z;
}


int DevFluka::Traccia ( Particella *pr )
{double tratok=0;

  // particella stabile 
  
  if(pr->Get_Ido()==1 && pr->Get_Ctau()==-1.0 &&nowpos==-1)     
  {
  //tracing in the device "buco".....
  // define the interaction point inside the device,..
  // even if the particle is stable. Mark the particle as decay
//    cout<<"\n Ev0 "<< Eventi_Fatti<< " "<<pr->Get_Nome()<<" f "<<pr->Get_Fato() <<" z "<<setw(8)<<pr->X.z<<"  pathmax "<<setw(8)<< pr->path_max
//    <<" done "<<setw(8)<<pr->path_done<<" pathok "<<setw(8)<<pr->pathok<<" last "<<setw(8)<<pr->last_path<<flush;
   
   
    if(pr->path_max==Zona)
    {
   tratok= Pran()* Lato_z;
    pr->path_max=Cface.z+ tratok;
    }
    
    if(pr->pathok > (  pr->path_max- pr->path_done ) )
        {
            pr->pathok= ( pr->path_max-pr->path_done+0.000001 );
           pr->Set_Fato(Decay,-2);   
        
        }

//    cout<<"\n Ev1 "<< Eventi_Fatti<< " "<<pr->Get_Nome()<<" f "<<pr->Get_Fato()<<" z "<<setw(8)<<pr->X.z<<"  pathmax "<<setw(8)<< pr->path_max
//    <<" done "<<setw(8)<<pr->path_done<<" pathok "<<setw(8)<<pr->pathok<<" last "<<setw(8)<<pr->last_path<<flush; 
 
  }
    
int ok = Device::Traccia(pr);

 return ok;
}



void DevFluka::Prgeom()
{
//----------------- --------
  Gout<<"==================  R i v e l a t o r e =================";
  Gout<<"\n --> < "<<nome <<" >  [ "<<fun
  <<" ]  Decay region Dev";

  Device::Prgeom();

// ridefinizione dei buffer per contenere gli hits se troppi

}