High efficiency gaseous tracking detector for cosmic muon radiography

A tracking detector system has been constructed with an innovative approach to the classical multi-wire proportional chamber concept, using contemporary technologies. The detectors, covering an area of 0.58 square meters each, are optimized for the application of muon radiography. The main features are high (>99.5%) and uniform detection efficiency, 9 mm FWHM position resolution, filling gas consumption below 2 liters per hour for the non toxic, non flammable argon and carbon dioxide mixture. These parameters, along with the simplicity of the construction and the tolerance for mechanical effects, make the detectors to be a viable option for a large area muography observation system.


Introduction
Cosmic muons are highly penetrating particles, which are created by interaction of high energy cosmic particles impinging on the upper atmosphere of Earth. Muons are decay products of pions, which make up the majority of the particles reaching the Earth surface. Production mechanisms and interaction properties are extensively discussed in the literature [1]. The mean energy of muons is around 2 GeV, and such an average muon would be able to cross 3 -4 meters of rock. This considerable penetrating power gives rise to the practical applicability of cosmic muon radiography, or "muography" in short.
There is a broad range of applications of muography. Very large objects can be imaged by measuring muon absorption, such as in the classical measurement by Alvarez et al [2]. Notable, high impact examples are volcano imaging [3,4,5] and mapping the core of nuclear reactors [6]. Underground measurements have smaller background due to the natural shielding, and can be exploited during ore mining [7] or natural cave search [8,9,10,11]. If the object of interest is not exceedingly large, measuring the muon scattering [12,13,14,15] can provide information about the internal structure.
When it comes to particle detectors for muon radiography, one of the key considerations is the flux of the signal particles. On the surface, the rate is in the order of one hundred particles per square meter per second -a relatively low value compared to high energy physics (HEP) experimentation. On the other hand, if a large object is to be imaged, the flux of the muons can be drastically reduced further [3]: this implies that large size detectors are mandatory for most meaningful measurements. An additional important point is the suppression of background against the relevant (transmitted) muon signal: this background turns out to be mostly low energy particles from various sources [16,17].
From the above it follows that the requirements for muography detectors differ from those generally expected from high energy physics detectors (for the latter, see textbooks such as Ref. [18]). The key points may be summarized as follows: • Large size at sensible cost. Due to the low and fixed cosmic rate, the statistics becomes the ultimate limiting factor of the measurements on large objects.
• Long term operation under ambient outdoor environmental conditions. Considerable temperature, pressure and humidity variations need to be tolerated.
• Tracking capability. In typical HEP conditions, the detector size is usually constrained by available space. On the contrary, most applications of muography require a specific angular resolution, which by increasing the measurable track length, allows for weaker position resolution. However, clear identification of the straight muon tracks is mandatory to suppress low energy background.
• Rate capability. The cosmic particle rate is low on the Earth surface, therefore none of the issues apply which characterize "high rate" detectors in HEP. Eliminating those design elements which are meant to improve rate capability potentially results in considerable simplification of the structure.
• Sustainable and safe operation. Given the long continuous measurement times, low maintenance and low running costs are desirable, whereas outdoor installation requires safe and environmentally friendly operation.
Various detector technologies have been considered for muon radiography. The three notable groups are scintillators, emulsion and gaseous detectors. So far scintillator-based trackers proved to be the most successful, applied for vulcanology [19,20,21] and mining [7]. Emulsions present very high position resolution [22], and hence the detector layers require small space [23,24,16]. The main drawback of this technology is that the information can be made available only after the film development, making real time measurements impossible. Gaseous detectors have the possibility to reach relatively large surfaces at reasonable cost, and therefore these tend to be the choice for the outer layers of most HEP experiments. Various groups proposed viable solutions with Resistive Plate Chambers [25], drift chambers [13,14] and GEMs [15].
Our objective was to conceive a tracking detector system which is optimized for muography from the very beginning of the design. Based on our own experience with various gaseous detectors, an innovative approach to the classical Multi Wire Proportional Chambers (MWPC-s) [18] was chosen. In smaller sizes, MWPC-s with simplified geometry may be used for educational purposes [28].
In order to address the "muography-requirements" above, design features include consistent material choice (Section 2) for improved thermal behaviour, an internal support structure with causes little efficiency loss (Section 7) and improves mechanical stability (8). As working gas, the mixture of argon and carbon dioxide was studied in detail (Section 10), being the most obvious choice from cost, environmental and safety considerations.
The paper presents the detector construction, followed by relevant performance studies, for a sensitive area of about 0.58 m 2 . The key guiding principle was the simplification of the design, without compromise on detection efficiency and tracking capability: our studies below attempt to give a comprehensive performance report.

