Esempio n. 1
0
def createProductionStep(name, type, inputQuery=None, outputQuery=None):
    # create a production step
    prodStep = ProductionStep()
    prodStep.Name = name
    prodStep.Type = type
    prodStep.Inputquery = inputQuery
    prodStep.Outputquery = outputQuery
    return prodStep
Esempio n. 2
0
def build_evndisp_step(DL0_data_set, nsb=1):
    ''' Define a new EventDisplay analysis production step

    @return ProductionStep object
    '''
    if nsb == 1:
        DIRAC.gLogger.notice('NSB1x Analysis')
        DL0_data_set_NSB = DL0_data_set
    elif nsb == 5:
        DIRAC.gLogger.notice('NSB5x Analysis')
        DL0_data_set_NSB = DL0_data_set.replace('NSB1x', 'NSB5x')

    prod_step_2 = ProductionStep()
    prod_step_2.Name = 'Analysis_'+DL0_data_set_NSB.replace('AdvancedBaseline_', '').replace('DL0', 'DL1')
    prod_step_2.Type = 'DataReprocessing'  # This corresponds to the Transformation Type
    prod_step_2.Inputquery = get_dataset_MQ(DL0_data_set_NSB)
    prod_step_2.Outputquery = get_dataset_MQ(DL0_data_set_NSB.replace('DL0', 'DL1'))

    # Here define the job description to be associated to the second ProductionStep
    job2 = EvnDispProd5Job()
    job2.setName('Prod5_EvnDisp')
    # output
    job2.setOutputSandbox(['*Log.txt'])
    # refine output meta data if needed
    output_meta_data = copy(prod_step_2.Outputquery)
    job2.set_meta_data(output_meta_data)
    job2.set_file_meta_data(nsb=output_meta_data['nsb']['='])

    # check if La Palma else use default that is Paranal
    if output_meta_data['site'] == 'LaPalma':
        job2.prefix = "CTA.prod5N"
        job2.layout_list = 'BL-0LSTs05MSTs-MSTF BL-0LSTs05MSTs-MSTN \
                            BL-4LSTs00MSTs-MSTN BL-4LSTs05MSTs-MSTF \
                            BL-4LSTs05MSTs-MSTN BL-4LSTs09MSTs-MSTF \
                            BL-4LSTs09MSTs-MSTN BL-4LSTs15MSTs-MSTF \
                            BL-4LSTs15MSTs-MSTN'
        DIRAC.gLogger.notice('LaPalma layouts:\n',job2.layout_list.split())
    elif output_meta_data['site'] == 'Paranal':
        DIRAC.gLogger.notice('Paranal layouts:\n',job2.layout_list.split())

    job2.ts_task_id = '@{JOB_ID}'  # dynamic
    job2.setupWorkflow(debug=False)
    job2.setType('EvnDisp3')  # mandatory *here*
    prod_step_2.Body = job2.workflow.toXML()
    prod_step_2.GroupSize = 5

    return prod_step_2
Esempio n. 3
0
def build_evndisp_step(DL0_data_set, nsb=1, name_tag=''):
    ''' Define a new EventDisplay analysis production step

    @return ProductionStep object
    '''
    if nsb == 1:
        DIRAC.gLogger.notice('NSB1x Analysis')
        DL0_data_set_NSB = DL0_data_set
    elif nsb == 5:
        DIRAC.gLogger.notice('NSB5x Analysis')
        DL0_data_set_NSB = DL0_data_set.replace('NSB1x', 'NSB5x')

    prod_step_2 = ProductionStep()
    prod_step_2.Name = 'Analysis_'+DL0_data_set_NSB.replace('AdvancedBaseline_', '').replace('DL0', 'DL1')
    prod_step_2.Name += '%s' % name_tag
    prod_step_2.Type = 'DataReprocessing'  # This corresponds to the Transformation Type
    prod_step_2.Inputquery = get_dataset_MQ(DL0_data_set_NSB)
    prod_step_2.Outputquery = get_dataset_MQ(DL0_data_set_NSB.replace('DL0', 'DL1'))

    # Here define the job description to be associated to the second ProductionStep
    job2 = EvnDispProd5SingJob(cpuTime=259200.)
    job2.version = 'eventdisplay-cta-dl1-prod5.v03'
    job2.setName('Prod5b_EvnDisp_Singularity')
    # output
    job2.setOutputSandbox(['*Log.txt'])
    # refine output meta data if needed
    output_meta_data = copy(prod_step_2.Outputquery)
    job2.set_meta_data(output_meta_data)
    job2.set_file_meta_data(nsb=output_meta_data['nsb']['='])

    job2.ts_task_id = '@{JOB_ID}'  # dynamic
    job2.group_size = 5   # for the input files verification
    job2.setupWorkflow(debug=False)
    job2.setType('EvnDisp3')  # mandatory *here*
    prod_step_2.Body = job2.workflow.toXML()
    prod_step_2.GroupSize = 5  # have to match the above group size?
    # return ProductionStep object
    return prod_step_2
