CoincidenceFinder.hh

Go to the documentation of this file.

#include <stdlib.h>
#include <stdio.h>
#include <vector>
#include <list>
#include <utility>
#include <string>

#include "GammaEnergyChart.hh"
#include "FileEventParser.hh"
#include "EventHit.hh"
using namespace std;

#ifndef COINCIDENCE_FINDER_HH
#define COINCIDENCE_FINDER_HH //!<Inclusion guard
class CoincidenceFinder
{
public:
  CoincidenceFinder(FileEventParser * fParser 
                    ,double maxWait 
                    , int nTubes = 40 
                    , float detEps = 1E-2 
                    , float minGe = 1E-2 
                    , float maxGe = 1E20 
                    , unsigned int multiImplant = 1 
                    , float decayCut = 1E20 
                    ); 
  
  ~CoincidenceFinder(); 

  pair<EventHit*, EventHit*> getNextCoincidence(); 
  
  bool hasMoreCoincidences(); 

  int getFailureReason(int reason ); 

  int getNoHits(); 

  int getNoImplantations(); 

  int getNoDecays(); 

  int getNoOther(); 

  int getNoCorrectImplantations(); 

  int getNoDoubleHit(); 

  int getNoGeHit();

  int getNoNonAssociatedDecays(); 

  int getNoCorrConnected(); 

  int getNoErrConnected(); 

  int getNoDestroyedByOther(); 

  int getNoCorrCoincidence(); 

  int getNoFoundCoincidences(); 

  double getTimeSpan(); 

  vector<int> getDetectorCount(); 
  
private:
  int noHits; 

  int failureReasons[8]; 
  
  int noImplantations; 
  int noDecays; 
  int noOther; 

  int noCorrectImplantations; 
  int noDoubleHit; 

  int noGeHit; 
  int noNonAssociatedDecays; 
  int noCorrConnected; 
  int noErrConnected; 
  int noCorrCoincidence; 
  
  int noDestroyedByOther; 
  int noFoundCoincidences; 

  double timeSpan; 

  unsigned int multipleImplant; 
  
  float decayCutTime; 


  bool myMoreCoincidences; 

  vector<pair<list< EventHit*>, double> > myScintillatorTubes; //Contains the EventHits and the time a scintillator tube becomes availible for implantation again.

  FileEventParser * myParser; 

  void FindNextCoincidence(); 

  list<pair<EventHit*, EventHit*> > nextCoincidence; 

  void DoImplantation(EventHit * myHit 
                      ); 
  int DoDecay(EventHit * myHit 
               ); 
  void DoOther(EventHit * myHit 
               ); 

  double maxDecayWait; 
  int numberOfTubes; 
  float detectorEpsilon; 
  float minGeEnergy; 
  float maxGeEnergy; 
  vector<int> detectorCount; 
};

#endif