CbmRoot
Loading...
Searching...
No Matches
CbmAnaConversionCutSettings.h
Go to the documentation of this file.
1/* Copyright (C) 2015-2016 Fakultaet fuer Mathematik und Naturwissenschaften, Bergische Universitaet Wuppertal, Wuppertal
2 SPDX-License-Identifier: GPL-3.0-only
3 Authors: Sascha Reinecke [committer], Florian Uhlig */
4
11#ifndef CBM_ANA_CONVERSION_CUT_SETTINGS
12#define CBM_ANA_CONVERSION_CUT_SETTINGS
13
14#include "CbmMCTrack.h"
15
16#include "TLorentzVector.h"
17#include "TMath.h"
18
19#define M2E 2.6112004954086e-7
20
22public:
23 /*
24 * Calculate cut for opening angle between electrons
25 */
26 static Double_t CalcOpeningAngleCut(Double_t pt)
27 {
28 Double_t cut = 0;
29 //cut = 1.5 - 0.5 * pt;
30 cut = 1.8 - 0.6 * pt;
31 return cut;
32 }
33
34
35 static Double_t CalcOpeningAngleCutAlt1(Double_t pt)
36 {
37 Double_t cut = 0;
38 cut = 1.5 - 0.5 * pt;
39 return cut;
40 }
41
42
43 static Double_t CalcOpeningAngleCutAlt2(Double_t /*pt*/)
44 {
45 Double_t cut = 0;
46 cut = 1.8;
47 return cut;
48 }
49
50
51 static Double_t CalcOpeningAngleCutAlt3(Double_t /*pt*/)
52 {
53 Double_t cut = 0;
54 cut = 1.5;
55 return cut;
56 }
57
58
59 static Double_t CalcOpeningAngleCutAlt4(Double_t pt)
60 {
61 Double_t cut = 0;
62 if (pt <= 1) cut = 1.5 - 0.5 * pt;
63 if (pt > 1) cut = 1.0;
64 return cut;
65 }
66
67
68 static Double_t CalcOpeningAngleCutAlt5(Double_t pt)
69 {
70 Double_t cut = 0;
71 if (pt <= 1) cut = 1.0 + 1.0 * pt;
72 if (pt > 1) cut = 2.0;
73 return cut;
74 }
75
76
77 /*
78 * Calculate cut for chi for reconstructed momenta
79 */
80 static Double_t CalcChiCut(Double_t pt)
81 {
82 Double_t cut = 0;
83 if (pt < 0.4) { cut = 31. - 70. * pt; }
84 if (pt >= 0.4) { cut = 3; }
85 return cut;
86 }
87
88
89 static Double_t CalcChiCutAlt1(Double_t /*pt*/)
90 {
91 Double_t cut = 0;
92 cut = 3;
93 return cut;
94 }
95
96
97 static Double_t CalcChiCutAlt2(Double_t pt)
98 {
99 Double_t cut = 0;
100 if (pt < 0.4) { cut = 15. - 30. * pt; }
101 if (pt >= 0.4) { cut = 3; }
102 return cut;
103 }
104
105
106 static Double_t CalcChiCutAlt3(Double_t /*pt*/)
107 {
108 Double_t cut = 0;
109 cut = 4;
110 return cut;
111 }
112};
113
114#endif
static Double_t CalcOpeningAngleCutAlt5(Double_t pt)
static Double_t CalcOpeningAngleCutAlt3(Double_t)
static Double_t CalcChiCutAlt3(Double_t)
static Double_t CalcChiCutAlt2(Double_t pt)
static Double_t CalcOpeningAngleCutAlt1(Double_t pt)
static Double_t CalcOpeningAngleCutAlt2(Double_t)
static Double_t CalcOpeningAngleCutAlt4(Double_t pt)
static Double_t CalcOpeningAngleCut(Double_t pt)
static Double_t CalcChiCutAlt1(Double_t)
static Double_t CalcChiCut(Double_t pt)