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 print '**********************************************' print 'Generating stopped %s decays' % flavor print 'stopped particle:', stoppedId