|
32 | 32 |
|
33 | 33 | #include "CCDB/BasicCCDBManager.h" |
34 | 34 | #include "CCDB/CcdbApi.h" |
| 35 | +#include "DataFormatsTPC/BetheBlochAleph.h" |
35 | 36 | #include "Framework/ASoA.h" |
36 | 37 | #include "Framework/ASoAHelpers.h" |
37 | 38 | #include "Framework/AnalysisDataModel.h" |
|
44 | 45 | #include "ReconstructionDataFormats/DCA.h" |
45 | 46 | #include "ReconstructionDataFormats/PID.h" |
46 | 47 | #include "ReconstructionDataFormats/Track.h" |
47 | | -#include "DataFormatsTPC/BetheBlochAleph.h" |
48 | 48 |
|
49 | 49 | #include "TGrid.h" |
50 | 50 | #include <TList.h> |
@@ -141,7 +141,7 @@ struct AntinucleiInJets { |
141 | 141 | Configurable<bool> setMCDefaultItsParams{"setMCDefaultItsParams", false, "set MC default parameters"}; |
142 | 142 | Configurable<bool> cfgCompensatePIDinTracking{"cfgCompensatePIDinTracking", false, "If true, divide tpcInnerParam by the electric charge"}; |
143 | 143 | Configurable<std::array<double, 5>> cfgBetheBlochParams{"cfgBetheBlochParams", {0.6539, 1.591, 0.8225, 2.363, 0.09}, "TPC Bethe-Bloch parameterisation for He3"}; |
144 | | - |
| 144 | + |
145 | 145 | // CCDB manager service for accessing condition data |
146 | 146 | Service<o2::ccdb::BasicCCDBManager> ccdb; |
147 | 147 |
|
@@ -249,7 +249,7 @@ struct AntinucleiInJets { |
249 | 249 |
|
250 | 250 | // nsigmaITS for antiproton candidates |
251 | 251 | registryData.add("antiproton_nsigma_its_data", "antiproton_nsigma_its_data", HistType::kTH2F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{ITS}"}}); |
252 | | - |
| 252 | + |
253 | 253 | // custom nsigma for He3 - needed for 24 pp data |
254 | 254 | registryData.add("tpcsignal_data", "Specific energy loss", HistType::kTH2F, {{600, -6., 6., "#it{p} (GeV/#it{c})"}, {1400, 0, 1400, "d#it{E} / d#it{X} (a. u.)"}}); |
255 | 255 | registryData.add("antihelium3_jet_tpc_custom", "antihelium3_jet_tpc_custom", HistType::kTH2F, {{nbins, 3 * min, 3 * max, "#it{p}_{T} (GeV/#it{c})"}, {400, -20.0, 20.0, "n#sigma_{TPC} custom"}}); |
@@ -796,7 +796,6 @@ struct AntinucleiInJets { |
796 | 796 | bool heliumPID = track.pidForTracking() == o2::track::PID::Helium3; |
797 | 797 | float correctedTpcInnerParam = (heliumPID && cfgCompensatePIDinTracking) ? track.tpcInnerParam() / 2 : track.tpcInnerParam(); |
798 | 798 | registryData.fill(HIST("tpcsignal_data"), correctedTpcInnerParam * track.sign(), track.tpcSignal()); |
799 | | - |
800 | 799 | } |
801 | 800 |
|
802 | 801 | // Reject empty events |
@@ -907,7 +906,7 @@ struct AntinucleiInJets { |
907 | 906 | } |
908 | 907 | if (passedItsPidHel) { |
909 | 908 | registryData.fill(HIST("antihelium3_jet_tpc"), 2.0 * pt, nsigmaTPCHe); |
910 | | - //custom nsigma He3 based on bethe bloch fit of TPC signal |
| 909 | + // custom nsigma He3 based on bethe bloch fit of TPC signal |
911 | 910 | double tpcSignal = track.tpcSignal(); |
912 | 911 | double expectedSignalHe3 = tpc::BetheBlochAleph(static_cast<double>(track.tpcInnerParam() * 2. / o2::constants::physics::MassHelium3), cfgBetheBlochParams.value[0], cfgBetheBlochParams.value[1], cfgBetheBlochParams.value[2], cfgBetheBlochParams.value[3], cfgBetheBlochParams.value[4]); |
913 | 912 | double sigmaHe3 = expectedSignalHe3 * cfgBetheBlochParams.value[4]; |
@@ -1003,7 +1002,7 @@ struct AntinucleiInJets { |
1003 | 1002 | } |
1004 | 1003 | if (passedItsPidHel) { |
1005 | 1004 | registryData.fill(HIST("antihelium3_ue_tpc"), 2.0 * pt, nsigmaTPCHe); |
1006 | | - //custom nsigma He3 based on bethe bloch fit of TPC signal |
| 1005 | + // custom nsigma He3 based on bethe bloch fit of TPC signal |
1007 | 1006 | double tpcSignal = track.tpcSignal(); |
1008 | 1007 | double expectedSignalHe3 = tpc::BetheBlochAleph(static_cast<double>(track.tpcInnerParam() * 2. / o2::constants::physics::MassHelium3), cfgBetheBlochParams.value[0], cfgBetheBlochParams.value[1], cfgBetheBlochParams.value[2], cfgBetheBlochParams.value[3], cfgBetheBlochParams.value[4]); |
1009 | 1008 | double sigmaHe3 = expectedSignalHe3 * cfgBetheBlochParams.value[4]; |
|
0 commit comments