Skip to content

Commit 5c3b91b

Browse files
Add files via upload
Added option to reject the inside or the outside of the detector ring of FT0
1 parent e63b492 commit 5c3b91b

File tree

1 file changed

+12
-2
lines changed

1 file changed

+12
-2
lines changed

PWGCF/TwoParticleCorrelations/Tasks/longRangeDihadronCor.cxx

Lines changed: 12 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -238,6 +238,16 @@ struct LongRangeDihadronCor {
238238
kTOF,
239239
kITS
240240
};
241+
enum DetectorChannels {
242+
kFT0AInnerRingMin = 0,
243+
kFT0AInnerRingMax = 31,
244+
kFT0AOuterRingMin = 32,
245+
kFT0AOuterRingMax = 95,
246+
kFT0CInnerRingMin = 96,
247+
kFT0CInnerRingMax = 143,
248+
kFT0COuterRingMin = 144,
249+
kFT0COuterRingMax = 207
250+
};
241251
std::array<float, 6> tofNsigmaCut;
242252
std::array<float, 6> itsNsigmaCut;
243253
std::array<float, 6> tpcNsigmaCut;
@@ -641,15 +651,15 @@ struct LongRangeDihadronCor {
641651
id = ft0.channelC()[iCh];
642652
id = id + Ft0IndexA;
643653
ampl = ft0.amplitudeC()[iCh];
644-
if ((cfgRejectInsideDetectors && ((id >= 31 && id <= 96) || (id >= 143 && id <= 208))) || (cfgRejectOutsideDetectors && ((id >= 0 && id <= 30) || (id >= 97 && id <= 142))))
654+
if ((cfgRejectInsideDetectors && ((id >= kFT0AInnerRingMin && id <= kFT0AInnerRingMax) || (id >= kFT0CInnerRingMin && id <= kFT0CInnerRingMax))) || (cfgRejectOutsideDetectors && ((id >= kFT0AOuterRingMin && id <= kFT0AOuterRingMax) || (id >= kFT0COuterRingMin && id <= kFT0COuterRingMax))))
645655
ampl = 0.;
646656
registry.fill(HIST("FT0Amp"), id, ampl);
647657
ampl = ampl / cstFT0RelGain[iCh];
648658
registry.fill(HIST("FT0AmpCorrect"), id, ampl);
649659
} else if (fitType == kFT0A) {
650660
id = ft0.channelA()[iCh];
651661
ampl = ft0.amplitudeA()[iCh];
652-
if ((cfgRejectInsideDetectors && ((id >= 31 && id <= 96) || (id >= 143 && id <= 208))) || (cfgRejectOutsideDetectors && ((id >= 0 && id <= 30) || (id >= 97 && id <= 142))))
662+
if ((cfgRejectInsideDetectors && ((id >= kFT0AInnerRingMin && id <= kFT0AInnerRingMax) || (id >= kFT0CInnerRingMin && id <= kFT0CInnerRingMax))) || (cfgRejectOutsideDetectors && ((id >= kFT0AOuterRingMin && id <= kFT0AOuterRingMax) || (id >= kFT0COuterRingMin && id <= kFT0COuterRingMax))))
653663
ampl = 0.;
654664
registry.fill(HIST("FT0Amp"), id, ampl);
655665
ampl = ampl / cstFT0RelGain[iCh];

0 commit comments

Comments
 (0)