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);
   _badsum.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++ ){
       if(_currenttime == 0) _badsum[j] = 0; else {_badsum[j] += _sumLT[j]; _gapNumb++;}
        _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. ;
}

Float_t LiveTimeCollection::GetBadsumLiveTime(unsigned int actNumber) {
   if(actNumber >= _sum.size()) return -1; 
    return _badsum[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;
      if (_actionsPositions[i] < (int) firstFailed){
         if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) || (_currenttime == 0)){
           for(unsigned int j = 0; j < _sumLT.size(); j++ ) {
             if(_currenttime == 0) _badsum[j] = 0; else {_badsum[j] += _sumLT[j]; _gapNumb++;}
             _sumLT[j] = 0;
           }
         }
         if( che ){
           inc = event->GetTrigLevel2()->dltime[0];
           Inc = inc;
        }
         _sum[i] += (_sumLT[i]+inc);
         event->GetTrigLevel2()->dltime[0] += _sumLT[i] + inc;
         che = false;
        _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() {
  SmartCollection::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"<<"\t"<<GetBadsumLiveTime(i)<<"\t"<<_gapNumb<<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	pam-mep	1.2	_badsum.push_back(0);
18	pam-mep	1.1	}
19
20			void LiveTimeCollection::OnBad(PamLevel2 *event, int actNumber, unsigned int inc) {
21			if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) \|\| (_currenttime == 0)){
22	pam-mep	1.2	for(unsigned int j = 0; j < _sumLT.size(); j++ ){
23			if(_currenttime == 0) _badsum[j] = 0; else {_badsum[j] += _sumLT[j]; _gapNumb++;}
24			_sumLT[j] = 0;
25			}
26	pam-mep	1.1	}
27			_currenttime = event->GetOrbitalInfo()->absTime;
28			unsigned int Inc = 0;
29			if(inc == 0) Inc = event->GetTrigLevel2()->dltime[0]; else Inc = inc;
30			for(unsigned int j = actNumber; j <= _actions.size()-1; j++) _sumLT[j] += Inc;
31			}
32
33			Float_t LiveTimeCollection::GetSumLiveTime(unsigned int actNumber) {
34			if(actNumber >= _sum.size()) return -1;
35			return _sum[actNumber] * 0.16 /1000. ;
36			}
37
38	pam-mep	1.2	Float_t LiveTimeCollection::GetBadsumLiveTime(unsigned int actNumber) {
39			if(actNumber >= _sum.size()) return -1;
40			return _badsum[actNumber] * 0.16 /1000. ;
41			}
42
43	pam-mep	1.1	int LiveTimeCollection::ApplyCut(PamLevel2 *event) {
44
45			_nEv++;
46
47			// Execute the actions placed before the cuts
48			unsigned int iBeforeCuts = 0;
49			if (_actions.size() > 0) {
50			while (_actionsPositions[iBeforeCuts] == -1) {
51			_actions[iBeforeCuts]->OnGood(event);
52			iBeforeCuts++;
53			if (iBeforeCuts == _actions.size())
54			break;
55			}
56			}
57
58			// Apply the cuts
59			if (_cuts.size() == 0) {
60			_nGood++;
61			OnGood(event);
62			return CUTOK;
63			}
64
65			unsigned int firstFailed = _cuts.size();
66			for (unsigned int iCut = 0; iCut < _cuts.size(); iCut++) {
67			if (_cuts[iCut]->ApplyCut(event) != CUTOK && firstFailed == _cuts.size()) {
68			firstFailed = iCut;
69			break;
70			}
71			}
72
73			Int_t actNumber = 0;
74			unsigned int Inc = 0;
75
76			//Do actions
77			if (_actions.size() > 0) {
78			Bool_t che = true;
79			Int_t inc = 0;
80			unsigned int lastPosition = _cuts.size();
81			for (unsigned int i = iBeforeCuts; i < _actions.size(); i++) {
82	pam-mep	1.2	actNumber = i;
83	pam-mep	1.1	if (_actionsPositions[i] < (int) firstFailed){
84			if((event->GetOrbitalInfo()->absTime - _currenttime > 1000.) \|\| (_currenttime == 0)){
85	pam-mep	1.2	for(unsigned int j = 0; j < _sumLT.size(); j++ ) {
86			if(_currenttime == 0) _badsum[j] = 0; else {_badsum[j] += _sumLT[j]; _gapNumb++;}
87			_sumLT[j] = 0;
88			}
89	pam-mep	1.1	}
90			if( che ){
91			inc = event->GetTrigLevel2()->dltime[0];
92			Inc = inc;
93	pam-mep	1.2	}
94			_sum[i] += (_sumLT[i]+inc);
95			event->GetTrigLevel2()->dltime[0] += _sumLT[i] + inc;
96			che = false;
97	pam-mep	1.1	_actions[i]->OnGood(event);
98			_sumLT[i] = 0;
99			_currenttime = event->GetOrbitalInfo()->absTime;
100			}
101			else { // A cut has failed
102			if (lastPosition == _cuts.size()) // Record the position of the end of the bunch
103			lastPosition = _actionsPositions[i];
104			if (_actionsPositions[i] > (int) lastPosition) // Don't do actions at the end of successive bunches
105			break;
106			_actions[i]->OnBad(event, i);
107			}
108			}
109			}
110
111			if (firstFailed == _cuts.size()) {
112			_nGood++;
113			OnGood(event);
114			return CUTOK;
115			}
116			else {
117			OnBad(event, actNumber, Inc);
118			return firstFailed;
119			}
120			}
121
122			void LiveTimeCollection::Finalize() {
123	pam-mep	1.2	SmartCollection::Finalize();
124	pam-mep	1.1	ofstream livetimereport;
125			livetimereport.open (_livetimereport);
126			livetimereport << "Action name\t\tsum live time"<<endl;
127	pam-mep	1.2	for (UInt_t i = 0; i < _sum.size(); i++) livetimereport << _actions[i]->GetName()<<"\t"<<GetSumLiveTime(i)<<"sec"<<"\t"<<GetBadsumLiveTime(i)<<"\t"<<_gapNumb<<endl;
128	pam-mep	1.1	livetimereport.close();
129			}