|
FlyoDoc_2011 Pisa 2011 by GmP --- 011
|
#include <devfluka.h>
Public Member Functions | |
| DevFluka () | |
| int | Traccia (Particella *) |
| ================ T R A C C I A ==================================== | |
| virtual void | Prgeom () |
Protected Attributes | |
| double | Lato_z |
| double | zintgas |
Definition at line 24 of file devfluka.h.
| DevFluka::DevFluka | ( | ) |
The "buco", of the Fluka device defines a cylindric space that corresponds to the decay region. The beam particles entering the region interact with the gas and in Trasport the tracing is stopped with the production of a cluster.
Definition at line 23 of file devfluka.cpp.
{
//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 | ) | [virtual] |
================ T R A C C I A ====================================
Traccia le particelle nell'apparato per un tratto come definito in Get_camm(). Nel tracciare la particella si tiene conto della possibilita' della pair production e del bremsstrhalung.
Reimplemented from Device.
Definition at line 33 of file devfluka.cpp.
{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;
}