JiveXML::TruthMuonTrackRetriever Class Reference

Retrieves the muons information from the TrackRecordCollection. More...

#include <TruthMuonTrackRetriever.h>

Inheritance diagram for JiveXML::TruthMuonTrackRetriever:

Collaboration diagram for JiveXML::TruthMuonTrackRetriever:

Public Member Functions
	TruthMuonTrackRetriever (const std::string &type, const std::string &name, const IInterface *parent)
	Standard Constructor.
virtual StatusCode	retrieve (ToolHandle< IFormatTool > &FormatTool)
	Retrieve all the data.
virtual std::string	dataTypeName () const
	Return the name of the data type.
StatusCode	initialize ()
	Default AthAlgTool methods.

Private Attributes
const std::string	m_typeName
	The data type that is generated by this retriever.
std::vector< std::string >	m_TrackRecCollNames
	A list of StoreGate names to probe in this order for the muon record collecton.

Detailed Description

Retrieves the muons information from the TrackRecordCollection.

• Properties
  • TrackRecCollNames = ["MuonEntryLayer","MuonEntryRecord"] :

• Retrieved Data
  • code : the PDG ID of the particle
  • id : the particle barcode
  • pt : transverse momentum
  • eta, phi : \(\eta\) and \(\phi\) of the momentum vector
  • rhoVertex,phiVertex,zVertex : position of the production vertex in \(\rho\), \(\phi\) and \(z\)

Definition at line 38 of file TruthMuonTrackRetriever.h.

Constructor & Destructor Documentation

TruthMuonTrackRetriever()

JiveXML::TruthMuonTrackRetriever::TruthMuonTrackRetriever	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Standard Constructor.

This is the standard AthAlgTool constructor.

Parameters

type	AlgTool type name
name	AlgTool instance name
parent	AlgTools parent owning this tool

Definition at line 28 of file TruthMuonTrackRetriever.cxx.

                                                                                                                        :
     base_class(type, name, parent),
     m_typeName("SMTr") {
     //Fill the list with the names of track record collections to try
     m_TrackRecCollNames.push_back("MuonEntryLayer");
     m_TrackRecCollNames.push_back("MuonEntryRecord");
 
     //Declare the properties
     declareProperty("TrackRecCollNames",m_TrackRecCollNames,"List of track records collections to try in this order - only one is retrieved");
 
   }

Member Function Documentation

dataTypeName()

virtual std::string JiveXML::TruthMuonTrackRetriever::dataTypeName ( ) const

inlinevirtual

Return the name of the data type.

Definition at line 49 of file TruthMuonTrackRetriever.h.

{ return m_typeName; };

initialize()

StatusCode JiveXML::TruthMuonTrackRetriever::initialize

(

)

Default AthAlgTool methods.

Initialize before event loop.

Definition at line 44 of file TruthMuonTrackRetriever.cxx.

                                                {
    
    //Nothing to be done here
    return StatusCode::SUCCESS;
  }

retrieve()

StatusCode JiveXML::TruthMuonTrackRetriever::retrieve ( ToolHandle< IFormatTool > & FormatTool )

virtual

Retrieve all the data.

Loop over all true particles, find the muons and get their basic parameters.

Parameters

FormatTool

the tool that will create formated output from the DataMap

