ATLAS Offline Software
Public Member Functions | Private Attributes | List of all members
DerivationFramework::TruthCollectionMaker Class Reference

#include <TruthCollectionMaker.h>

Inheritance diagram for DerivationFramework::TruthCollectionMaker:
Collaboration diagram for DerivationFramework::TruthCollectionMaker:

Public Member Functions

 TruthCollectionMaker (const std::string &t, const std::string &n, const IInterface *p)
 
 ~TruthCollectionMaker ()
 
StatusCode initialize ()
 
StatusCode finalize ()
 
virtual StatusCode addBranches () const
 

Private Attributes

std::atomic< unsigned int > m_ntotpart
 
std::atomic< unsigned int > m_npasspart
 
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_particlesKey {this, "ParticlesKey", "TruthParticles", "ReadHandleKey for input TruthParticleContainer"}
 
SG::WriteHandleKey< xAOD::TruthParticleContainer > m_outputParticlesKey {this, "NewCollectionName", "OutputTruthCollection", "WriteHandleKey for new TruthParticleContainer"}
 
Gaudi::Property< std::string > m_partString {this, "ParticleSelectionString", "", "ExpressionEvaluation string for particle selection"}
 
Gaudi::Property< bool > m_do_compress {this, "Do_Compress", false, "Removes particles with the same pdgId in a decay chain (but keeps first and last)"}
 
Gaudi::Property< bool > m_do_sherpa {this, "Do_Sherpa", false, "Checks if there are truth W bosons in the current record. If not, tries to combine W daughters to create one"}
 