Construction
The detector uses a relatively large, 22 mm gas gap between two grounded cathode planes. The gas amplification takes place on an anode (sense, SW) wire plane nearly half way between the cathodes. In the same plane, field shaping wires (FW) are fixed as well parallel to the anode wires, as shown in Figure 1. There is an other wire plane, for the so called "pick-up wires" (PW), which run perpendicularly to the anode wires. The concept is similar to the CCC case [26,27]: the anode wires are used as trigger and amplitude measurement, whereas the field shaping wires and the pick-up wires give the position information. The pick-up wire plane is 2 mm above one of the cathodes, as shown in Fig. 1. The advantage of using such pick-up wires is cost consideration and simplification of the structure: an alternative would be using printed circuit boards (PCB-s) with strips, however PCB-s at large sizes tend to be expensive and their surfaces are not sufficiently flat. The anode wires for the specific detectors presented in this paper are 22 and 25 micrometers thick Au-plated tungsten, whereas both FW and PW are brass wires of 100 micrometers diameter. The two different anode wire diameters presented no practical difference except for a slight change in operational voltage.
In one detector panel, there are 64 field shaping wires (hence 63 sense wires), with wire spacing of 12 mm, covering a 768 mm long region. The pick-up wires cover the same length with 4 mm spacing (hence there are 192 of them), which are connected in groups of 3 (that is, each of the 64 pick-up wire groups are 12 mm wide). All the wires are approximately 770 mm long. The front view outline is shown in Fig. 2.  From the electronic connection point of view, the 64 field wires and the 64 pick-up wire groups are directly connected to the readout (see in Sec. 4), whereas the 63 sense wires are all interconnected. From the sense wire signal, two output signals are formed after pre-amplification and shaping: one is a TTL-compatible discriminated trigger signal, and the other is a direct analog output for amplitude measurement.  The detector is constructed out of fiberglass (glass fiber reinforced epoxy). Specifically the sidewalls are G-10, whereas the cathode planes and other printed circuit components (soldering points, electronics support) are standard FR-4. This way, all parts of the detector feature a consistent thermal behavior. The cathode planes are 1.5 mm thick FR-4 layers, fully covered with copper on both sides. The complete detector is glued with epoxy: the units are thus not repairable in case of broken wires (however, no such case was observed so far after a year of operation). The total weight of a detector panel is about 6 kg including electronics.
Given the relatively large area, the distance between the cathode planes need to be kept at the nominal 22 mm with reasonable precision. In order to achieve this, internal support pillars were installed in a particular way between the wires, as shown in Fig. 1. The pillars have a rectangular 3mm by 4mm cross section around the anode wire plane. These pillars, made out of insulating material (PLA) using 3D printing technology, do not touch any of the field and sense wires. The pillars were placed closer to the field wires, with around 0.5 mm clearance, which leaves approximately 2.5 mm distance from the sense wire. One can expect that near a support pillar, the field lines still terminate on the sense wires, therefore the detector remains sensitive to the traversing particles -a feature which proves to be correct according to the studies presented in Sec. 7. For our particular detectors, 6 such pillars have been installed, according to Fig. 2. The possibility of installing many of such pillars without introducing a large dead zone in the detector is of key importance for mechanical stability.
The detectors were operated with mixture of argon and carbon dioxide, in proportion of 82:18, at various flow rates ranging from 1 to 10 liters per hour. Detailed studies on gas behavior is presented in Section 10.

