/[PAMELA software]/PamCut/CollectionActions/RefitTrackAction/RefitTrackAction.cpp
ViewVC logotype

Diff of /PamCut/CollectionActions/RefitTrackAction/RefitTrackAction.cpp

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

revision 1.1.2.8 by pam-fi, Wed Nov 30 16:59:03 2011 UTC revision 1.2 by pam-fi, Thu Jul 8 14:00:14 2010 UTC
# Line 48  void RefitTrackAction::OnGood(PamLevel2 Line 48  void RefitTrackAction::OnGood(PamLevel2
48      _currEvent = event->GetReadEntry();      _currEvent = event->GetReadEntry();
49      // Save the original cluster status and measured positions      // Save the original cluster status and measured positions
50      for (int i = 0; i < 6; i++) {      for (int i = 0; i < 6; i++) {
51        _layerXBackup[i] = trkTrack->xgood[i];        _layerXBackup[i] = trkTrack->xgood[0];
52        _layerYBackup[i] = trkTrack->ygood[i];        _layerYBackup[i] = trkTrack->ygood[0];
53        _xmBackup[i] = trkTrack->xm[i];        _xmBackup[i] = trkTrack->xm[i];
54        _ymBackup[i] = trkTrack->ym[i];        _ymBackup[i] = trkTrack->ym[i];
55      }      }
# Line 60  void RefitTrackAction::OnGood(PamLevel2 Line 60  void RefitTrackAction::OnGood(PamLevel2
60      for (int i = 0; i < 6; i++) {      for (int i = 0; i < 6; i++) {
61        trkTrack->xm[i] = _xmBackup[i];        trkTrack->xm[i] = _xmBackup[i];
62        trkTrack->ym[i] = _ymBackup[i];        trkTrack->ym[i] = _ymBackup[i];
63        trkTrack->xgood[i] = _layerXBackup[i];        _layerXBackup[i] = trkTrack->xgood[0];
64        trkTrack->ygood[i] = _layerYBackup[i];        _layerYBackup[i] = trkTrack->ygood[0];
65      }      }
66    
67    }    }
68    
69    // Apply the gaussian spread    // Apply the gaussian spread
70    
71    if (_spreadX != 0.) {    for (Int_t ip = 0; ip < 6; ip++) {
72      for (Int_t ip = 0; ip < 6; ip++) {      // gaussian distributed random number Gauss(mean,sigma) (in cm)
73        // gaussian distributed random number Gauss(mean,sigma) (in cm)      trkTrack->xm[ip] += _r3->Gaus(0., _spreadX);
74        trkTrack->xm[ip] += _r3->Gaus(0., _spreadX);      trkTrack->ym[ip] += _r3->Gaus(0., _spreadY);
     }  
   }  
   if (_spreadY != 0.) {  
     for (Int_t ip = 0; ip < 6; ip++) {  
       trkTrack->ym[ip] += _r3->Gaus(0., _spreadY);  
     }  
75    }    }
76    
77    // Mask the layers    // Mask the layers
Legend:
Removed from v.1.1.2.8  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23