Collections/LiveTimeCollection/LiveTimeCollection.cpp

/*
 * LiveTimeCollection.cpp
 *
 *  Created on: 30-marth-2010
 *      Author: Vitaly Malakhov
 */

/*! @file LiveTimeCollection.cpp The LiveTimeCollection class implementation file */

#include "LiveTimeCollection.h"

void LiveTimeCollection::AddAction(CollectionAction *action) {
   _actions.push_back(action);
   _actionsPositions.push_back(GetSize() - 1);
   _sumLT.push_back(0);
   _sum.push_back(0);
}

void LiveTimeCollection::OnBad(PamLevel2 *event, int actNumber, unsigned int inc) {
   if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) || (_currenttime == 0)){
     for(unsigned int j = 0; j < _sumLT.size(); j++ ) _sumLT[j] = 0;
   } 
   _currenttime = event->GetOrbitalInfo()->absTime;
   unsigned int Inc = 0;
   if(inc == 0) Inc = event->GetTrigLevel2()->dltime[0]; else Inc = inc;
   for(unsigned int j = actNumber; j <= _actions.size()-1; j++) _sumLT[j] += Inc;
}

Float_t LiveTimeCollection::GetSumLiveTime(unsigned int actNumber) {
   if(actNumber >= _sum.size()) return -1; 
    return _sum[actNumber] * 0.16 /1000. ;
}

int LiveTimeCollection::ApplyCut(PamLevel2 *event) {

  _nEv++;

  // Execute the actions placed before the cuts
  unsigned int iBeforeCuts = 0;
  if (_actions.size() > 0) {
    while (_actionsPositions[iBeforeCuts] == -1) {
      _actions[iBeforeCuts]->OnGood(event);
      iBeforeCuts++;
      if (iBeforeCuts == _actions.size())
        break;
    }
  }

  // Apply the cuts
  if (_cuts.size() == 0) {
    _nGood++;
    OnGood(event);
    return CUTOK;
  }

  unsigned int firstFailed = _cuts.size();
  for (unsigned int iCut = 0; iCut < _cuts.size(); iCut++) {
    if (_cuts[iCut]->ApplyCut(event) != CUTOK && firstFailed == _cuts.size()) {
      firstFailed = iCut;
      break;
    }
  }

  Int_t actNumber = 0;
  unsigned int Inc = 0;

  //Do actions
  if (_actions.size() > 0) {
    Bool_t che = true;
    Int_t inc = 0;
    unsigned int lastPosition = _cuts.size();
    for (unsigned int i = iBeforeCuts; i < _actions.size(); i++) {
        actNumber = i-1;
      if (_actionsPositions[i] < (int) firstFailed){
         if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) || (_currenttime == 0)){
           for(unsigned int j = 0; j < _sumLT.size(); j++ ) _sumLT[j] = 0;
         }
         if( che ){
           inc = event->GetTrigLevel2()->dltime[0];
           Inc = inc;
           _sum[i] += _sumLT[i]+inc;
           event->GetTrigLevel2()->dltime[0] += _sumLT[i] ;
           che = false;
         }else _sum[i] += (_sumLT[i]+inc);
        _actions[i]->OnGood(event);
         _sumLT[i] = 0;
         _currenttime = event->GetOrbitalInfo()->absTime;
        }
      else { // A cut has failed
        if (lastPosition == _cuts.size()) // Record the position of the end of the bunch
          lastPosition = _actionsPositions[i];
        if (_actionsPositions[i] > (int) lastPosition) // Don't do actions at the end of successive bunches
          break;
        _actions[i]->OnBad(event, i);
      }
    }
  }

  if (firstFailed == _cuts.size()) {
    _nGood++;
    OnGood(event);
    return CUTOK;
  }
  else {
    OnBad(event, actNumber, Inc);
    return firstFailed;
  }
}

