FlyoDoc_2011 Pisa 2011 by GmP --- 011

flyopunta/src/src_dev48/devmwc.cpp

00001 /***************************************************************************
00002                           devmwc.cpp  -  description
00003                              -------------------
00004     begin                : Sat May 26 2001
00005     copyright            : (C) 2001 by Giuseppe Pierazzini
00006     email                : peppe@unipi.it
00007   ***************************************************************************
00008   *                                                                         *
00009   *   NA48  simulation program.                                             *
00010   *                                                                         *
00011   ***************************************************************************/
00012 #include "parm.h"
00013 #include "devmwc.h"
00014 #include "particella.h"
00015 
00016 //using namespace std;
00017 
00018 void DevMwc::Prgeom()
00019 {
00020         Gout<<"==================  R i v e l a t o r e =================";
00021         Gout<<"\n --> < "<<nome <<" >  [ "<<fun
00022         <<" ]  Wire Chamer  == Corona ";
00023         Device::Prgeom();
00024 
00025 }
00026 //=======-----------------------------------
00027 void DevMwc::DataSmear()
00028 {
00029 //  calcolo degli errori
00030 // Opera sulle variabili di tipo (2)
00031 // vedi Device::Fill_Xdevhit()
00032         int i=0;
00033         double norma;
00034         sepm=120.;
00035         while ( i<=nhit )
00036         {
00037                 x=M_Hits[i].Xdev.x;
00038                 y=M_Hits[i].Xdev.y;
00039 
00040 // errori di 120 micron.
00041                 xc=M_Hits[i].Xdev.x=Pgauss ( 0.012,x );
00042                 yc=M_Hits[i].Xdev.y=Pgauss ( 0.012,y );
00043                 M_Hits[i].Xlab=Devc2Lab ( M_Hits[i].Xdev );
00044 
00046                 /*  e il tempo? se si vuole attivare con errore fare qualche cosa di simile
00047                     tm_visto[i]=tm_visto[i] + errore
00048                 */
00050                 norma=xc*xc+yc*yc;
00051 // distanza minima dal centro
00052                 if ( norma<sepm ) sepm=norma;
00053 
00054                 i++;
00055         }
00056 }
 All Classes Namespaces Files Functions Variables