fastsim::ParticleFilter Class Reference

(Kinematic) cuts on the particles that are propagated. More...

#include <ParticleFilter.h>

Public Member Functions
bool	accepts (const Particle &particle) const
	Check all if all criteria are fullfilled. More...
bool	acceptsEn (const Particle &particle) const
	Kinematic cuts on the particle. More...
bool	acceptsVtx (const math::XYZTLorentzVector &originVertexPosition) const
	Vertex within tracker volume. More...
	ParticleFilter (const edm::ParameterSet &cfg)
	Default Constructor. More...

Private Attributes
double	chargedPtMin2_
	Minimum pT^2 of a charged particle. More...
double	cos2ThetaMax_
	Particles must have abs(eta) < etaMax if close to beampipe. More...
double	EMin_
	Minimum energy of a particle. More...
double	protonEMin_
	Allow ALL protons with energy > protonEMin. More...
std::vector< int >	skipParticles_
	List of invisible particles (neutrinos are excluded by default) More...
double	vertexRMax2_
	Radius^2 of tracker volume. More...
double	vertexZMax_
	Z of tracker volume. More...

Detailed Description

(Kinematic) cuts on the particles that are propagated.

All other particles are skipped.

Definition at line 30 of file ParticleFilter.h.

Constructor & Destructor Documentation

fastsim::ParticleFilter::ParticleFilter

(

const edm::ParameterSet &

cfg

)

Default Constructor.

Definition at line 6 of file ParticleFilter.cc.

References particleFlowSimParticle_cfi::chargedPtMin, chargedPtMin2_, cos2ThetaMax_, EMin_, ALCARECOTkAlBeamHalo_cff::etaMax, edm::ParameterSet::getParameter(), protonEMin_, skipParticles_, vertexRMax2_, and vertexZMax_.

{
    // Charged particles must have pt greater than chargedPtMin [GeV]
    double chargedPtMin  = cfg.getParameter<double>("chargedPtMin");
    chargedPtMin2_  = chargedPtMin*chargedPtMin;
    
    // (All) particles must have energy greater than EMin [GeV]
    EMin_   = cfg.getParameter<double>("EMin");
    
    // Allow *ALL* protons with energy > protonEMin
    protonEMin_   = cfg.getParameter<double>("protonEMin");

    // List of invisible particles if extension needed
    // Predefined: Neutrinos, Neutralino_1
    skipParticles_   = cfg.getParameter<std::vector<int>>("invisibleParticles");
    
    // Particles must have abs(eta) < etaMax (if close enough to 0,0,0)
    double etaMax = cfg.getParameter<double>("etaMax");
    cos2ThetaMax_ = (vdt::fast_exp(2.*etaMax)-1.) / (vdt::fast_exp(2.*etaMax)+1.);
    cos2ThetaMax_ *= cos2ThetaMax_;

    // Particles must have vertex inside the tracker
    vertexRMax2_ = 129.0*129.0; 
    vertexZMax_ = 303.353;
}

Member Function Documentation

bool fastsim::ParticleFilter::accepts

(

const Particle &

particle

)

const

Check all if all criteria are fullfilled.

• Particle is invisible (neutrinos by default, list can be extended)
• Kinematic cuts (calls acceptsEN(...))
• Vertex within tracker volume (calls acceptsVtx(...))

  See also

  acceptsEn(const Particle & particle)

  acceptsVtx(const math::XYZTLorentzVector & originVertexPosition)

Definition at line 32 of file ParticleFilter.cc.

References funct::abs(), acceptsEn(), acceptsVtx(), cos2ThetaMax_, fastsim::Particle::momentum(), fastsim::Particle::pdgId(), fastsim::Particle::position(), protonEMin_, and skipParticles_.

Referenced by fastsim::ParticleManager::nextParticle().

