Skip to content

FT0: align trigger definition in Digitizer with data trigger logic #15156

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
spulawsk wants to merge 1 commit into
base: dev
Choose a base branch
from
+167 −6
Open

FT0: align trigger definition in Digitizer with data trigger logic #15156

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 3 additions & 2 deletions Detectors/FIT/FT0/base/include/FT0Base/FT0DigParam.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -31,8 +31,8 @@ struct FT0DigParam : o2::conf::ConfigurableParamHelper<FT0DigParam> {
float mAmpRecordUp = 15; // to [ns]
float hitTimeOffsetA = 0; ///< hit time offset on the A side [ns]
float hitTimeOffsetC = 0; ///< hit time offset on the C side [ns]
int mtrg_central_trh = 600.; // channels
int mtrg_semicentral_trh = 300.; // channels
int mtrg_central_trh = 1433.; // Tclu
int mtrg_semicentral_trh = 35.; // Tclu units = TCM charge level unit = 16 ADCunits

float mMip_in_V = 7; // MIP to mV
float mPe_in_mip = 0.004; // invserse Np.e. in MIP 1./250.
Expand All @@ -43,6 +43,7 @@ struct FT0DigParam : o2::conf::ConfigurableParamHelper<FT0DigParam> {
float mNoiseVar = 0.1; // noise level
float mNoisePeriod = 1 / 0.9; // GHz low frequency noise period;
short mTime_trg_gate = 153; // #channels as in TCM as in Pilot beams ('OR gate' setting in TCM tab in ControlServer)
short mTime_trg_vertex_gate = 100; // #channels as in TCM as in Pilot beams ('OR gate' setting in TCM tab in ControlServer)
float mAmpThresholdForReco = 5; // only channels with amplitude higher will participate in calibration and collision time: 0.3 MIP
short mTimeThresholdForReco = 1000; // only channels with time below will participate in calibration and collision time

Expand Down
168 changes: 164 additions & 4 deletions Detectors/FIT/FT0/simulation/src/Digitizer.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,13 @@
#include "CommonConstants/PhysicsConstants.h"
#include "CommonDataFormat/InteractionRecord.h"

#include "DataFormatsFT0/LookUpTable.h"
#include "FT0Base/Constants.h"
#include <map>
#include <array>
#include <regex>
#include <string>

#include "TMath.h"
#include "TRandom.h"
#include <TH1F.h>
Expand All @@ -34,6 +41,66 @@
// signal shape function
template <typename Float>
Float signalForm_i(Float x)
{
float p0, p1, p2, p3, p4, p5, p6, p7;
p0 = 1.30853;

Check failure on line 46 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Trailing spaces 

Remove the trailing spaces at the end of the line.
p1 = 0.516807;
p2 = 3.36714;
p3 = -1.01206;
p4 = 1.42832;
p5 = 1.1589;
p6 = 1.22019;
p7 = 0.426818;

Double_t val = 0;

if (x > p3) {
Double_t x1 = x - p3;
Double_t arg1 = (log(x1) - p0) / p1;
val += p2 * (1.0 / (x1 * p1 * sqrt(2 * TMath::Pi()))) * exp(-0.5 * arg1 * arg1);
}

if (x > p7) {
Double_t x2 = x - p7;
Double_t arg2 = (log(x2) - p4) / p5;
val += p6 * (1.0 / (x2 * p5 * sqrt(2 * TMath::Pi()))) * exp(-0.5 * arg2 * arg2);
}

return val;
};

// integrated signal shape function
inline float signalForm_integral(float x)
{
float p0, p1, p2, p3, p4, p5, p6, p7;
p0 = 1.30853;

Check failure on line 76 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Trailing spaces 

Remove the trailing spaces at the end of the line.
p1 = 0.516807;
p2 = 3.36714;
p3 = -1.01206;
p4 = 1.42832;
p5 = 1.1589;
p6 = 1.22019;
p7 = 0.426818;
Double_t val = 0;

if (x > p3) {
Double_t x1 = x - p3;
Double_t z1 = (log(x1) - p0) / (sqrt(2) * p1);
val += p2 * 0.5 * (1 + TMath::Erf(z1)); // norm1 * CDF1
}

if (x > p7) {
Double_t x2 = x - p7;
Double_t z2 = (log(x2) - p4) / (sqrt(2) * p5);
val += p6 * 0.5 * (1 + TMath::Erf(z2)); // norm2 * CDF2
}

return val;
};
/*

Check failure on line 100 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Tab characters found 

Indent code using spaces instead of tabs.
// signal shape function
template <typename Float>
Float signalForm_i(Float x)
{
using namespace std;
Float const a = -0.45458;
Expand All @@ -52,7 +119,7 @@
}
return -(exp(b * x) / b - exp(a * x) / a) / 7.8446501;
};

*/
// SIMD version of the integrated signal shape function
inline Vc::float_v signalForm_integralVc(Vc::float_v x)
{
Expand Down Expand Up @@ -249,8 +316,64 @@
if (bc.hits.empty()) {
return;
}
// Initialize mapping channelID -> PM hash and PM side (A/C) using FT0 LUT
static bool pmLutInitialized = false;
static std::array<uint8_t, o2::ft0::Constants::sNCHANNELS_PM> mChID2PMhash{};
static std::map<uint8_t, bool> mMapPMhash2isAside; // hashed PM -> is A side

if (!pmLutInitialized) {
std::map<std::string, uint8_t> mapFEE2hash; // module name -> hashed PM id
uint8_t tcmHash = 0;

const auto& lut = o2::ft0::SingleLUT::Instance().getVecMetadataFEE();
auto lutSorted = lut;
std::sort(lutSorted.begin(), lutSorted.end(),
[](const auto& first, const auto& second) { return first.mModuleName < second.mModuleName; });

uint8_t binPos = 0;
for (const auto& lutEntry : lutSorted) {
const auto& moduleName = lutEntry.mModuleName;
const auto& moduleType = lutEntry.mModuleType;
const auto& strChID = lutEntry.mChannelID;

auto [it, inserted] = mapFEE2hash.insert({moduleName, binPos});
if (inserted) {
if (moduleName.find("PMA") != std::string::npos) {
mMapPMhash2isAside.insert({binPos, true});
} else if (moduleName.find("PMC") != std::string::npos) {
mMapPMhash2isAside.insert({binPos, false});
}
++binPos;
}

if (std::regex_match(strChID, std::regex("^[0-9]{1,3}$"))) {
int chID = std::stoi(strChID);
if (chID < o2::ft0::Constants::sNCHANNELS_PM) {
mChID2PMhash[chID] = mapFEE2hash[moduleName];
} else {
LOG(fatal) << "Incorrect LUT entry: chID " << strChID << " | " << moduleName;
}
} else if (moduleType != "TCM") {
LOG(fatal) << "Non-TCM module w/o numerical chID: chID " << strChID << " | " << moduleName;
} else { // TCM
tcmHash = mapFEE2hash[moduleName];
}
}

pmLutInitialized = true;
}

int n_hit_A = 0, n_hit_C = 0, mean_time_A = 0, mean_time_C = 0;
int summ_ampl_A = 0, summ_ampl_C = 0;
int sum_A_ampl = 0, sum_C_ampl = 0;
int nPMTs = mGeometry.NCellsA * 4 + mGeometry.NCellsC * 4;
std::vector<int> sum_ampl_ipmt(nPMTs, 0);
// Per-PM summed charge (like in digits2trgFT0)
std::map<uint8_t, int> mapPMhash2sumAmpl;
for (const auto& entry : mMapPMhash2isAside) {
mapPMhash2sumAmpl.insert({entry.first, 0});
}

int vertex_time;
const auto& params = FT0DigParam::Instance();
int first = digitsCh.size(), nStored = 0;
Expand Down Expand Up @@ -297,6 +420,10 @@
if (is_time_in_signal_gate) {
chain |= (1 << o2::ft0::ChannelData::EEventDataBit::kIsCFDinADCgate);
chain |= (1 << o2::ft0::ChannelData::EEventDataBit::kIsEventInTVDC);
// Sum channel charge per PM (similar logic as in digits2trgFT0)
if (ipmt < o2::ft0::Constants::sNCHANNELS_PM) {
mapPMhash2sumAmpl[mChID2PMhash[static_cast<uint8_t>(ipmt)]] += static_cast<int>(amp);
}
}
digitsCh.emplace_back(ipmt, smeared_time, int(amp), chain);
nStored++;
Expand All @@ -308,6 +435,8 @@
continue;
}

sum_ampl_ipmt[ipmt] += amp;

if (is_A_side) {
n_hit_A++;
summ_ampl_A += amp;
Expand All @@ -318,17 +447,48 @@
mean_time_C += smeared_time;
}
}

Check failure on line 450 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Trailing spaces 

Remove the trailing spaces at the end of the line.
for (size_t i = 0; i < sum_ampl_ipmt.size(); i++) {
sum_ampl_ipmt[i] = sum_ampl_ipmt[i] >> 3;

Check failure on line 452 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Tab characters found 

Indent code using spaces instead of tabs.

Check failure on line 452 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Trailing spaces 

Remove the trailing spaces at the end of the line.
if (i < 4 * mGeometry.NCellsA) {
sum_A_ampl += sum_ampl_ipmt[i];

Check failure on line 454 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Tab characters found 

Indent code using spaces instead of tabs.

Check failure on line 454 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Trailing spaces 

Remove the trailing spaces at the end of the line.
} else {

Check failure on line 455 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Tab characters found 

Indent code using spaces instead of tabs.
sum_C_ampl += sum_ampl_ipmt[i];

Check failure on line 456 in Detectors/FIT/FT0/simulation/src/Digitizer.cxx

View workflow job for this annotation

GitHub Actions / PR formatting / whitespace

Trailing spaces 

Remove the trailing spaces at the end of the line.
}
}

// Sum over PMs (using per-PM map) for debug/monitoring
int sum_PM_ampl_debug = 0;
int sum_PM_ampl_A_debug = 0;
int sum_PM_ampl_C_debug = 0;
for (const auto& entry : mapPMhash2sumAmpl) {
int pmAmpl = (entry.second >> 3);
sum_PM_ampl_debug += pmAmpl;
auto itSide = mMapPMhash2isAside.find(entry.first);
if (itSide != mMapPMhash2isAside.end()) {
if (itSide->second) {
sum_PM_ampl_A_debug += pmAmpl;
} else {
sum_PM_ampl_C_debug += pmAmpl;
}
}
}
LOG(debug) << "Sum PM amplitude (LUT-based): total=" << sum_PM_ampl_debug
<< " A-side=" << sum_PM_ampl_A_debug
<< " C-side=" << sum_PM_ampl_C_debug;


Bool_t is_A, is_C, isVertex, is_Central, is_SemiCentral = 0;
is_A = n_hit_A > 0;
is_C = n_hit_C > 0;
is_Central = summ_ampl_A + summ_ampl_C >= params.mtrg_central_trh;
is_SemiCentral = summ_ampl_A + summ_ampl_C >= params.mtrg_semicentral_trh;
is_Central = sum_PM_ampl_A_debug + sum_PM_ampl_C_debug >= 2*params.mtrg_central_trh;
is_SemiCentral = sum_PM_ampl_A_debug + sum_PM_ampl_C_debug >= 2*params.mtrg_semicentral_trh && !is_Central;
uint32_t amplA = is_A ? summ_ampl_A * 0.125 : -5000; // sum amplitude A side / 8 (hardware)
uint32_t amplC = is_C ? summ_ampl_C * 0.125 : -5000; // sum amplitude C side / 8 (hardware)
int timeA = is_A ? mean_time_A / n_hit_A : -5000; // average time A side
int timeC = is_C ? mean_time_C / n_hit_C : -5000; // average time C side
vertex_time = (timeC - timeA) * 0.5;
isVertex = is_A && is_C && (vertex_time > -params.mTime_trg_gate && vertex_time < params.mTime_trg_gate);
isVertex = is_A && is_C && (vertex_time > -params.mTime_trg_vertex_gate && vertex_time < params.mTime_trg_vertex_gate);
LOG(debug) << " A " << is_A << " timeA " << timeA << " mean_time_A " << mean_time_A << " n_hit_A " << n_hit_A << " C " << is_C << " timeC " << timeC << " mean_time_C " << mean_time_C << " n_hit_C " << n_hit_C << " vertex_time " << vertex_time;
Triggers triggers;
bool isLaser = false;
Expand Down
Loading