####################################################33
##Now define the productions
pslic = ProductionJob()
pslic.setInputDataQuery(meta)
res = pslic.append(slic)
if not res['OK']:
print res['Message']
exit(1)
pslic.addFinalization(True, True, True, True)
pslic.setProdType("MCSimulation")
pslic.setProdGroup(prodGroup)
pslic.setWorkflowName(workflowName + '_sim_sid_cdr')
pslic.setCPUTime(cpuLimit)
pslic.setOutputSE("CERN-SRM")
pslic.setDescription('Simulating ' + workflowDescription)
res = pslic.createProduction()
if not res['OK']:
print res['Message']
exit(1)
res = pslic.finalizeProd()
if not res['OK']:
print res['Message']
exit(1)
meta = pslic.getMetadata() ##This is needed to link Sim and Rec
###Reconstruction w/o overlay
plcsim = ProductionJob()
####################################################33
##Now define the productions
pslic = ProductionJob()
pslic.setInputDataQuery(meta)
res = pslic.append(slic)
if not res['OK']:
print res['Message']
exit(1)
pslic.addFinalization(True,True,True,True)
pslic.setProdType("MCSimulation")
pslic.setProdGroup(prodGroup)
pslic.setWorkflowName(workflowName+'_sim_sid_cdr')
pslic.setCPUTime(cpuLimit)
pslic.setOutputSE("CERN-SRM")
pslic.setDescription('Simulating '+workflowDescription)
res = pslic.createProduction()
if not res['OK']:
print res['Message']
exit(1)
res = pslic.finalizeProd()
if not res['OK']:
print res['Message']
exit(1)
meta = pslic.getMetadata() ##This is needed to link Sim and Rec
###Reconstruction w/o overlay
psl.addMetadataToFinalFiles({"BeamParticle1":proddict['pname1'], "BeamParticle2":proddict['pname2'],
"EPA_B1":proddict['epa_b1'], "EPA_B2":proddict['epa_b2']})
res = psl.finalizeProd()
if not res['OK']:
print res['Message']
exit(1)
#As before: get the metadata for this production to input into the next
meta = psl.getMetadata()
if activesplit and meta:
#######################
## Split the input files.
psplit = ProductionJob()
psplit.setCPUTime(30000)
psplit.setLogLevel("verbose")
psplit.setProdType('Split')
psplit.setDestination("LCG.CERN.ch")
res = psplit.setInputDataQuery(meta)
if not res['OK']:
print res['Message']
exit(1)
psplit.setOutputSE("CERN-SRM")
wname = process+"_"+str(energy)+"_split"
wname += prod_name
psplit.setWorkflowName(wname)
psplit.setProdGroup(analysis+"_"+str(energy))
#Add the application
res = psplit.append(split)
