daqMC.h

Go to the documentation of this file.

#ifndef daqMC_h
#define daqMC_h


#include "HitMatrix.h"
#include "PhotonSource.h"
#include "sipmMC.h"

#include <TRandom3.h>
#include <TGraph.h>
#include <TH1.h>
#include <TF1.h>
#include <TGraph.h>
#include <TGraphErrors.h>
#include <TLegend.h>
#include <stdlib.h>
#include <vector>
#include <TROOT.h>
#include <string>
#include <RQ_OBJECT.h>

using std::vector;
using std::cout;
using std::endl;

struct GResonseCurve 
{
        TGraphErrors *response;
        TGraphErrors *resolution;
        TGraphErrors *responsePE;
        TGraphErrors *resolutionPE;
        TGraphErrors *responseDR;
        TGraphErrors *resolutionDR;
        TGraphErrors *responseAP;
        TGraphErrors *resolutionAP;
        TGraphErrors *responseCT;
        TGraphErrors *resolutionCT;
        TGraphErrors *responseENF;
        TGraphErrors *resolutionENF;
        TGraphErrors *responseEN;
        TGraphErrors *resolutionEN;
};

class daqMC{
        RQ_OBJECT("daqMC")

        public:

                daqMC();
                ~daqMC();

                void                    SetSiPM( sipmMC *mySiPM){ sipm = mySiPM; };
                void                    SetPhotonSource( PhotonSource *myPhotonSource ){ photonSource = myPhotonSource; };
                void                    Statistic( int N );                                                                     
                TGraph*                 Scope();                                                                                
                TH1D*                   QDCSpectrum( int N );                                                                   
                TGraph*                 ThreshScan( double gate, double tstart, double tstop, double tstep );                   
                TH1D*                   TDCSpectrum( int N );                                                                   
                GResonseCurve           DynamicRange( int N, double Ngamma_max, double Ngamma_step );                           
                double                  QDC( double charge );                                                                   
                void                    SetQDCChannels( int nChannels ){ h_QDC->SetBins(nChannels,0,nChannels); };
                void            SetPedestal( double Pedestal ){ pedestal = Pedestal; };
                void                    SetDiscriMinTime( double mintime ){ discriMinTime = mintime; };
                void                    SetDiscriWidth( double width ){ discriWidth = width; };

                TGraph*                 Discriminator( TGraph* waveform, double threshold );                                    
                //   void               ScanCT( int );                                                                          /**<Single pixel scan.*/

                void                    Progress( int ); //*SIGNAL*

                GResonseCurve           responseCurve;                                                                          
                TGraph                  *g_threshScan;                                                                          
                TH1D                    *h_pe;                                                                                  
                TH1D                    *h_dr;                                                                                  
                TH1D                    *h_xt;                                                                                  
                TH1D                    *h_ap;                                                                                  
                TH1D                    *h_QDC;                                                                                 
                TH1D                    *h_TDC;                                                                                 
                //   TGraph             *g_wf;                                                                                  /**<Oscilloscope waveform*/


                bool                    cancel;                                                                                 
                double                  progress;                                                                               
        private:

                bool                    Check();

                double                  pedestal;

                int                     nTDC_channels;

                TRandom3                r;
                HitMatrix*              hitMatrix;
                vector<double>          hit;
                sipmMC*                 sipm;
                PhotonSource*           photonSource;

                double                  discriMinTime;
                double                  discriWidth;

                int                     plast;

                TLegend*                legend;
                TLegend*                legend2;

};

#endif