Uploaded image for project: 'Data Management'
  1. Data Management
  2. DM-21481

Validate gains from AuxTel sensor

    Details

    • Type: Story
    • Status: Done
    • Resolution: Done
    • Fix Version/s: None
    • Component/s: None
    • Labels:
      None

      Description

      Use the gains from ticket DM-18052, and use them to correct a flatfield (or whatever we have that's closest to one in the configuration that this data was taken in) and check to see how good the continuity across amp boundaries. 

        Attachments

          Issue Links

            Activity

            Hide
            plazas Andrés Alejandro Plazas Malagón added a comment -

            Using (lsst_distrib 19.0.0+65 w_2020_10 current setup):

            measurePhotonTransferCurve.py /project/shared/auxTel/ --rerun plazas/auxtelPTC2020MAR10 --id detector=0 --visit-pairs 2020022100013,2020022100014 2020022100015,2020022100016 2020022100017,2020022100018 2020022100019,2020022100020 2020022100021,2020022100022 2020022100023,2020022100024 2020022100025,2020022100026 2020022100027,2020022100028 2020022100029,2020022100030 2020022100031,2020022100032 2020022100033,2020022100034 2020022100035,2020022100036 -c makePlots=True ptcFitType=POLYNOMIAL polynomialFitDegree=3 -c isr.doBias=False isr.doDark=False isr.doSuspect=False isr.doLinearize=False

            The following gains were produced:

            gain_pickle_file = "/project/shared/auxTel/rerun/plazas/auxtelPTC2020MAR10/calibrations/ptc/ptcDataset-det000.pkl"
            gain_file = open(gain_pickle_file, 'rb')
            gain_data = pkl.load(gain_file)
            print (gain_data.gain)

            {'C10': 1.0348085853531022, 'C11': 1.02347459717329, 'C12': 1.0378724084326139, 'C13': 1.0358289608332496, 'C14': 1.012234515720022, 'C15': 1.0060929875013593, 'C16': 1.037898801107405, 'C17': 1.081185663200865, 'C07': 1.1520424164621434, 'C06': 1.1012153572252976, 'C05': 1.113872220848673, 'C04': 1.0939783643655419, 'C03': 1.0698829435732895, 'C02': 1.0906280187126696, 'C01': 1.104278016310778, 'C00': 1.1038149461060722}

            The, i applied the gains to image "dataId = dict(dayObs="2020-02-20", seqnum=433)" (expId 2020022000433 in '/project/shared/auxTel'):

            code (modified versioon of coode by Craig Lage: /project/cslage/ComCam/notebooks/Ingest_ComCam_06Mar20.ipynb):

            bias = butler.get('bias', dataId)

            isrConfig = IsrTask.ConfigClass()
            isrConfig.doLinearize = False
            isrConfig.doBias = True
            isrConfig.doApplyGains = True
            isrConfig.doFlat = False
            isrConfig.doDark = False
            isrConfig.doFringe = False
            isrConfig.doDefect = False
            isrConfig.doAddDistortionModel = False
            isrConfig.doWrite = False
            isrTask = IsrTask(config=isrConfig)

            gain_pickle_file = "/project/shared/auxTel/rerun/plazas/auxtelPTC2020MAR10/calibrations/ptc/ptcDataset-det000.pkl"
            gain_file = open(gain_pickle_file, 'rb')
            gain_data = pkl.load(gain_file).gain
            raw = butler.get('raw', dataId)
            ccd = raw.getDetector()
            builder = ccd.rebuild()
            for amp in builder:
            print (amp.getName())
            amp.setGain(gain_data[amp.getName()])
            raw.setDetector(builder.finish())
            ccd = raw.getDetector()
            postIsr = isrTask.run(raw, bias=bias).exposure

            Show
            plazas Andrés Alejandro Plazas Malagón added a comment - Using (lsst_distrib 19.0.0+65 w_2020_10 current setup): measurePhotonTransferCurve.py /project/shared/auxTel/ --rerun plazas/auxtelPTC2020MAR10 --id detector=0 --visit-pairs 2020022100013,2020022100014 2020022100015,2020022100016 2020022100017,2020022100018 2020022100019,2020022100020 2020022100021,2020022100022 2020022100023,2020022100024 2020022100025,2020022100026 2020022100027,2020022100028 2020022100029,2020022100030 2020022100031,2020022100032 2020022100033,2020022100034 2020022100035,2020022100036 -c makePlots=True ptcFitType=POLYNOMIAL polynomialFitDegree=3 -c isr.doBias=False isr.doDark=False isr.doSuspect=False isr.doLinearize=False The following gains were produced: gain_pickle_file = "/project/shared/auxTel/rerun/plazas/auxtelPTC2020MAR10/calibrations/ptc/ptcDataset-det000.pkl" gain_file = open(gain_pickle_file, 'rb') gain_data = pkl.load(gain_file) print (gain_data.gain) {'C10': 1.0348085853531022, 'C11': 1.02347459717329, 'C12': 1.0378724084326139, 'C13': 1.0358289608332496, 'C14': 1.012234515720022, 'C15': 1.0060929875013593, 'C16': 1.037898801107405, 'C17': 1.081185663200865, 'C07': 1.1520424164621434, 'C06': 1.1012153572252976, 'C05': 1.113872220848673, 'C04': 1.0939783643655419, 'C03': 1.0698829435732895, 'C02': 1.0906280187126696, 'C01': 1.104278016310778, 'C00': 1.1038149461060722} The, i applied the gains to image "dataId = dict(dayObs="2020-02-20", seqnum=433)" (expId 2020022000433 in '/project/shared/auxTel'): code (modified versioon of coode by Craig Lage: /project/cslage/ComCam/notebooks/Ingest_ComCam_06Mar20.ipynb): bias = butler.get('bias', dataId) isrConfig = IsrTask.ConfigClass() isrConfig.doLinearize = False isrConfig.doBias = True isrConfig.doApplyGains = True isrConfig.doFlat = False isrConfig.doDark = False isrConfig.doFringe = False isrConfig.doDefect = False isrConfig.doAddDistortionModel = False isrConfig.doWrite = False isrTask = IsrTask(config=isrConfig) gain_pickle_file = "/project/shared/auxTel/rerun/plazas/auxtelPTC2020MAR10/calibrations/ptc/ptcDataset-det000.pkl" gain_file = open(gain_pickle_file, 'rb') gain_data = pkl.load(gain_file).gain raw = butler.get('raw', dataId) ccd = raw.getDetector() builder = ccd.rebuild() for amp in builder: print (amp.getName()) amp.setGain(gain_data [amp.getName()] ) raw.setDetector(builder.finish()) ccd = raw.getDetector() postIsr = isrTask.run(raw, bias=bias).exposure

              People

              • Assignee:
                plazas Andrés Alejandro Plazas Malagón
                Reporter:
                plazas Andrés Alejandro Plazas Malagón
                Watchers:
                Andrés Alejandro Plazas Malagón, John Swinbank, Merlin Fisher-Levine
              • Votes:
                0 Vote for this issue
                Watchers:
                3 Start watching this issue

                Dates

                • Created:
                  Updated:
                  Resolved:

                  Summary Panel