Skip to content

Navigation Menu

Sign in
Appearance settings

Search code, repositories, users, issues, pull requests...

Provide feedback

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly
Appearance settings

[FDD] Now propagate avg time for all the BCs in Rec #10545

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

Merged
shahor02 merged 3 commits into from
Jan 23, 2023
Merged

[FDD] Now propagate avg time for all the BCs in Rec #10545

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
45a3781
[FV0] Changed long to long64
Arvindiiti Sep 29, 2022
fcca577
[FDD] Now propagate avg time for all the BCs in Rec (Requested by Evg…
Arvindiiti Jan 10, 2023
c47d4dd
Update Digitizer.cxx
arvindkhuntia Jan 11, 2023
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
20 changes: 13 additions & 7 deletions 20 Detectors/FIT/FDD/reconstruction/src/Reconstructor.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,7 +37,7 @@ void Reconstructor::process(o2::fdd::Digit const& digitBC, gsl::span<const o2::f
bool inTime = inChData[ich].getFlag(ChannelData::EEventDataBit::kIsEventInTVDC);
bool inAdcGate = inChData[ich].getFlag(ChannelData::EEventDataBit::kIsCFDinADCgate);
if (inAdcGate) {
outChData.emplace_back(inChData[ich].mPMNumber, (inChData[ich].mTime) * timePerTDC,
outChData.emplace_back((int)inChData[ich].mPMNumber, (inChData[ich].mTime) * timePerTDC,
(double)inChData[ich].mChargeADC, inChData[ich].mFEEBits);
nStored++;
}
Expand All @@ -49,9 +49,9 @@ void Reconstructor::process(o2::fdd::Digit const& digitBC, gsl::span<const o2::f
}
Float_t timeErr = 1;
if (adc > 3) {
timeErr = 1 / adc;
timeErr = 1. / adc;
}
if (outChData[ich].mPMNumber < 8) {
if ((int)inChData[ich].mPMNumber < 8) {
nInTimeC += inTime;
timeFDC += time / (timeErr * timeErr);
weightFDC += 1. / (timeErr * timeErr);
Expand All @@ -63,18 +63,24 @@ void Reconstructor::process(o2::fdd::Digit const& digitBC, gsl::span<const o2::f
}
const int nsToPs = 1e3;
std::array<int, 2> mCollisionTime = {o2::fdd::RecPoint::sDummyCollissionTime, o2::fdd::RecPoint::sDummyCollissionTime};
/// Avg time for each side, only if one channel satisfy the TVDC condition
/// Avg time for each side, only if one channel satisfy the TVDC condition (if not, also avg time is propagated for background study using AO2D)
if (nInTimeA > 0) {
mCollisionTime[o2::fdd::RecPoint::TimeA] = (weightFDA > 1) ? round(timeFDA / weightFDA * nsToPs)
: o2::fdd::RecPoint::sDummyCollissionTime;
} else {
if (weightFDA > 0) {
mCollisionTime[o2::fdd::RecPoint::TimeA] = round(timeFDA / weightFDA * nsToPs);
}
}
if (nInTimeC > 0) {
mCollisionTime[o2::fdd::RecPoint::TimeC] = (weightFDC > 1) ? round(timeFDC / weightFDC * nsToPs)
: o2::fdd::RecPoint::sDummyCollissionTime;
} else {
if (weightFDC > 0) {
mCollisionTime[o2::fdd::RecPoint::TimeC] = round(timeFDC / weightFDC * nsToPs);
}
}
if (nStored != 0) {
RecPoints.emplace_back(mCollisionTime, firstEntry, nStored, digitBC.getIntRecord(), digitBC.mTriggers);
}
RecPoints.emplace_back(mCollisionTime, firstEntry, nStored, digitBC.getIntRecord(), digitBC.mTriggers);
}
//________________________________________________________
void Reconstructor::finish()
Expand Down
Morty Proxy This is a proxified and sanitized view of the page, visit original site.