<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>October</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>HIT FINDER AND TRACK RECONSTRUCTION ALGORITHMS IN THE MULTI-WIRE PROPORTIONAL CHAMBERS OF THE BM@N EXPERIMENT</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>S.P. Merts</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>S.A. Nemnyugin</string-name>
          <email>s.nemnyugin@spbu.ru</email>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>V.A. Roudnev</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>M.M. Stepanova</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>D.P. Usov</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
          <xref ref-type="aff" rid="aff2">2</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Joint Institute for Nuclear Research</institution>
          ,
          <addr-line>Joliot-Curie, 6, 141980 Dubna, Moscow region, Russian Federation</addr-line>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Saint Petersburg State University</institution>
          ,
          <addr-line>7-9 Universitetskaya emb., Saint Petersburg, 199034, Russian Federation</addr-line>
        </aff>
        <aff id="aff2">
          <label>2</label>
          <institution>Sergei Merts</institution>
          ,
          <addr-line>Sergei Nemnyugin, Vladimir Roudnev, Margarita Stepanova, Daniil Usov</addr-line>
        </aff>
      </contrib-group>
      <pub-date>
        <year>2019</year>
      </pub-date>
      <volume>4</volume>
      <issue>2019</issue>
      <fpage>397</fpage>
      <lpage>401</lpage>
      <abstract>
        <p>The BM@N experiment is a part of the NICA accelerator complex at the Joint Institute for Nuclear Research, Dubna. The main goal of the research program on this detector is to study the interactions of relativistic heavy ion beams with fixed targets. BmnRoot software is used to process both experimental and simulated data. Not all the detectors, however, are included in the reconstruction chain. One of them is the Multi-Wire Proportional Chamber (MWPC). Here we present the description of the MWPC detector geometry and the algorithms for digitizing, hit finding and track reconstruction. The results of simulation are given.</p>
      </abstract>
      <kwd-group>
        <kwd>BM@N</kwd>
        <kwd>NICA</kwd>
        <kwd>Multi-Wire Proportional Chamber</kwd>
        <kwd>MWPC</kwd>
        <kwd>track reconstruction</kwd>
        <kwd>hit finder</kwd>
        <kwd>digitizer</kwd>
        <kwd>BmnRoot</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <sec id="sec-1-1">
        <title>NICA (Nuclotron-based Ion Collider fAсility) [1] is a new accelerator complex designed at</title>
        <p>
          the Joint Institute for Nuclear Research (Dubna, Russian Federation) to study properties of dense
baryonic matter. BM@N experiment (Baryonic Matter at Nuclotron) [
          <xref ref-type="bibr" rid="ref2">2</xref>
          ] has the purpose of studying
the collisions of elementary particles and ions with a fixed target. The BM@N setup consists of the set
of coordinate detectors for reconstruction of charge particle trajectories, time-of-flight detectors for
particle identification and calorimeters for energy measurements (Fig. 1).
        </p>
      </sec>
      <sec id="sec-1-2">
        <title>BmnRoot framework [3] is used for the BM@N experiment data processing. It provides a</title>
        <p>
          powerful tool for detector performance studies, event simulation and development of algorithms for
event reconstruction and physics analysis of experimental data registered by the BM@N facility. The
BmnRoot is implemented in the programming language C++ and based on the ROOT [
          <xref ref-type="bibr" rid="ref4">4</xref>
          ] environment
and the object-oriented framework FairRoot [
          <xref ref-type="bibr" rid="ref5">5</xref>
          ].
        </p>
        <p>The detector inclusion in the simulation and reconstruction chain can be divided into several
stages:



creating the detector geometry to describe particle interactions with the detector;
detector digitization – transformation of simulated data into detector signals;
development of particle track reconstruction algorithms.</p>
      </sec>
    </sec>
    <sec id="sec-2">
      <title>2. Multi-Wire Proportional Chamber</title>
      <sec id="sec-2-1">
        <title>Multi-Wire Proportional Chambers are used to determine the trajectory of a particle beam.</title>
        <p>The detector is a chamber filled with a mixture of gases, in which anode wires are located between the
cathode planes (Fig. 2). Upon entry, a charged particle causes ionization of the gas, and the generated
electrons under the action of the electric field of the chamber trigger the wire closest to the particle’s
trajectory.</p>
      </sec>
      <sec id="sec-2-2">
        <title>In the BM@N experiment, the MWPC detector incorporates six flat planes consisting of 96</title>
        <p>wires each. Each plane is rotated by 60° relative to the previous one. The distance between the wires