Tracking detector design
A tracking system has been set up, using six individual detector layers, arranged each with 20 cm spacing (that is, 1 m distance between the first and the last layer). The detectors were mounted vertically, standing on the side without connectors (bottom side in Fig. 2), that is, the anode wires were horizontal. The detectors are counted from 1 to 6 for the rest of the paper. For most of the measurements, the detectors shared the same anode high voltage. The gas line was connecting the layers in series: the gas entered to detector 1, and was vented to the atmosphere from detector 6.
The setting with vertical detectors resembles most to the application in volcanology, where the object is imaged from the "side", and where, because of the reduced muon flux close to the horizon, the system is most sensitive to background. Such dedicated issues are however outside of the scope of the present paper: for the tracking system we did not use absorbers, which are otherwise mandatory for such arrangements [5], and will be subject to later detailed studies. The image of the actual system with vertical detectors is shown in Fig. 3, including two additional non-functional layers. For various performance studies, we applied a similar arrangement, but with horizontal detectors inside the laboratory: that is, capturing the most abundant nearly vertical muons, high statistics mappings of the detector parameters became possible.
The six detector layers are sufficient to identify clear particle tracks, and allow for the actual detector performance studies. We have used tracking algorithms from our own earlier systems [10]. From the tracking information one can tune the detector geometry (by requiring straight tracks on average), which is needed for a more precise track finding. The tracking system has shown a consistently excellent performance, enabling a reliable characterization of the individual detector layers using cosmic muons.

Readout and data acquisition
The detector readout was adapted from our own earlier system [11], with a single Raspberry Pi microcomputer as the core of the data acquisition [29]. An event sequence starts with the arrival of the trigger signals: if three coincident triggers appear (from any of the 6 layers), the readout will wait for a pre-defined delay time, and then issue a "master" trigger pulse. At the instance of this master trigger, the position information is recorded from discriminated signals of all the channels of the field wires and the pick-up wires. Also the master trigger time defines when the sense wire amplitude measurement, via 10 bit ADC-s, takes place.
The total data from the tracking system is moderate: there are 6×2×64 = 768 bits from the position information, and including the converted ADC values as well as recording the trigger pattern, the raw data is still around 100 Byte per event. This explains why a single Raspberry Pi can handle the DAQ task.
The input signals from the wires are protected from external electronics noise as much as possible, and since the grounded cathodes act as a sort of Faraday cage, a reasonable noise performance has been observed. The trigger level from each of the layers are set approximately to 2 × 10 5 electrons (32 fC) on input; for the field shaping and pick-up wires the discrimination threshold was also set to 2 × 10 5 electrons (measured as charge on the sense wires).

Efficiency of cosmic muon detection
A key figure of merit of a detector is the tracking detection efficiency, that is, the probability that a valid particle hit (position) is actually detected. This we define in the following way for the rest of the paper: considering a specific detector layer, a particle trajectory is formed from the 5 detectors excluding the one under study. If this 5-point pattern is a unique, straight track from both directions (see Section 6), one checks if there is a hit at the expected position (within 5 sigma of the position resolution) in the detector under study. The probability of finding the hit is defined as the "tracking efficiency". In a similar way, one can check if the specific detector layer provided a trigger signal for the event: this defines is the "trigger efficiency". Trigger and tracking efficiencies are found to be well above 90% for the operational conditions. First, the tracking and trigger efficiency has been determined as a function of the sense wire voltage. A clear plateau is observed, shown in Fig. 4. The efficiency reaches around 95% (>99%) at 1650 V (1700 V) for the detectors with 25 micron anode wires. The plateau continues to 1800 V, where stable operation has still been observed.
The efficiency values shown in Fig. 4 are consistent with the expectation from the signal amplitude spectrum, shown in Figure 5 for a sense wire voltage of 1700 V for the first tracking layer. The spectrum of all tracks with at least 5 points are overlayed with the spectrum of 6-point tracks. The difference between these two sets, in case of full efficiency, is just those tracks which miss the first layer and hence result in a 5-point track. The approximate detection level is 120 ADC units which has also been indicated with an arrow. The horizontal scale of Fig. 5 is expressed in ADC units, which have been calibrated by direct charge injection: 1 ADC unit corresponds to 1.8 ×10 3 e, or 0.3 fC.
The true amplification gain, that is, the multiplication factor of the avalanches, has been estimated measuring the mean amplitude on the sense wires. Taking into account that the mean number of electrons over the 2 cm of the drift gap is around 200 in the applied gas mixture [18], one finds that the detectors with 25 micron anode wires at 1700 V (full efficiency) has a gain of about 7 ×10 3 , a conveniently low figure.
The trigger signal is created at the point when the input signal reaches the trigger threshold. There is a certain time delay needed after this to reach the highest amplitude on the field wire and pick-up wire discriminators, and hence to reach optimal tracking efficiency. The tracking efficiency as a function of this trigger delay is shown in Figure 6, with a maximum around 2-3 µs consistent with the peaking of the signal. The data has been taken keeping detectors 1,2,5,6 at 1700 V and reducing the voltage on the middle detectors, 3 and 4, to the values indicated (1600 V and 1550 V). This procedure ensures reasonably high trigger efficiency (from four detectors) and a measurable efficiency (that is, well below 100%) for the inner layers. The setting for all other measurements use a fixed 2 µs delay.

