ATLAS Offline Software
Thin_vtxDuplicates.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
3 */
4 
6 // Thin_vtxDuplicates.cxx
8 // Matteo Bedognetti (matteo.bedognetti@cern.ch)
9 //Based on Thin_vtxTrk.cxx, by
10 
11 
12 
13 #include "Thin_vtxDuplicates.h"
16 #include <vector>
17 #include <string>
18 #include <algorithm> // for the sort function
19 #include <iomanip>
21 // Constructor
22 DerivationFramework::Thin_vtxDuplicates::Thin_vtxDuplicates(const std::string& t, const std::string& n, const IInterface* p ) :
23  base_class(t,n,p),
24  // m_acceptanceR(-1.),
25  m_noFlags(false),
26  m_nVtxTot(0),
27  m_nVtxPass(0)
28 {
29  declareProperty("VertexContainerName" , m_vertexContainerNames);
30  declareProperty("PassFlags" , m_passFlags);
31  //declareProperty("AcceptanceRadius" , m_acceptanceR);
32  declareProperty("IgnoreFlags" , m_noFlags);
33  //declareProperty("ApplyAndForTracks" , m_trackAnd = false);
34  //declareProperty("ThinTracks" , m_thinTracks = true);
35 }
36 
37 // Destructor
39 
40 // Athena initialize and finalize
42 {
43  // Decide which collections need to be checked for ID TrackParticles
44  ATH_MSG_VERBOSE("initialize() ...");
45  ATH_CHECK(m_vertexContainerNames.initialize(m_streamName));
46 
47 
48  if (m_passFlags.empty()) {
49  ATH_MSG_FATAL("No pass flags provided for thinning.");
50  return StatusCode::FAILURE;
51  } else {
52  for(auto itr = m_passFlags.cbegin(); itr!=m_passFlags.cend(); ++itr) {
53  ATH_MSG_INFO("Vertices must pass the \"" << itr->key() << "\" selection");
54  }
55  }
56 
57  for(auto &key : m_passFlags){
58  key = m_vertexContainerNames.key() + '.' + key.key();
59  }
60  ATH_CHECK(m_passFlags.initialize());
61  return StatusCode::SUCCESS;
62 }
63 
65 {
66  ATH_MSG_VERBOSE("finalize() ...");
67  ATH_MSG_INFO("Processed "<< m_nVtxTot <<" vertices, "<< m_nVtxPass<< " were retained ");
68 
69  return StatusCode::SUCCESS;
70 }
71 
72 // The thinning itself
74 {
75  // retieve vertex
76  const EventContext& ctx = Gaudi::Hive::currentContext();
77  SG::ThinningHandle< xAOD::VertexContainer > vertexContainer(m_vertexContainerNames, ctx);
78  std::vector<bool> vtxMask(vertexContainer->size(), true); // default: keep all vertices
79  int vtxTot = 0;
80  int nVtxPass = 0;
81  // loop over vertices
82  int k = 0;
83  std::vector<SG::ReadDecorHandle<xAOD::VertexContainer, Char_t>> handles;
84  handles.reserve(m_passFlags.size());
85  for(const auto &key : m_passFlags){
86  handles.emplace_back(key, ctx);
87  if(!handles.back().isPresent()) return StatusCode::FAILURE;
88  }
89  for(auto vtxItr = vertexContainer->cbegin(); vtxItr!=vertexContainer->cend(); ++vtxItr, ++k) {
90  const xAOD::Vertex* vtx = *vtxItr;
91  // check if the vertex passed the required selections criteria (is run when the vertex is already excluded, because the counter needs the info)
92  bool passed = false;
93  if(m_noFlags){passed = true; vtxTot++; }
94  else{
95  for(auto &flagAcc : handles) {
96  if(flagAcc(*vtx) != 0) {
97  passed = true;
98  vtxTot++;//Have to count the ones which are accepted to start with
99  break;
100  }
101  } // end of loop over flags
102  }
103 
104  // Skip if it has already been identified as duplicate
105  if(vtxMask[k] == false)continue; //After the flag-check to have the total-passed work correctly
106 
107  if(!passed)vtxMask[k]= false;
108 
109  if(passed) {
110  // vertex passed the selection
111  nVtxPass++;
112 
113  // determine the sum of the tracks at vertex as centre for the cone
114  std::vector<const xAOD::TrackParticle*> presentVertex, compareVertex;
115 
116  //Fill in the present vertex, for later comparison against other vertices
117  presentVertex.clear();
118  for(uint j=0; j<vtx->nTrackParticles(); ++j) {
119  presentVertex.push_back(vtx->trackParticle(j));
120  }
121  sort( presentVertex.begin(), presentVertex.end() ); //Sort the trackparticles BY POINTER ADDRESS
122 
123  //Loop over the remaining vertices and remove them if needed
124  int loop_k = k+1;
125  for(auto vtxLoopItr = vtxItr+1; vtxLoopItr!=vertexContainer->cend(); vtxLoopItr++, loop_k++){
126 
127  const xAOD::Vertex* loop_vtx = *vtxLoopItr;
128 
129  //Vertices are distinct if have different size
130  if(vtx->nTrackParticles() != loop_vtx->nTrackParticles())continue;
131 
132  //If the vertex is still active load and compare
133  if(vtxMask[loop_k]){
134 
135  compareVertex.clear();
136  for(uint j=0; j<loop_vtx->nTrackParticles(); ++j) {
137  compareVertex.push_back(loop_vtx->trackParticle(j));
138  }
139 
140  std::sort( compareVertex.begin(), compareVertex.end());
141 
142  vtxMask[loop_k] = false;
143 
144  ATH_MSG_DEBUG("Compared tracks: ");
145  ATH_MSG_DEBUG(std::setw(14)<<compareVertex[0]<<std::setw(14) << compareVertex[1]<<std::setw(14)<<compareVertex[2]);
146  ATH_MSG_DEBUG(std::setw(14)<<presentVertex[0]<<std::setw(14) << presentVertex[1]<<std::setw(14)<<presentVertex[2]);
147 
148  for(uint j=0; j<loop_vtx->nTrackParticles(); ++j) {
149  if( compareVertex[j] != presentVertex[j] ){vtxMask[loop_k] = true; break;}
150  }
151  ATH_MSG_DEBUG("Verdict:"<<(vtxMask[loop_k]? "keep": "erase") );
152  }
153 
154  } // Endo of extra loop over remaining vertices
155 
156  } // if( passed )
157  } // end of loop over vertices
158 
159  // Execute the thinning service based on the vtxMask.
160  vertexContainer.keep(vtxMask);
161 
162  m_nVtxTot.fetch_add( vtxTot, std::memory_order_relaxed);
163  m_nVtxPass.fetch_add( nVtxPass, std::memory_order_relaxed);
164 
165 
166 
167  return StatusCode::SUCCESS;
168 }
169 
DerivationFramework::Thin_vtxDuplicates::Thin_vtxDuplicates
Thin_vtxDuplicates(const std::string &t, const std::string &n, const IInterface *p)
Definition: Thin_vtxDuplicates.cxx:22
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
xAOD::Vertex_v1::nTrackParticles
size_t nTrackParticles() const
Get the number of tracks associated with this vertex.
Definition: Vertex_v1.cxx:270
TrigCompositeUtils::passed
bool passed(DecisionID id, const DecisionIDContainer &idSet)
checks if required decision ID is in the set of IDs in the container
Definition: TrigCompositeUtilsRoot.cxx:117
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
ThinningHandle.h
Handle for requesting thinning for a data object.
DerivationFramework::Thin_vtxDuplicates::m_passFlags
SG::ReadDecorHandleKeyArray< xAOD::VertexContainer > m_passFlags
Definition: Thin_vtxDuplicates.h:32
DerivationFramework::Thin_vtxDuplicates::initialize
virtual StatusCode initialize()
Definition: Thin_vtxDuplicates.cxx:41
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
SG::ThinningHandle
Handle for requesting thinning for a data object.
Definition: ThinningHandle.h:84
uint
unsigned int uint
Definition: LArOFPhaseFill.cxx:20
SG::ThinningHandleBase::keep
void keep(size_t ndx)
Mark that index ndx in the container should be kept (not thinned away).
Definition: ThinningHandleBase.cxx:75
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
Thin_vtxDuplicates.h
DerivationFramework::Thin_vtxDuplicates::m_vertexContainerNames
SG::ThinningHandleKey< xAOD::VertexContainer > m_vertexContainerNames
Definition: Thin_vtxDuplicates.h:31
beamspotman.n
n
Definition: beamspotman.py:729
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
BPhysHypoHelper.h
: B-physics xAOD helpers.
DerivationFramework::Thin_vtxDuplicates::finalize
virtual StatusCode finalize()
Definition: Thin_vtxDuplicates.cxx:64
DerivationFramework::Thin_vtxDuplicates::~Thin_vtxDuplicates
~Thin_vtxDuplicates()
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
xAOD::Vertex_v1::trackParticle
const TrackParticle * trackParticle(size_t i) const
Get the pointer to a given track that was used in vertex reco.
Definition: Vertex_v1.cxx:249
DerivationFramework::Thin_vtxDuplicates::doThinning
virtual StatusCode doThinning() const
Definition: Thin_vtxDuplicates.cxx:73
DerivationFramework::Thin_vtxDuplicates::m_noFlags
bool m_noFlags
Definition: Thin_vtxDuplicates.h:27
xAOD::Vertex_v1
Class describing a Vertex.
Definition: Vertex_v1.h:42
TrackParticleContainer.h
fitman.k
k
Definition: fitman.py:528
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37