/[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.37 by pam-fi, Wed Jun 6 09:36:07 2007 UTC revision 1.38 by pam-fi, Fri Aug 17 13:25:14 2007 UTC
# Line 397  Float_t TrkTrack::GetLnLX(){ Line 397  Float_t TrkTrack::GetLnLX(){
397              lnl += (tailx[ip]+1.) * log( (tailx[ip]*pow(resx[ip],2.) + pow(xv[ip]-xm[ip],2.)) / (tailx[ip]*pow(resx[ip],2)) );              lnl += (tailx[ip]+1.) * log( (tailx[ip]*pow(resx[ip],2.) + pow(xv[ip]-xm[ip],2.)) / (tailx[ip]*pow(resx[ip],2)) );
398      if(GetNX()>3)lnl=lnl/(GetNX()-3);      if(GetNX()>3)lnl=lnl/(GetNX()-3);
399      else lnl=0;      else lnl=0;
400      if(lnl==0)cout << " Float_t  TrkTrack::GetLnLX() -- WARNING -- value not defined "<<lnl<<endl;      if(lnl==0){
401            cout << " Float_t  TrkTrack::GetLnLX() -- WARNING -- value not defined "<<lnl<<endl;
402            Dump();
403        }
404      return lnl;      return lnl;
405            
406  }  }
# Line 411  Float_t TrkTrack::GetLnLY(){ Line 414  Float_t TrkTrack::GetLnLY(){
414              lnl += (taily[ip]+1.) * log( (taily[ip]*pow(resy[ip],2.) + pow(yv[ip]-ym[ip],2.)) / (taily[ip]*pow(resy[ip],2)) );              lnl += (taily[ip]+1.) * log( (taily[ip]*pow(resy[ip],2.) + pow(yv[ip]-ym[ip],2.)) / (taily[ip]*pow(resy[ip],2)) );
415      if(GetNY()>2)lnl=lnl/(GetNY()-2);      if(GetNY()>2)lnl=lnl/(GetNY()-2);
416      else lnl=0;      else lnl=0;
417      if(lnl==0)cout << " Float_t  TrkTrack::GetLnLY() -- WARNING -- value not defined "<<lnl<<endl;      if(lnl==0){
418            cout << " Float_t  TrkTrack::GetLnLY() -- WARNING -- value not defined "<<lnl<<endl;
419            Dump();
420        }
421      return lnl;      return lnl;
422            
423  }  }
Legend:
Removed from v.1.37  
changed lines
  Added in v.1.38
  ViewVC Help
Powered by ViewVC 1.1.23