Position resolution
With the 6-layer tracking detector setup, the position resolution of the layers can be quantified. Figure 7 shows the distance between the measured hit in one of the middle layers and the extrapolation from the 5-point track (that is, the U sw = + 1600 V: 3 U sw = + 1600 V: 4 U sw = + 1550 V: 3 U sw = + 1550 V: 4 Figure 6: Signal timing curves: the tracking efficiency of the two inner layers as a function of the delay relative to the arrival of the coincident trigger (mostly from the four other detectors, operated using high efficiency settings, see text). The measurement is consistent with the observed signal peaking time of 2 µs, the value used as the "standard" delay setting. fitted track excluding the specific layer). One finds that the position resolution (RMS) is σ FW =3.73 mm in the field wire, and σ PW =3.92 mm in the pick-up wire direction, that is, around 9 mm FWHM. This is roughly consistent with the 12 mm spacing expectation (12mm/ √ 12 = 3.5mm). The tracking algorithm incorporates this position resolution to find a set of points on 5 or 6 of the chambers which form a straight line from both directions.

Dead zones due to support pillars
The detector, as seen above, is highly efficient except for the support pillars, which deserve a closer look. In the tracking system configuration with horizontal detectors, the local efficiency has been determined, as shown in Figure  8. For this measurement, nearly vertical tracks were selected, with less than 200 mrad relative to the zenith. The pixel size in this map is one half of the readout unit (6 mm), and one immediately notes the surprisingly small reduction. (The measured efficiency does not go to zero due to the fact that the pillars are smaller than the tracking detector position resolution). Taking a 3 cm wide strip over the pillars (indicated in the Figure by two dashed lines), the mean efficiency is indicated in Figure 9. To first order, one can conclude that the detector is sensitive for any muons which enters the gas volume around the pillars, or in other words, the presence of the support pillars does not distort the field structure to the extent to considerably reduce efficiency. There seems to be some fluctuation on the apparent loss caused by the pillars, which may have to do with the details of the actual positioning.

Mechanical stability
The mechanical stability and required tolerances of the detector is a key question for practical outdoor applications. In the studies below, some of the mechanical effects on the detectors are quantified in relation to its performance. Despite the simplicity of the structure, it turned out that most of the relevant mechanical variations (e.g. difference between vertical or horizontal position) give insignificant changes, therefore we decided to expose the detectors to considerable mechanical effects, outside of the normally expected conditions.
The gain map, that is, the mean signal amplitude over the detector surface under reference conditions is shown in the left panel of Figure 10. In the first study, the detectors were ran with an increased inner gas pressure of 1.5 mbar above the ambient atmosphere. This results in a slight bulging; the inner support pillars in this case act against the forces on the cathode layers. The corresponding gain map is shown in the middle panel of Figure 10, demonstrating that the gain drop under this condition is below 10% (and hence, under normal operating conditions, one expects practically no efficiency loss).
The second study involved placing a weight directly on the detector. In this case, the detector was supported at the four corners, and two weights, 1kg each, was placed at the positions indicated in Figure 10. One position was between the pillars, the other was approximately on one of the pillars. The gain increase (due to the reduction of the distance between the cathode layers) is apparent in the first case (by about 15%), whereas at the position of the pillar, the change is smaller. There is no effect observable outside of the weight positions (e.g. no global detector bending).  One can conclude that the mechanical structure, despite the simplicity, is very robust: the internal support pillars ensure the constant distance between cathode layers, and hence allow for acceptable low gain variation.

