In This Package:

MCTruthManager.h

Go to the documentation of this file.

// $Id: MCTruthManager.h,v 1.5 2008/10/20 08:23:11 robbep Exp $
#ifndef COMPONENTS_MCTRUTHMANAGER_H 
#define COMPONENTS_MCTRUTHMANAGER_H 1

#include <utility>
// Include files
//#include "Event/HepMCEvent.h"
#include "HepMC/GenEvent.h"

namespace LHCb {
  class MCParticle ;
};

class MCTruthManager 
{
public: 

  // Offset value for intermediate dummy and split particles
  enum MCTruthBarCode {
    DummyBarCode =  20000000,
    SplitBarCode = 100000000       
  };
  
  enum MCTruthPDGID {
    DummyPDGID = -99000000 
  };

  virtual ~MCTruthManager( ); 

  static MCTruthManager* GetInstance();

  void NewEvent();
  HepMC::GenEvent* GetCurrentEvent() {
    return event;
  };

  void PrintEvent();
 
  void AddParticle(HepMC::FourVector& fourMom, 
                   HepMC::FourVector& prodPos, HepMC::FourVector& endPos, 
                   int pdg, int trID, int parentID, bool directParent,
                   int creatorID, LHCb::MCParticle * motherMCP , 
                   bool oscillated = false );
 
  const std::vector<int>& GetPrimaryBarcodes();

  int GetCreatorID(int barcode);

  const std::vector<int>& GetOscillatedBarcodes();

  LHCb::MCParticle * GetMotherMCParticle( const int barcode ) ;

protected:
  MCTruthManager( );


private:
  HepMC::GenEvent* event;

  // vector containing barcodes of primary particles (not having any mother)
  std::vector<int> primarybarcodes;

  // map containing barcode of vertices and process identifier
  std::map<int,int> creators;

  // vector containing bar code of oscillated particles
  std::vector<int> oscillated;

  // map containing number of 'segmentations' for each particle
  // (i.e. number of additional vertices introduced in order to attach secondary particles
  // which were created 'in-flight', for instance bremstrahlung gammas, etc)
  std::map<int,int> segmentations;

  std::map< int , LHCb::MCParticle * > m_mcparticles ;
  
  // recursive printing of the tree
  void printTree(HepMC::GenParticle*,std::string);

};
#endif // COMPONENTS_MCTRUTHMANAGER_H

---