forked from mcoquet642/AliceO2
-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathdynamicColumns.cxx
More file actions
65 lines (58 loc) · 2.17 KB
/
dynamicColumns.cxx
File metadata and controls
65 lines (58 loc) · 2.17 KB
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
// Copyright CERN and copyright holders of ALICE O2. This software is
// distributed under the terms of the GNU General Public License v3 (GPL
// Version 3), copied verbatim in the file "COPYING".
//
// See http://alice-o2.web.cern.ch/license for full licensing information.
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
#include "Framework/runDataProcessing.h"
#include "Framework/AnalysisTask.h"
#include "Framework/AnalysisDataModel.h"
namespace o2::aod
{
namespace etaphi
{
DECLARE_SOA_COLUMN(Tgl, tgl, float, "fTgl");
DECLARE_SOA_COLUMN(Snp, snp, float, "fSnp");
DECLARE_SOA_COLUMN(Alpha, alpha, float, "fAlpha");
DECLARE_SOA_DYNAMIC_COLUMN(Eta, eta, [](float tgl) { return log(tan(0.25 * M_PI - 0.5 * atan(tgl))); });
DECLARE_SOA_DYNAMIC_COLUMN(Phi, phi, [](float snp, float alpha) { return asin(snp) + alpha + M_PI; });
} // namespace etaphi
DECLARE_SOA_TABLE(EtaPhis, "RN2", "ETAPHI",
etaphi::Tgl, etaphi::Snp, etaphi::Alpha,
etaphi::Eta<etaphi::Tgl>,
etaphi::Phi<etaphi::Snp, etaphi::Alpha>);
} // namespace o2::aod
using namespace o2;
using namespace o2::framework;
// This is a very simple example showing how to iterate over tracks
// and create a new collection for them.
// FIXME: this should really inherit from AnalysisTask but
// we need GCC 7.4+ for that
struct ATask {
Produces<aod::EtaPhis> etaphi;
void process(aod::Tracks const& tracks)
{
for (auto& track : tracks) {
etaphi(track.tgl(), track.snp(), track.alpha());
}
}
};
struct BTask {
void process(aod::EtaPhis const& etaPhis)
{
for (auto& etaPhi : etaPhis) {
auto phi = asin(etaPhi.snp()) + etaPhi.alpha() + M_PI;
auto eta = log(tan(0.25 * M_PI - 0.5 * atan(etaPhi.tgl())));
LOGF(ERROR, "(%f, %f, %f, %f)", etaPhi.eta(), etaPhi.phi(), eta - etaPhi.eta(), phi - etaPhi.phi());
}
}
};
WorkflowSpec defineDataProcessing(ConfigContext const&)
{
return WorkflowSpec{
adaptAnalysisTask<ATask>("produce-track-copy"),
adaptAnalysisTask<BTask>("check-eta-phi")};
}