/[PAMELA software]/PamCut/CollectionActions/Histo2DActions/TofDedxVsBetaHistoAction/TofDedxVsBetaHistoAction.cpp
ViewVC logotype

Diff of /PamCut/CollectionActions/Histo2DActions/TofDedxVsBetaHistoAction/TofDedxVsBetaHistoAction.cpp

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

revision 1.1 by pam-fi, Fri Sep 25 15:36:35 2009 UTC revision 1.1.2.2 by pam-fi, Tue Jan 3 08:45:27 2012 UTC
# Line 11  Line 11 
11    
12  TofDedxVsBetaHistoAction::TofDedxVsBetaHistoAction(const char *actionName, unsigned int layers, bool standAlone,  TofDedxVsBetaHistoAction::TofDedxVsBetaHistoAction(const char *actionName, unsigned int layers, bool standAlone,
13      TString outFileBase, TString mode, bool outRoot, bool outText, TString title) :      TString outFileBase, TString mode, bool outRoot, bool outText, TString title) :
14    Histo2DAction<Int_t> (actionName, title, outFileBase, mode, outRoot, outText), _layers(layers), _nLayers(0),      Histo2DAction<Int_t>(actionName, title, outFileBase, mode, outRoot, outText), _layers(layers), _nLayers(0), _badEvents(
15        _badEvents(0), _standAlone(standAlone) {          0), _standAlone(standAlone) {
16    
17      bool buildTitle = false;
18    if (title == "") {    if (title == "") {
19        buildTitle = true;
20      title = "TOF dE/dx (";      title = "TOF dE/dx (";
21      if ((_layers & S11) == S11) {    }
22      if ((_layers & S11) == S11) {
23        if (buildTitle)
24        title += "S11";        title += "S11";
25        _nLayers++;      _nLayers++;
26      }    }
27      if ((_layers & S12) == S12) {    if ((_layers & S12) == S12) {
28        if (buildTitle) {
29        if (_nLayers > 0)        if (_nLayers > 0)
30          title += ",";          title += ",";
31        title += "S12";        title += "S12";
       _nLayers++;  
32      }      }
33      if ((_layers & S21) == S21) {      _nLayers++;
34      }
35      if ((_layers & S21) == S21) {
36        if (buildTitle) {
37        if (_nLayers > 0)        if (_nLayers > 0)
38          title += ",";          title += ",";
39        title += "S21";        title += "S21";
       _nLayers++;  
40      }      }
41      if ((_layers & S22) == S22) {      _nLayers++;
42      }
43      if ((_layers & S22) == S22) {
44        if (buildTitle) {
45        if (_nLayers > 0)        if (_nLayers > 0)
46          title += ",";          title += ",";
47        title += "S22";        title += "S22";
       _nLayers++;  
48      }      }
49      if ((_layers & S31) == S31) {      _nLayers++;
50      }
51      if ((_layers & S31) == S31) {
52        if (buildTitle) {
53        if (_nLayers > 0)        if (_nLayers > 0)
54          title += ",";          title += ",";
55        title += "S31";        title += "S31";
       _nLayers++;  
56      }      }
57      if ((_layers & S32) == S32) {      _nLayers++;
58      }
59      if ((_layers & S32) == S32) {
60        if (buildTitle) {
61        if (_nLayers > 0)        if (_nLayers > 0)
62          title += ",";          title += ",";
63        title += "S32";        title += "S32";
       _nLayers++;  
64      }      }
65        _nLayers++;
66      }
67      if (buildTitle)
68      title += ") Vs Beta";      title += ") Vs Beta";
69    
70      SetTitle(title);    SetTitle(title);
   }  
71    
72  }  }
73    
# Line 64  void TofDedxVsBetaHistoAction::OnGood(Pa Line 78  void TofDedxVsBetaHistoAction::OnGood(Pa
78    static unsigned int badLayers;    static unsigned int badLayers;
79    
80    static int trkSeqNo;    static int trkSeqNo;
81      static int idx;
82    if (_standAlone)    if (_standAlone)
83      trkSeqNo = 0;      trkSeqNo = -1;
84    else    else
85      for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {      trkSeqNo = event->GetTrack(0)->GetTrkTrack()->GetSeqNo();
86        if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo())  
87          trkSeqNo = i;    for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) {
88      }      if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == trkSeqNo)
89          idx = i;
90      }
91    
92    dEdx = 0.;    dEdx = 0.;
93    badLayers = 0;    badLayers = 0;
94    if ((_layers & S11) == S11) {    if ((_layers & S11) == S11) {
95      dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 0, 100);      dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 0, 100);
96      if (dEdxLayer > 0. && dEdxLayer < 4090)      if (dEdxLayer > 0. && dEdxLayer < 4090)
97        dEdx += dEdxLayer;        dEdx += dEdxLayer;
98      else      else
99        badLayers++;        badLayers++;
100    }    }
101    if ((_layers & S12) == S12) {    if ((_layers & S12) == S12) {
102      dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 1, 100);      dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 1, 100);
103      if (dEdxLayer > 0. && dEdxLayer < 4090)      if (dEdxLayer > 0. && dEdxLayer < 4090)
104        dEdx += dEdxLayer;        dEdx += dEdxLayer;
105      else      else
106        badLayers++;        badLayers++;
107    }    }
108    if ((_layers & S21) == S21) {    if ((_layers & S21) == S21) {
109      dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 2, 100);      dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 2, 100);
110      if (dEdxLayer > 0. && dEdxLayer < 4090)      if (dEdxLayer > 0. && dEdxLayer < 4090)
111        dEdx += dEdxLayer;        dEdx += dEdxLayer;
112      else      else
113        badLayers++;        badLayers++;
114    }    }
115    if ((_layers & S22) == S22) {    if ((_layers & S22) == S22) {
116      dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 3, 100);      dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 3, 100);
117      if (dEdxLayer > 0. && dEdxLayer < 4090)      if (dEdxLayer > 0. && dEdxLayer < 4090)
118        dEdx += dEdxLayer;        dEdx += dEdxLayer;
119      else      else
120        badLayers++;        badLayers++;
121    }    }
122    if ((_layers & S31) == S31) {    if ((_layers & S31) == S31) {
123      dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 4, 100);      dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 4, 100);
124      if (dEdxLayer > 0. && dEdxLayer < 4090)      if (dEdxLayer > 0. && dEdxLayer < 4090)
125        dEdx += dEdxLayer;        dEdx += dEdxLayer;
126      else      else
127        badLayers++;        badLayers++;
128    }    }
129    if ((_layers & S32) == S32) {    if ((_layers & S32) == S32) {
130      dEdxLayer = event->GetToFLevel2()->GetdEdx(trkSeqNo, 5, 100);      dEdxLayer = event->GetToFLevel2()->GetdEdx(idx, 5, 100);
131      if (dEdxLayer > 0. && dEdxLayer < 4090)      if (dEdxLayer > 0. && dEdxLayer < 4090)
132        dEdx += dEdxLayer;        dEdx += dEdxLayer;
133      else      else
# Line 118  void TofDedxVsBetaHistoAction::OnGood(Pa Line 135  void TofDedxVsBetaHistoAction::OnGood(Pa
135    }    }
136    if (badLayers == 0) {    if (badLayers == 0) {
137      dEdx /= _nLayers - badLayers;      dEdx /= _nLayers - badLayers;
138      Fill(event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12], dEdx);      Fill(event->GetToFLevel2()->GetToFTrkVar(idx)->beta[12], dEdx);
139    }    }
140    else    else
141      _badEvents++;      _badEvents++;
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.1.2.2
  ViewVC Help
Powered by ViewVC 1.1.23