Esempio n. 4
0
    def test_SeqProduction(self):

        # Define the first step of the production
        prodStep1 = ProductionStep()
        prodStep1.Name = 'Sim_prog'
        prodStep1.Type = 'MCSimulation'
        outputquery = {
            'zenith': {
                'in': [20, 40]
            },
            'particle': 'gamma',
            'tel_sim_prog': 'simtel',
            'outputType': {
                'in': ['Data', 'Log']
            }
        }
        prodStep1.Outputquery = outputquery

        # Add the step to the production
        res = self.prodClient.addProductionStep(prodStep1)
        self.assertTrue(res['OK'])

        # Define the second step of the production
        prodStep2 = ProductionStep()
        prodStep2.Name = 'Reco_prog'
        prodStep2.Type = 'DataProcessing'
        prodStep2.ParentStep = prodStep1

        inputquery = {
            'zenith': 20,
            'particle': 'gamma',
            'tel_sim_prog': 'simtel',
            'outputType': 'Data'
        }
        outputquery = {
            'zenith': 20,
            'particle': 'gamma',
            'analysis_prog': 'evndisp',
            'data_level': 1,
            'outputType': {
                'in': ['Data', 'Log']
            }
        }

        prodStep2.Inputquery = inputquery
        prodStep2.Outputquery = outputquery

        # Add the step to the production
        res = self.prodClient.addProductionStep(prodStep2)
        self.assertTrue(res['OK'])

        # Define the third step of the production
        prodStep3 = ProductionStep()
        prodStep3.Name = 'Analyis_prog'
        prodStep3.Type = 'DataProcessing'
        prodStep3.ParentStep = prodStep2

        inputquery = {
            'zenith': 20,
            'particle': 'gamma',
            'analysis_prog': 'evndisp',
            'data_level': 1,
            'outputType': 'Data'
        }
        outputquery = {
            'zenith': 20,
            'particle': 'gamma',
            'analysis_prog': 'evndisp',
            'data_level': 2,
            'outputType': {
                'in': ['Data', 'Log']
            }
        }

        prodStep3.Inputquery = inputquery
        prodStep3.Outputquery = outputquery

        # Add the step to the production
        res = self.prodClient.addProductionStep(prodStep3)
        self.assertTrue(res['OK'])

        # Get the production description
        prodDesc = self.prodClient.prodDescription

        # Create the production
        prodName = 'SeqProd'
        res = self.prodClient.addProduction(prodName, json.dumps(prodDesc))
        self.assertTrue(res['OK'])

        # Start the production, i.e. instatiate the transformation steps
        res = self.prodClient.startProduction(prodName)
        self.assertTrue(res['OK'])

        # Get the transformations of the production
        res = self.prodClient.getProduction(prodName)
        self.assertTrue(res['OK'])
        prodID = res['Value']['ProductionID']

        res = self.prodClient.getProductionTransformations(prodID)
        self.assertTrue(res['OK'])
        self.assertEqual(len(res['Value']), 3)

        # Delete the production
        res = self.prodClient.deleteProduction(prodName)
        self.assertTrue(res['OK'])
Esempio n. 5
0
    def test_SeqProduction(self):

        # Define the first step of the production
        prodStep1 = ProductionStep()
        prodStep1.Name = "Sim_prog"
        prodStep1.Type = "MCSimulation"
        outputquery = {
            "zenith": {
                "in": [20, 40]
            },
            "particle": "gamma",
            "tel_sim_prog": "simtel",
            "outputType": {
                "in": ["Data", "Log"]
            },
        }
        prodStep1.Outputquery = outputquery

        # Add the step to the production
        res = self.prodClient.addProductionStep(prodStep1)
        self.assertTrue(res["OK"])

        # Define the second step of the production
        prodStep2 = ProductionStep()
        prodStep2.Name = "Reco_prog"
        prodStep2.Type = "DataProcessing"
        prodStep2.ParentStep = prodStep1

        inputquery = {
            "zenith": 20,
            "particle": "gamma",
            "tel_sim_prog": "simtel",
            "outputType": "Data"
        }
        outputquery = {
            "zenith": 20,
            "particle": "gamma",
            "analysis_prog": "evndisp",
            "data_level": 1,
            "outputType": {
                "in": ["Data", "Log"]
            },
        }

        prodStep2.Inputquery = inputquery
        prodStep2.Outputquery = outputquery

        # Add the step to the production
        res = self.prodClient.addProductionStep(prodStep2)
        self.assertTrue(res["OK"])

        # Define the third step of the production
        prodStep3 = ProductionStep()
        prodStep3.Name = "Analyis_prog"
        prodStep3.Type = "DataProcessing"
        prodStep3.ParentStep = prodStep2

        inputquery = {
            "zenith": 20,
            "particle": "gamma",
            "analysis_prog": "evndisp",
            "data_level": 1,
            "outputType": "Data",
        }
        outputquery = {
            "zenith": 20,
            "particle": "gamma",
            "analysis_prog": "evndisp",
            "data_level": 2,
            "outputType": {
                "in": ["Data", "Log"]
            },
        }

        prodStep3.Inputquery = inputquery
        prodStep3.Outputquery = outputquery

        # Add the step to the production
        res = self.prodClient.addProductionStep(prodStep3)
        self.assertTrue(res["OK"])

        # Get the production description
        prodDesc = self.prodClient.prodDescription

        # Create the production
        prodName = f"SeqProd{random.randint(0, 100_000)}"
        res = self.prodClient.addProduction(prodName, json.dumps(prodDesc))
        self.assertTrue(res["OK"])

        # Start the production, i.e. instatiate the transformation steps
        res = self.prodClient.startProduction(prodName)
        self.assertTrue(res["OK"])

        # Get the transformations of the production
        res = self.prodClient.getProduction(prodName)
        self.assertTrue(res["OK"])
        prodID = res["Value"]["ProductionID"]

        res = self.prodClient.getProductionTransformations(prodID)
        self.assertTrue(res["OK"])
        self.assertEqual(len(res["Value"]), 3)

        # Delete the production
        res = self.prodClient.deleteProduction(prodName)
        self.assertTrue(res["OK"])