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 60 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

useful debug

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 !!!
 
  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;
    }  
  }
  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());
}