def getStdhepSplit(): """ Get some stdhep split """ from ILCDIRAC.Interfaces.API.NewInterface.Applications import StdHepSplit stdhepsplit = StdHepSplit() stdhepsplit.setVersion("V2") stdhepsplit.setNumberOfEventsPerFile(5) stdhepsplit.setOutputFile("teststdhepsplit.stdhep") return stdhepsplit
def getStdhepSplit(): """ Get some stdhep split """ from ILCDIRAC.Interfaces.API.NewInterface.Applications import StdHepSplit stdhepsplit = StdHepSplit() stdhepsplit.setVersion("V2") stdhepsplit.setNumberOfEventsPerFile(5) stdhepsplit.setOutputFile("teststdhepsplit.stdhep") stdhepsplit.setMaxRead(10) return stdhepsplit
def createSplitApplication( eventsPerJob, eventsPerBaseFile, splitType='stdhep' ): """ create Split application """ from ILCDIRAC.Interfaces.API.NewInterface.Applications import StdHepSplit, SLCIOSplit if splitType.lower() == 'stdhep': stdhepsplit = StdHepSplit() stdhepsplit.setVersion("V3") stdhepsplit.setNumberOfEventsPerFile( eventsPerJob ) stdhepsplit.datatype = 'gen' stdhepsplit.setMaxRead( eventsPerBaseFile ) return stdhepsplit if splitType.lower() == 'lcio': split = SLCIOSplit() split.setNumberOfEventsPerFile( eventsPerJob ) return stdhepsplit raise NotImplementedError( 'unknown splitType: %s ' % splitType )
def createSplitApplication( eventsPerJob, eventsPerBaseFile, splitType='stdhep' ): """ create Split application """ from ILCDIRAC.Interfaces.API.NewInterface.Applications import StdHepSplit, SLCIOSplit if splitType.lower() == 'stdhep': stdhepsplit = StdHepSplit() stdhepsplit.setVersion("V3") stdhepsplit.setNumberOfEventsPerFile( eventsPerJob ) stdhepsplit.datatype = 'gen' stdhepsplit.setMaxRead( eventsPerBaseFile ) return stdhepsplit if splitType.lower() == 'lcio': split = SLCIOSplit() split.setNumberOfEventsPerFile( eventsPerJob ) return stdhepsplit raise NotImplementedError( 'unknown splitType: %s ' % splitType )
#DoSplit activesplit = False nbevtsperfile = 200 #Do Reco sid_rec = True #Do Reco with Overlay sid_rec_ov = False n_events = 100 #How many events per job. (Is also obtained from the FC in principle) ##Split stdhepsplit = StdHepSplit() stdhepsplit.setVersion("V2") stdhepsplit.setNumberOfEventsPerFile(nbevtsperfilestdhep) ##Simulation SID slic = SLIC() slic.setVersion('v2r9p8')#This will change only once slic.setSteeringFile('defaultClicCrossingAngle.mac')#This will change only once slic.setDetectorModel(detector_model) slic.setNumberOfEvents(n_events) ##Split split = SLCIOSplit() split.setNumberOfEventsPerFile(nbevtsperfile) ## SID Reco w/o overlay lcsim_prepandora = LCSIM()
if res in [ 'n', 'N' ]: sys.exit(0) # Do Split activesplitstdhep = True # Do Sim sid_sim = True # Do Replication replicateFiles = True ## Split stdhepsplit = StdHepSplit() stdhepsplit.setVersion( stdhepSplitVersion ) stdhepsplit.setNumberOfEventsPerFile( nEvtsPerJob ) ## Simulation slic = SLIC() slic.setVersion( slicVersion ) slic.setSteeringFile( slicMacro ) slic.setDetectorModel( detectorModel ) slic.setNumberOfEvents( nEvtsPerJob ) ############################################ #What is below WILL NEVER NEED TO BE TOUCHED #(I'm not kidding, if you touch and break, not my problem) # # Define production step splitting the stdhep files if activesplitstdhep and meta: pstdhepsplit = SIDProductionJob()
print "matchToInput_mokka (input) = " + matchToInput_mokka print "matchToInput_marlin (input) = " + matchToInput_marlin print "Basepath for sim&rec (output) = " + basepath print "Diskpath for stdhepsplit&dst (output) = " + diskpath print "Outtput SE : stdhepsplit(%s), SUM(%s), psplit(%s), REC&DST(%s)" % (SE_stdhepsplit, SE_sim, SE_psplit, SE_rec) print "####################################################################################" ############################################################################################## ## [PART2] Below is not to be touched ############################################################################################## ###### Whatever is below is not to be touched... Or at least only when something changes ##Split stdhepsplit = StdHepSplit() stdhepsplit.setVersion("V2") stdhepsplit.setNumberOfEventsPerFile(nbevtsperfilestdhep) ##Simulation ILD mo = Mokka() mo.setVersion(MokkaVer) ###SET HERE YOUR MOKKA VERSION, the software will come from the ILDConfig mo.setDetectorModel(detectorModel) mo.setSteeringFile("bbudsc_3evt.steer") ### Do not include '.tgz' mo.setDbSlice(dbslice) ##Simulation ILD ddsim = None if UseDD4hepGeometry: ddsim = DDSim() ddsim.setVersion(DDSimVer) ###SET HERE YOUR MOKKA VERSION, the software will come from the ILDConfig ddsim.setDetectorModel(detectorModel)
if res in ['n', 'N']: sys.exit(0) # Do Split activesplitstdhep = True # Do Sim sid_sim = True # Do Replication replicateFiles = True ## Split stdhepsplit = StdHepSplit() stdhepsplit.setVersion(stdhepSplitVersion) stdhepsplit.setNumberOfEventsPerFile(nEvtsPerJob) ## Simulation slic = SLIC() slic.setVersion(slicVersion) slic.setSteeringFile(slicMacro) slic.setDetectorModel(detectorModel) slic.setNumberOfEvents(nEvtsPerJob) ############################################ #What is below WILL NEVER NEED TO BE TOUCHED #(I'm not kidding, if you touch and break, not my problem) # # Define production step splitting the stdhep files if activesplitstdhep and meta: pstdhepsplit = SIDProductionJob()