Skip to content

lstautorta.High_Level_analysis

Functions:

Name Description
get_current_obs_info

Retrieve current observation info from the shared observation file.
get_night_timestamp

Get current timestamp (YYYYMMDD) with respect to the night.
main

Initialize current_run_is and quit variables and get the current timestamp with Start_time = now()
now

Returns the current timestamp in seconds, relative to the Unix epoch.
wait_for_directory

Wait for the directory to exist

get_current_obs_info 

get_current_obs_info(db_hostname='lst101', obs_target_query_timeout_s=10)

Retrieve current observation info from the shared observation file. Waits until the file is available.

Source code in src/lstautorta/High_Level_analysis.py 
33
34
35
36
37
38
39
40
41
42
def get_current_obs_info(
    db_hostname: str = "lst101",
    obs_target_query_timeout_s: int = 10,
):
    """
    Retrieve current observation info from the shared observation file.
    Waits until the file is available.
    """
    logging.info("Waiting for shared observation file at %s", DEFAULT_SHARED_OBS_PATH)
    return wait_for_shared_obs(DEFAULT_SHARED_OBS_PATH)

get_night_timestamp 

get_night_timestamp(today)

Get current timestamp (YYYYMMDD) with respect to the night. It will consider a same night everything run between 08h00 until 07h59 of the next day.

Ex: An observation launched at 04h00 (local time) will refer to the previous day.

:param today: datetime object datetime object with current local time

:return: timestamp: str

Source code in src/lstautorta/High_Level_analysis.py 
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
def get_night_timestamp(today):
    """
    Get current timestamp (YYYYMMDD) with respect to the night.
    It will consider a same night everything run between 08h00 until 07h59 of the next day.

    Ex: An observation launched at 04h00 (local time) will refer to the previous day.

    :param today: datetime object
        datetime object with current local time

    :return:
        timestamp: str
    """
    # TODO check that today is a datetime object ?
    if today.hour < 8:  # script launched between
        yesterday = today - datetime.timedelta(days=1)
        timestamp = f"{yesterday.year:04d}{yesterday.month:02d}{yesterday.day:02d}"
    else:
        timestamp = f"{today.year:04d}{today.month:02d}{today.day:02d}"

    return timestamp

main 

main()

Initialize current_run_is and quit variables and get the current timestamp with Start_time = now() Initialize the logging system by setting the log file name and logging level Change the current working directory to /fefs/onsite/pipeline/rta/data Find the current reservations for nodes using the scontrol command and parse the output to find the individual node names Modify a configuration file with sed commands, replacing placeholders with the node names found in step 4 Check if the ib0 network interface is in connected mode for each of the nodes and log the results