withinin one plane is dw = 2.5 mm, and the inter-plane distance within in a chamber is 10 mm. In a
Cartesian coordinate system the z axis is perpendicular to the planes. The coordinate system of the
detector is shown in (Fig. 3). In this representation, the planes have the sign "+" or "−" depending on
the increasing or decreasing numbers of the wires along the coordinate axis.</p>
      </sec>
      <sec id="sec-2-3">
        <title>The detectors are located before the magnet, so straight lines can approximate particle trajectories.</title>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>3. Detector geometry</title>
      <sec id="sec-3-1">
        <title>Simulations require a description of the detector geometry, its proportions, materials and</title>
        <p>location. The ROOT geometry package is used for this purpose. It is a tool to build, browse and
visualize detector geometries. Previous version of MWPC geometry described only the active parts
involved in the interactions, but the detector also contains frames and some additional materials where
new particles may appear, and this should be taken into account (Fig. 4).
Proceedings of the 27th International Symposium Nuclear Electronics and Computing (NEC’2019)</p>
      </sec>
      <sec id="sec-3-2">
        <title>We have added aluminum frames around the detector and layers of copper and fiberglass on both sides of the stations to the new version of the geometry. The comparison between the detector models is presented in (Fig.5).</title>
      </sec>
      <sec id="sec-3-3">
        <title>a) Old geometry</title>
      </sec>
      <sec id="sec-3-4">
        <title>b) New geometry Figure 5. Comparison of detector models</title>
      </sec>
    </sec>
    <sec id="sec-4">
      <title>4. Monte Carlo data digitization and hit reconstruction</title>
      <p>The main purpose of the reconstruction procedure is to process the experimental data. For
further tests of track reconstruction, however, processing the simulated data is also important. It
should
be
brought to the same form
with the
experimental data.</p>
      <sec id="sec-4-1">
        <title>For this</title>
        <p>purpose the
BmnMwpcDigitizer class has been developed and incorporated into BmnRoot framework. The
coordinate data obtained during the simulation must be digitized, that is, presented as a discrete set of
wire numbers and detector planes. Based on the geometry of the detector, this can be done by
rounding to integers the following expressions for different planes:   + = 
+ 47.5,   + = 
47.5,   − = 47.5 −</p>
        <p>,   − = 47.5 −
 ,   − = 47.5 −
 ,   + = 47. +</p>
        <p>The next step is creating three-dimensional reconstructed points - hits, on which tracking will
be built. The BmnMwpcHitMaker class has been developed and incorporated into BmnRoot
framework. Tracks are built based on three hits obtained from six values of the detector wire numbers,
 .
 
 
 
+
one from each plane.</p>
      </sec>
      <sec id="sec-4-2">
        <title>The intersections of wires in planes 0-1, 2-3 or 4-5 are converted into a hit, the z coordinate of which is taken as the average between the coordinates of the planes, and the remaining coordinates are calculated based on the detector geometry. Hit 1</title>
        <p>= (  + − 47.5) 
 = (  + − 47.5) 
 =
2 − 
√3
Hit 2
 = (47.5 −   −) 
 = (47.5 −   −
) 
 =
 − 2
√3
Hit 3
 = (47.5 −   −
) 
 = (  + − 47.5) 
 =  + 
 =
2 − 
√3</p>
      </sec>
    </sec>
    <sec id="sec-5">
      <title>5. Track reconstruction algorithm</title>
      <sec id="sec-5-1">
        <title>The particle tracks are built on three hits. This is being performed by the BmnMwpcTracking</title>
        <p>classthat we have developed. At first, the data read from the BmnMwpcHits branch is sorted into three
planes corresponding to different z coordinates. Further, all the hit combinations from different planes
are approximated by straight lines using the least squares method and the best ones are selected from
them by the chi-squared test.</p>
      </sec>
      <sec id="sec-5-2">
        <title>The following functional is being minimized:</title>
        <p>∑

 2
 2 → 
where   is the difference between the hit coordinates and the approximated coordinates in the  -th
plane,  =   ⁄√12 = 0.072 cm is the standard deviation.</p>
      </sec>
      <sec id="sec-5-3">
        <title>If more than 60% of the hits in a track refer to the same simulated track, then the</title>
        <p>reconstruction is considered correct. The efficiency of the algorithm for track reconstruction for 1, 2,
and 3 simulated tracks in an event as a function of the polar angle  is presented in (Fig. 6).</p>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>6. Conclusion</title>
      <sec id="sec-6-1">
        <title>The realistic geometry of the MWPC detector that includes the surrounding materials is</title>
        <p>described. It turned out that the material around the detector has a small effect on the appearance of
