Skip to content

ITS Calib QC: from rise time to time of arrival #2534

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
Barthelemy merged 1 commit into from
Jun 26, 2025
Merged

ITS Calib QC: from rise time to time of arrival #2534

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
8 changes: 4 additions & 4 deletions Modules/ITS/include/ITS/ITSThresholdCalibrationTask.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -61,8 +61,8 @@ class ITSThresholdCalibrationTask : public TaskInterface
float status;
float Tot; // time over threshold
float TotRms; // time over threshold rms
float Rt; // rise time
float RtRms; // rise time rms
float ToA; // time of arrival
float ToARms; // time of arrival rms
};
struct CalibrationResStructPixel {
int Layer;
Expand Down Expand Up @@ -135,8 +135,8 @@ class ITSThresholdCalibrationTask : public TaskInterface
TH1F *hCalibrationThrNoiseLayer[NLayer], *hCalibrationThrNoiseRMSLayer[NLayer];

// Histograms for pulse-length scan
TH2F *hTimeOverThreshold[3], *hTimeOverThresholdRms[3], *hRiseTime[3], *hRiseTimeRms[3];
TH1F *hTimeOverThresholdLayer[NLayer], *hTimeOverThresholdRmsLayer[NLayer], *hRiseTimeLayer[NLayer], *hRiseTimeRmsLayer[NLayer];
TH2F *hTimeOverThreshold[3], *hTimeOverThresholdRms[3], *hToA[3], *hToARms[3];
TH1F *hTimeOverThresholdLayer[NLayer], *hTimeOverThresholdRmsLayer[NLayer], *hToALayer[NLayer], *hToARmsLayer[NLayer];

o2::itsmft::ChipMappingITS mp;
};
Expand Down
64 changes: 32 additions & 32 deletions Modules/ITS/src/ITSThresholdCalibrationTask.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -60,8 +60,8 @@ ITSThresholdCalibrationTask::~ITSThresholdCalibrationTask()
if (iScan == TOT) {
delete hTimeOverThreshold[iBarrel];
delete hTimeOverThresholdRms[iBarrel];
delete hRiseTime[iBarrel];
delete hRiseTimeRms[iBarrel];
delete hToA[iBarrel];
delete hToARms[iBarrel];
}
if (iScan > 3) {
delete hCalibrationDColChipAverage[iBarrel];
Expand All @@ -74,8 +74,8 @@ ITSThresholdCalibrationTask::~ITSThresholdCalibrationTask()
for (int iLayer = 0; iLayer < NLayer; iLayer++) {
delete hTimeOverThresholdLayer[iLayer];
delete hTimeOverThresholdRmsLayer[iLayer];
delete hRiseTimeLayer[iLayer];
delete hRiseTimeRmsLayer[iLayer];
delete hToALayer[iLayer];
delete hToARmsLayer[iLayer];
}
}
}
Expand Down Expand Up @@ -211,13 +211,13 @@ void ITSThresholdCalibrationTask::doAnalysisTHR(string inString, int iScan)
// fill 2D plots for the pulse length scan
hTimeOverThreshold[iBarrel]->SetBinContent(currentChip, currentStave, result.Tot);
hTimeOverThresholdRms[iBarrel]->SetBinContent(currentChip, currentStave, result.TotRms);
hRiseTime[iBarrel]->SetBinContent(currentChip, currentStave, result.Rt);
hRiseTimeRms[iBarrel]->SetBinContent(currentChip, currentStave, result.RtRms);
hToA[iBarrel]->SetBinContent(currentChip, currentStave, result.ToA);
hToARms[iBarrel]->SetBinContent(currentChip, currentStave, result.ToARms);
// fill 1D plots for the pulse length scan
hTimeOverThresholdLayer[result.Layer]->Fill(result.Tot);
hTimeOverThresholdRmsLayer[result.Layer]->Fill(result.TotRms);
hRiseTimeLayer[result.Layer]->Fill(result.Rt);
hRiseTimeRmsLayer[result.Layer]->Fill(result.RtRms);
hToALayer[result.Layer]->Fill(result.ToA);
hToARmsLayer[result.Layer]->Fill(result.ToARms);
}
}
}
Expand Down Expand Up @@ -339,10 +339,10 @@ ITSThresholdCalibrationTask::CalibrationResStructTHR ITSThresholdCalibrationTask
result.Tot = std::stof(data) / 1000; // to get micro-seconds
} else if (name == "TotRms") {
result.TotRms = std::stof(data) / 1000; // to get micro-seconds
} else if (name == "Rt") {
result.Rt = std::stof(data); // this will be in nano-seconds automatically
} else if (name == "RtRms") {
result.RtRms = std::stof(data); // this will be in nano-seconds automatically
} else if (name == "ToA") {
result.ToA = std::stof(data); // this will be in nano-seconds automatically
} else if (name == "ToARms") {
result.ToARms = std::stof(data); // this will be in nano-seconds automatically
}
}
}
Expand Down Expand Up @@ -391,8 +391,8 @@ void ITSThresholdCalibrationTask::reset()
if (iScan == TOT) {
hTimeOverThreshold[iBarrel]->Reset();
hTimeOverThresholdRms[iBarrel]->Reset();
hRiseTime[iBarrel]->Reset();
hRiseTimeRms[iBarrel]->Reset();
hToA[iBarrel]->Reset();
hToARms[iBarrel]->Reset();
}
if (iScan > 3) {
hCalibrationDColChipAverage[iBarrel]->Reset();
Expand All @@ -404,8 +404,8 @@ void ITSThresholdCalibrationTask::reset()
for (int iLayer = 0; iLayer < NLayer; iLayer++) {
hTimeOverThresholdLayer[iLayer]->Reset();
hTimeOverThresholdRmsLayer[iLayer]->Reset();
hRiseTimeLayer[iLayer]->Reset();
hRiseTimeRmsLayer[iLayer]->Reset();
hToALayer[iLayer]->Reset();
hToARmsLayer[iLayer]->Reset();
}
}
}
Expand Down Expand Up @@ -507,20 +507,20 @@ void ITSThresholdCalibrationTask::createAllHistos()
addObject(hTimeOverThresholdRms[iBarrel]);