Definition at line 54 of file TruthMuonTrackRetriever.cxx.

                                                                                  {
 
    ATH_MSG_DEBUG( "Retrieving " << dataTypeName() ); 
 
    //Try to retrieve the track record collection
    const TrackRecordCollection* TrackRecordColl = NULL ;
    //Loop over all the collections and try a retrieve (more efficenct than
    //contain-retrieve combination)
    for (const auto & CollNameItr : m_TrackRecCollNames ) {
      //be verbose
      ATH_MSG_DEBUG( "Trying to retrieve " << CollNameItr );
      //try to retrieve
      if ( !evtStore()->contains<TrackRecordCollection>( CollNameItr )){ continue; } // skip if not in SG
      if (evtStore()->retrieve(TrackRecordColl, CollNameItr).isSuccess()) break ;
    }
 
    //If we didnt' get any, return
    if (TrackRecordColl == NULL ) {
      ATH_MSG_WARNING( "Unable to retrieve any track collection from " << m_TrackRecCollNames );
      return StatusCode::RECOVERABLE;
    }
 
    //Reserve space for the output
    DataVect pt; pt.reserve(TrackRecordColl->size());
    DataVect phi; phi.reserve(TrackRecordColl->size());
    DataVect eta; eta.reserve(TrackRecordColl->size());
    DataVect rhoVertex; rhoVertex.reserve(TrackRecordColl->size());
    DataVect phiVertex; phiVertex.reserve(TrackRecordColl->size());
    DataVect zVertex; zVertex.reserve(TrackRecordColl->size());
    DataVect code; code.reserve(TrackRecordColl->size());
    DataVect id; id.reserve(TrackRecordColl->size()); 
 
    //Now loop over the collection and retrieve data
    for (const auto & record : *TrackRecordColl ) {
 
      //Get the pdg code
      int pdgCode = record.GetPDGCode();
 
      //Only accept muons
      if (abs(pdgCode) != 13) {
        ATH_MSG_DEBUG( "Reject non-muon track with PDG ID " << pdgCode );
        continue;
      }
 
      //Get vertex and momentum
      HepGeom::Point3D<double> vertex = record.GetPosition();
      HepGeom::Vector3D<double> momentum = record.GetMomentum();
 
      //And store output
      pt.emplace_back( momentum.perp()/CLHEP::GeV );
      phi.emplace_back( momentum.phi() < 0 ? momentum.phi() + 2*M_PI : momentum.phi() );
      eta.emplace_back( momentum.pseudoRapidity() );
      rhoVertex.emplace_back( vertex.perp()*CLHEP::mm/CLHEP::cm );
      phiVertex.emplace_back( vertex.phi() < 0 ? vertex.phi() + 2*M_PI : vertex.phi() );
      zVertex.emplace_back( vertex.z()*CLHEP::mm/CLHEP::cm );
      code.emplace_back( pdgCode );
      id.emplace_back( HepMC::barcode(record) ); // FIXME barcode-based
    }
    
    //Finall add everything to the datamap
    DataMap dataMap;
    const auto nEntries = pt.size();
    dataMap["pt"] = std::move(pt);
    dataMap["phi"] = std::move(phi);
    dataMap["eta"] = std::move(eta);
    dataMap["rhoVertex"] = std::move(rhoVertex);
    dataMap["phiVertex"] = std::move(phiVertex);
    dataMap["zVertex"] = std::move(zVertex);
    dataMap["code"] = std::move(code);
    dataMap["id"] = std::move(id);
    
    //some summary
    ATH_MSG_DEBUG( dataTypeName() << ": "<< nEntries );
 
    //forward data to formating tool
    //return FormatTool->AddToEvent(dataTypeName(), (*CollNameItr), &dataMap);
    std::string emptyStr="";
    return FormatTool->AddToEvent(dataTypeName(), emptyStr, &dataMap);
  }

Member Data Documentation

m_TrackRecCollNames

std::vector<std::string> JiveXML::TruthMuonTrackRetriever::m_TrackRecCollNames

private

A list of StoreGate names to probe in this order for the muon record collecton.

Definition at line 59 of file TruthMuonTrackRetriever.h.

m_typeName

const std::string JiveXML::TruthMuonTrackRetriever::m_typeName

private

The data type that is generated by this retriever.

Definition at line 56 of file TruthMuonTrackRetriever.h.

---

The documentation for this class was generated from the following files:

• TruthMuonTrackRetriever.h
• TruthMuonTrackRetriever.cxx