AnalysisConfigMT_Ntuple.cxx File Reference

#include "TrigInDetAnalysisExample/AnalysisConfigMT_Ntuple.h"
#include "xAODTruth/TruthParticleContainer.h"
#include "xAODTruth/TruthVertexContainer.h"
#include "xAODEventInfo/EventInfo.h"
#include "TrigAnalysisHelpers/FeatureRequestDescriptor.h"
#include "TrigInDetAnalysis/Filter_AcceptAll.h"
#include "TrigInDetAnalysisUtils/Filter_etaPT.h"
#include "TrigInDetAnalysisUtils/Filter_RoiSelector.h"
#include "TrigInDetAnalysisUtils/Filters.h"
#include "TrigInDetAnalysisUtils/TIDAVertexBuilder.h"
#include "TrigInDetAnalysis/TIDDirectory.h"
#include "TrigInDetAnalysisUtils/TIDARoiDescriptorBuilder.h"

Go to the source code of this file.

Functions
std::string	date ()
	sadly, includes a return at the end
HepMC::ConstGenParticlePtr	fromParent (int pdg_id, HepMC::ConstGenParticlePtr p, bool printout=false)
template<class T>
void	remove_duplicates (std::vector< T > &vec)

Detailed Description

Author

mark sutton

Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration

Definition in file AnalysisConfigMT_Ntuple.cxx.

Function Documentation

date()

std::string date

(

)

sadly, includes a return at the end

Definition at line 58 of file hcg.cxx.

                 { 
  time_t t;
  time(&t);
  std::string a = ctime(&t);
  std::string b = "";
  for ( unsigned i=0 ; i<a.size()-1 ; i++ ) b+=a[i];
  return b;
}

fromParent()

HepMC::ConstGenParticlePtr fromParent	(	int	pdg_id,
		HepMC::ConstGenParticlePtr	p,
		bool	printout = false )

recursive stopping conditions

recursive stopping conditions

recursive stopping conditions

Definition at line 78 of file AnalysisConfig_Ntuple.cxx.

                                                                                                   { 
 
  if ( p==0 ) return 0;
  if (std::abs(p->pdg_id())==11 || std::abs(p->pdg_id())==13 ) return 0; //don't want light leptons from tau decays
  if ( std::abs(p->pdg_id())==pdg_id ) return p;   
    
  auto vertex = p->production_vertex();
  if ( !vertex) return 0; // has no production vertex !!!
 
#ifdef HEPMC3
  if ( vertex->particles_in().size() < 1 ) return 0;  

  //  if ( printout ) { 
  //     TruthParticle t(p);
  //     std::cout << "particle " << *p << "  " << t.pdgId() << "\tparent " << p << std::endl;
  //  }
  
  for ( auto in: vertex->particles_in()) {
    auto parent = fromParent( pdg_id, in, printout );
    TruthParticle t(in);
    if ( parent && std::abs(parent->pdg_id())==pdg_id) { 
      return parent;
    }  
  }
#else  
  if ( vertex->particles_in_size() < 1 ) return 0;  
 
  HepMC::GenVertex::particles_in_const_iterator in  = vertex->particles_in_const_begin();
  HepMC::GenVertex::particles_in_const_iterator end = vertex->particles_in_const_end();
  while ( in!=end ) {
    const HepMC::GenParticle* parent = fromParent( pdg_id, *in, printout );
    TruthParticle t(*in);
    // if ( printout ) std::cout << "\tvalue for particle " << *in << "  " << t.pdgId() << "\tparent " << parent << std::endl;
    if ( parent && std::abs(parent->pdg_id())==pdg_id) { 
      //if ( printout ) std::cout << "found tau! - in parents" << std::endl; 
      return parent;
    }   
    in++;
  }
#endif
  
  return 0;
}

remove_duplicates()

template<class T>

void remove_duplicates

(

std::vector< T > &

vec

)

Definition at line 37 of file AnalysisConfigMT_Ntuple.cxx.

                                          {
  std::sort(vec.begin(), vec.end());
  vec.erase(std::unique(vec.begin(), vec.end()), vec.end());
}