QTriggerThreshold.cc

Go to the documentation of this file.

#include "QTriggerThreshold.hh"
#include "QReader.hh"
#include <vector>
 
ClassImp(Diana::QTriggerThreshold);
using std::cout;
using std::endl;
 
Q_BEGIN_NAMESPACE
 
void QTriggerThreshold::GetParametersFromCFG(){
    if(reader==0){
        cout<<"NO READER PASSED TO TRIGGER"<<endl;
        exit(-1);
    }
    
    SetCFGParameters(reader->GetDouble(("ThresholdTrigLVL"+triggerLabel).c_str(),600.),reader->GetBool(("SubtractBaseline"+triggerLabel).c_str(),false));
}
 
void QTriggerThreshold::TriggerStream(const vector<double>& stream){
    if (!read_parameters){
            GetParametersFromCFG();
    }
    //empty triggers
    triggers.clear();
    triggerHeights.clear();
    
    if(samplingFrequency ==0 or windowsize==0)cout<<"Sampling frequecy "<< samplingFrequency<<"\tWindow Size = "<<windowsize<<endl;    
    
    size_t first_trig=stream.size()+1;
    size_t prev_trig=first_trig;
    double height =0;
    int sample;
    double baseline=0;
    int samps=0;
    filteredStream=stream;
    if(subtract_baseline){
       for(size_t i =0;(double) i<0.1*(double) windowsize and i<stream.size();i++){
           samps++;
           baseline +=stream[i];
       }
       baseline/=samps;
    }else{
       samps=1;
    }
    for(size_t i=0; i<stream.size();i++){
        if(stream[i]-baseline/samps>=triggerLevel){//If sample triggers
            if(first_trig==stream.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 stream[i]<triggerLevel) or(i == stream.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 stream
            if(subtract_baseline){
                baseline = stream[i];
                samps=1;
            }
            if(triggers.size()==0){//if first trigger
                for(size_t j =first_trig; j<=prev_trig;j++){
                    if(j==first_trig){
                        height=stream[j];
                        sample=(int)j;
                    }else{
                        if(height<stream[j]){
                            height=stream[j];
                            sample=(int)j;
                        }
                    }
                }
                size_t min = first_trig -1;
                if(min<0){
                    min = i;
                }
                if(min >= stream.size()){
                    min = first_trig;
                }
                triggers.push_back(sample);
                triggerHeights.push_back(height-stream[min]);
            }else{
                if(subtract_baseline){
                        baseline += stream[i];
                        samps++;
                }
                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=stream[j];
                            sample=(int)j;
                        }else{
                            if(height<stream[j]){
                                height=stream[j];
                                sample=(int)j;
                            }
                        }
                    }
                    triggers.push_back(sample);
                    size_t min = first_trig -1;
                    if(min<0){
                        min = i;
                    }
                    if(min >= stream.size()){
                        min = first_trig;
                    }
                    triggerHeights.push_back(height-stream[min]);
                }
            }
        }
    }
}
 
Q_END_NAMESPACE