PndHypGeGammaAna.h

Go to the documentation of this file.

/******************************************************

Analysis Task created by M.Steinen steinen@kph.uni-mainz.de
Analysis of Gamma Simulation with hypGe detectors
*******************************************************/ 


#ifndef PNDHYPGEGAMMAANA_H
#define PNDHYPGEGAMMAANA_H 1

#include <FairTask.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TClonesArray.h>
#include <TRandom.h>
#include <fstream>
#include <set>
//#include "PndHypGeGammaAnaStorage.h"
#include "PndHypGePoint.h"
#include "../../pnddata/HypData/PndHypPoint.h"
#include "PndMCTrack.h"
#include <TSpectrum.h>
#include <TF1.h>
#include "TGeoManager.h"


class PndHypGeGammaAna: public FairTask
{
        public:
                PndHypGeGammaAna();
                PndHypGeGammaAna(TString TxtFileName_Ext, Double_t GammaEnergy_Ext,Int_t nEvents = 0,Int_t nPeaks_Ext=1,Int_t PeakToLook_Ext = 1);
                ~PndHypGeGammaAna();
        
                virtual InitStatus Init();
                virtual void Exec(Option_t* opt);
                void Finish();
                void SetParContainers();
                void SetTarget(Bool_t  TargetAvailable = 1);
                void SetOmegaQuadrupolMode(Bool_t useOmegaQuadrupolMode_ext = 1);
                void SetPeakWidtchStrechFactor(Double_t PeakWidtchStrechFactor_ext = 1);

        protected:
        
                Double_t                                                                        Compton(Double_t E, Double_t Th);
                Double_t                                                                        invCompton(Double_t E, Double_t Eprime);
                Double_t                                                                        SmearHit(Double_t Energy);
                Double_t                                                                        CalculatePeakWidth(Double_t Energy);
                void                                                                                    FitSpectrum();
                void                                                                                    FittingOutput();
                void                                                                                    FindPeaks();
                void                                                                                    GetPeakContent();
                void                                                                                    ActivateComptonCalculations(Bool_t isActivated = 1);
                void                                                                                    CalculateCompton();
                void                                                                                    WriteHistogramsToFile();
                void                                                                                    HistogramCosmetics(TH1D* histo, TString XTitle, Double_t XTitleSize, Double_t XTitleOffset, Double_t XLabelSize, TString YTitle, Double_t YTitleSize, Double_t YTitleOffset, Double_t YLabelSize);
                void                                                                                    HistogramCosmetics2D(TH2D* histo, TString XTitle, Double_t XTitleSize, Double_t XTitleOffset, Double_t XLabelSize, TString YTitle, Double_t YTitleSize, Double_t YTitleOffset, Double_t YLabelSize);

                bool                                                                                    InTargetAbsorbedGammaWouldHitGermanium();
                TString                                                                         TxtFileName;
                ofstream                                                                        TxtFile;
                
                Bool_t                                                                          useOmegaQuadrupolMode;

                Double_t                                                                        ResConstPar;
                Double_t                                                                        ResGradientPar;
                Double_t                                                                        PeakWidtchStrechFactor;
                Double_t                                                                        PeakWidth;

                Double_t                                                                        bufferEnergy;
                
                TH1D*                                                                                   hGamEnergy;
                TH1D*                                                                                   hGamEnergyTargetBackground;
                TH1D*                                                                                   hNumberOfHits;
                TH1D*                                                                                   hThetaCheck;
                TH1D*                                                                                   hThetaCheckPrimary;
                TH1D*                                                                                   hAbsorption;
                TH2D*                                                                                   h2AbsorptionDistanceAngle;
                
                TClonesArray*                                                   fMcTr;
                TClonesArray*                                                   fMc;
                TClonesArray*                                                   fGe;
                TClonesArray*                                                   fGeAl;
                TClonesArray*                                                   fSilicon;
                TClonesArray*                                                   fTargetOther;

                PndHypGePoint                                           *fHitGe;;
                PndHypPoint                                                     *fHitSi;
                PndHypPoint                                                     *fHitTargetOther;
                PndMCTrack                                                              *fMCGam;
                PndMCTrack                                                              *fMCGamPrimary;

                Long_t                                                                          EvtCount;
                Long_t                                                                          NumberOfEvents;
                Double_t                                                                        GammaEnergy;
                Double_t                                                                        iHistoUpperLimit;
                TString                                                                         fName;
                TRandom*                                                                        Rng;
                
                //PndHypGeGammaAnaStorage*      fStorage;
                
                FairRootManager*                                        ioman;
                // Exec
                std::set<int>                                           SetOfCrystalHit;
                TVector3                                                                        fStartVertex;
                TString                                                                         fVertexVolumeName;

                //Analysis
                Int_t                                                                           nPeaks;
                Int_t                                                                           PeakToLook;
                Double_t                                                                        PeakX;
                Double_t                                                                        PeakY; 
                Double_t                                                                        DPeakY; 
                Double_t                                                                        SumPeak;
                Double_t                                                                        DSumPeak;
                Bool_t                                                                          DoComptonCalculations;
                Double_t                                                                        lowComptonAngle;
                Double_t                                                                        highComptonAngle;
                Double_t                                                                        SumCompton;
                Double_t                                                                        DSumCompton;
                Double_t                                                                        PeakToCompton;
                Double_t                                                                        DPeakToCompton;
                
                TF1*                                                                                    GausBG;
                
                Float_t                                                                         *xpeaks;        
                Float_t                                                                         *ypeaks;        

                Bool_t                                                                          fTargetAvailable;



        ClassDef(PndHypGeGammaAna,1);
};

#endif /* PNDHYPGEGAMMAANA_H */