/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkLevel2.cpp
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/TrkLevel2.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.38 by pam-fi, Fri Aug 17 13:25:14 2007 UTC revision 1.39 by pam-fi, Wed Aug 22 07:03:45 2007 UTC
# Line 421  Float_t TrkTrack::GetLnLY(){ Line 421  Float_t TrkTrack::GetLnLY(){
421      return lnl;      return lnl;
422            
423  }  }
424    /**
425     * Returns the effective angle, relative to the sensor, on each plane.
426     * @param ip plane (0-5)
427     * @param iv view (0=x 1=y)
428     */
429    Float_t TrkTrack::GetEffectiveAngle(int ip, int iv){
430    
431        if(ip<0 || ip>5){
432            cout << "Float_t TrkTrack::GetEffectiveAngle(int "<<ip<<", int "<<iv<<") ==> wrong input"<<endl;
433            return 0.;
434        }
435    
436        float v[3]={xv[ip],yv[ip],zv[ip]};
437        //-----------------------------------------
438        // effective angle (relative to the sensor)
439        //-----------------------------------------
440        float axv_geo  = axv[ip];
441        float muhall_h = 297.61; //cm**2/Vs
442        float BY = TrkParams::GetBY(v);
443        float axv_eff = 0;
444        if(ip==5) axv_geo = -1*axv_geo;
445        if(ip==5) BY      = -1*BY;
446        axv_eff = 180.*atan( tan(axv_geo*acos(-1.)/180.) + muhall_h * BY * 0.0001)/acos(-1.);
447        //-----------------------------------------
448        // effective angle (relative to the sensor)
449        //-----------------------------------------
450        float ayv_geo = ayv[ip];
451        float muhall_e = 1258.18; //cm**2/Vs
452        float BX = TrkParams::GetBX(v);
453        float ayv_eff = 0;
454        ayv_eff = 180.*atan( tan(ayv_geo*acos(-1.)/180.) + muhall_e * BX * 0.0001)/acos(-1.);
455      
456        if     (iv==0)return axv_eff;
457        else if(iv==1)return ayv_eff;
458        else{
459            cout << "Float_t TrkTrack::GetEffectiveAngle(int "<<ip<<", int "<<iv<<") ==> wrong input"<<endl;
460            return 0.;
461        }
462      
463    };
464    
465  //--------------------------------------  //--------------------------------------
466  //  //
467  //  //
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.39
  ViewVC Help
Powered by ViewVC 1.1.23