if opt == "stoppedId": stoppedId = int(value)
if opt == "randomize": randomize = int (value)
if not flavor or not nEvents or not iMass or not iNeutralinoMass or not outputFile:
print '======> Missing parameters! <======'
print doc
sys.exit(1)
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
if opt == "globalTag": globalTag = value
if opt == "rhadronId": rhadronId = int(value)
if opt == "randomize": randomize = int (value)
if not nEvents or not iGluinoMass or not iNeutralinoMass or not outputFile:
print '======> Missing parameters! <======'
print doc
sys.exit(1)
gluinoMass = float (iGluinoMass)
neutralinoMass = float (iNeutralinoMass)
if not stopPointsFile:
stopPointsFile = 'StoppedHSCP/Simulation/data/stopped_rhadrons_gluino%s.txt' % iGluinoMass
decayFile = 'Rhadron_1jet_%s_%s.dat' % (iGluinoMass, iNeutralinoMass)
hscp.makeSingleJetRhadronDecay (decayFile, rhadronId, gluinoMass, neutralinoMass)
firstRun = 1
if randomize: firstRun = randomize
print '**********************************************'
print 'Generating stopped R-hadrons decays'
print 'R-hadron:', rhadronId
print 'Total events:', nEvents
print 'Output file:', outputFile
print 'Gluino mass:', gluinoMass
print 'Neutralinoino mass:', neutralinoMass
print 'Global Tag:', globalTag
print 'Stopped points:', stopPointsFile
print 'Decay description file:', decayFile
print 'randomize job:', randomize
print '**********************************************'