new particles in it. Classes for digitizing the simulated data and finding reconstructed hits on which
tracking is built is developed and incorporated into BmnRoot framework. Track reconstruction
algorithm and an algorithm for merging tracks in a two chamber system have been developed.</p>
      </sec>
    </sec>
    <sec id="sec-7">
      <title>7. Acknowledgement</title>
      <sec id="sec-7-1">
        <title>This work is supported by Russian Foundation for Basic Research grant 18-02-40104 mega.</title>
        <p>We are also grateful to the Physics Educational Center of the Research Park of the Saint-Petersburg</p>
      </sec>
      <sec id="sec-7-2">
        <title>State University for support of educational projects related to the subject of the present study.</title>
      </sec>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>V.D.</given-names>
            <surname>Kekelidze</surname>
          </string-name>
          .
          <article-title>NICA project at JINR: status and prospects</article-title>
          .
          <source>Journal of Instrumentation</source>
          , vol.
          <volume>12</volume>
          ,
          <string-name>
            <surname>Iss</surname>
          </string-name>
          .6,
          <issue>2017</issue>
          , pp.
          <fpage>C06012</fpage>
          . https://doi.org/ 10.1088/
          <fpage>1748</fpage>
          -0221/12/06/C06012.
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <given-names>D.</given-names>
            <surname>Baranov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Kapishin</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Mamontova</surname>
          </string-name>
          , G.Pokatashkin,
          <string-name>
            <given-names>I.</given-names>
            <surname>Rufanov</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.</given-names>
            <surname>Vasendina</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Zinchenko. The BM@N Experiment at</surname>
          </string-name>
          <string-name>
            <surname>JINR</surname>
          </string-name>
          :
          <article-title>Status and Physics Program</article-title>
          .
          <source>KnE Energy</source>
          , vol
          <volume>3</volume>
          ,
          <source>Iss. 1</source>
          ,
          <issue>2018</issue>
          , pp.
          <fpage>291</fpage>
          -
          <lpage>296</lpage>
          . https://doi.org/10.18502/ken.v3i1.
          <fpage>1757</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>K.</given-names>
            <surname>Gertsenberger</surname>
          </string-name>
          ,
          <string-name>
            <given-names>S.</given-names>
            <surname>Merts</surname>
          </string-name>
          ,
          <string-name>
            <given-names>O.</given-names>
            <surname>Rogachevsky</surname>
          </string-name>
          ,
          <string-name>
            <given-names>A.</given-names>
            <surname>Zinchenko</surname>
          </string-name>
          .
          <article-title>Simulation and analysis software for the NICA experiments</article-title>
          .
          <source>Eur. Phys. J. A</source>
          .
          <year>2016</year>
          . Vol.
          <volume>52</volume>
          . p.
          <volume>214</volume>
          . https://doi.org/10.1140/epja/i2016- 16214-y.
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          <article-title>[4] CERN</article-title>
          .ROOT. https://root.cern.ch/.
          <source>(Accessed 9 Nov</source>
          <year>2019</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>M.</given-names>
            <surname>Al-Turany</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Bertini</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Karabowicz</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Kresan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>P.</given-names>
            <surname>Malzacher</surname>
          </string-name>
          ,
          <string-name>
            <given-names>T.</given-names>
            <surname>Stockmanns</surname>
          </string-name>
          ,
          <string-name>
            <given-names>F.</given-names>
            <surname>Uhlig</surname>
          </string-name>
          .
          <article-title>TheFairRoot framework</article-title>
          .
          <source>Journal of Physics: Conference Series</source>
          ,
          <year>2012</year>
          , vol.
          <volume>396</volume>
          ,
          <string-name>
            <surname>Iss</surname>
          </string-name>
          .
          <volume>2</volume>
          :
          <fpage>022001</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>V.V.</given-names>
            <surname>Lenivenko</surname>
          </string-name>
          ,
          <string-name>
            <given-names>V.V.</given-names>
            <surname>Palichik</surname>
          </string-name>
          .
          <article-title>Reconstruction of Charged-Particle Trajectories in Multiwire Proportional Chambers at the BM@</article-title>
          N Experiment // Physics of Particles and Nuclei Letters, vol.
          <volume>15</volume>
          ,
          <string-name>
            <surname>Issue</surname>
            <given-names>6</given-names>
          </string-name>
          ,
          <year>November 2018</year>
          , pp
          <fpage>637</fpage>
          -
          <lpage>649</lpage>
          .
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>