AnalysisConfig_Ntuple.cxx File Reference

#include <cstdio>
#include <sys/time.h>
#include "McParticleEvent/TruthParticleContainer.h"
#include "GeneratorObjects/McEventCollection.h"
#include "AtlasHepMC/GenEvent.h"
#include "AtlasHepMC/GenVertex.h"
#include "AtlasHepMC/GenParticle.h"
#include "EventInfo/EventInfo.h"
#include "EventInfo/EventID.h"
#include "xAODEventInfo/EventInfo.h"
#include "TrigInDetAnalysis/TIDDirectory.h"
#include "TrigInDetAnalysisUtils/TIDARoiDescriptorBuilder.h"
#include "TrigInDetAnalysis/Filter_AcceptAll.h"
#include "TrigInDetAnalysisUtils/Filter_etaPT.h"
#include "TrigInDetAnalysisUtils/Filter_RoiSelector.h"
#include "TrigInDetAnalysisUtils/Filters.h"
#include "TrigInDetAnalysisExample/AnalysisConfig_Ntuple.h"
#include "TrkParameters/TrackParameters.h"
#include "TrkTrack/TrackCollection.h"
#include "TrkTrack/Track.h"
#include "VxVertex/VxContainer.h"
#include "muonEvent/MuonContainer.h"
#include "egammaEvent/ElectronContainer.h"
#include "tauEvent/TauJetContainer.h"
#include "TrigSteeringEvent/HLTResult.h"
#include "TrigDecisionTool/ExpertMethods.h"
#include "TrigSteeringEvent/TrigRoiDescriptorCollection.h"
#include "xAODTracking/TrackParticle.h"
#include "xAODTracking/TrackParticleContainer.h"

Go to the source code of this file.

Macros
#define	endmsg   endmsg

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-2025 CERN for the benefit of the ATLAS collaboration

Definition in file AnalysisConfig_Ntuple.cxx.

Macro Definition Documentation

endmsg

#define endmsg   endmsg

Definition at line 63 of file AnalysisConfig_Ntuple.cxx.

Function Documentation

date()

std::string date

(

)

sadly, includes a return at the end

Definition at line 67 of file AnalysisConfig_Ntuple.cxx.

                 { 
  time_t t;
  time(&t);
  char buf[26];
  std::string mtime = ctime_r(&t, buf);
  mtime.erase( std::remove(mtime.begin(), mtime.end(), '\n'), mtime.end() );
  return mtime;
}

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 126 of file AnalysisConfig_Ntuple.cxx.

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