-
Notifications
You must be signed in to change notification settings - Fork 494
Expand file tree
/
Copy pathSAMPAProcessing.cxx
More file actions
60 lines (53 loc) · 2.09 KB
/
SAMPAProcessing.cxx
File metadata and controls
60 lines (53 loc) · 2.09 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
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// 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.
/// \file SAMPAProcessing.cxx
/// \brief Implementation of the SAMPA response
/// \author Andi Mathis, TU München, andreas.mathis@ph.tum.de
#include "TPCSimulation/SAMPAProcessing.h"
#include "TPCBaseRecSim/CDBInterface.h"
#include <fstream>
#include <iostream>
#include <sstream>
#include <string>
#include "Framework/Logger.h"
using namespace o2::tpc;
SAMPAProcessing::SAMPAProcessing() : mRandomNoiseRing()
{
updateParameters();
}
void SAMPAProcessing::updateParameters(float vdrift)
{
mGasParam = &(ParameterGas::Instance());
mDetParam = &(ParameterDetector::Instance());
mEleParam = &(ParameterElectronics::Instance());
auto& cdb = CDBInterface::instance();
mPedestalMap = &(cdb.getPedestals());
mPedestalMapCRU = &(cdb.getPedestalsCRU());
mNoiseMap = &(cdb.getNoise());
mZeroSuppression = &(cdb.getZeroSuppressionThreshold());
mVDrift = vdrift > 0 ? vdrift : mGasParam->DriftV;
}
void SAMPAProcessing::getShapedSignal(float ADCsignal, float driftTime, std::vector<float>& signalArray) const
{
const float timeBinTime = getTimeBinTime(driftTime);
const float offset = driftTime - timeBinTime;
for (float bin = 0; bin < mEleParam->NShapedPoints; bin += Vc::float_v::Size) {
Vc::float_v binvector;
for (int i = 0; i < Vc::float_v::Size; ++i) {
binvector[i] = bin + i;
}
Vc::float_v time = timeBinTime + binvector * mEleParam->ZbinWidth;
Vc::float_v signal = getGamma4(time, Vc::float_v(timeBinTime + offset), Vc::float_v(ADCsignal));
for (int i = 0; i < Vc::float_v::Size; ++i) {
signalArray[bin + i] = signal[i];
}
}
}