PandaRoot
PndHypGeCOSYBackgroundAna.h
Go to the documentation of this file.
1 //****************************************************************************
2 //* This file is part of PandaRoot. *
3 //* *
4 //* PandaRoot is distributed under the terms of the *
5 //* GNU General Public License (GPL) version 3, *
6 //* copied verbatim in the file "LICENSE". *
7 //* *
8 //* Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot *
9 //* The copyright holders are listed in the file "COPYRIGHTHOLDERS". *
10 //* The authors are listed in the file "AUTHORS". *
11 //****************************************************************************
12 
13 /******************************************************
14 
15 Analysis Task created by M.Steinen steinen@kph.uni-mainz.de
16 Analysis of Gamma Simulation with hypGe detectors
17 *******************************************************/
18 
19 #ifndef PNDHYPGECOSYBACKGROUNDANA_H
20 #define PNDHYPGECOSYBACKGROUNDANA_H 1
21 
22 #include <FairTask.h>
23 #include <TH1D.h>
24 #include <TH2D.h>
25 #include <TVector3.h>
26 #include <TString.h>
27 #include <TClonesArray.h>
28 #include <TRandom.h>
29 #include <fstream>
30 #include <TGeoManager.h>
31 #include <TROOT.h>
32 //#include "PndHypGeCOSYBackgroundAnaStorage.h"
33 #include <TSpectrum.h>
34 #include <TF1.h>
35 #include <TCanvas.h>
36 
37 class PndHypGeCOSYBackgroundAna : public FairTask {
38  public:
39  PndHypGeCOSYBackgroundAna();
40  PndHypGeCOSYBackgroundAna(Int_t nEvents);
41  ~PndHypGeCOSYBackgroundAna();
42 
43  virtual InitStatus Init();
44  virtual void Exec(Option_t *opt);
45  void Finish();
46 
47  protected:
48  TGeoManager *fgeom;
49 
50  TClonesArray *fMcTr;
51  TClonesArray *fHyp;
52  TClonesArray *fHypGe;
53 
54  TH1D *hNHits;
55  TH1D *hCrystalHit;
56  TH1D *hNeutronOrigin;
57  TH2D *hCrystalOrigin;
58  TH1D *hNeutronEkin;
59  TH2D *hNeutronEkinOrigin;
60  TH1D *hNeutronEnergyLoss;
61 
62  TH1D *hAllParticlesGe;
63  TH1D *hAllParticlesCrystal1;
64  TH1D *hAllParticlesCrystal2;
65  TH1D *hAllParticlesPiezo;
66  TH1D *hAllParticlesSiPm1;
67  TH1D *hAllParticlesSiPm2;
68 
69  TH2D *hEkinAllParticles;
70  TH2D *hEkinAllParticlesCrystal1;
71  TH2D *hEkinAllParticlesCrystal2;
72  TH2D *hEkinAllParticlesPiezo;
73  TH2D *hEkinAllParticlesSiPm1;
74  TH2D *hEkinAllParticlesSiPm2;
75 
76  TH2D *hGammaEkinOrigin;
77 
78  TCanvas *cNHits;
79  TCanvas *cCrystalHit;
80  TCanvas *cNeutronOrigin;
81  TCanvas *cCrystalOrigin;
82  TCanvas *cNeutronEkin;
83  TCanvas *cNeutronEkinOrigin;
84  TCanvas *cNeutronEnergyLoss;
85  TCanvas *cAllParticlesGe;
86  TCanvas *cAllParticlesCrystal1;
87  TCanvas *cAllParticlesCrystal2;
88  TCanvas *cEkinAllParticles;
89  TCanvas *cGammaEkinOrigin;
90 
91  TCanvas *cAllParticlesPiezo;
92  TCanvas *cAllParticlesSiPm1;
93  TCanvas *cAllParticlesSiPm2;
94  TCanvas *cEkinAllParticlesSplit;
95 
96  TVector3 StartVertex;
97  TString VertexVolumeName;
98 
99  TVector3 NeutronMomentum;
100  Double_t NeutronEkin;
101  TVector3 ParticleMomentum;
102  Double_t ParticleEkin;
103 
104  Long_t EvtCount;
105  Long_t NumberOfEvents;
106 
107  TString fName;
108 
109  // Analysis
110 
111  ClassDef(PndHypGeCOSYBackgroundAna, 1);
112 };
113 
114 #endif /* PNDHYPGECOSYBACKGROUNDANA_H */
PndHypGeCOSYBackgroundAna::ClassDef
ClassDef(PndHypGeCOSYBackgroundAna, 1)
PndHypGeCOSYBackgroundAna::Init
virtual InitStatus Init()
PndHypGeCOSYBackgroundAna::Exec
virtual void Exec(Option_t *opt)