11 #include <tracking/ckf/pxd/findlets/CKFToPXDFindlet.h>
13 #include <tracking/ckf/general/findlets/StateCreator.icc.h>
14 #include <tracking/ckf/general/findlets/CKFRelationCreator.icc.h>
15 #include <tracking/ckf/general/findlets/StateCreatorWithReversal.icc.h>
16 #include <tracking/ckf/general/findlets/TreeSearcher.icc.h>
17 #include <tracking/ckf/general/findlets/OverlapResolver.icc.h>
18 #include <tracking/ckf/general/findlets/SpacePointTagger.icc.h>
19 #include <tracking/ckf/general/findlets/ResultStorer.icc.h>
21 #include <tracking/ckf/pxd/entities/CKFToPXDResult.h>
22 #include <tracking/ckf/pxd/entities/CKFToPXDState.h>
24 #include <tracking/dataobjects/RecoTrack.h>
25 #include <tracking/spacePointCreation/SpacePoint.h>
26 #include <pxd/dataobjects/PXDCluster.h>
28 #include <framework/core/ModuleParamList.h>
31 using namespace TrackFindingCDC;
37 addProcessingSignalListener(&m_dataHandler);
38 addProcessingSignalListener(&m_hitsLoader);
39 addProcessingSignalListener(&m_stateCreatorFromTracks);
40 addProcessingSignalListener(&m_stateCreatorFromHits);
41 addProcessingSignalListener(&m_relationCreator);
42 addProcessingSignalListener(&m_treeSearchFindlet);
43 addProcessingSignalListener(&m_overlapResolver);
44 addProcessingSignalListener(&m_spacePointTagger);
45 addProcessingSignalListener(&m_resultStorer);
50 Super::exposeParameters(moduleParamList, prefix);
52 m_dataHandler.exposeParameters(moduleParamList, prefix);
53 m_hitsLoader.exposeParameters(moduleParamList, prefix);
54 m_stateCreatorFromTracks.exposeParameters(moduleParamList, prefix);
55 m_stateCreatorFromHits.exposeParameters(moduleParamList, prefix);
56 m_relationCreator.exposeParameters(moduleParamList, prefix);
57 m_treeSearchFindlet.exposeParameters(moduleParamList, prefix);
58 m_overlapResolver.exposeParameters(moduleParamList, prefix);
59 m_spacePointTagger.exposeParameters(moduleParamList, prefix);
60 m_resultStorer.exposeParameters(moduleParamList, prefix);
62 moduleParamList->
addParameter(
"minimalHitRequirement", m_param_minimalHitRequirement,
63 "Minimal Hit requirement for the results (counted in space points)",
64 m_param_minimalHitRequirement);
65 moduleParamList->
addParameter(
"onlyUseTracksWithSVD", m_param_onlyUseTracksWithSVD,
66 "Only use tracks which have an SVD hit associated.",
67 m_param_onlyUseTracksWithSVD);
68 moduleParamList->
addParameter(
"reverseSeedState", m_param_reverseSeed,
"Reverse the seed.", m_param_reverseSeed);
71 moduleParamList->
getParameter<std::string>(
"advanceHighFilter").setDefaultValue(
"advance");
72 moduleParamList->
getParameter<std::string>(
"updateHighFilter").setDefaultValue(
"fit");
74 moduleParamList->
getParameter<std::string>(
"firstHighFilter").setDefaultValue(
"mva_with_direction_check");
75 moduleParamList->
getParameter<std::string>(
"secondHighFilter").setDefaultValue(
"mva");
76 moduleParamList->
getParameter<std::string>(
"thirdHighFilter").setDefaultValue(
"mva");
78 moduleParamList->
getParameter<
bool>(
"useAssignedHits").setDefaultValue(
false);
80 moduleParamList->
getParameter<std::string>(
"hitFilter").setDefaultValue(
"sensor");
81 moduleParamList->
getParameter<std::string>(
"seedFilter").setDefaultValue(
"sensor");
82 moduleParamList->
getParameter<std::string>(
"preSeedFilter").setDefaultValue(
"loose");
83 moduleParamList->
getParameter<std::string>(
"preHitFilter").setDefaultValue(
"loose");
85 moduleParamList->
getParameter<std::string>(
"hitsSpacePointsStoreArrayName").setDefaultValue(
"PXDSpacePoints");
87 moduleParamList->
getParameter<std::string>(
"filter").setDefaultValue(
"mva");
94 m_recoTracksVector.clear();
95 m_spacePointVector.clear();
102 m_filteredResults.clear();
107 m_dataHandler.apply(m_recoTracksVector);
108 m_hitsLoader.apply(m_spacePointVector);
110 if (m_spacePointVector.empty() or m_recoTracksVector.empty()) {
115 const auto notFromPXD = [](
const SpacePoint * spacePoint) {
118 TrackFindingCDC::erase_remove_if(m_spacePointVector, notFromPXD);
120 if (m_param_onlyUseTracksWithSVD) {
121 const auto hasNoSVD = [
this](
const RecoTrack * recoTrack) {
122 const auto& svdHitList = recoTrack->getSortedSVDHitList();
123 if (svdHitList.empty())
return true;
125 return m_param_reverseSeed ? svdHitList.back()->getSensorID().getLayerNumber() > 4
126 : svdHitList.front()->getSensorID().getLayerNumber() > 4;
128 TrackFindingCDC::erase_remove_if(m_recoTracksVector, hasNoSVD);
131 B2DEBUG(50,
"Now have " << m_spacePointVector.size() <<
" hits.");
133 m_stateCreatorFromTracks.apply(m_recoTracksVector, m_seedStates);
134 m_stateCreatorFromHits.apply(m_spacePointVector, m_states);
135 m_relationCreator.apply(m_seedStates, m_states, m_relations);
137 B2DEBUG(50,
"Created " << m_relations.size() <<
" relations.");
139 m_treeSearchFindlet.apply(m_seedStates, m_states, m_relations, m_results);
141 B2DEBUG(50,
"Having found " << m_results.size() <<
" results before overlap check");
144 return result.getHits().size() < m_param_minimalHitRequirement;
146 TrackFindingCDC::erase_remove_if(m_results, hasLowHitNumber);
148 m_overlapResolver.apply(m_results, m_filteredResults);
150 B2DEBUG(50,
"Having found " << m_filteredResults.size() <<
" results");
152 m_resultStorer.apply(m_filteredResults);
153 m_spacePointTagger.apply(m_filteredResults, m_spacePointVector);