cameraMJD = 59580.0
obs = ObservationMetaData(pointingRA=RA,
pointingDec=Dec,
rotSkyPos=cameraRotation,
mjd=mjdTime)
# Instantiate the subsystems
phoSimCommu = PhosimCommu()
skySim = SkySim()
# Instantiate the telescope
tele = TeleFacade(phoSimCommu=phoSimCommu)
# Set the subsystem directory
tele.setSubSysConfigFile(phosimDir=phosimDir)
# Update the telescope degree of freedom
# dofInUm = np.zeros(50)
# Camera piston in um
# dofInUm[5] = 1000
# Camera dx in um
# dofInUm[6] = 500
# Set the telescope DOF
# tele.setDofInUm(dofInUm)
# Add the star on WFS
sensorName = [
# Instantiate the subsystems
M1M3 = M1M3Sim()
M2 = M2Sim()
phoSimCommu = PhosimCommu()
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")