Skip to content

Commit 41109c7

Browse files
authored
PWGJE : Add amplitude cut for cell time histograms in EMCAL task [EMCAL-733] (AliceO2Group#1614)
Exclude low-E cells which have a poor time resolution from time distribution histograms. Amplitude cut is settable.
1 parent 7f88118 commit 41109c7

File tree

1 file changed

+7
-4
lines changed

1 file changed

+7
-4
lines changed

PWGJE/Tasks/emccellmonitor.cxx

+7-4
Original file line numberDiff line numberDiff line change
@@ -49,6 +49,7 @@ struct CellMonitor {
4949
o2::framework::Configurable<double> mMinCellAmplitude{"minCellAmplitude", 0., "Minimum cell amplitude for histograms."};
5050
o2::framework::Configurable<double> mMinCellTimeMain{"minCellTimeMain", -50, "Min. cell time of main bunch selection"};
5151
o2::framework::Configurable<double> mMaxCellTimeMain{"maxCellTimeMain", 100, "Max. cell time of main bunch selection"};
52+
o2::framework::Configurable<double> mMinCellAmplitudeTimeHists{"minCellAmplitudeTimeHists", 0, "Min. cell amplitude used for time distribution"};
5253
o2::framework::Configurable<std::string> mVetoBCID{"vetoBCID", "", "BC ID to be excluded"};
5354
o2::framework::Configurable<std::string> mSelectBCID{"selectBCID", "all", "BC ID to be included"};
5455

@@ -157,13 +158,15 @@ struct CellMonitor {
157158
continue;
158159
mHistManager.fill(HIST("cellAmplitudeCut"), cell.amplitude(), cell.cellNumber());
159160
mHistManager.fill(HIST("cellFrequency"), cell.cellNumber());
160-
mHistManager.fill(HIST("cellTime"), cell.time(), cell.cellNumber());
161-
if (cell.time() > mMinCellTimeMain && cell.time() < mMaxCellTimeMain) {
162-
mHistManager.fill(HIST("cellTimeMain"), cell.time(), cell.cellNumber());
161+
if (cell.amplitude() >= mMinCellAmplitudeTimeHists) {
162+
mHistManager.fill(HIST("cellTime"), cell.time(), cell.cellNumber());
163+
if (cell.time() > mMinCellTimeMain && cell.time() < mMaxCellTimeMain) {
164+
mHistManager.fill(HIST("cellTimeMain"), cell.time(), cell.cellNumber());
165+
}
166+
mHistManager.fill(HIST("celTimeBC"), cellIR.bc, cell.time());
163167
}
164168

165169
mHistManager.fill(HIST("cellAmplitudeBC"), cellIR.bc, cell.amplitude());
166-
mHistManager.fill(HIST("celTimeBC"), cellIR.bc, cell.time());
167170

168171
// Get Cell index in eta-phi of sm
169172
auto [supermodule, module, phiInModule, etaInModule] = mGeometry->GetCellIndex(cell.cellNumber());

0 commit comments

Comments
 (0)