ParticleManager.h

Go to the documentation of this file.

#ifndef FASTSIM_PARTICLEMANAGER_H
#define FASTSIM_PARTICLEMANAGER_H

#include "DataFormats/Math/interface/LorentzVector.h"
#include "HepMC/GenEvent.h"
#include <vector>
#include <memory>

#include "SimDataFormats/Track/interface/SimTrack.h"
#include "SimDataFormats/Vertex/interface/SimVertex.h"


// Author: L. Vanelderen, S. Kurz
// Date: 29 May 2017


namespace HepPDT
{
    class ParticleDataTable;
}

class RandomEngineAndDistribution;

namespace fastsim {
    class Particle;
    class ParticleFilter;
    class SimplifiedGeometry;


    class ParticleManager
    {
        public:

        ParticleManager(
            const HepMC::GenEvent & genEvent,
            const HepPDT::ParticleDataTable & particleDataTable,
            double beamPipeRadius,
            double deltaRchargedMother,
            const ParticleFilter & particleFilter,
            std::vector<SimTrack> & simTracks,
            std::vector<SimVertex> & simVertices,
            bool fixLongLivedBug);
        
        ~ParticleManager();


        std::unique_ptr<Particle> nextParticle(const RandomEngineAndDistribution & random);
        

        void addSecondaries(
            const math::XYZTLorentzVector & vertexPosition,
            int motherSimTrackId,
            std::vector<std::unique_ptr<Particle> > & secondaries,
            const SimplifiedGeometry * layer = nullptr);

        const SimVertex getSimVertex(unsigned i) { return simVertices_->at(i); }

        const SimTrack getSimTrack(unsigned i) { return simTracks_->at(i); }


        unsigned addEndVertex(const Particle * particle);

        private:

        unsigned addSimVertex(
            const math::XYZTLorentzVector & position,
            int motherIndex);
        

        unsigned addSimTrack(const Particle * particle);
        void exoticRelativesChecker(const HepMC::GenVertex* originVertex, int& hasExoticAssociation, int ngendepth);


        std::unique_ptr<Particle> nextGenParticle();

        // data members
        const HepMC::GenEvent * const genEvent_;  
        HepMC::GenEvent::particle_const_iterator genParticleIterator_;  
        const HepMC::GenEvent::particle_const_iterator genParticleEnd_;  
        int genParticleIndex_;  
        const HepPDT::ParticleDataTable * const particleDataTable_;  
        const double beamPipeRadius2_;  
        const double deltaRchargedMother_;  
        const ParticleFilter * const particleFilter_;  
        std::vector<SimTrack> * simTracks_;  
        std::vector<SimVertex> * simVertices_;  
        double momentumUnitConversionFactor_;  
        double lengthUnitConversionFactor_;  
        double lengthUnitConversionFactor2_;  
        double timeUnitConversionFactor_;  
        std::vector<std::unique_ptr<Particle> > particleBuffer_;  
        bool fixLongLivedBug_;
    };
}

inline bool isExotic(bool fixLongLivedBug, int pdgid_) {
  unsigned int pdgid = std::abs(pdgid_);
  return ((pdgid >= 1000000 && pdgid < 4000000 && pdgid != 3000022) ||  // SUSY, R-hadron, and technicolor particles
          pdgid == 17 ||                                                // 4th generation lepton
          pdgid == 34 ||                                                // W-prime
          pdgid == 37 ||                                                // charged Higgs  
          (pdgid == 39 && fixLongLivedBug));                           // bulk graviton
 
}

#endif