stoppedMass = float (iMass)
neutralinoMass = float (iNeutralinoMass)
decayFile = 'decays_%s_%s_%s.dat' % (flavor, iMass, iNeutralinoMass)
if flavor == "gluino1jet":
if not stoppedId: stoppedId=1092224
hscp.makeSingleJetRhadronDecay (decayFile, stoppedId, stoppedMass, neutralinoMass)
elif flavor == "gluino2jet":
if not stoppedId: stoppedId=1092224
hscp.makeDoubleJetRhadronDecay (decayFile, stoppedId, stoppedMass, neutralinoMass)
elif flavor == "stop":
if not stoppedId: stoppedId=1006223
hscp.makeStopHadronDecay (decayFile, stoppedId, stoppedMass, neutralinoMass)
elif flavor == "stau":
hscp.makeStauDecay (decayFile, stauMass, neutralinoMass)
else:
print 'Unknown flavor "%s". Expected: "gluino1jet", "gluino2jet", "stop", "stau"' % flavor
sys.exit(1)
firstRun = 1
if randomize: firstRun = randomize
print '**********************************************'
print 'Generating stopped %s decays' % flavor
print 'stopped particle:', stoppedId
print 'Total events:', nEvents
print 'Output file:', outputFile
print 'stopped mass:', iMass
print 'Neutralinoino mass:', neutralinoMass
print 'Global Tag:', globalTag
