Skip to content
Snippets Groups Projects
Commit 11c2836b authored by Christian Ohm's avatar Christian Ohm Committed by Graeme Stewart
Browse files

Adding patch from Jordi Duarte Campderros fixing issue discussed in ATLASG-711...

Adding patch from Jordi Duarte Campderros fixing issue discussed in ATLASG-711 (VrtSecInclusive-00-04-20)
parent 5a645ddd
No related branches found
No related tags found
No related merge requests found
......@@ -16,8 +16,6 @@ namespace VKalVrtAthena {
//
// get intersection point of track with various surfaces
//
MsgStream mlog( messageService(), name() );
//--------------------
// main loop
......
......@@ -26,8 +26,7 @@ namespace VKalVrtAthena {
{
MsgStream mlog( messageService(), name() );
mlog << MSG::DEBUG << " > SelGoodTrkParticle: begin" << endreq;
ATH_MSG_DEBUG( " > SelGoodTrkParticle: begin" );
xAOD::TrackParticleContainer *m_selectedBaseTracks ( nullptr );
ATH_CHECK( evtStore()->retrieve(m_selectedBaseTracks, "VrtSecInclusive_SelectedTrackParticles") );
......@@ -206,10 +205,6 @@ namespace VKalVrtAthena {
trk->auxdecor<bool>("is_selected") = true;
trk->auxdecor<unsigned long>("selected_trk_id") = a_trk->index();
// Also, need to put the element link both sides
ElementLink<xAOD::TrackParticleContainer> trackElementLink( *m_selectedBaseTracks, a_trk->index() );
trk->auxdecor<ElementLink<xAOD::TrackParticleContainer> >("selectedTrackLink") = trackElementLink;
ElementLink<xAOD::TrackParticleContainer> trackElementLink_orig( *trackParticleContainer, trk->index() );
a_trk->auxdata<ElementLink<xAOD::TrackParticleContainer> >("recoTrackLink") = trackElementLink_orig;
......@@ -233,10 +228,10 @@ namespace VKalVrtAthena {
} // end of track loop
if( m_FillNtuple ) {
mlog << MSG::DEBUG << " > SelGoodTrkParticle: Size of Inp Reco / w/barcode = "
<< trackParticleContainer->size() << ","
<< m_ntupleVars->get< vector<int> >( "RecoTrk_barcode" ).size()
<< endreq;
ATH_MSG_DEBUG( " > SelGoodTrkParticle: Size of Inp Reco / w/barcode = "
<< trackParticleContainer->size() << ","
<< m_ntupleVars->get< vector<int> >( "RecoTrk_barcode" ).size()
);
}
return StatusCode::SUCCESS;
......
......@@ -719,8 +719,9 @@ namespace VKalVrtAthena {
vertex->auxdata<char>( "truth_vtx_pattern" ) = truth_vtx_pattern;
if( noTruthVertex != truth_vtx_pattern ) {
ElementLink<xAOD::TruthVertexContainer> vtx_link;
if( noTruthVertex != truth_vtx_pattern ) {
// Retrieve the truth vertex container for element link
ATH_MSG_VERBOSE( "categorizeVertexTruthTopology(): Retrieve the truth vertex container for element link" );
const xAOD::TruthVertexContainer* truthVertexContainer ( nullptr );
......@@ -730,15 +731,15 @@ namespace VKalVrtAthena {
ATH_MSG_VERBOSE( "categorizeVertexTruthTopology(): create the element link" );
const auto *theVertex = get<0>( tmp_tuple );
if( theVertex ) {
ElementLink<xAOD::TruthVertexContainer> vtx_link( *truthVertexContainer, theVertex->index() );
// Add the truth vertex element link to the reco vertex
ATH_MSG_VERBOSE( "categorizeVertexTruthTopology(): Add the truth vertex element link to the reco vertex" );
vertex->auxdata<ElementLink<xAOD::TruthVertexContainer> >( "truth_vtx_link" ) = vtx_link;
// Add the truth vertex element link to the reco vertex
vtx_link.toIndexedElement(*truthVertexContainer,theVertex->index());
ATH_MSG_VERBOSE( "categorizeVertexTruthTopology(): Add the truth vertex element link to the reco vertex" );
}
}
}
// [JDC] a ElementLink decorator should be filled every event
// although using a null link
vertex->auxdata<ElementLink<xAOD::TruthVertexContainer> >( "truth_vtx_link" ) = vtx_link;
return StatusCode::SUCCESS;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment