forked from AliceO2Group/O2Physics
-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
8c3ba2a
commit ce5e4ae
Showing
4 changed files
with
19 additions
and
19 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -9,8 +9,8 @@ | |
// granted to it by virtue of its status as an Intergovernmental Organization | ||
// or submit itself to any jurisdiction. | ||
|
||
// Task to produce a SV indices table joinable to the jet tables and 3/2-prong SV for hf jet tagging | ||
// | ||
/// \file secondaryVertexReconstruction.cxx | ||
/// \brief Task to produce a SV indices table joinable to the jet tables and 3/2-prong SV for hf jet tagging | ||
/// \author Hadi Hassan <[email protected]> | ||
|
||
#include <type_traits> | ||
|
@@ -272,12 +272,12 @@ struct SecondaryVertexReconstruction { | |
|
||
// fill histograms | ||
if (fillHistograms) { | ||
double DecayLengthNormalised = RecoDecay::distance(std::array{primaryVertex.getX(), primaryVertex.getY(), primaryVertex.getZ()}, std::array{secondaryVertex[0], secondaryVertex[1], secondaryVertex[2]}) / errorDecayLength; | ||
double DecayLengthXYNormalised = RecoDecay::distanceXY(std::array{primaryVertex.getX(), primaryVertex.getY()}, std::array{secondaryVertex[0], secondaryVertex[1]}) / errorDecayLengthXY; | ||
double decayLengthNormalised = RecoDecay::distance(std::array{primaryVertex.getX(), primaryVertex.getY(), primaryVertex.getZ()}, std::array{secondaryVertex[0], secondaryVertex[1], secondaryVertex[2]}) / errorDecayLength; | ||
double decayLengthXYNormalised = RecoDecay::distanceXY(std::array{primaryVertex.getX(), primaryVertex.getY()}, std::array{secondaryVertex[0], secondaryVertex[1]}) / errorDecayLengthXY; | ||
|
||
registry.fill(HIST("hMassNProngs"), massSV, numProngs); | ||
registry.fill(HIST("hLxySNProngs"), DecayLengthXYNormalised, numProngs); | ||
registry.fill(HIST("hLSNProngs"), DecayLengthNormalised, numProngs); | ||
registry.fill(HIST("hLxySNProngs"), decayLengthXYNormalised, numProngs); | ||
registry.fill(HIST("hLSNProngs"), decayLengthNormalised, numProngs); | ||
registry.fill(HIST("hFeNProngs"), energySV / analysisJet.energy() > 1. ? 0.99 : energySV / analysisJet.energy(), numProngs); | ||
} | ||
|
||
|
@@ -306,7 +306,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processData3Prongs(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets, JetTracksData const& tracks, OriginalTracks const& /*tracks*/, aod::BCsWithTimestamps const& /*bcWithTimeStamps*/) | ||
{ | ||
for (auto& jet : jets) { | ||
for (const auto& jet : jets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<3, false>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df3); | ||
sv3prongIndicesTableData(svIndices); | ||
|
@@ -316,7 +316,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processData3ProngsExternalMagneticField(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets, JetTracksData const& tracks, OriginalTracks const& /*tracks*/) | ||
{ | ||
for (auto& jet : jets) { | ||
for (const auto& jet : jets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<3, true>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df3); | ||
sv3prongIndicesTableData(svIndices); | ||
|
@@ -326,7 +326,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processData2Prongs(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets, JetTracksData const& tracks, OriginalTracks const& /*tracks*/, aod::BCsWithTimestamps const& /*bcWithTimeStamps*/) | ||
{ | ||
for (auto& jet : jets) { | ||
for (const auto& jet : jets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<2, false>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df2); | ||
sv2prongIndicesTableData(svIndices); | ||
|
@@ -336,7 +336,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processData2ProngsExternalMagneticField(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedJets, aod::ChargedJetConstituents> const& jets, JetTracksData const& tracks, OriginalTracks const& /*tracks*/) | ||
{ | ||
for (auto& jet : jets) { | ||
for (const auto& jet : jets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<2, true>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df2); | ||
sv2prongIndicesTableData(svIndices); | ||
|
@@ -346,7 +346,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processMCD3Prongs(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents> const& mcdjets, JetTracksMCDwPIs const& tracks, OriginalTracks const& /*tracks*/, aod::BCsWithTimestamps const& /*bcWithTimeStamps*/) | ||
{ | ||
for (auto& jet : mcdjets) { | ||
for (const auto& jet : mcdjets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<3, false>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df3); | ||
sv3prongIndicesTableMCD(svIndices); | ||
|
@@ -356,7 +356,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processMCD3ProngsExternalMagneticField(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents> const& mcdjets, JetTracksMCDwPIs const& tracks, OriginalTracks const& /*tracks*/) | ||
{ | ||
for (auto& jet : mcdjets) { | ||
for (const auto& jet : mcdjets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<3, true>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df3); | ||
sv3prongIndicesTableMCD(svIndices); | ||
|
@@ -366,7 +366,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processMCD2Prongs(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents> const& mcdjets, JetTracksMCDwPIs const& tracks, OriginalTracks const& /*tracks*/, aod::BCsWithTimestamps const& /*bcWithTimeStamps*/) | ||
{ | ||
for (auto& jet : mcdjets) { | ||
for (const auto& jet : mcdjets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<2, false>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df2); | ||
sv2prongIndicesTableMCD(svIndices); | ||
|
@@ -376,7 +376,7 @@ struct SecondaryVertexReconstruction { | |
|
||
void processMCD2ProngsExternalMagneticField(JetCollisionwPIs::iterator const& collision, aod::Collisions const& /*realColl*/, soa::Join<aod::ChargedMCDetectorLevelJets, aod::ChargedMCDetectorLevelJetConstituents> const& mcdjets, JetTracksMCDwPIs const& tracks, OriginalTracks const& /*tracks*/) | ||
{ | ||
for (auto& jet : mcdjets) { | ||
for (const auto& jet : mcdjets) { | ||
std::vector<int> svIndices; | ||
runCreatorNProng<2, true>(collision.template collision_as<aod::Collisions>(), jet, tracks, svIndices, df2); | ||
sv2prongIndicesTableMCD(svIndices); | ||
|
@@ -387,5 +387,5 @@ struct SecondaryVertexReconstruction { | |
|
||
WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) | ||
{ | ||
return WorkflowSpec{adaptAnalysisTask<SecondaryVertexReconstruction>(cfgc, TaskName{"jet-sv-reconstruction-charged"})}; | ||
return WorkflowSpec{adaptAnalysisTask<SecondaryVertexReconstruction>(cfgc, TaskName{"jet-sv-reconstruction-charged"})}; // o2-linter: disable=name/o2-task | ||
} |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters