skySim = SkySim() # Instantiate the telescope # There is no camera correction for the comcam at this moment tele = TeleFacade(M1M3=M1M3, M2=M2, phoSimCommu=phoSimCommu) # Set the configuration file path tele.setConfigFile(configFilePath) # Set the subsystem directory tele.setSubSysConfigFile(M1M3dataDir=M1M3dataDir, M2dataDir=M2dataDir, phosimDir=phosimDir) # Set the comcam camera tele.setInstName("comcam10") # Add the star skySim.addStarByRaDecInDeg(0, 0, 0, 17) skySim.addStarByRaDecInDeg(1, 0.0033, 0.0033, 18) # Set the degree of freedom for intra- and extra-focal images # Double check the direction with Bo instFileNameList = {"-1": "starExtra.inst", "1": "starIntra.inst"} logFileNameList = {"-1": "starExtraPhoSim.log", "1": "starIntraPhoSim.log"} outputImgDirIntra = os.path.join(outputImgDir, "Intra") outputImgDirExtra = os.path.join(outputImgDir, "Extra") outputImgDirList = {"-1": outputImgDirExtra, "1": outputImgDirIntra} argStringList = [] for ii in (-1, 1):
rotSkyPos=cameraRotation, mjd=mjdTime) # Instantiate the subsystems phoSimCommu = PhosimCommu() skySim = SkySim() # Instantiate the telescope # There is no camera correction for the comcam at this moment tele = TeleFacade(phoSimCommu=phoSimCommu) # Set the subsystem directory tele.setSubSysConfigFile(phosimDir=phosimDir) # Set the comcam camera tele.setInstName("lsst") # Add the star sensorName = ["R22_S11", "R22_S10"] starMag = [15, 16.5, 15, 15] xInpixelInCam = [500, 500, 800, 950] yInPixelInCam = [1000, 1080, 1300, 1500] starId = 23 for sensor in sensorName: for ii in range(len(xInpixelInCam)): skySim.addStarByChipPos(camera, obs, sensor, starId, xInpixelInCam[ii], yInPixelInCam[ii], starMag[ii], folderPath2FocalPlane) starId += 1