Source code in src/lstautorta/High_Level_analysis.py 
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
def main():
    """
    Initialize current_run_is and quit variables and get the current timestamp with Start_time = now()
    Initialize the logging system by setting the log file name and logging level
    Change the current working directory to /fefs/onsite/pipeline/rta/data
    Find the current reservations for nodes using the scontrol command and parse the output to find the individual node names
    Modify a configuration file with sed commands, replacing placeholders with the node names found in step 4
    Check if the ib0 network interface is in connected mode for each of the nodes and log the results
    """
    current_run_is = -1
    quit = False

    Start_time = now()

    # today = str(datetime.datetime.now().year) + str(datetime.datetime.now().month) + str(datetime.datetime.now().day)
    today = get_night_timestamp(datetime.datetime.now())

    wait_for_directory("/fefs/onsite/pipeline/rta/data/" + today_to_directory(today))
    logging.basicConfig(
        filename="/fefs/onsite/pipeline/rta/data/" + today_to_directory(today) + "/log_High_Level_" + today + ".txt",
        level=logging.INFO,
    )
    logging.info("Start RTA High Level for the day " + today)

    RTA_ready = True

    loop_id = 0
    data_directory = "/fefs/onsite/pipeline/rta/data/"
    while quit == False and RTA_ready:
        time.sleep(10)
        obs_info = get_current_obs_info()
        if obs_info is None or obs_info.obs_id is None:
            logging.info("None")
            current_run_is = -1
            continue
        if current_run_is == -1:
            current_run_is = obs_info.obs_id
        if current_run_is != obs_info.obs_id:
            time.sleep(1)
            obs_info = get_current_obs_info()
            if obs_info is None or obs_info.obs_id is None:
                logging.info("None")
                current_run_is = -1
                continue
            current_run_is = obs_info.obs_id
            if current_run_is != -1:
                logging.info("Start RTA check for run " + str(current_run_is))
                time.sleep(2)
                source_ra, source_dec = get_source_ra_dec(obs_info)
                logging.info(
                    # ~ "srun --reservation=rta_one_node ./Theta_square.py -da "
                    "./Theta_square.py -da "
                    + today_to_directory(today)
                    + " -r "
                    + str(current_run_is)
                    + " -RA "
                    + str(source_ra)
                    + " -DEC "
                    + str(source_dec)
                )
                os.system(
                    # ~ "srun --reservation=rta_one_node ./Theta_square.py -da "
                    "./Theta_square.py -da "
                    + today_to_directory(today)
                    + " -r "
                    + str(current_run_is)
                    + " -RA "
                    + str(source_ra)
                    + " -DEC "
                    + str(source_dec)
                )
                logging.info(
                    # ~ "srun --reservation=rta_one_node ./Spectra.py -da "
                    "./Spectra.py -da "
                    + today_to_directory(today)
                    + " -r "
                    + str(current_run_is)
                    + " -RA "
                    + str(source_ra)
                    + " -DEC "
                    + str(source_dec)
                )
                os.system(
                    # ~ "srun --reservation=rta_one_node ./Spectra.py -da "
                    "./Spectra.py -da "
                    + today_to_directory(today)
                    + " -r "
                    + str(current_run_is)
                    + " -RA "
                    + str(source_ra)
                    + " -DEC "
                    + str(source_dec)
                )
                logging.info(
                    # ~ "srun --reservation=rta_one_node ./rta_var.py --dl3_folder "
                    "./rta_var.py --dl3_folder "
                    + data_directory
                    + today_to_directory(today)
                    + "/"
                    + str(current_run_is)
                    + "/DL3 -ra "
                    + str(source_ra)
                    + " -dec "
                    + str(source_dec)
                    + "  --output_dir "
                    + data_directory
                    + today_to_directory(today)
                    + "/"
                    + str(current_run_is)
                    + "/plots --date "
                    + today
                )
                os.system(
                    # ~ "srun --reservation=rta_one_node ./rta_var.py --dl3_folder "
                    "./rta_var.py --dl3_folder "
                    + data_directory
                    + today_to_directory(today)
                    + "/"
                    + str(current_run_is)
                    + "/DL3 -ra "
                    + str(source_ra)
                    + " -dec "
                    + str(source_dec)
                    + "  --output_dir "
                    + data_directory
                    + today_to_directory(today)
                    + "/"
                    + str(current_run_is)
                    + "/plots --date "
                    + today
                )

        if loop_id % 1 == 0:
            logging.info("Start RTA check for run " + str(current_run_is))
            time.sleep(2)
            logging.info(datetime.datetime.now())
            source_ra, source_dec = get_source_ra_dec(obs_info)
            logging.info(
                "./Theta_square.py -da "
                + today_to_directory(today)
                + " -r "
                + str(current_run_is)
                + " -RA "
                + str(source_ra)
                + " -DEC "
                + str(source_dec)
            )
            os.system(
                "./Theta_square.py -da "
                + today_to_directory(today)
                + " -r "
                + str(current_run_is)
                + " -RA "
                + str(source_ra)
                + " -DEC "
                + str(source_dec)
            )
            logging.info(
                "./Spectra.py -da "
                + today_to_directory(today)
                + " -r "
                + str(current_run_is)
                + " -RA "
                + str(source_ra)
                + " -DEC "
                + str(source_dec)
            )
            os.system(
                "./Spectra.py -da "
                + today_to_directory(today)
                + " -r "
                + str(current_run_is)
                + " -RA "
                + str(source_ra)
                + " -DEC "
                + str(source_dec)
            )
            logging.info(
                "./rta_var.py --dl3_folder  "
                + data_directory
                + today_to_directory(today)
                + "/"
                + str(current_run_is)
                + "/DL3 -ra "
                + str(source_ra)
                + " -dec "
                + str(source_dec)
                + "  --output_dir "
                + data_directory
                + today_to_directory(today)
                + "/"
                + str(current_run_is)
                + "/plots --date "
                + today
            )
            os.system(
                "./rta_var.py --dl3_folder  "
                + data_directory
                + today_to_directory(today)
                + "/"
                + str(current_run_is)
                + "/DL3 -ra "
                + str(source_ra)
                + " -dec "
                + str(source_dec)
                + "  --output_dir "
                + data_directory
                + today_to_directory(today)
                + "/"
                + str(current_run_is)
                + "/plots --date "
                + today
            )

        loop_id = loop_id + 1
        if (datetime.datetime.timestamp(datetime.datetime.now()) - Start_time) > 14 * 3600:
            quit = True
            logging.info("End of the night, Stop the RTA check")

now 

now()

Returns the current timestamp in seconds, relative to the Unix epoch.

:return: float

Source code in src/lstautorta/High_Level_analysis.py 
24
25
26
27
28
29
30
def now():
    """
    Returns the current timestamp in seconds, relative to the Unix epoch.

    :return: float
    """
    return datetime.datetime.timestamp(datetime.datetime.now())

wait_for_directory 

wait_for_directory(directory, interval=5)

Wait for the directory to exist

:param directory: directory path :param interval: Time interval between verification

Source code in src/lstautorta/High_Level_analysis.py 
11
12
13
14
15
16
17
18
19
20
21
def wait_for_directory(directory, interval=5):
    """
    Wait for the directory to exist

    :param directory: directory path
    :param interval: Time interval between verification
    """
    while not os.path.exists(directory):
        print(f"{directory} does not exist yet. let's wait {interval} seconds...")
        time.sleep(interval)
    print(f"The directory {directory} finally exist !!!")