hephysics.hepsim

Class PromcUtil

java.lang.Object

hephysics.hepsim.PromcUtil

public class PromcUtil
extends Object

Utils classes to fill particles from ProMC data format of the HepSim repository with Monte Carlo predictions.

Constructor Summary

Constructors

Constructor and Description
PromcUtil()

Method Summary

Modifier and Type	Method and Description
static LParticle[]	getLParticleArray(ProMCHeaderFile.ProMCHeader header, ProMC.ProMCEvent.Particles pa, int status, double pTmin, double etaMax) Fill particles from ProMC and apply into array after some cuts if available.
static List<LParticle>	getLParticleList(ProMCHeaderFile.ProMCHeader header, ProMC.ProMCEvent.Particles pa, int status, double pTmin, double etaMax) Fill particles from ProMC and apply some cuts if available.
static HEParticle	getParticle(ProMC.ProMCEvent.Particles pa, int energy_unit, int length_unit, int j) For a given momentum unit, create a Lorentz particle object at the position j.
static List<ParticleD>	getParticleDList(double unit, proio.model.Mc.VarintPackedParticles pa, double pTmin, double etaMax) Fill particles from ProIO and apply some cuts if available.
static List<ParticleD>	getParticleDList(double unit, proio.model.Mc.VarintPackedParticles pa, int status, double pTmin, double etaMax) Fill particles from ProIO and apply some cuts if available.
static List<ParticleD>	getParticleDList(ProMCHeaderFile.ProMCHeader header, ProMC.ProMCEvent.Particles pa, double pTmin, double etaMax) Fill particles from ProMC from a typical NLO program, and apply some cuts if available.
static List<ParticleD>	getParticleDList(ProMCHeaderFile.ProMCHeader header, ProMC.ProMCEvent.Particles pa, int status, double pTmin, double etaMax) Fill particles from ProMC and apply some cuts if available.
static List<ParticleF>	getParticleFList(ProMCHeaderFile.ProMCHeader header, ProMC.ProMCEvent.Particles pa, double pTmin, double etaMax) Fill particles from ProMC from a typical NLO program, and apply some cuts if available.
static List<ParticleF>	getParticleFList(ProMCHeaderFile.ProMCHeader header, ProMC.ProMCEvent.Particles pa, int status, double pTmin, double etaMax) Fill particles from ProMC and apply some cuts if available.
static LParticle	getPosition(ProMC.ProMCEvent.Particles pa, int unit, int j) For a given space/time unit, create a particle with postions X, Y,Z,Time.

Methods inherited from class java.lang.Object

equals, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Constructor Detail

PromcUtil

public PromcUtil()

Method Detail

getParticle

public static HEParticle getParticle(ProMC.ProMCEvent.Particles pa,
                                     int energy_unit,
                                     int length_unit,
                                     int j)

For a given momentum unit, create a Lorentz particle object at the position j. It keeps 4-momentum, postion, PDG ID, status etc.

Parameters:

pa - particle record from a ProMC file

energy_unit - momentum / energy unit used to create a 4-momentum from a varint (int64)

length_unit - length unit used to create (x,y,z,t) from a varint (int64)

j - position of the particle in the ProMC record

Returns:

a HEP Lorentz particle of the record.

getPosition

public static LParticle getPosition(ProMC.ProMCEvent.Particles pa,
                                    int unit,
                                    int j)

For a given space/time unit, create a particle with postions X, Y,Z,Time.

Parameters:

pa - particle record from a ProMC file.

unit - space/time lenght unit used to convert a varint to a value.

j - position of particle in the record

Returns:

particle position

getLParticleList

public static List<LParticle> getLParticleList(ProMCHeaderFile.ProMCHeader header,
                                               ProMC.ProMCEvent.Particles pa,
                                               int status,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProMC and apply some cuts if available. Also you can apply some basic cuts.

Parameters:

header - header from ProMC file

pa - particle record from ProMC file

status - status code for particles to be accepted. For final state, use status==1

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getLParticleArray

public static LParticle[] getLParticleArray(ProMCHeaderFile.ProMCHeader header,
                                            ProMC.ProMCEvent.Particles pa,
                                            int status,
                                            double pTmin,
                                            double etaMax)

Fill particles from ProMC and apply into array after some cuts if available. Also you can apply some basic cuts.

Parameters:

header - header from ProMC file

pa - particle record from ProMC file

status - status code for particles to be accepted. For final state, use status==1

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getParticleFList

public static List<ParticleF> getParticleFList(ProMCHeaderFile.ProMCHeader header,
                                               ProMC.ProMCEvent.Particles pa,
                                               int status,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProMC and apply some cuts if available. Use float precisions for particle definitions. Also you can apply some basic cuts.

Parameters:

header - header from ProMC file

pa - particle record from ProMC file

status - status code for particles to be accepted. For final state, use status==1

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getParticleFList

public static List<ParticleF> getParticleFList(ProMCHeaderFile.ProMCHeader header,
                                               ProMC.ProMCEvent.Particles pa,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProMC from a typical NLO program, and apply some cuts if available. Use float precisions for particle definitions. Also you can apply some basic cuts.

Parameters:

header - header from ProMC file

pa - particle record from ProMC file

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getParticleDList

public static List<ParticleD> getParticleDList(ProMCHeaderFile.ProMCHeader header,
                                               ProMC.ProMCEvent.Particles pa,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProMC from a typical NLO program, and apply some cuts if available. Use double precisions for particle definitions. Also you can apply some basic cuts.

Parameters:

header - header from ProMC file

pa - particle record from ProMC file

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getParticleDList

public static List<ParticleD> getParticleDList(double unit,
                                               proio.model.Mc.VarintPackedParticles pa,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProIO and apply some cuts if available. Use double precisions for particle definitions. Also you can apply some basic cuts.

Parameters:

unit - energy unit used to convert varint type of px,py,pz,e to floats

pa - particle record from ProIO file

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getParticleDList

public static List<ParticleD> getParticleDList(ProMCHeaderFile.ProMCHeader header,
                                               ProMC.ProMCEvent.Particles pa,
                                               int status,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProMC and apply some cuts if available. Use double precision particle definition. Also you can apply some basic cuts.

Parameters:

header - header from ProMC file

pa - particle record from ProMC file

status - status code for particles to be accepted. For final state, use status==1

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.

getParticleDList

public static List<ParticleD> getParticleDList(double unit,
                                               proio.model.Mc.VarintPackedParticles pa,
                                               int status,
                                               double pTmin,
                                               double etaMax)

Fill particles from ProIO and apply some cuts if available. Use double precision particle definition. Also you can apply some basic cuts.

Parameters:

unit - energy unit for varint to float conversions

pa - particle record from ProIO file

status - status code for particles to be accepted. For final state, use status==1

pTmin - minimum pT cut for filled particles. If pTmin less than, do not apply any pT cut

etaMax - maximum Eta (pseudorapidity) cut on filled particles. If etaMax greater than 1000, no any Eta cut is applied.