Rate capability
In the standard setting, the tracking system has been operated in outdoor conditions on ground level, which is the maximum available cosmic muon rate. In order to quantify if a higher local rate is tolerable without efficiency loss, a 90 S r beta source has been placed close to one of the layers. The native trigger count rate of the specific detector layer increased by 3.5 kHz (from around 1 kHz without source in horizontal position) from this concentrated irradiation covering about 20 cm 2 . Correspondingly, the number of hits appearing in triggered events (events in the sense of the tracking system) increases locally over the homogeneous distribution dominated by cosmic events, as shown in Fig. 11. Figure 11: (Color online) Distribution of all hits identified in the detector layer irradiated with the source. The local rate increase is apparent relative to the otherwise homogeneous distribution of cosmic-induced hits.
The tracking efficiency has been determined along strips of 36 mm width (indicated by the dashed lines in Figure  11) in both directions. The data was taken at an anode voltage of 1700 V, that is, at efficiency of 99.5 %. According to Figure 12, there is no efficiency loss observed for the source region (shown with the peaked continuous line).  The count rate by cosmic muons is in the order of 0.1 Hz/cm 2 , whereas in our experiment the induced local rate is above 100 Hz/cm 2 , corresponding to a three orders of magnitude increase. From the result, one can conclude that rate capability limitations, if any, are not relevant for the presented detector design in cosmic radiography applications.

Required gas flow and tolerance on gas purity
The detector runs with industrial grade argon and carbon dioxide mixture in 82:18 proportion (Ferroline C18). The total volume of a single detector panel is around 13 liters, which implies long gas equilibration times for typical operational gas flow rates of 1-3 liters per hour.
In order to understand the behavior of the gas system, we have performed two studies. One of these involved direct injection of ambient air into the gas line, whereas the other was part of the realistic outdoor studies presented in Section 12. In the first case, the gas flow was set to 1 liter per hour, and the detectors were connected in series to the gas line (sequentially from 1 to 6). Using a gas bypass between detectors 1 and 2, 166 cm 3 ambient air has been injected into the gas line, which then entered the second detector and later gradually moved along the subsequent panels. This way, the signal amplitude change relative to the first detector quantifies the effect of air (oxygen) contamination. Figure 13 shows the mean anode wire signal amplitude in the layers as a function of time t, normalized to the mean signal in detector 1 as reference before the gas injection (at t = 0 in the Figure). The pattern is very clear, with a sharp decrease in the second detector, followed with an exponential purging.  Assuming complete mixing of the air bubble in the gas volume, indeed an exponential relaxation is expected for the mean signal amplitude A, described by where the time constant T is fixed by the ratio of the detector volume and the flow rate and found to be T = 12.5 hours. The only free parameter C, the relative gain drop, is determined to be C= 34 %. The function is shown in Figure  13 with the continuous line. In fact, the assumption of "complete mixing" predicts the gain drop in the subsequent layers, without any free parameters: for k = 2, 3, 4, 5, 6, where k is the number of the detector layer, and with k = 2 we get back the pure exponential behavior. Figure 13 confirms this behavior, indicated with the lines: the air bubble walks along the detectors.
The gain drop of 34% corresponds to 1.3 % (13000 ppm) air contamination; or assuming that only the oxygen content of air counts, 2500 ppm of O 2 contamination level. One can conclude that the gas distribution inside the detector system is fully predictable, with the fixed time constant T determined as the detector volume filling time. Oxygen contamination of 1000 ppm, leading to less than 15% gain drop (and hence no efficiency loss) thus seems to be tolerable, which is indeed a relaxed requirement for gas quality.