// Rise time 2D plot
hRiseTime[iBarrel] = new TH2F(Form("RiseTime%s", sBarrelType[iBarrel].Data()), Form("Rise time for %s (in ns)", sBarrelType[iBarrel].Data()), nChips[iBarrel], -0.5, nChips[iBarrel] - 0.5, nStaves[iBarrel], -0.5, nStaves[iBarrel] - 0.5);
hRiseTime[iBarrel]->SetStats(0);
hToA[iBarrel] = new TH2F(Form("ToA%s", sBarrelType[iBarrel].Data()), Form("ToA for %s (in ns)", sBarrelType[iBarrel].Data()), nChips[iBarrel], -0.5, nChips[iBarrel] - 0.5, nStaves[iBarrel], -0.5, nStaves[iBarrel] - 0.5);
hToA[iBarrel]->SetStats(0);
if (iBarrel != 0)
formatAxes(hRiseTime[iBarrel], "Chip", "", 1, 1.10);
formatLayers(hRiseTime[iBarrel], iBarrel);
addObject(hRiseTime[iBarrel]);
formatAxes(hToA[iBarrel], "Chip", "", 1, 1.10);
formatLayers(hToA[iBarrel], iBarrel);
addObject(hToA[iBarrel]);

// Rise time RMS 2D plot
hRiseTimeRms[iBarrel] = new TH2F(Form("RiseTimeRms%s", sBarrelType[iBarrel].Data()), Form("Rise time RMS for %s (in ns)", sBarrelType[iBarrel].Data()), nChips[iBarrel], -0.5, nChips[iBarrel] - 0.5, nStaves[iBarrel], -0.5, nStaves[iBarrel] - 0.5);
hRiseTimeRms[iBarrel]->SetStats(0);
hToARms[iBarrel] = new TH2F(Form("ToARms%s", sBarrelType[iBarrel].Data()), Form("ToA RMS for %s (in ns)", sBarrelType[iBarrel].Data()), nChips[iBarrel], -0.5, nChips[iBarrel] - 0.5, nStaves[iBarrel], -0.5, nStaves[iBarrel] - 0.5);
hToARms[iBarrel]->SetStats(0);
if (iBarrel != 0)
formatAxes(hRiseTimeRms[iBarrel], "Chip", "", 1, 1.10);
formatLayers(hRiseTimeRms[iBarrel], iBarrel);
addObject(hRiseTimeRms[iBarrel]);
formatAxes(hToARms[iBarrel], "Chip", "", 1, 1.10);
formatLayers(hToARms[iBarrel], iBarrel);
addObject(hToARms[iBarrel]);
}
}

Expand Down Expand Up @@ -549,14 +549,14 @@ void ITSThresholdCalibrationTask::createAllHistos()
addObject(hTimeOverThresholdRmsLayer[iLayer]);

// Rise time for each layer
hRiseTimeLayer[iLayer] = new TH1F(Form("RiseTimeLayer%d", iLayer), Form("Rise time distribution for Layer%d", iLayer), 1600, 0, 800);
formatAxes(hRiseTimeLayer[iLayer], "Rise time (ns) ", "Counts", 1, 1.10);
addObject(hRiseTimeLayer[iLayer]);
hToALayer[iLayer] = new TH1F(Form("ToALayer%d", iLayer), Form("ToA distribution for Layer%d", iLayer), 1600, 0, 800);
formatAxes(hToALayer[iLayer], "ToA (ns) ", "Counts", 1, 1.10);
addObject(hToALayer[iLayer]);

// Rise time RMS for each layer
hRiseTimeRmsLayer[iLayer] = new TH1F(Form("RiseTimeRmsLayer%d", iLayer), Form("Rise time Rms distribution for Layer%d", iLayer), 1000, 0, 500);
formatAxes(hRiseTimeRmsLayer[iLayer], "Rise time Rms (ns) ", "Counts", 1, 1.10);
addObject(hRiseTimeRmsLayer[iLayer]);
hToARmsLayer[iLayer] = new TH1F(Form("ToARmsLayer%d", iLayer), Form("ToA Rms distribution for Layer%d", iLayer), 1000, 0, 500);
formatAxes(hToARmsLayer[iLayer], "ToA Rms (ns) ", "Counts", 1, 1.10);
addObject(hToARmsLayer[iLayer]);
}
}

Expand Down