Gaudi::Property< bool > m_keep_navigation_info {this, "KeepNavigationInfo", true, "m_do_sherpa currently only works for W+jets"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_linkDecoratorKey {this, "originalTruthParticle", m_outputParticlesKey, "originalTruthParticle", "Name of the decoration linking to the original truth particle"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_originDecoratorKey {this, "classifierParticleOrigin", m_outputParticlesKey, "classifierParticleOrigin", "Name of the decoration which records the particle origin as determined by the MCTruthClassifier"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_typeDecoratorKey {this, "classifierParticleType", m_outputParticlesKey, "classifierParticleType", "Name of the decoration which records the particle type as determined by the MCTruthClassifier"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_outcomeDecoratorKey {this, "classifierParticleOutCome", m_outputParticlesKey, "classifierParticleOutCome", "Name of the decoration which records the particle outcome as determined by the MCTruthClassifier"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_classificationDecoratorKey {this, "Classification", m_outputParticlesKey, "Classification", "Name of the decoration which records the particle classification as determined by the MCTruthClassifier"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_motherIDDecoratorKey {this, "motherID", m_outputParticlesKey, "motherID", "Name of the decoration which records the ID of the particle's mother"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_daughterIDDecoratorKey {this, "daughterID", m_outputParticlesKey, "daughterID", "Name of the decoration which records the ID of the particle's daughter"}
 
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_hadronOriginDecoratorKey {this, "TopHadronOriginFlag", m_outputParticlesKey, "TopHadronOriginFlag", "Name of the decoration which records the origin of hadrons from top decays"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_originReadDecorKey {this, "inputClassifierParticleOrigin", m_particlesKey, "classifierParticleOrigin", "Particle origin"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_typeReadDecorKey {this, "inputClassifierParticleType", m_particlesKey, "classifierParticleType", "Particle type"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_outcomeReadDecorKey {this, "inputClassifierParticleOutCome", m_particlesKey, "classifierParticleOutCome", "Particle outcome"}
 
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_classificationReadDecorKey {this, "inputClassification", m_particlesKey, "Classification", "Classification code"}
 
ServiceHandle< StoreGateSvc > m_metaStore
 Handle on the metadata store for init. More...
 

Detailed Description

Definition at line 29 of file TruthCollectionMaker.h.

Constructor & Destructor Documentation

◆ TruthCollectionMaker()

DerivationFramework::TruthCollectionMaker::TruthCollectionMaker ( const std::string &  t,
const std::string &  n,
const IInterface *  p 
)

Definition at line 29 of file TruthCollectionMaker.cxx.

32  : base_class(t,n,p)
33  , m_ntotpart(0)
34  , m_npasspart(0)
35  , m_metaStore( "MetaDataStore", n )
36 {
37  declareProperty("MetaDataStore", m_metaStore );
38 }

◆ ~TruthCollectionMaker()

DerivationFramework::TruthCollectionMaker::~TruthCollectionMaker ( )

Definition at line 41 of file TruthCollectionMaker.cxx.

41  {
42 }

Member Function Documentation

◆ addBranches()

StatusCode DerivationFramework::TruthCollectionMaker::addBranches ( ) const
virtual

Definition at line 95 of file TruthCollectionMaker.cxx.

96 {
97  // Event context for AthenaMT
98  const EventContext& ctx = Gaudi::Hive::currentContext();
99 
100  // Set up for some metadata handling
101  // TODO: this isn't MT compliant. This information should go into the config level and avoid meta store
102  static const bool is_sherpa = [this]() {
103  bool is_sherpa = false;
104  if (m_metaStore->contains<xAOD::TruthMetaDataContainer>("TruthMetaData")){
105  // Note that I'd like to get this out of metadata in general, but it seems that the
106  // metadata isn't fully available in initialize, and since this is a const function
107  // I can only do the retrieve every event, rather than lazy-initializing, since this
108  // metadata ought not change during a run
109  const xAOD::TruthMetaDataContainer* truthMetaData(nullptr);
110  // Shamelessly stolen from the file meta data tool
111  if (m_metaStore->retrieve(truthMetaData).isSuccess() && !truthMetaData->empty()){
112  // Let's just be super sure...
113  const std::string gens = truthMetaData->at(0)->generators();
114  is_sherpa = (gens.find("sherpa")==std::string::npos &&
115  gens.find("Sherpa")==std::string::npos &&
116  gens.find("SHERPA")==std::string::npos) ? false : true;
117  } // Seems to be the only sure way...
118  else {
119  ATH_MSG_WARNING("Found xAODTruthMetaDataContainer empty! Configuring to be NOT Sherpa.");
120  }
121  ATH_MSG_INFO("From metadata configured: Sherpa? " << is_sherpa);
122  } else {
123  ATH_MSG_WARNING("Could not find metadata container in storegate; assuming NOT Sherpa");
124  }
125  return is_sherpa;
126  }();
127 
128  // Retrieve truth collections
129  SG::ReadHandle<xAOD::TruthParticleContainer> truthParticles(m_particlesKey,ctx);
130  if (!truthParticles.isValid()) {
131  ATH_MSG_ERROR("Couldn't retrieve TruthParticle collection with name " << m_particlesKey);
132  return StatusCode::FAILURE;
133  }
134 
135  // Create the new particle containers and WriteHandles
136  SG::WriteHandle<xAOD::TruthParticleContainer> newParticlesWriteHandle(m_outputParticlesKey, ctx);
137  ATH_CHECK(newParticlesWriteHandle.record(std::make_unique<xAOD::TruthParticleContainer>(),
138  std::make_unique<xAOD::TruthParticleAuxContainer>()));
139  ATH_MSG_DEBUG( "Recorded new TruthParticleContainer with key: " << (m_outputParticlesKey.key()));
140 
141  // Set up a mask with the same entries as the full collections
142  unsigned int nParticles = truthParticles->size();
143  m_ntotpart += nParticles;
144 
145  // Set up decor readers
146  SG::ReadDecorHandle<xAOD::TruthParticleContainer, unsigned int > originReadDecor(m_originReadDecorKey, ctx);
147  SG::ReadDecorHandle<xAOD::TruthParticleContainer, unsigned int > typeReadDecor(m_typeReadDecorKey, ctx);
148  SG::ReadDecorHandle<xAOD::TruthParticleContainer, unsigned int > outcomeReadDecor(m_outcomeReadDecorKey, ctx);
149  SG::ReadDecorHandle<xAOD::TruthParticleContainer, unsigned int > classificationReadDecor(m_classificationReadDecorKey, ctx);
150 
151  // Set up decorators
152  SG::WriteDecorHandle<xAOD::TruthParticleContainer, ElementLink<xAOD::TruthParticleContainer> > linkDecorator(m_linkDecoratorKey, ctx);
153  SG::WriteDecorHandle<xAOD::TruthParticleContainer, unsigned int > originDecorator(m_originDecoratorKey, ctx);
154  SG::WriteDecorHandle<xAOD::TruthParticleContainer, unsigned int > typeDecorator(m_typeDecoratorKey, ctx);
155  SG::WriteDecorHandle<xAOD::TruthParticleContainer, unsigned int > outcomeDecorator(m_outcomeDecoratorKey, ctx);
156  SG::WriteDecorHandle<xAOD::TruthParticleContainer, unsigned int > classificationDecorator(m_classificationDecoratorKey, ctx);
157  SG::WriteDecorHandle< xAOD::TruthParticleContainer, int > motherIDDecorator(m_motherIDDecoratorKey, ctx);
158  SG::WriteDecorHandle< xAOD::TruthParticleContainer, int > daughterIDDecorator(m_daughterIDDecoratorKey, ctx);
159  SG::WriteDecorHandle< xAOD::TruthParticleContainer, int > hadronOriginDecorator(m_hadronOriginDecoratorKey, ctx);
160 
161  // Execute the text parsers and update the mask
162  if (!m_partString.empty()) {
163  std::vector<int> entries = m_parser->evaluateAsVector();
164  unsigned int nEntries = entries.size();
165  // check the sizes are compatible
166  if (nParticles != nEntries ) {
167  ATH_MSG_ERROR("Sizes incompatible! Are you sure your selection string used TruthParticles?");
168  return StatusCode::FAILURE;
169  } else {
170  // add relevant particles to new collection
171 
172  //---------------
173  //This is some code to *add* new particles. Probably a good idea to break this off as a sub-function, but I'll let James C decide where that should go.
174  //---------------
175 
176  //Let's check if we want to build W/Z bosons
177  bool SherpaW = false;
178  bool SherpaZ = false;
179  if (m_partString.value().find("24") != std::string::npos && m_do_sherpa && is_sherpa) {
180  SherpaW = true;
181  }
182  if (m_partString.value().find("23") != std::string::npos && m_do_sherpa && is_sherpa) {
183  SherpaZ = true;
184  }
185  if ((SherpaW or SherpaZ) && is_sherpa){
186  if (std::accumulate(entries.begin(),entries.end(),0) > 0){ //We actually have some W and Z bosons in there.
187  SherpaW = false;
188  SherpaZ = false;
189  }
190  }
191 
192  if ((SherpaW || SherpaZ) && is_sherpa){
193  // Currently only handles un-ambiguous cases
194  std::vector<const xAOD::TruthParticle*> status_nonPhysical;
195  for (unsigned int i=0; i<nParticles; ++i) {
196  // Nullptr check
197  if (!truthParticles->at(i)) continue;
198  // Only collect leptons
199  if (!MC::isSMLepton(truthParticles->at(i))) continue;
200  // Gather by status
201  if (!MC::isPhysical(truthParticles->at(i))) status_nonPhysical.push_back( truthParticles->at(i) );
202  } // Done with loop over truth particles
203  // Boson cases that we can actually deal with -- generically up to VVV
204  if ((status_nonPhysical.size()==2 || status_nonPhysical.size()==4 || status_nonPhysical.size()==6) && (SherpaZ || SherpaW)){
205  // Basic boson pairing...
206  int gens[3] = {0,0,0};
207  for (size_t i=0;i<status_nonPhysical.size();++i){
208  if (status_nonPhysical[i]->absPdgId()<13) gens[0]++;
209  else if (status_nonPhysical[i]->absPdgId()<15) gens[1]++;
210  else gens[2]++;
211  } // Loop over status_nonPhysical particles
212  // Should only have even numbers per generation. Any number greater than 2 or ==1 and we're dead
213  if (gens[0]>2 || gens[0]==1 || gens[1]>2 || gens[1]==1 || gens[2]>2 || gens[2]==1){
214  // In agreeing with Sherpa authors, these are Q-M ambiguous states. Do not let users be evil.
215  ATH_MSG_VERBOSE("Too many leptons of one generation. Cannot make bosons. Give up");
216  return StatusCode::SUCCESS;
217  }
218  std::vector<const xAOD::TruthParticle*> boson;
219  for (size_t i=0;i<status_nonPhysical.size();++i){
220  if (status_nonPhysical[i]->absPdgId()<13) boson.push_back(status_nonPhysical[i]);
221  else if (gens[0]==0 && status_nonPhysical[i]->absPdgId()<15) boson.push_back(status_nonPhysical[i]);
222  else if (gens[0]==0 && gens[1]==0) boson.push_back(status_nonPhysical[i]);
223  if (boson.size()==2){
224  // Make a boson! Just have to figure out _which_ boson!
225  int pdg_id=0;
226  // Easy case: Z boson
227  if (boson[0]->pdgId()==-boson[1]->pdgId()) pdg_id=23;
228  else if (std::abs(boson[0]->pdgId()+boson[1]->pdgId())!=1){
229  // No idea what you were
230  ATH_MSG_WARNING("Do not know how to interpret as a boson: " << boson[0]->pdgId() << " " << boson[1]->pdgId());
231  } else {
232  // W boson
233  pdg_id=24*(boson[0]->pdgId()+boson[1]->pdgId());
234  }
235  if ( (SherpaW && MC::isW(pdg_id)) ||
236  (SherpaZ && MC::isZ(pdg_id)) ){
237  // Make a Z or a W
238  xAOD::TruthParticle* xTruthParticle = new xAOD::TruthParticle();
239  newParticlesWriteHandle->push_back( xTruthParticle );
240  if (m_keep_navigation_info){
241  if (boson[0]->hasProdVtx()) {
242  if ((boson[0]->prodVtx()->nIncomingParticles() > 0) && (boson[0]->prodVtx()->incomingParticle(0)!=nullptr)) {
243  motherIDDecorator(*xTruthParticle) = boson[0]->prodVtx()->incomingParticle(0)->pdgId();
244  } else {motherIDDecorator(*xTruthParticle) = 0;}
245  } else {motherIDDecorator(*xTruthParticle) = 0;}
246  if (boson[0]->hasDecayVtx()) {
247  if ((boson[0]->decayVtx()->nOutgoingParticles() > 0) && (boson[0]->decayVtx()->outgoingParticle(0)!=nullptr)) {
248  daughterIDDecorator(*xTruthParticle) = boson[0]->decayVtx()->outgoingParticle(0)->pdgId();
249  } else {daughterIDDecorator(*xTruthParticle) = 0;}
250  } else {daughterIDDecorator(*xTruthParticle) = 0;}
251  }
252  // Set with what makes sense here
253  xTruthParticle->setPdgId(pdg_id);
254  // Set dummy values
255  xTruthParticle->setBarcode(HepMC::INVALID_PARTICLE_ID);
256  xTruthParticle->setStatus(3);
257  // Use the sum of the momenta
258  xAOD::IParticle::FourMom_t new_mom = boson[0]->p4()+boson[1]->p4();
259  xTruthParticle->setM(new_mom.M());
260  xTruthParticle->setPx(new_mom.Px());
261  xTruthParticle->setPy(new_mom.Py());
262  xTruthParticle->setPz(new_mom.Pz());
263  xTruthParticle->setE(new_mom.E());
264  }
265  // Now clear the vectors
266  boson.clear();
267  // And move to the next generation
268  if (gens[0]==0 && gens[1]==0) gens[2]=0;
269  else if (gens[0]==0) gens[1]=0;
270  else gens[0]=0;
271  } // Done making a boson
272  } // Done looping over particles
273  }
274  if (status_nonPhysical.size()==1 || status_nonPhysical.size()==3 || status_nonPhysical.size()==5 || status_nonPhysical.size()>6){
275  ATH_MSG_WARNING(status_nonPhysical.size() << " leptons found in the Sherpa event record. Not sure how to deal with this.");
276  }
277  return StatusCode::SUCCESS;
278  }
279  for (unsigned int i=0; i<nParticles; ++i) {
280  ElementLink<xAOD::TruthParticleContainer> eltp(*truthParticles,i);
281  if (entries[i]==1) {
282  //In TRUTH3, we want to remove all particles but the first and last in a decay chain. This is off in TRUTH1. The first and last particles in the decay chain are decorated as such.
283 
284  const xAOD::TruthParticle* theParticle = (*truthParticles)[i];
285  if (m_do_compress){
286  bool same_as_mother = false;
287  bool same_as_daughter = false;
288  if (theParticle->hasProdVtx()){
289  if ((theParticle->prodVtx()->nIncomingParticles() > 0) && (theParticle->prodVtx()->incomingParticle(0)!=nullptr)) {
290  if (theParticle->prodVtx()->incomingParticle(0)->pdgId() == theParticle->pdgId()){
291  same_as_mother = true;
292  }
293  }
294  }
295  if (theParticle->hasDecayVtx()){
296  if ((theParticle->decayVtx()->nOutgoingParticles() > 0) && (theParticle->decayVtx()->outgoingParticle(0)!=nullptr)) {
297  if (theParticle->decayVtx()->outgoingParticle(0)->pdgId() == theParticle->pdgId()){
298  same_as_daughter = true;
299  }
300  }
301  }
302  if (same_as_mother && same_as_daughter){
303  entries[i]=0;
304  continue;
305  }
306  }
307  xAOD::TruthParticle* xTruthParticle = new xAOD::TruthParticle();
308  newParticlesWriteHandle->push_back( xTruthParticle );
309  if (m_keep_navigation_info){
310  if (theParticle->hasProdVtx()) {
311  if ((theParticle->prodVtx()->nIncomingParticles() > 0) && (theParticle->prodVtx()->incomingParticle(0)!=nullptr)) {
312  motherIDDecorator(*xTruthParticle) = theParticle->prodVtx()->incomingParticle(0)->pdgId();
313  } else {motherIDDecorator(*xTruthParticle) = 0;}
314  } else {motherIDDecorator(*xTruthParticle) = 0;}
315  if (theParticle->hasDecayVtx()) {
316  if ((theParticle->decayVtx()->nOutgoingParticles() > 0) && (theParticle->decayVtx()->outgoingParticle(0)!=nullptr)) {
317  daughterIDDecorator(*xTruthParticle) = theParticle->decayVtx()->outgoingParticle(0)->pdgId();
318  } else {daughterIDDecorator(*xTruthParticle) = 0;}
319  } else {daughterIDDecorator(*xTruthParticle) = 0;}
320  }
321  // Fill with numerical content
322  *xTruthParticle=*theParticle;
323  // Copy over the decorations if they are available
324  typeDecorator(*xTruthParticle) =
325  typeReadDecor.withDefault(*theParticle, 0);
326 
327  originDecorator(*xTruthParticle) =
328  originReadDecor.withDefault(*theParticle, 0);
329 
330  outcomeDecorator(*xTruthParticle) =
331  outcomeReadDecor.withDefault(*theParticle, 0);
332 
333  classificationDecorator(*xTruthParticle) =
334  classificationReadDecor.withDefault(*theParticle, 0);
335 
336  if (m_outputParticlesKey.key()=="TruthHFHadrons"){
337  static const SG::ConstAccessor<int> TopHadronOriginFlagAcc("TopHadronOriginFlag");
338  hadronOriginDecorator(*xTruthParticle) =
339  TopHadronOriginFlagAcc.withDefault (*theParticle, 0);
340  }
341 
342  if(m_keep_navigation_info) linkDecorator(*xTruthParticle) = eltp;
343  }
344  }
345  }
346  // Count the mask
347  for (unsigned int i=0; i<nParticles; ++i) if (entries[i]) ++m_npasspart;
348  }
349 
350  return StatusCode::SUCCESS;
351 }

◆ finalize()

StatusCode DerivationFramework::TruthCollectionMaker::finalize ( )

Definition at line 85 of file TruthCollectionMaker.cxx.

86 {
87  ATH_MSG_VERBOSE("finalize() ...");
88  //ATH_MSG_INFO("Processed "<< m_ntotvtx <<" truth vertices, "<< m_npassvtx << " were retained ");
89  ATH_MSG_INFO("Processed "<< m_ntotpart <<" truth particles, "<< m_npasspart << " were retained ");
90  ATH_CHECK( finalizeParser() );
91  return StatusCode::SUCCESS;
92 }

◆ initialize()

StatusCode DerivationFramework::TruthCollectionMaker::initialize ( )

Definition at line 45 of file TruthCollectionMaker.cxx.

46 {
47  ATH_MSG_VERBOSE("initialize() ...");
48  // Input truth particles
49  ATH_CHECK( m_particlesKey.initialize() );
50  ATH_MSG_INFO("Using " << m_particlesKey.key() << " as the input truth container key");
51 
52  // Output (new) truth particles
53  ATH_CHECK(m_outputParticlesKey.initialize());
54  ATH_MSG_INFO("New truth particles container key: " << m_outputParticlesKey.key() );
55 
56  if (m_partString.empty()) {
57  ATH_MSG_FATAL("No selection string provided");
58  return StatusCode::FAILURE;
59  } else {ATH_MSG_INFO("Truth particle selection string: " << m_partString );}
60 
61  // Set up the text-parsing machinery for thinning the truth directly according to user cuts
62  if (!m_partString.empty()) {
63  ATH_CHECK( initializeParser(m_partString) );
64  }
65 
66  // Initialise read handle keys
67  ATH_CHECK(m_originReadDecorKey.initialize());
68  ATH_CHECK(m_typeReadDecorKey.initialize());
69  ATH_CHECK(m_outcomeReadDecorKey.initialize());
70  ATH_CHECK(m_classificationReadDecorKey.initialize());
71 
72  // Initialise decorator handle keys
73  ATH_CHECK(m_linkDecoratorKey.initialize());
74  ATH_CHECK(m_originDecoratorKey.initialize());
75  ATH_CHECK(m_typeDecoratorKey.initialize());
76  ATH_CHECK(m_outcomeDecoratorKey.initialize());
77  ATH_CHECK(m_classificationDecoratorKey.initialize());
78  ATH_CHECK(m_motherIDDecoratorKey.initialize());
79  ATH_CHECK(m_daughterIDDecoratorKey.initialize());
80  ATH_CHECK(m_hadronOriginDecoratorKey.initialize());
81 
82  return StatusCode::SUCCESS;
83 }

Member Data Documentation

◆ m_classificationDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_classificationDecoratorKey {this, "Classification", m_outputParticlesKey, "Classification", "Name of the decoration which records the particle classification as determined by the MCTruthClassifier"}
private

Definition at line 62 of file TruthCollectionMaker.h.

◆ m_classificationReadDecorKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_classificationReadDecorKey {this, "inputClassification", m_particlesKey, "Classification", "Classification code"}
private

Definition at line 79 of file TruthCollectionMaker.h.

◆ m_daughterIDDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_daughterIDDecoratorKey {this, "daughterID", m_outputParticlesKey, "daughterID", "Name of the decoration which records the ID of the particle's daughter"}
private

Definition at line 66 of file TruthCollectionMaker.h.

◆ m_do_compress

Gaudi::Property<bool> DerivationFramework::TruthCollectionMaker::m_do_compress {this, "Do_Compress", false, "Removes particles with the same pdgId in a decay chain (but keeps first and last)"}
private

Definition at line 47 of file TruthCollectionMaker.h.

◆ m_do_sherpa

Gaudi::Property<bool> DerivationFramework::TruthCollectionMaker::m_do_sherpa {this, "Do_Sherpa", false, "Checks if there are truth W bosons in the current record. If not, tries to combine W daughters to create one"}
private

Definition at line 49 of file TruthCollectionMaker.h.

◆ m_hadronOriginDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_hadronOriginDecoratorKey {this, "TopHadronOriginFlag", m_outputParticlesKey, "TopHadronOriginFlag", "Name of the decoration which records the origin of hadrons from top decays"}
private

Definition at line 68 of file TruthCollectionMaker.h.

◆ m_keep_navigation_info

Gaudi::Property<bool> DerivationFramework::TruthCollectionMaker::m_keep_navigation_info {this, "KeepNavigationInfo", true, "m_do_sherpa currently only works for W+jets"}
private

Definition at line 51 of file TruthCollectionMaker.h.

◆ m_linkDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_linkDecoratorKey {this, "originalTruthParticle", m_outputParticlesKey, "originalTruthParticle", "Name of the decoration linking to the original truth particle"}
private

Definition at line 54 of file TruthCollectionMaker.h.

◆ m_metaStore

ServiceHandle<StoreGateSvc> DerivationFramework::TruthCollectionMaker::m_metaStore
private

Handle on the metadata store for init.

Definition at line 82 of file TruthCollectionMaker.h.

◆ m_motherIDDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_motherIDDecoratorKey {this, "motherID", m_outputParticlesKey, "motherID", "Name of the decoration which records the ID of the particle's mother"}
private

Definition at line 64 of file TruthCollectionMaker.h.

◆ m_npasspart

std::atomic<unsigned int> DerivationFramework::TruthCollectionMaker::m_npasspart
private

Definition at line 38 of file TruthCollectionMaker.h.

◆ m_ntotpart

std::atomic<unsigned int> DerivationFramework::TruthCollectionMaker::m_ntotpart
mutableprivate

Definition at line 38 of file TruthCollectionMaker.h.

◆ m_originDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_originDecoratorKey {this, "classifierParticleOrigin", m_outputParticlesKey, "classifierParticleOrigin", "Name of the decoration which records the particle origin as determined by the MCTruthClassifier"}
private

Definition at line 56 of file TruthCollectionMaker.h.

◆ m_originReadDecorKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_originReadDecorKey {this, "inputClassifierParticleOrigin", m_particlesKey, "classifierParticleOrigin", "Particle origin"}
private

Definition at line 73 of file TruthCollectionMaker.h.

◆ m_outcomeDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_outcomeDecoratorKey {this, "classifierParticleOutCome", m_outputParticlesKey, "classifierParticleOutCome", "Name of the decoration which records the particle outcome as determined by the MCTruthClassifier"}
private

Definition at line 60 of file TruthCollectionMaker.h.

◆ m_outcomeReadDecorKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_outcomeReadDecorKey {this, "inputClassifierParticleOutCome", m_particlesKey, "classifierParticleOutCome", "Particle outcome"}
private

Definition at line 77 of file TruthCollectionMaker.h.

◆ m_outputParticlesKey

SG::WriteHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_outputParticlesKey {this, "NewCollectionName", "OutputTruthCollection", "WriteHandleKey for new TruthParticleContainer"}
private

Definition at line 42 of file TruthCollectionMaker.h.

◆ m_particlesKey

SG::ReadHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_particlesKey {this, "ParticlesKey", "TruthParticles", "ReadHandleKey for input TruthParticleContainer"}
private

Definition at line 40 of file TruthCollectionMaker.h.

◆ m_partString

Gaudi::Property<std::string> DerivationFramework::TruthCollectionMaker::m_partString {this, "ParticleSelectionString", "", "ExpressionEvaluation string for particle selection"}
private

Definition at line 45 of file TruthCollectionMaker.h.

◆ m_typeDecoratorKey

SG::WriteDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_typeDecoratorKey {this, "classifierParticleType", m_outputParticlesKey, "classifierParticleType", "Name of the decoration which records the particle type as determined by the MCTruthClassifier"}
private

Definition at line 58 of file TruthCollectionMaker.h.

◆ m_typeReadDecorKey

SG::ReadDecorHandleKey<xAOD::TruthParticleContainer> DerivationFramework::TruthCollectionMaker::m_typeReadDecorKey {this, "inputClassifierParticleType", m_particlesKey, "classifierParticleType", "Particle type"}
private

Definition at line 75 of file TruthCollectionMaker.h.

The documentation for this class was generated from the following files:
xAOD::TruthVertex_v1::nOutgoingParticles
size_t nOutgoingParticles() const
Get the number of outgoing particles.
xAOD::TruthParticle_v1::setStatus
void setStatus(int value)
Set status code.
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
DerivationFramework::TruthCollectionMaker::m_daughterIDDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_daughterIDDecoratorKey
Definition: TruthCollectionMaker.h:67
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
SG::ReadHandle< xAOD::TruthParticleContainer >
xAOD::TruthParticle_v1::setE
void setE(float value)
Set the energy of the particle.
Definition: TruthParticle_v1.cxx:235
xAOD::TruthParticle_v1::setBarcode
void setBarcode(int value)
Set barcode.
DerivationFramework::TruthCollectionMaker::m_do_sherpa
Gaudi::Property< bool > m_do_sherpa
Definition: TruthCollectionMaker.h:50
DerivationFramework::TruthCollectionMaker::m_outcomeDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_outcomeDecoratorKey
Definition: TruthCollectionMaker.h:61
DerivationFramework::TruthCollectionMaker::m_do_compress
Gaudi::Property< bool > m_do_compress
Definition: TruthCollectionMaker.h:48
HepMC::INVALID_PARTICLE_ID
constexpr int INVALID_PARTICLE_ID
Definition: MagicNumbers.h:57
DerivationFramework::TruthCollectionMaker::m_typeReadDecorKey
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_typeReadDecorKey
Definition: TruthCollectionMaker.h:76
DerivationFramework::TruthCollectionMaker::m_particlesKey
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_particlesKey
Definition: TruthCollectionMaker.h:41
xAOD::TruthParticle_v1::setPx
void setPx(float value)
Set the x component of the particle's momentum.
SG::ConstAccessor< int >
DerivationFramework::TruthCollectionMaker::m_hadronOriginDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_hadronOriginDecoratorKey
Definition: TruthCollectionMaker.h:69
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
DerivationFramework::TruthCollectionMaker::m_npasspart
std::atomic< unsigned int > m_npasspart
Definition: TruthCollectionMaker.h:38
isSMLepton
bool isSMLepton(const T &p)
Definition: AtlasPID.h:183
xAOD::IParticle::FourMom_t
TLorentzVector FourMom_t
Definition of the 4-momentum type.
Definition: Event/xAOD/xAODBase/xAODBase/IParticle.h:69
MC::isPhysical
bool isPhysical(const T &p)
Identify if the particle is physical, i.e. is stable or decayed.
Definition: HepMCHelpers.h:51
xAOD::TruthParticle_v1::hasDecayVtx
bool hasDecayVtx() const
Check for a decay vertex on this particle.
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:210
xAOD::TruthParticle_v1::setM
void setM(float value)
Also store the mass.
Definition: TruthParticle_v1.cxx:241
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
SG::ReadDecorHandle
Handle class for reading a decoration on an object.
Definition: StoreGate/StoreGate/ReadDecorHandle.h:94
lumiFormat.i
int i
Definition: lumiFormat.py:85
SG::ReadDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
beamspotman.n
n
Definition: beamspotman.py:731
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
isZ
bool isZ(const T &p)
Definition: AtlasPID.h:352
xAOD::TruthParticle_v1
Class describing a truth particle in the MC record.
Definition: TruthParticle_v1.h:37
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:100
xAOD::TruthParticle
TruthParticle_v1 TruthParticle
Typedef to implementation.
Definition: Event/xAOD/xAODTruth/xAODTruth/TruthParticle.h:15
xAOD::TruthParticle_v1::hasProdVtx
bool hasProdVtx() const
Check for a production vertex on this particle.
Definition: TruthParticle_v1.cxx:74
xAOD::TruthParticle_v1::setPy
void setPy(float value)
Set the y component of the particle's momentum.
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
DerivationFramework::TruthCollectionMaker::m_classificationDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_classificationDecoratorKey
Definition: TruthCollectionMaker.h:63
xAOD::TruthVertex_v1::incomingParticle
const TruthParticle_v1 * incomingParticle(size_t index) const
Get one of the incoming particles.
Definition: TruthVertex_v1.cxx:69
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
DerivationFramework::TruthCollectionMaker::m_typeDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_typeDecoratorKey
Definition: TruthCollectionMaker.h:59
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
DerivationFramework::TruthCollectionMaker::m_originReadDecorKey
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_originReadDecorKey
Definition: TruthCollectionMaker.h:74
DerivationFramework::TruthCollectionMaker::m_ntotpart
std::atomic< unsigned int > m_ntotpart
Definition: TruthCollectionMaker.h:38
xAOD::TruthParticle_v1::decayVtx
const TruthVertex_v1 * decayVtx() const
The decay vertex of this particle.
xAOD::TruthParticle_v1::prodVtx
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Definition: TruthParticle_v1.cxx:80
DerivationFramework::TruthCollectionMaker::m_linkDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_linkDecoratorKey
Definition: TruthCollectionMaker.h:55
xAOD::TruthParticle_v1::setPdgId
void setPdgId(int pid)
Set PDG ID code.
DerivationFramework::TruthCollectionMaker::m_keep_navigation_info
Gaudi::Property< bool > m_keep_navigation_info
Definition: TruthCollectionMaker.h:52
DerivationFramework::TruthCollectionMaker::m_outputParticlesKey
SG::WriteHandleKey< xAOD::TruthParticleContainer > m_outputParticlesKey
Definition: TruthCollectionMaker.h:43
SG::WriteHandle
Definition: StoreGate/StoreGate/WriteHandle.h:73
isW
bool isW(const T &p)
Definition: AtlasPID.h:355
runIDAlign.accumulate
accumulate
Update flags based on parser line args.
Definition: runIDAlign.py:63
SG::WriteDecorHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
xAOD::TruthVertex_v1::nIncomingParticles
size_t nIncomingParticles() const
Get the number of incoming particles.
Definition: TruthVertex_v1.cxx:47
xAOD::TruthParticle_v1::setPz
void setPz(float value)
Set the z component of the particle's momentum.
DerivationFramework::TruthCollectionMaker::m_outcomeReadDecorKey
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_outcomeReadDecorKey
Definition: TruthCollectionMaker.h:78
DerivationFramework::TruthCollectionMaker::m_originDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_originDecoratorKey
Definition: TruthCollectionMaker.h:57
entries
double entries
Definition: listroot.cxx:49
DerivationFramework::TruthCollectionMaker::m_metaStore
ServiceHandle< StoreGateSvc > m_metaStore
Handle on the metadata store for init.
Definition: TruthCollectionMaker.h:82
DerivationFramework::TruthCollectionMaker::m_motherIDDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_motherIDDecoratorKey
Definition: TruthCollectionMaker.h:65
xAOD::TruthParticle_v1::pdgId
int pdgId() const
PDG ID code.
dqBeamSpot.nEntries
int nEntries
Definition: dqBeamSpot.py:73
xAOD::TruthVertex_v1::outgoingParticle
const TruthParticle_v1 * outgoingParticle(size_t index) const
Get one of the outgoing particles.
Definition: TruthVertex_v1.cxx:119
DerivationFramework::TruthCollectionMaker::m_classificationReadDecorKey
SG::ReadDecorHandleKey< xAOD::TruthParticleContainer > m_classificationReadDecorKey
Definition: TruthCollectionMaker.h:80
DerivationFramework::TruthCollectionMaker::m_partString
Gaudi::Property< std::string > m_partString
Definition: TruthCollectionMaker.h:46
Generated on Fri Apr 18 2025 22:02:05 for ATLAS Offline Software by   doxygen 1.8.18