Background and noise performance
The instrumental background induced by the detector system is of vital importance from the point of view of practical muon radiography. Vaguely, we can call "background" anything which looks like a muon that traversed the object under study, while it actually did not. The most relevant component of such backgrounds is given by (low energy) particles which actually reach the detector [17]. This component can only be eliminated by a careful detector design and will be subject of future detailed studies; however, at this point we aim at quantifying instrumental artifacts (background which may not even be connected to any particle depositing charge in the detector). One possibility is to apply a "random" (uncorrelated oscillator) trigger and observe patterns resembling trajectories. In the first place, it has been verified that the pure "electronics" noise is exceedingly small: with the high voltage off, simply there are no triggers, and there are no observable trajectories. The probability of any of the single channels (on FW or PW) firing is in the order of 10 −6 , that is, the chance to have 5 or 6 of them in a line to fake a track is negligible. In order to have a more relevant background estimate, the tracking system was operated in a fully efficient configuration (U S W = 1700V) with vertical detectors, and was running for a day with oscillator trigger at 110 Hz. This data contains reconstructed tracks, and allows to estimate the "background" flux f BG . The true particle flux f can be calculated with standard procedures [9], and is proportional to the number of observed tracks per event: where the number of tracks N are measured over time of t, whereas the number of triggered events N trig divided by t is the trigger rate R trig . N/N trig is the mean track number per event. Most triggered events contain only short tracks (due to the weak trigger condition), therefore the background tracks N BG add up to the "true" track set. The probability of random tracks, per event (taken with the oscillator trigger), can be used to estimate the background flux f BG : where N osc is the event number taken with the oscillator (random) trigger. The background flux f BG has been evaluated as a function of elevation angle and is shown in Figure 14. There is a clear minimum at the horizon (elevation angle 0), reaching a background flux value of around 2 × 10 −4 m −2 sr −1 s −1 . The background flux happens to be around 4 orders of magnitude below the true measured flux f (also indicated in the Figure). One has to note that this measured flux has nothing to do with the flux of cosmic muons (which should reduce to far lower values near the horizon): the measurement includes a considerable and undetermined low energy background.
The key conclusion of the observations above is that there is no apparent uncorrelated track background beyond that expected from the truly random appearance of otherwise physical tracks; such tracks are expected to be suppressed along with the low energy background component in true muography tracking systems. Measured Background Figure 14: The measured background flux f BG , using random (oscillator) trigger. It is found to be nearly four orders of magnitude below the "triggered" particle flux f (upper curve), and has a similar shape (indicated with the line) with a clear minimum at the horizon.

Experience under ambient outdoor conditions
For testing under ambient outdoor conditions, the tracking system with vertical detectors have been installed on a support frame and placed on the flat roof of a building at the Wigner Research Centre campus. To allow for considerable environmental influence, the system was only lightly covered with a box, thus receiving strong sunlight in the early afternoon. In Figure 15, the outside temperature is shown on the top panel, featuring more than 25 • C daily temperature variations. During the 6 days of continuous data taking in our example, the system sustained a light raining on the 5th day, otherwise the weather was clear.
The change of the temperature by 30 • C corresponds to nearly 10% change in absolute temperature, and hence in gas density. The system thus can not run at constant anode voltage, but needs to be compensated to first order: in our system, the anode HV was reduced by 1 V for a 1 • C temperature increase (from the nominal 1700 V at 25 • C). Though the compensation is far from perfect (possibly containing non linear terms), it has been verified that this ensures only a slight, 10% mean amplitude variation (demonstrated on the middle panel).
The key figure of merit for the detector is the tracking efficiency, shown in the lower panel of Figure 15. The measured tracking efficiency is reliably high (>99.5%) and stable for the whole period (shown only for two detector layers but found consistently for all panels).
The detector was running with a constant, 2 liters per hour gas flow rate. This presents an interesting issue during decreasing temperatures. The total tracking system gas volume was about 100 liters. Assuming that the temperature decreases by 6 • C in one hour, the 2% relative change in absolute temperature causes the total gas volume to reduce accordingly, that is, by 2 liters. In this case, actually air may be sucked in from the outside atmosphere to the last detector. In order to avoid this, a long buffer tube has been added after the last detector, with sufficiently large (in our case 1 liters) total volume.

Summary and conclusions
A new approach to the classical MWPC-s has led to a highly efficient, low weight, mechanically and environmentally tolerant tracking detector, which by its simple and cost efficient construction, can be an ideal part of large size  Figure 15: Time dependence of relevant operational parameters during 6 days of data taking in outdoor conditions. Temperature variations up to 25 • C has been sustained (top panel), with not more than 10% mean amplitude changes (due to partial compensation of the anode HV). The tracking efficiency has been consistently high, shown in the lower panel, during the whole period (note the scale on the lower panel starts at 99.5%). muon radiography detectors. An internal support pillar structure has been shown to present a negligible dead zone.
The tracking system has been demonstrated to tolerate broad temperature variations maintaining very high detection efficiency (>99.5%). The position resolution of better than 4 mm RMS (9mm FWHM) is compatible with the segmentation. Gas flow rate using non-flammable, non-toxic Ar and CO 2 mixture of 1-2 liters per hour is sufficient for the operation, with a well predictable gas flow pattern. Future complete muography tracking systems will be completed with absorber (scatterer) layers and are intended to be applied for large size object imaging, such as the internal structure of volcanos.