{
    int absPdgId = abs(particle.pdgId());

    // skip invisible particles
    if(absPdgId == 12 || absPdgId == 14 || absPdgId == 16 || absPdgId == 1000022)
    {
       return false;
    }
    // keep all high-energy protons
    else if(absPdgId == 2212 && particle.momentum().E() >= protonEMin_)
    {
        return true;
    }
    
    // cut on eta if the origin vertex is close to the beam
    else if(particle.position().Perp2() < 25. && particle.momentum().Pz()*particle.momentum().Pz()/particle.momentum().P2() > cos2ThetaMax_)
    {
       return false;
    }

    // possible to extend list of invisible particles
    for(unsigned InvIdx = 0; InvIdx < skipParticles_.size(); InvIdx++){
        if(absPdgId == abs(skipParticles_.at(InvIdx)))
        {
            return false;
        }
    }

    // particles must have vertex in volume of tracker
    return acceptsVtx(particle.position()) && acceptsEn(particle);
    //return (acceptsVtx(particle.position()) || particle.momentum().Pz()*particle.momentum().Pz()/particle.momentum().P2() > (vdt::fast_exp(2.*3.0)-1.) / (vdt::fast_exp(2.*3.0)+1.)*(vdt::fast_exp(2.*3.0)-1.) / (vdt::fast_exp(2.*3.0)+1.)) && acceptsEn(particle);
}

bool fastsim::ParticleFilter::acceptsEn

(

const Particle &

particle

)

const

Kinematic cuts on the particle.

Definition at line 66 of file ParticleFilter.cc.

References funct::abs(), fastsim::Particle::charge(), chargedPtMin2_, EMin_, fastsim::Particle::momentum(), fastsim::Particle::pdgId(), and protonEMin_.

Referenced by accepts().

{
    int absPdgId = abs(particle.pdgId());

    // keep all high-energy protons
    if(absPdgId == 2212 && particle.momentum().E() >= protonEMin_)
    {
        return true;
    }
    
    // cut on the energy
    else if(particle.momentum().E() < EMin_)
    {
        return false;
    }
    
    // cut on pt of charged particles
    else if(particle.charge()!=0 && particle.momentum().Perp2()<chargedPtMin2_)
    {
        return false;
    }
    
    return true;
}

bool fastsim::ParticleFilter::acceptsVtx

(

const math::XYZTLorentzVector &

originVertexPosition

)

const

Vertex within tracker volume.

Parameters

originVertexPosition

Position of origin vertex.

Definition at line 92 of file ParticleFilter.cc.

References vertexRMax2_, and vertexZMax_.

Referenced by accepts(), and fastsim::ParticleManager::addSecondaries().

{
    // origin vertex is within the tracker volume
    return (originVertex.Perp2() < vertexRMax2_ && fabs(originVertex.Z()) < vertexZMax_);
}

Member Data Documentation

double fastsim::ParticleFilter::chargedPtMin2_

private

Minimum pT^2 of a charged particle.

Definition at line 56 of file ParticleFilter.h.

Referenced by acceptsEn(), and ParticleFilter().

double fastsim::ParticleFilter::cos2ThetaMax_

private

Particles must have abs(eta) < etaMax if close to beampipe.

Definition at line 59 of file ParticleFilter.h.

Referenced by accepts(), and ParticleFilter().

double fastsim::ParticleFilter::EMin_

private

Minimum energy of a particle.

Definition at line 57 of file ParticleFilter.h.

Referenced by acceptsEn(), and ParticleFilter().

double fastsim::ParticleFilter::protonEMin_

private

Allow ALL protons with energy > protonEMin.

Definition at line 58 of file ParticleFilter.h.

Referenced by accepts(), acceptsEn(), and ParticleFilter().

std::vector<int> fastsim::ParticleFilter::skipParticles_

private

List of invisible particles (neutrinos are excluded by default)

Definition at line 62 of file ParticleFilter.h.

Referenced by accepts(), and ParticleFilter().

double fastsim::ParticleFilter::vertexRMax2_

private

Radius^2 of tracker volume.

Definition at line 60 of file ParticleFilter.h.

Referenced by acceptsVtx(), and ParticleFilter().

double fastsim::ParticleFilter::vertexZMax_

private

Z of tracker volume.

Definition at line 61 of file ParticleFilter.h.

Referenced by acceptsVtx(), and ParticleFilter().