FlyoDoc_2011 Pisa 2011 by GmP --- 011
|
00001 /*************************************************************************** 00002 devquad.cpp - description 00003 ------------------- 00004 begin : Thu Jan 9 2003 00005 copyright : (C) 2003 by Giuseppe Pierazzini 00006 email : peppe@unipi.it 00007 *************************************************************************** 00008 * * 00009 * NA48 simulation program. * 00010 * * 00011 ***************************************************************************/ 00012 00013 00014 #include "evento.h" 00015 #include "parm.h" 00016 #include "devquad.h" 00017 00018 //using namespace std; 00019 DevQuad::DevQuad(){ 00020 devtype=TypDevQuad; 00021 devclass="Quadrupol"; 00022 // quadropulo: campo quadrupolare ==> campo = 1 00023 campo=1; 00024 if(Dbref<=0.0) 00025 { Dbref=1.; 00026 Gout<<"\n In Quad "<<nome<<" Distance for field normalization in 'vfield' card not defined...forced to 1 cm."; 00027 00028 } 00029 Bref/=Dbref; 00030 } 00031 DevQuad::~DevQuad(){ 00032 } 00033 //---- 00034 00035 gvet& DevQuad::Get_Field(gvet &P) 00036 { 00037 // Il campo e' descritto da curve equipotenziali iperboliche. 00038 // l'intensita' e' normalizzata a B ad 1 cm dal centro. 00039 // l' intensita' cresce con la distanza dal centro. 00040 // P e' il raggio vettore. 00041 // Field.setvn( P[0]*Bref,-P[1]*Bref,0.); 00042 Field.setvn( P[1]*Bref,P[0]*Bref,0.); 00043 /* 00044 Puls.print("Puls"); 00045 P.print( evento_.Event, nome); 00046 Field.print( evento_.Event, "B"); 00047 */ 00048 return Field; 00049 } 00050