Example #1
0
#!/usr/bin/env python

import os
import sys
import SDHCALSim as sim

import math

os.environ["SIMEXE"] = '/home/libo/workplace/SDHCALSim/bin/SDHCALSim'

for iTask in range(10, 20):
    params = sim.Params()
    params.physicsList = 'QGSP_BERT'
    params.nEvent = 2000
    params.seed = iTask * 100 + 12345
    params.outputFileName = 'pi_' + str(iTask)

    particle = sim.Particle()
    particle.particleName = 'pi-'
    particle.positionOption = 'gaus'
    particle.positionX = 0
    particle.positionY = 0
    particle.positionZ = -20
    particle.sigmaPos = 0

    particle.momentumOption = 'gaus'
    particle.momentumPhi = 0
    particle.momentumTheta = 0
    particle.sigmaMomentum = 0
    particle.energyDistribution = 'uniform'
    particle.minEnergy = 1
Example #2
0
import os
import sys
import SDHCALSim as sim

import math

if __name__ == '__main__':

    #os.environ["SIMEXE"] = '/path/to/SDHCALSim/bin/SDHCALSim'
    os.environ["SIMEXE"] = '/Users/libo/workplace/SDHCALSim/bin/SDHCALSim'

    particle = 'pi-'
    energy = 10

    params = sim.Params()
    params.physicsList = 'QGSP_BERT'
    params.nEvent = 3
    params.seed = 0
    params.outputFileName = 'SDHCALSim'

    params.oldConfig = False
    params.killNeutrons = False

    charged = sim.Particle()
    charged.particleName = particle
    charged.energy = energy

    charged.momentumOption = 'gaus'
    #charged.sigmaMomentum = 0.15
    charged.sigmaMomentum = 0.
Example #3
0
#!/usr/bin/env python

import os
import sys
import SDHCALSim

import math

if __name__ == '__main__':

    os.environ["SIMEXE"] = '/home/garillot/SDHCALSim2/bin/SDHCALSim'

    params = SDHCALSim.Params()

    params.particle = "pi-"
    params.physicsList = "FTF_BIC"

    #params.energy = sys.argv[1]
    params.energyDistribution = "uniform"
    params.minEnergy = 1
    params.maxEnergy = 20

    params.nEvent = 1000
    params.seed = 0

    params.momentumOption = "gaus"
    params.sigmaMomentum = 0.20

    params.positionOption = "uniform"
    params.positionX = 0
    params.positionY = 0
Example #4
0
#!/usr/bin/env python

import os
import sys
import SDHCALSim as sim

import math

os.environ[
    "SIMEXE"] = '/Users/libo/workplace/SDHCALSimTools_libo/SDHCALSim/bin/SDHCALSim'

for iEnergy in range(10, 100, 10):
    params = sim.Params()
    params.physicsList = 'QGSP_BERT'
    params.nEvent = 2
    params.seed = iEnergy * 100 + 654321
    params.outputFileName = 'pi_3_' + str(iEnergy) + 'GeV'

    particle = sim.Particle()
    particle.particleName = 'pi-'
    particle.positionOption = 'gaus'
    particle.positionX = 0
    particle.positionY = 0
    particle.positionZ = -20
    particle.sigmaPos = 0

    particle.momentumOption = 'gaus'
    particle.momentumPhi = 0
    particle.momentumTheta = 0
    particle.sigmaMomentum = 0
    particle.energy = iEnergy