FlyoDoc_2011 Pisa 2011 by GmP --- 011
|
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 }