00004 :
00005
00006 FLAVOR = process.generator.hscpFlavor.value()
00007 MASS_POINT = process.generator.massPoint.value()
00008 SLHA_FILE = process.generator.slhaFile.value()
00009 PROCESS_FILE = process.generator.processFile.value()
00010 PARTICLE_FILE = process.generator.particleFile.value()
00011
00012
00013 process.load("SimG4Core.CustomPhysics.CustomPhysics_cfi")
00014 process.customPhysicsSetup.particlesDef = PARTICLE_FILE
00015
00016 process.g4SimHits.Watchers = cms.VPSet (
00017 cms.PSet(
00018 type = cms.string('RHStopTracer'),
00019 RHStopTracer = cms.PSet(
00020 verbose = cms.untracked.bool (False),
00021 traceParticle = cms.string ("(~|tau1).*"),
00022 stopRegularParticles = cms.untracked.bool (False)
00023 )
00024 )
00025 )
00026
00027
00028 if FLAVOR=="gluino" or FLAVOR=="stop":
00029 process.customPhysicsSetup.processesDef = PROCESS_FILE
00030 process.g4SimHits.Physics = cms.PSet(
00031 process.customPhysicsSetup,
00032 DummyEMPhysics = cms.bool(True),
00033 G4BremsstrahlungThreshold = cms.double(0.5),
00034 DefaultCutValue = cms.double(1.),
00035 CutsPerRegion = cms.bool(True),
00036 Verbosity = cms.untracked.int32(0),
00037 type = cms.string('SimG4Core/Physics/CustomPhysics'),
00038 EMPhysics = cms.untracked.bool(True),
00039 HadPhysics = cms.untracked.bool(True),
00040 FlagBERT = cms.untracked.bool(False),
00041 FlagCHIPS = cms.untracked.bool(False),
00042 FlagFTF = cms.untracked.bool(False),
00043 FlagGlauber = cms.untracked.bool(False),
00044 FlagHP = cms.untracked.bool(False),
00045 GFlash = cms.PSet(
00046 GflashHistogram = cms.bool(False),
00047 GflashEMShowerModel = cms.bool(False),
00048 GflashHadronPhysics = cms.string('QGSP_BERT_EMV'),
00049 GflashHadronShowerModel = cms.bool(False)
00050 )
00051 )
00052 process.g4SimHits.G4Commands = cms.vstring('/tracking/verbose 1')
00053
00054 elif FLAVOR =="stau":
00055 process.g4SimHits.Physics = cms.PSet(
00056 process.customPhysicsSetup,
00057 DummyEMPhysics = cms.bool(True),
00058 G4BremsstrahlungThreshold = cms.double(0.5),
00059 DefaultCutValue = cms.double(1.),
00060 CutsPerRegion = cms.bool(True),
00061 Verbosity = cms.untracked.int32(0),
00062 type = cms.string('SimG4Core/Physics/CustomPhysics'),
00063 )
00064 process.g4SimHits.G4Commands = cms.vstring('/tracking/verbose 1')
00065
00066 else:
00067 print "Wrong flavor %s. Only accepted are gluino, stau, stop." % FLAVOR
00068
00069 return process
00070