MFlagCoincidence.cc

Go to the documentation of this file.

#include "MFlagCoincidence.hh"
#include "QEvent.hh"
#include "QRawEvent.hh"
#include "QFiltersData.hh"
 
 
REGISTER_MODULE(MFlagCoincidence)
 
MFlagCoincidence::MFlagCoincidence(QSequence* s) : QModule("FlagCoincidence",s)
{
    // called once in a sequence
}
 
MFlagCoincidence::~MFlagCoincidence()
{
    // called once in a sequence
}
 
void MFlagCoincidence::Init()
{
   fIteration = GetIteration();
         
        if(fIteration == 1)
        {               
                fWindow = GetDouble("Window", 100.);
                fWindow *= 1000000.0; // Convert to ns
                
                fFlagName = GetString("FlagName", "coincidence");
                Info("Coincidence flag will be stored in AuxData variable %s",
                          fFlagName.c_str());
                
                fMultName = GetString("MultiplicityName", "multiplicity");
                Info("Multiplicity will be stored in AuxData variable %s", fMultName.c_str());
        }
        else
        {
                fnCoincidences = fCoincidentEvents.size();
                
                if(fnCoincidences == 0)
                        Panic("No coincidences!");
                
                Debug("%d elements in coincidence multiset.", fnCoincidences);
        }
}
 
QEvent* MFlagCoincidence::Do(QEvent* ev)
{       
        
        fCurrentEvent = ev->GetReadNumber();
        
        if(fIteration == 1)
        {
                if(!(ev->GetFilters().GetPassed())) return ev; 
        
                const QRawEvent& raw = ev->GetRawEvent();       
                QTime t;
                t = raw.GetTime(); 
                fCurrentTime = t.GetFromStartRunNs();
                bool isSignal = raw.GetPulse().GetIsSignal();
                        
                //Only include real events -- no heaters or noise acquisitions
                if(isSignal)
                {
                        fCoincidentEvents.insert(fCurrentEvent);
                        int i = fPreviousEvents.size()-1;
                        if(i>=0)
                        {
                                while((fCurrentTime-fPreviousTimes[i]) < fWindow)
                                {               
                                        fCoincidentEvents.insert(fPreviousEvents[i]);
                                        fCoincidentEvents.insert(fCurrentEvent);
                                        i--;
                                        if(i<0) break;
                                }
                        }
                        fPreviousEvents.push_back(fCurrentEvent);
                        fPreviousTimes.push_back(fCurrentTime);
                }
        }
        else
        {
                int multiplicity = fCoincidentEvents.count(fCurrentEvent);
                ev->AuxData().SetInt(fMultName, multiplicity, "save");  
                
                if(multiplicity>1)
                        ev->AuxData().SetInt(fFlagName, 1, "save");
                else
                        ev->AuxData().SetInt(fFlagName, 0, "save");
        }
        return ev;
}
 
void MFlagCoincidence::Done()
{
        if(fIteration == 1)
                SetRunAgain(true);
        else
                SetRunAgain(false);
}