/[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.4 by pam-fi, Mon Jun 7 13:33:15 2010 UTC revision 1.1.2.6 by pam-fi, Mon Jun 7 14:30:20 2010 UTC
# Line 33  void RefitTrackAction::OnGood(PamLevel2 Line 33  void RefitTrackAction::OnGood(PamLevel2
33    if (!trkTrack)    if (!trkTrack)
34      return;      return;
35    
36      // Check if good clusters have meaningful energy releases
37      for (UInt_t ip = 0; ip < 6; ip++) {
38        if (trkTrack->XGood(ip) == 1 && trkTrack->GetDEDX(ip, 0) <= 0.)
39          trkTrack->xgood[ip] = 0; // Set it to bad
40        if (trkTrack->YGood(ip) == 1 && trkTrack->GetDEDX(ip, 1) <= 0.)
41          trkTrack->ygood[ip] = 0; // Set it to bad
42      }
43    
44    if (_currEvent != event->GetReadEntry()) {    if (_currEvent != event->GetReadEntry()) {
45      // Saves the original track data if the current event (recorded in the static member _currEvent) is      // Saves the original track data if the current event (recorded in the static member _currEvent) is
46      // different from the actually processed event. This way, multiple instances of the class      // different from the actually processed event. This way, multiple instances of the class
# Line 60  void RefitTrackAction::OnGood(PamLevel2 Line 68  void RefitTrackAction::OnGood(PamLevel2
68    
69    // Apply the gaussian spread    // Apply the gaussian spread
70    
   cout << GetName() << "  " << trkTrack->xm[0];  
   cout.flush();  
71    for (Int_t ip = 0; ip < 6; ip++) {    for (Int_t ip = 0; ip < 6; ip++) {
72      // gaussian distributed random number Gauss(mean,sigma) (in cm)      // gaussian distributed random number Gauss(mean,sigma) (in cm)
73      trkTrack->xm[ip] += _r3->Gaus(0., _spreadX);      trkTrack->xm[ip] += _r3->Gaus(0., _spreadX);
74      trkTrack->ym[ip] += _r3->Gaus(0., _spreadY);      trkTrack->ym[ip] += _r3->Gaus(0., _spreadY);
75    }    }
76    cout << GetName() << "  " << trkTrack->xm[0] << endl;  
77    // Mask the layers    // Mask the layers
78    if (_maskedViews) {    if (_maskedViews) {
79      if ((_maskedViews & T1X) == T1X)      if ((_maskedViews & T1X) == T1X)
Legend:
Removed from v.1.1.2.4  
changed lines
  Added in v.1.1.2.6
  ViewVC Help
Powered by ViewVC 1.1.23