void LiveTimeCollection::Finalize() {
  ofstream livetimereport;
  livetimereport.open (_livetimereport);
  livetimereport << "Action name\t\tsum live time"<<endl;
  for (UInt_t i = 0; i < _sum.size(); i++) livetimereport << _actions[i]->GetName()<<"\t"<<GetSumLiveTime(i)<<"sec"<<endl;
  livetimereport.close();
}
1	pam-mep	1.1	/*
2			* LiveTimeCollection.cpp
3			*
4			* Created on: 30-marth-2010
5			* Author: Vitaly Malakhov
6			*/
7
8			/! @file LiveTimeCollection.cpp The LiveTimeCollection class implementation file /
9
10			#include "LiveTimeCollection.h"
11
12			void LiveTimeCollection::AddAction(CollectionAction *action) {
13			_actions.push_back(action);
14			_actionsPositions.push_back(GetSize() - 1);
15			_sumLT.push_back(0);
16			_sum.push_back(0);
17			}
18
19			void LiveTimeCollection::OnBad(PamLevel2 *event, int actNumber, unsigned int inc) {
20			if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) \|\| (_currenttime == 0)){
21			for(unsigned int j = 0; j < _sumLT.size(); j++ ) _sumLT[j] = 0;
22			}
23			_currenttime = event->GetOrbitalInfo()->absTime;
24			unsigned int Inc = 0;
25			if(inc == 0) Inc = event->GetTrigLevel2()->dltime[0]; else Inc = inc;
26			for(unsigned int j = actNumber; j <= _actions.size()-1; j++) _sumLT[j] += Inc;
27			}
28
29			Float_t LiveTimeCollection::GetSumLiveTime(unsigned int actNumber) {
30			if(actNumber >= _sum.size()) return -1;
31			return _sum[actNumber] * 0.16 /1000. ;
32			}
33
34			int LiveTimeCollection::ApplyCut(PamLevel2 *event) {
35
36			_nEv++;
37
38			// Execute the actions placed before the cuts
39			unsigned int iBeforeCuts = 0;
40			if (_actions.size() > 0) {
41			while (_actionsPositions[iBeforeCuts] == -1) {
42			_actions[iBeforeCuts]->OnGood(event);
43			iBeforeCuts++;
44			if (iBeforeCuts == _actions.size())
45			break;
46			}
47			}
48
49			// Apply the cuts
50			if (_cuts.size() == 0) {
51			_nGood++;
52			OnGood(event);
53			return CUTOK;
54			}
55
56			unsigned int firstFailed = _cuts.size();
57			for (unsigned int iCut = 0; iCut < _cuts.size(); iCut++) {
58			if (_cuts[iCut]->ApplyCut(event) != CUTOK && firstFailed == _cuts.size()) {
59			firstFailed = iCut;
60			break;
61			}
62			}
63
64			Int_t actNumber = 0;
65			unsigned int Inc = 0;
66
67			//Do actions
68			if (_actions.size() > 0) {
69			Bool_t che = true;
70			Int_t inc = 0;
71			unsigned int lastPosition = _cuts.size();
72			for (unsigned int i = iBeforeCuts; i < _actions.size(); i++) {
73			actNumber = i-1;
74			if (_actionsPositions[i] < (int) firstFailed){
75			if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) \|\| (_currenttime == 0)){
76			for(unsigned int j = 0; j < _sumLT.size(); j++ ) _sumLT[j] = 0;
77			}
78			if( che ){
79			inc = event->GetTrigLevel2()->dltime[0];
80			Inc = inc;
81			_sum[i] += _sumLT[i]+inc;
82			event->GetTrigLevel2()->dltime[0] += _sumLT[i] ;
83			che = false;
84			}else _sum[i] += (_sumLT[i]+inc);
85			_actions[i]->OnGood(event);
86			_sumLT[i] = 0;
87			_currenttime = event->GetOrbitalInfo()->absTime;
88			}
89			else { // A cut has failed
90			if (lastPosition == _cuts.size()) // Record the position of the end of the bunch
91			lastPosition = _actionsPositions[i];
92			if (_actionsPositions[i] > (int) lastPosition) // Don't do actions at the end of successive bunches
93			break;
94			_actions[i]->OnBad(event, i);
95			}
96			}
97			}
98
99			if (firstFailed == _cuts.size()) {
100			_nGood++;
101			OnGood(event);
102			return CUTOK;
103			}
104			else {
105			OnBad(event, actNumber, Inc);
106			return firstFailed;
107			}
108			}
109
110			void LiveTimeCollection::Finalize() {
111			ofstream livetimereport;
112			livetimereport.open (_livetimereport);
113			livetimereport << "Action name\t\tsum live time"<<endl;
114			for (UInt_t i = 0; i < _sum.size(); i++) livetimereport << _actions[i]->GetName()<<"\t"<<GetSumLiveTime(i)<<"sec"<<endl;
115			livetimereport.close();
116			}