Transient detections and other real-time data processing from wide-field chambers MASTER-VWF

At present time Robotic observatory making is of current importance. Having a large field of view and being able to point at anywhere, Robotic astronomical systems are indispensable when they looking for transients like grb, supernovae explosions, novae etc, as it's impossible in these cases to foresee what you should point you telescope at and when. In work are described prompt GRB observations received on wide-field chambers MASTER-VWF, and also methods of the images analysis and transients classifications applied in real-time data processing in this experiment. For 7 months of operation 6 synchronous observations of gamma-ray burst had been made by MASTER VWF in Kislovodsk and Irkutsk. In all cases a high upper limits have been received (see tabl \ref {tab_grbwf} and fig. \ref {allgrb}).


Very wide field telescopes.
Currently, several dozen robotic telescopes are operated worldwide, each of them usually meant for certain dedicated tasks. A special case of these instruments are super-wide field telescopes with fields of view ranging from hundreds to thousands square degrees. It is impossible to make a very wide field system with a large aperture using now available CCDs and optical systems, and that is why the apertures of very wide field cameras are limited to 5-15 cm. All very wide field cameras are designed to record various optical transients. Most of the optical transients are of circumterraneous origin -meteors and satellites, however, very wide-field cameras may also provide valuable data on the prompt emission of cosmological grbs. A striking example is the incredibly successful observation of GRB080319B that was performed within the framework of the Russian-Italian experiment TORTORA [5] and Polish experiment "Pi of the Sky" [11]. All transients are highly variable and therefore snapshots must be taken at intervals of several seconds or even less. Hence solutions capable of operating without dead time (such as the time needed for ccd readout), like MASTER VWF [1], TORTORA and others (see table 1), have an important advantage over other systems.
The operation of any automatic astronomical facility is supported by a dedicated software package, which also performs primary image processing and data storage. Which software to use for an automatic astronomical system depends on the particular task the facility is meant to perform, and multi-purpose solutions are therefore unsuitable in such cases. Classical universal techniques are even less suited for super wide field observations, because of the problems due to the field curvature and very large data amount.

