NTuple2EventTree/Converter.hh at main · GRIFFINCollaboration/NTuple2EventTree · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
#ifndef __CONVERTER_HH
#define __CONVERTER_HH

#include <vector>

#include "TChain.h"
#include "TFile.h"
#include "TTree.h"
#include "TRandom3.h"
#include "TVector3.h"

#include "TRunInfo.h"
#include "TChannel.h"
#include "TFragment.h"
#include "TGriffin.h"
#include "TGriffinBgo.h"
#include "TSceptar.h"
#include "TPaces.h"
#include "TLaBr.h"
#include "TDescant.h"

#include "Settings.hh"

class Converter {
public:
	Converter(std::vector<std::string>& inputFileNames, const int& runNumber, const int& subRunNumber, const TRunInfo* runInfo, Settings* settings, bool writeFragmentTree);
	~Converter();

	bool Run();

private:
	int  Cfd(EDigitizer);
	bool AboveThreshold(double, int);
	bool InsideTimeWindow();
	bool DescantNeutronDiscrimination();
	void FillDetectors();

	void PrintStatistics();

	Settings* fSettings;
	TChain fChain;
	TFile* fFragmentFile;
	TFile* fAnalysisFile;
	TTree fEventTree;
	TFragment* fFragment;
	std::map<uint32_t, TFragment> fFragments;
	bool fWriteFragmentTree;
	TTree fFragmentTree;
	int fFragmentTreeEntries;
	int fRunNumber;
	int fSubRunNumber;
	const TRunInfo* fRunInfo;
	int fKValue;
	TRandom3 fRandom;

	//branches of input tree/chain
	Int_t fEventNumber;
	Int_t fTrackID;
	Int_t fParentID;
	Int_t fStepNumber;
	Int_t fParticleType;
	Int_t fProcessType;
	Int_t fSystemID;
	Int_t fCryNumber;
	Int_t fDetNumber;
	Double_t fDepEnergy;
	Double_t fPosx;
	Double_t fPosy;
	Double_t fPosz;
	Double_t fTime;

	//branches of output tree
	// GRIFFIN
	TGriffin* fGriffin;

	// BGO
	TGriffinBgo* fGriffinBgo;

	// LaBr
	TLaBr* fLaBr;

	// Sceptar
	TSceptar* fSceptar;

	// Descant
	TDescant* fDescant;

	// Paces
	TPaces* fPaces;
};
#endif