QTriggerTTL.cc

Go to the documentation of this file.

#include "QTriggerTTL.hh"
#include "QReader.hh"
#include <vector>
 
ClassImp(Diana::QTriggerTTL);
using std::cout;
using std::endl;
 
Q_BEGIN_NAMESPACE
 
void QTriggerTTL::GetParametersFromCFG(){
    if(reader==0){
        cout<<"NO READER PASSED TO TRIGGER"<<endl;
        exit(-1);
    }
    
    SetCFGParameters(reader->GetDouble(("TTLThreshold"+triggerLabel).c_str(),600.),reader->GetBool(("TTLUseDiff"+triggerLabel).c_str(),true));
}
 
void QTriggerTTL::TriggerStream(const vector<double>& stream){
    //if parameters not read from cfg
    if (!read_parameters){
            GetParametersFromCFG();
    }
    //empty triggers
    triggers.clear();
    triggerHeights.clear();
    
    //check sampling frequency and windowsize
    if(samplingFrequency ==0 or windowsize==0)cout<<"Sampling frequecy "<< samplingFrequency<<"\tWindow Size = "<<windowsize<<endl;
    
    //start triggering
    size_t first_trig=stream.size()+1;
    size_t prev_trig=first_trig;
    double height =0;
    int sample;
    filteredStream.resize(stream.size());
    /*
    double mean = 0;
    for(size_t i=0; i<stream.size();i++){
        mean +=stream[i];
    }
    mean/=stream.size();
    for(size_t i=0; i<stream.size();i++){
        filteredStream[i] = stream[i]-mean;
    }
    */
    if(fuseDiff==true){
        for(size_t i=1; i<stream.size();i++){
                filteredStream[i] = stream[i]-stream[i-1];
        }
 
        filteredStream[0]=filteredStream[1];
    }else{
        filteredStream=stream;
    }
    for(size_t i=0; i<filteredStream.size();i++){
        if(filteredStream[i]>=triggerLevel){//If sample triggers
            if(first_trig==filteredStream.size()+1){//first trigger point in chunk
                first_trig = i;
                prev_trig = i;
            }else if (i==prev_trig+1){//same peak
                prev_trig=i;
            }else if(i>prev_trig+1){//new peak
                first_trig=i;
                prev_trig=i;
            }
        } 
        
        if((i==prev_trig+1 and filteredStream[i]<triggerLevel) or(i == filteredStream.size()-1 and i ==prev_trig)){//as soon as the peak ends or if the peak is ongoing and we are at the end of the filteredStream
            if(triggers.size()==0){//if first trigger
                for(size_t j =first_trig; j<=prev_trig;j++){
                    if(j==first_trig){
                        height=filteredStream[j];
                        sample=(int)j;
                    }else{
                        if(height<filteredStream[j]){
                            height=filteredStream[j];
                            sample=(int)j;
                        }
                    }
                }
                size_t min = first_trig -5;
                if(min<0){
                    min = i;
                }
                if(min >= filteredStream.size()){
                    min = first_trig;
                }
                triggers.push_back(sample);
                triggerHeights.push_back(stream[sample]-stream[min]);
            }else{
                if((double)((int)first_trig-triggers[triggers.size()-1]) >= (double)debounce){//don't trigger if trigger is too near the previous one
                    for(size_t j =first_trig; j<=prev_trig;j++){
                        if(j==first_trig){
                            height=filteredStream[j];
                            sample=(int)j;
                        }else{
                            if(height<filteredStream[j]){
                                height=filteredStream[j];
                                sample=(int)j;
                            }
                        }
                    }
                    triggers.push_back(sample);
                    size_t min = first_trig -5;
                    if(min<0){
                        min = i;
                    }
                    if(min >= stream.size()){
                        min = first_trig;
                    }
                    triggerHeights.push_back(stream[sample]-stream[min]);
                }
            }
        }
    }
}
 
Q_END_NAMESPACE