MASTER -Mobile Astronomical System of the TElescope Robot
MASTER was the first robotic telescope of this kind constructed in Russia. It was built and tested by students, post-graduate students and lecturers of SAI MSU in 2002 (http://observ.pereplet.ru). MASTER Table 1: Comparative characteristics of super-wide field systems Project FOV Exposure Limiting magnitude WIDGET [10] 44 o x 44 o x 3 5 sec 10 m RAPTOR 40 o x 40 o 60 sec 12 m Pi of the Sky [11] 22 o x 22 o 10 sec 11.5 m Yatsugatake Cameras 85x70 8 5 m FAVOR [6] 16x24 0. 13 11.5 m TORTORA [7] 24x32 0. 13 10.5 m MASTER-VWF(1) [1] 18 has evolved into a broad network of telescopes spread across the entire country. A detailed description of the MASTER facility can be found in [2]. Below we discuss the basic properties of MASTER wield field cameras (MASTER-VWF).
The main task of the MASTER-VWF experiment is to perform continuous all-sky monitoring with the aim to detect all objects absent in the available astronomical catalogues. This includes, in particular: • Prompt detection of grb emission synchronously with gamma-ray space observatories.
• Search for star-like transients of unknown origin and orphan bursts.
• Detection of meteors and determination of their basic parameters: brightness (photometry), velocity (astrometry), and altitude of combustion in the atmosphere (triangulation). • Detection of satellites and space debris and determination of their basic parameters -astrometry, velocity, brightness, and altitude above the earth's surface.
For these purposes, a completely robotic MASTER VWF4 facility with four wide-field cameras and a total field of view of 4000 square degrees has been deployed at the MSU Caucasian Mountain Astronomical Observatory near Kislovodsk(2075m). Another facility -MASTER VWF2 (two cameras and 2000 square degrees field of view) -has been deployed near Irkutsk. In Kislovodsk cameras are set on parallactic mounts in pairs 702 meters apart (see fig. 2). Each mount has an automated dome and two fast and powerful (11 Megapixels) Prosilica GE 4000 CCD cameras with Nikkor 50mm (f /1.2) lenses. The camera is capable of continuous imaging with a speed of one frame in 0.2 to 60 seconds with no time gaps. The 702-m baseline allows the facility in Kislovodsk to determine the parallaxes of circumterraneous objects and reconstruct their altitudes from the results of astrometric reduction.
All telescopes operate in a fully autonomous mode: at night, if the weather conditions permit (a "Boltwood Cloud Sensor" provides the current weather report updated every three seconds), the cameras begin to automatically survey the sky. The frames taken are then processed and stored in a special data storage device. Because of the very large data flow (up to 700 Gb per night), the images can be kept (2) the main tower and shelter of the southern MASTER-VWF4 installation, and (3) the main tower of the MASTER-2 telescope [2] only for several days and then most of them have to be erased, except for some images with objects of special interest like meteors, candidate optical transients, and satellites. We also keep the images taken simultaneously with grbs. The survey is halted and the telescope dome is closed every time when the weather deteriorates or at dawn. The operation of the entire MASTER VWF4 facility is controlled by seven computers (see Fig. 3 ): four computers to control the operation of the Prosilica GE 4000 CCD camera and perform source extraction; two computers to control the operation of the dome and mount; one big quad-core server to process and store the data obtained and publish it in the Web.
The facility is a complex engineering system, which needs sophisticated and powerful software to control it, perform observations, and data processing. The data flow can be as high as several Terabytes/day (in the case of very short exposures on the order of 0.2s), requiring highly optimized software, powerful data storage facilities, and interfaces. With 5-second long exposures the limiting is 11.5-12 on each frame depending on weather conditions. Thus every image contains about 10-15 thousand objects, resulting in a data flow that is impossible to store and has to bee processed in real time.
A special software package has been developed for this purpose. It allows the coordinates of all objects and their photometric characteristics (profile, magnitude) to be determined in real time, and is also capable of finding and analysing optical transients. MASTER VWF can measure the coordinates of the objects with an error smaller than 10 arcsec. The system is supplied with a daily updated database of the ephemerides of artificial satellites and a program developed to identify the satellites in real time. The program is capable of finding new artificial satellites, to observe the destruction of known artificial satellites and not only to control space debris, but also observe the process of its formation. Sometimes unknown or lost objects are recorded. For example, on October, 20th, 2008 MASTER VWF recorded the debris of the American military satellite USA114Deb (see fig. 4).
A total of six synchronous observations of gamma-ray bursts have been made by MASTER VWF in Kislovodsk and Irkutsk during seven months of their operation. In all cases high upper limits have been determined (see tabl 4.3 and fig. 1).
In the following chapters we describe these and other results in more detail.
2 Primary reduction in wide fields: astrometry and photometry Let us consider the methods of astromeric and photometric reduction of a wide field in more detail. The identification of the field of view and calculation of the coordinate grid for the image is based on the table of Cartesian image-plane coordinates (X,Y) produced by program used for the extraction of objects (SExtractor) and the the Tycho-2 astrometric catalogue, which contains exact coordinates for most of the stars down to a limiting magnitude of 11 m .5. To perform this task, the program should meet high requirements in terms of speed of processing (10-15 seconds and less) for a field containing about 10000 stars; ensure high probability of identification (more 99 %); be applicable to various fields ranging in size from from narrow (tens of square minutes) to super-wide (thousand square degrees) fields; allow for a certain freedom in the telescope parameters, because the focal length and other parameters may change in the course of the year due to the variation of the outside temperature, and be capable of operating with approximate initial midfield coordinates provided that they are within the field of view. All processing in MASTER experiments is performed under Linux OS and hence no ready-made software was likely to be found. The current implementation of the "astrom" program is based on the algorithm of pairing similar triangles of bright stars in the image and in the catalogue. This program works perfectly with the most different images from those taken under poor weather conditions and containing a few dozen stars (typical situation for alert observations) to images in the Milky-Way fields taken during a good observing night and containing up to 100000 objects. The root-mean-square error of the measured coordinates ranges from 1 / 4 of arcsecond (scale 2.1"/pixel) for most of the images taken with the basic telescope of the MASTER system to about 7-8 arcseconds for very wide-field cameras (scale 36"/pixel) (see fig. 2).
The system may also operate in a "fast" astrometry mode to process consecutive frames of the same area taken with wide-field cameras. In this mode the software uses a priori information about the field of view based on the previous frames, allowing the operating time to be reduced to less than 0.2 seconds. There is also a dedicated mode for more accurate determination of the coordinates of special objects (gamma-ray bursts, supernovae, asteroids, meteors, and others) with the local coordinate adjustment computed only in the field of the object searched using 10-25 nearest stars with well accurate coordinates. Our experience shows that this mode provides a factor of two improvement of the accuracy of the inferred coordinates.
Our program for extracting objects from images is based on the SExtractor software package written for the TERAPIX [9] project. SExtractor is written in the C language and offers sufficiently high performance. We adapted this program to the task of reducing images taken with MASTER, and describe our version of the program in the following chapter. The main parameter used by the extraction program is the threshold signal-to-noise ratio. This value directly affects both the program execution time and the accuracy of the results obtained. We adopt S/N=2 for primary selection, whereas at the subsequent stages of processing we select from the initial list only the objects with higher S/N values.

Photometric calibration
Practically all wide-field systems operate without filters. MASTER VWF is not an exception. This poses the problem of the calibration of the photometric band of the measurements. It is clear that our instrumental photometric band is limited to a certain extent by the CCD sensitivity curve and lens transmission curve. Our task is therefore to determine instrumental magnitudes based on the magnitudes of the reference stars adopted from photometric catalogues.
To this end, colour-colour diagrams are analysed, which relate the difference between the instrumental magnitude and the magnitude in any of the modelled bands (or combinations thereof) to a chosen colour index for each detected star. The closer is the modelled band to the instrumental band, the more horizontal appears the the cloud of points on the diagram. Fig. 6 shows the corresponding colour-colour diagrams for the MASTER VWF-4 system consisting (in its photometric part) of a Nikkor 50 mm f/1.4 lens and Prosilica GE4000 CCD. The instrumental photometric band can be described fairly well by the V T Y CHO2 filter. The absolute accuracy of photometric measurements made with wide-field cameras is, on the average, of about ∆m abs ∼ 0.25 − 0.35 m . Note that the relative accuracy of measurements is substantially better with errors as low as ∆m rel 3 − 5%.

Automatic determination of the limiting magnitude of an image
Automatic determination of the limiting magnitude for the images taken is an important part of the operation of modern automatic sky surveys. An automatic survey must naturally include automated quality control. The automatically determined limiting magnitude may serve as some kind of objective . A more detailed analysis shows that the third-degree fit yields slightly better results. That is logical, because most of the known aberrations (e.g., distortion) depend on the third power of radial distance.
quality assessment for the images obtained. Furthermore, the determination of the limiting magnitude is also needed for gamma-ray burst observations. The task of the determination of the limiting magnitude can be addressed formally by setting the limiting magnitude of the frame equal to the magnitude of the faintest object. This approach, however, has a number of shortcomings. First, when the system is operated with a low extraction threshold most of the faintest objects are noise artifacts. If the sky is covered with clouds, no real faint objects can be seen and practically all objects at the visibility limit are chance fluctuations of cloud brightness. Moreover, because of the large number of such false objects some of them may be identified with real stars and hence the situation cannot be remedied by restricting the list to objects identified with catalogued stars. Second, this approach allows one to estimate the limiting magnitude only within a local area and the quantity inferred does not characterise the entire image and is therefore unsuitable for quality control.
These problems can be addressed as follows. All stars of the reference catalogue that are located within the field considered are subdivided into magnitude bins of fixed width. For each bin the fraction of identified stars to the total number of catalogued stars is computed, and the corresponding histogram is drawn (see fig. 7). One can see that the ratio decreases sharply when the average magnitude of the bin approaches the limiting magnitude. In practice, it is convenient to adopt a certain threshold ratio and define the limiting magnitude as the magnitude at which the curve reaches this threshold. The thresholds . colour-colour diagrams are used to calibrate instrumental magnitudes. On these plots m = 2.5lg(Finst) and the zero point is of no importance and cannot be determined at the this stage. Note that the width of the sequences is real and due to the fact that stars have different spectral types, and not due to measurement errors.
of σ lim = 25% and σ lim = 50% have been shown empirically to correspond to objects with S / N ≈ 9 − 11σ and S / N ≈ 14 − 15σ, respectively. Note that for robust operation of the algorithm the limiting magnitude of the reference catalogue should be fainter than that of the frame studied. Otherwise the method becomes inapplicable because the catalogue lacks lack stars of the given magnitude. However, the USNO-B1 catalogue is sufficient for most of the modern wide-field surveys, whereas very wide-field surveys can be reduced successfully even with the TYCHO-2 catalogue (and this is very convenient).
3 Automatic classification and analysis of astronomical transients.
3.1 Extraction of objects and data transfer to the server.
The search for and classification of transients starts once the primary reduction of the image is finished. All transient events in frames can be subdivided into two classes: starlike objects and strips, which differ substantially in the methods used for primary reduction. All rapidly moving objects are usually classified as strips: low-altitude satellites and meteors. All transients of stellar or extragalactic nature (novas and supernovas, orphan bursts, stellar flares) and geostationary satellites (see, for example, fig.18) are starlike. Starlike objects are much easier to process, because the program of object extraction and instrumental and σ lim = 25% thresholds. Averaging was made within 0.5 m and 0.1 m bins in the right and left histograms, respectively. The right (coarser binned) histogram is smoother, resulting in reduced probability of accidentally reaching the threshold at a bright magnitude. It is therefore convenient to determine the limit on the coarsely binned histogram and then refine it on the finer-binned curve. photometry (SExtractor, described above) has been developed to analyse such objects. Therefore all main parameters of starlike objects (flux and magnitude, FWHM, Cartesian frame coordinates (X, Y), and semiaxes (a, b)) are determined very reliably. Extraction of strip-shaped objects is fraught with serious problems. First, a fine adjustment of SExtractor's additional internal parameters is necessary. Second, a change of these parameters strongly affects the image processing time. Furthermore, strip-shaped objects contribute appreciably to the background within the small area suitable for identification of stars. Fine adjustment of all parameters allowed us nevertheless not only to identify strips in the image, but also achieve high speed of operation (half exposure time) (see fig.8). However, SExtractor algorithm cannot identify strips in the frame, to say nothing about performing their high-quality photometry. To this end, a special program had to be developed, which we describe below.
To reduce the server load, identification of objects in the image can be physically performed not only by the server, but also by the control computer of each particular camera (see fig. 3). Note that the resulting catalogue of detected objects should have higher priority for the transmission to the server compared to the images that are to be stored in the database.

Primary reduction of astronomical transients.
The next stage after identification of objects is their astrometric and photometric reduction, which we describe in detail in the previous section. We thus obtain, after identification, astrometric and photometric calibration, a catalogue of objects identified in the frame together with calibrated coordinates and photometry. All unidentified objects (identification is performed using both coordinate and photometric criteria: usually ∆S = 20 − 30" and ∆m = 1.5 m ) are viewed as possible optical transients. This approach allows us to identify stellar flares as transients and to a certain extent prevent spurious identifications with fainter and hence more numerous stars. To filter noise, objects that could not be identified with the catalogue are correlated with the objects of several (usually 2-3) previous sets. An object is considered to be a candidate transient if it could not be found either in the catalogue or in in the previous frames. All candidate transients are then passed through a number of special filters. For example, if the semimajor axis of the object, which should typically be of about b ∼ 0.7 − 2 pix, exceeds b ∼ 4 − 5 pixels, it must be a trail of a cloud with a probability of 95 %. If f whm < 1 pixel, we must be dealing with a hot pixel or a cosmic-ray hit. We also set for all transients the trigger threshold, which is always higher than the The left-hand figure shows a frame with objects extracted using standard SExtractor parameters (with the extraction threshold of S /N > 1.5) with no filters applied. Ax a result, the program identified all starlike objects, heaps of faint garbage, hot pixels, and simply random noise, i.e., virtually everything except the bolide. The right-hand frame shows the result of extraction made using the currently adopted refined parameters. Several bright elongated objects are extracted on the bolide image and the noise level is appreciably lower. extraction threshold. Currently, the trigger threshold is S / N > 10σ. We consider a candidate found on two or more consecutive frames to be a boda fide transient.
Preliminary classification of transients is performed after the objects are passed through all filters. First, all objects are searched for possible identifications in the catalogue of artificial satellites and all identified objects are included into the database of satellites. The remaining ones are subdivided into starlike and strip-shaped objects. The object is classified as a strip-shaped if a/b > 3, otherwise it is considered to be starlike (here a and b major and minor semiaxis). All starlike objects are then passed through the coincidence circuit, and strip-shaped objects are. Strips are provisionally subdivided into meteor candidates and unknown artificial satellites. If each of the three successive frames contains a candidate transient such that the three objects lie on a straight line, all of the them are considered to be the images of an artificial satellite. The remaining objects are considered to be meteor candidates.

Coincidence scheme of and determination of heights.
All the procedures described above are performed for each of the four cameras in real time mode. Real time processing ends here if the cameras operate in uncoupled (alert) mode. In the case of survey observations, when two cameras point to the same sky area, all starlike transients are passed through the coincidence scheme, and all strip-shaped objects have their heights determined (see fig. 9, 10). To improve the accuracy of identification, special transformation is applied to the frame taken with the Northern camera to convert it to the coordinate frame of the Southern camera. This transformation is based on the same principles  as astrometric reduction. It takes little time to compute (∼ 0.2 seconds) and improves the accuracy (and hence reduces the correlation radius and the number of chance identifications) by more than a factor of 2.5 compared to the accuracy of the results obtained in the case of identification by celestial coordinates. Should a real starlike transient be found, an alert can be immediately issued to larger (45-cm) MASTER telescopes (located near Moscow, in Kislovodsk and in the Urals) to study the event in more detail. However, this mode is currently in the stage of development, mostly because of the debris of high-orbital satellites, which often contaminate the frames (see fig. 4, 18).

Analysis of strips.
As we pointed out above, SExtractor computes only the approximate position of a strip, breaking it into several bright elongated objects. To perform complete photometry of a strip, determine the coordinates of its centre and ends and other parameters, a special program is used dedicated to the analysis of strips on astronomical images. To prevent server overload, this program is run only when cameras are idle (in the afternoon or during bad weather). The program starts with the coordinates of any location inside the strip determined at the previous stage. The working area of the given size (usually 512 on 512 pix, which is sufficient for the analyses of any satellite or meteor trails) is chosen around this initial point. The size of the working area is automatically increased when a very bright and long bolide is to be analysed.

Determination of background.
To determine the background, we subdivide the working area into subareas of the given size. For each subarea the median and standard deviation are computed twice: the first time over almost all pixels in the given subarea, and the second time more precisely, after discarding the pixels that deviate strongly (by more than 6-10σ) from the preliminary background estimate. We number all these subareas -there are a total of N zone of them (currently N zone = 256)-and compute the background and standard deviation in each subarea as where Bg (i,j) is the background at the point with coordinates (i,j); Med (k) is the median in k-th area used to compute the background; Σ is the standard deviation within the corresponding subaera, and r k is the distance between the point (i.j) and the centre of k-th area. Figures 11:3-4 show the corresponding background maps.

Extraction of strips.
To extract a strip, we first compute the map of bright pixels exceeding the background level by n σ ·Σ(i, j) or more. The optimum threshold has been shown empirically to be n σ = 1.5. The corresponding points are shown by black shade on the map (see fig 11:2).
A recursive algorithm of the search for connected area is then used starting at the centre found by SExtractor. Once the connected area is found, the main parameters of the object (the position of the centre, a, b, and the angle θ between the axis of the object and the X axis) are determined. The standard deviation σ line from the straight line drawn at angle θ through the centre of the object is computed if the object found can be classified as a strip (a/b > 3). The recursive algorithm of searching for connected area is then started again, this time within 2-3σ line of the straight line just mentioned. The algorithm is allowed to jump over one to two faint pixels within 2σ line of the straight line. Thus the second pass yields a completely extracted strip including faint edges in the case of meteors.
The program then attempts to find objects located on the extension of the strip line on the subsequent and previous frames. If such object can be found, the strip is classified as a satellite, otherwise, as a meteor. The brightness profiles along the strip are then computed both for satellites and meteors (see fig. 12 ). Portions of images are kept for all transients.

The importance of synchronous observations.
As is well known, optical grb glow can be roughly subdivided into two parts: prompt emission and afterglow. We do not discuss the afterglow phenomenon here, because it has already been thoroughly studied and it is not the target of wide field observations. Let us now consider prompt emission.
Prompt optical grb emission is the emission that appears simultaneously with gamma rays detected by various space-based gamma observatories (Hete2, Swift, Fermi). Unlike afterglow, it bears information about the object that gives birth to grb, and not about its environment. There are two ways to detect prompt optical grb emission: synchronous observations and very quick alert pointing. To date, about dozen [8] optical prompt grb have been made and most of them were alert based. That is to say, when prompted by a quick trigger from a gamma-ray observatory, a robotic optical telescope points at the object within several tens of seconds after the corresponding gamma-ray event is detected. If the grb is long enough, its prompt emission is still detectable by this time. Another, less common type of observations, are synchronous wide-field camera observations. Perhaps the most amazing example of the latter is the prompt GRB 080319B detection [5]. Synchronous observations are much more efficient than alert observations for at least two reasons. First, alert observations are subject to selection effects. The point is that an alert may result in the detection of prompt emission only if the grb lasts longer than 30-40 seconds, because of the delay due to the time it takes to perform data processing onboard the space telescope, transmit the signal to the ground, and point the optical telescope. Hence alerts cannot be used to detect prompt optical emission from short grbs. It is significant that this kind of emission has never been observed and moreover, it has never had its upper limit measured.

Energy calibrations of the MASTER-VWF4 photometric passband.
To make any conclusions about the brightness of objects not just in terms of relative magnitudes, but also in terms of energy units, we must calibrate the MASTER-VWF4 zero-point. Let us determine the energy flux (in the units if erg/cm 2 /s) produced by a zero-magnitude star (e.g., Vega). The MASTER-VWF4 band is close to the V band of the Tycho2 catalogue. Figure 13 shows the spectrum of Vega and the same spectrum folded with the lens passband and sensitivity curve of the CCD. Integration over the entire spectrum finally yields the flux produced within the MASTER-VW4 band by a zero-magnitude star

Synchronous GRB observations with MASTER VWF very wide-field cameras
.
MASTER VWF cameras perform nightly monitoring of the sky with the aim to detect optical transients. These cameras may operate in two different modes. In the first mode, which is meant for searching for optical transients, two cameras are pointed to the same sky area in order to implement the coincidence scheme (see fig. 9, 10). In the second mode the fields of view of the cameras are set apart as far as possible in order to provide a combined field of view of more than 4000 and 2000 deg 2 for the Kislovodsk and Irkutsk observatories, respectively. The second mode is dedicated for synchronous observations of gamma-ray bursts.

The GRB081102 and GRB090424 gamma-ray bursts recorded by Swift spacecraft
GRBs recorded by the Swift gamma-ray observatory satellite have very accurate coordinates, allowing the location of the burst origin on very wide-field cameras to be determined with the highest precision possible (i.e., to within less than one-pixel large error box).

Long gamma-ray burst GRB081102
(see fig. 14 ) was recorded by Swift gamma-ray space observatory on November, 2nd, 2008 17:44:39.5 UT [29]. Although it was a typical long (T 90 = 63sec) gamma-ray burst, the telegram from Swift was issued only after 15 minutes. This delay prevented early observations by automated telescopes. However, MASTER VWF-4 in Kislovodsk covered the entire error box at the edge of its field of view. MASTER VWF-4 operated in the survey mode with aligned cameras, allowing the grb error-box to be recorded simultaneously in two channels. This area has been observed for two hours before and seven hours after the grb with no time gaps [14].
Unfortunately, no optical counterpart could have been recorded from this burst (see fig. 14, 15). The lack of optical counterpart may have been due to extremely strong absorption toward the burst (4). However, a very high upper limit (for synchronous (!) observations with very wide-field cameras) of V grb081102 < 13.0 m has been obtained, which is the highest synchronous upper limit in history [13].
Let us now compare GRB081102 with GRB080319B from which bright prompt optical emission V = 5.3 m ( [5]) has been recorded. Consider the F opt /F γ ratio (where F opt and F γ are the optical and gammaray flux, respectively) for both bursts. The gamma-ray spectrum is known for all bursts discussed in this paper and therefore we reduce F γ to the flux in the 15 − 150keV energy interval for the sake of uniformity.
Let us estimate F grb081102 opt . The V-band extinction can be estimated by the following simple empirical formula [12] Fig. 15. GRB081102. Hash (red) from the Northern and Southern cameras and from the two cameras simultaneously inside the GRB081102 error-box and the signal from the V ref = 11.5 m reference star (blue) Fig. 16. A 6x6 deg 2 field around the GRB090424 Swift XRT position. One-minute exposure images taken before, during and after TGRB. Clouds can be seen on the images, resulting in the optical limit of V < 9.7 m on each frame. Animation available here: http://observ.pereplet.ru/images/GRB090424/grb film.html Given the known hydrogen column density toward GRB081102, N H = 4.9 * 10 21 cm −2 [31], we can use the inferred optical depth and Pogson's formula to estimate F grb081102 γ = 2.3 · 10 −6 erg/cm 2 [22] or τ grb081102 V = 5.2 · 10 −22 · 4.9 · 10 21 = 2.548 => δm = 2.8 m (4) 2.512 0−13.0+2.8 · 6.4 · 10 −6 erg/s/cm 2 · 63s 2.3 · 10 −6 erg/cm 2 = 1 83 The "less than" sign is used because only the upper limit could have been determined for grb081102. For GRB080319B the column density is N H = 9.2 * 10 20 cm −2 [5], i.e., five times lower than toward GRB081102 τ grb080319B 2.512 0−13.0+2.8 · 6.4 · 10 −6 erg/s/cm 2 · 50s 8.25 · 10 −5 erg/cm 2 = 1 22 It thus follows that the fraction of prompt optical emission for GRB081102 is at least four times smaller than for grb080319B.  [23]. At that time MASTER VWF-2 in Irkutsk operated in the alert mode (with cameras pointed apart) and GRB090424 happened to be at the centre of the field of view of one of the cameras. It was also a typical long (T 90 = 48sec) gamma-ray burst and UVOT detected afterglow with unfiltered magnitude of about 15.3 m 167 seconds after the burst [24]. MASTER VWF-2 observed the grb error-box without time gaps for 1.2 hours before, during and 1.5 hours after the gamma-ray burst with 1-second exposures. Unfortunately, weather conditions were far from favourable, resulting in a not so high upper limit of V 1sec < 8.0 m and V 60sec < 60.0 m for a coadd of 60 images. [16] and [17]. According to both Swift/BAT [23] data and FERMI observatory data converted to the 15 − 150keV energy interval [27], F grb090424 γ = 2.1e − 5, implying that Note that despite the rather modest upper limit the burst was very powerful and even such an upper limit restricts substantially the possible flux.

Gamma-ray bursts recorded by FERMI observatory
Unlike Swift, the FERMI observatory does not specialise in gamma-ray burst observations. Its onboard equipment allows the burst coordinates to be determined only to within several degrees or even worse. FERMI currently issues two to three times more alerts than Swift and is capable of recording harder (up to 10 MeV) radiation (Swift has an upper sensitivity limit of about 300 keV), and responds to short hard gamma-ray bursts, whose prompt emission has so far never been detected. Given the large error of the coordinate measurement, gamma-ray bursts detected by FERMI observatory may at present (and, very likely, also in the future) have their prompt emission observed only with wide-field cameras, such as MASTER-VWF.
Until now (July, 1 2009), a total of five FERMI GRB error boxes have overlapped with the MASTER-VWF field of view, and two more (GRB081215 and GRB090526) proved to be outside the field of view after a very substantial final final refinement of their coordinates.

GRB081130
was the first burst recorded by FERMI and synchronously observed by the MASTER VWF4 facility in Kislovodsk. Its coordinates and error-box have been refined repeatedly (with a characteristic scatter of ±10 degrees), but always fell within the field of view of our cameras. Figure 17 shows the final localisation of this gamma-ray burst. The error box of grb081130 has been observed without time gaps for several hours before, during, and for hours after the grb.
No proper candidates could have been found inside the error box. However, given the rather large uncertainty of the coordinates measured by FERMI, we extended the search area for transient candidates to the entire frame. As a result, a very interesting object resembling a grb optical counterpart has been found 14 degrees from the localisation centre synchronously with grb ( fig. 18). However, after cross identification with the catalogue of artificial satellites this transient was found to coincide with the high-orbital satellite, which occasionally produces such flashes. Thus a rather high upper limit of V grb081130B > 12 m .0 has been obtained again after the coadd of four images (T grb081130B = 15 − 20s taken in different channels [30] ) [20].
When converted to the 15−150keV energy interval, the gamma-ray fluence from GRB081130 is F grb081130B γ = 1.01 · 10 −6 erg/cm 2 . We then apply the same procedures as in the case of GRB081102 to infer F grb081102 opt = 3.82 · 10 −9 erg/cm 2 , implying a flux ratio of:

090305B
The gamma-ray burst GRB090305 (Fig. 19) was recorded by FERMI on March, 5th, 2009 14:12:09.33 UT [25]. MASTER VWF-4 in Kislovodsk operated in the alert mode (with cameras apart) and observed 80 % part of the GRB090305 error-box at the edge of the field of view of one of cameras. Despite of its rather short duration (T 90 = 2sec) [25], GRB090305 was nevertheless a rather long and bright gamma-ray burst. However, even the final position estimate has a large error-box (α = 10 h 20 m , δ = 68 • 06 ′ R 1sigma = 5.4 • )  and therefore no alert observations have been made. MASTER VWF-4 has observed the error-box of this GRB without time gaps for 3.5 hours before, during and for two hours after the with 1-second exposures.
No optical transients close to T GRB could be found and we obtain an upper-limit of V 1sec < 9.5 m [18].

GRB090320B
The MASTER VWF-4 system in Kislovodsk observed synchronously 80 % of the GRB090320B (see MASTER-VWF4 Kislovodsk GRB090320B follow-up observation. This is a 40x24 deg 2 full-frame image. The black circle shows the error box corresponding to the 3 sigma FERMI error limits plus the 3-degree systematic error (R = 12.6deg) [26]. Animation available here http://observ.pereplet.ru/images/GRB090320B/GRB090320B 60sec.gif the GRB090424 gamma-ray burst indicates that even the measured optical limit may constitute a result of real importance and that and this burst deserves separate analysis.
In this article we describe the methods of observation and analysis of images taken with very widefield cameras, and their application to synchronous observations of prompt gamma-ray burst emission in the MASTER-VWF experiment. We finally recall that the renowned astrophysicist Bohdan Paczynski [4] pointed out in his last paper the exceptional astrophysical significance of robotic telescope networks, which, like MASTER, are composed of wide field searching cameras and larger telescopes. Now his and our ideas at last become reality. This work was supported by the Ministry of Science of the Russian Federation (state contract N 02.740.11.0249 ). Fig. 21. MASTER-VWF4 follow-up observations of GRB090328B in Kislovodsk. This is a 40x24 deg 2 full-frame image. The black circle shows the error box corresponding to the 3 sigma FERMI error limits plus the 3-degree systematic error (R = 12.6deg) [21]. Animation is available here http://observ.pereplet.ru/images/GRB090320B/GRB090320B 60sec.gif