C Output information from fiTQun
*                                 IPX=1 ; GAMMA
*                                     2 ; ELECTRON
*                                     3 ; MUON
*                                     4 ; PION
*                                     5 ; KAON
*                                     6 ; PROTON
*                                     7 ; CONE GENERATOR

c  Fit information

      integer ifqver
      real fqproctime

      common /fitquninfo/ ifqver, fqproctime(20)

c  Time window information

      integer maxnpeak
      parameter(maxnpeak=10) ! also hardcoded in atmpd/src/analysis/official_ntuple/ntuple_module_fitqun.F

      real trgoff

      integer fqntwnd
      integer fqtwnd_iclstr
      integer fqtwnd_npeak
      real fqtwnd_prftt0
      real fqtwnd_prftpos
      real fqtwnd

      real fqtwnd_peakt0
      real fqtwnd_peakiness

      common /fitquntwnd/ trgoff, fqntwnd, fqtwnd_iclstr(maxse), fqtwnd_npeak(maxse), fqtwnd_prftt0(maxse)
     + , fqtwnd_prftpos(3,maxse), fqtwnd(2,maxse), fqtwnd_peakt0(maxnpeak,maxse), fqtwnd_peakiness(maxnpeak,maxse)

c  one ring fit information

      integer maxpx
      parameter(maxpx=7) ! also hardcoded in atmpd/src/analysis/official_ntuple/ntuple_module_fitqun.F

      integer fqnse
      integer fqitwnd
      integer fqipeak
      integer fqnhitpmt
      
      real fqtotq
      real fq0rtotmu
      real fq0rnll

      integer fqn50
      real fqq50

      integer fq1rpcflg
      real fq1rmom
      real fq1rt0
      real fq1rtotmu
      real fq1rnll
      real fq1rpos
      real fq1rdir
      real fq1rdconv
      real fq1reloss

      common /fitqun1r/ fqnse, fqitwnd(maxnpeak), fqipeak(maxnpeak), fqnhitpmt(maxnpeak), fqtotq(maxnpeak)
     + , fq0rtotmu(maxnpeak), fq0rnll(maxnpeak), fqn50(maxnpeak), fqq50(maxnpeak)
     + , fq1rpcflg(maxpx,maxnpeak), fq1rmom(maxpx,maxnpeak), fq1rt0(maxpx,maxnpeak)
     + , fq1rtotmu(maxpx,maxnpeak), fq1rnll(maxpx,maxnpeak), fq1rpos(3,maxpx,maxnpeak)
     + , fq1rdir(3,maxpx,maxnpeak), fq1rdconv(maxpx,maxnpeak), fq1reloss(maxpx,maxnpeak)

c  pi0 fit information

      integer fqpi0pcflg
      real fqpi0mom1
      real fqpi0mom2
      real fqpi0momtot
      real fqpi0dconv1
      real fqpi0dconv2
      real fqpi0t0
      real fqpi0totmu
      real fqpi0nll
      real fqpi0mass
      real fqpi0photangle
      real fqpi0pos
      real fqpi0dir1
      real fqpi0dir2
      real fqpi0dirtot

      common /fitqunpi0/ fqpi0pcflg(2)
     + , fqpi0mom1(2), fqpi0mom2(2), fqpi0momtot(2)
     + , fqpi0dconv1(2), fqpi0dconv2(2), fqpi0t0(2)
     + , fqpi0totmu(2), fqpi0nll(2)
     + , fqpi0mass(2), fqpi0photangle(2)
     + , fqpi0pos(3,2), fqpi0dir1(3,2), fqpi0dir2(3,2), fqpi0dirtot(3,2)


c  MR fit information

      integer fqmaxnmrfit
      parameter(fqmaxnmrfit=200) ! also hardcoded in atmpd/src/analysis/official_ntuple/ntuple_module_fitqun.F

      integer fqmaxnring
      parameter(fqmaxnring=6) ! also hardcoded in atmpd/src/analysis/official_ntuple/ntuple_module_fitqun.F

      integer fqnmrfit
      integer fqmrifit
      integer fqmrnring
      integer fqmrpcflg
      real fqmrnll
      real fqmrtotmu

      integer fqmrpid
      real fqmrmom
      real fqmrdconv
      real fqmreloss
      real fqmrt0
      real fqmrpos
      real fqmrdir

      common /fitqunmr/ fqnmrfit, fqmrifit(fqmaxnmrfit), fqmrnring(fqmaxnmrfit), fqmrpcflg(fqmaxnmrfit)
     + , fqmrnll(fqmaxnmrfit), fqmrtotmu(fqmaxnmrfit), fqmrpid(fqmaxnring,fqmaxnmrfit)
     + , fqmrmom(fqmaxnring,fqmaxnmrfit), fqmrdconv(fqmaxnring,fqmaxnmrfit), fqmreloss(fqmaxnring,fqmaxnmrfit)
     + , fqmrt0(fqmaxnring,fqmaxnmrfit), fqmrpos(3,fqmaxnring,fqmaxnmrfit), fqmrdir(3,fqmaxnring,fqmaxnmrfit)


c  MS fit information

      integer fqmsmaxring
C     ALSO HARDCODED IN NTUPLE MODULE
      parameter(fqmsmaxring=20) 

      integer fqmsmaxfit
C     ALSO HARDCODED IN NTUPLE MODULE
      parameter(fqmsmaxfit=5)

      integer fqmsnfit

c     for each ms fit
      integer fqmsnseg
      integer fqmspid
      integer fqmsifit
      integer fqmsimer
      integer fqmspcflg
      real fqmstotmu   
      real fqmsnll

c     for each segment
      real fqmsmom
      real fqmseloss
      real fqmst0
      real fqmspos
      real fqmsdir

      common /fitqunms/ fqmsnfit, fqmspcflg(fqmsmaxfit), fqmsnseg(fqmsmaxfit), fqmspid(fqmsmaxfit),
     +       fqmsifit(fqmsmaxfit), fqmsimer(fqmsmaxfit), fqmstotmu(fqmsmaxfit), fqmsnll(fqmsmaxfit),
     +       fqmsmom(fqmsmaxring,fqmsmaxfit), fqmseloss(fqmsmaxring,fqmsmaxfit),
     +       fqmst0(fqmsmaxring,fqmsmaxfit), fqmspos(3,fqmsmaxring,fqmsmaxfit),
     +       fqmsdir(3,fqmsmaxring,fqmsmaxfit)


c  2electron fit information
      
      real fq2ElecMom1
      real fq2ElecMom2
      real fq2Elect0
      real fq2Electotmu
      real fq2Elecnll
      real fq2Elecpos
      real fq2Elecdir1
      real fq2Elecdir2

      common /fitqun2Elec/ fq2ElecMom1, fq2ElecMom2, fq2Elect0
     + , fq2Electotmu, fq2Elecnll, fq2Elecpos(3), fq2Elecdir1(3)
     + , fq2Elecdir2(3)

      integer fqtestmax1r
      parameter(fqtestmax1r=30) ! also hardcoded in atmpd/src/analysis/official_ntuple/ntuple_module_fitqun.F

      integer fqtestmaxpi0
      parameter(fqtestmaxpi0=30) ! also hardcoded in atmpd/src/analysis/official_ntuple/ntuple_module_fitqun.F

      integer fqtestn1r
      integer fqtest1rstage
      integer fqtest1rse
      integer fqtest1rpid
      integer fqtest1rpcflg
      real fqtest1rmom
      real fqtest1rt0
      real fqtest1rtotmu
      real fqtest1rnll
      real fqtest1rpos
      real fqtest1rdir
      real fqtest1rdconv
      real fqtest1reloss

      common /fitqun1rtest/ fqtestn1r, fqtest1rstage(fqtestmax1r)
     + , fqtest1rse(fqtestmax1r), fqtest1rpid(fqtestmax1r)
     + , fqtest1rpcflg(fqtestmax1r)
     + , fqtest1rmom(fqtestmax1r), fqtest1rt0(fqtestmax1r)
     + , fqtest1rtotmu(fqtestmax1r), fqtest1rnll(fqtestmax1r)
     + , fqtest1rpos(3,fqtestmax1r), fqtest1rdir(3,fqtestmax1r)
     + , fqtest1rdconv(fqtestmax1r), fqtest1reloss(fqtestmax1r)

      integer fqtestnpi0
      integer fqtestpi0stage
      integer fqtestpi0pcflg
      real fqtestpi0mom1
      real fqtestpi0mom2
      real fqtestpi0momtot
      real fqtestpi0dconv1
      real fqtestpi0dconv2
      real fqtestpi0t0
      real fqtestpi0totmu
      real fqtestpi0nll
      real fqtestpi0mass
      real fqtestpi0photangle
      real fqtestpi0pos
      real fqtestpi0dir1
      real fqtestpi0dir2
      real fqtestpi0dirtot

      common /fitqunpi0test/ fqtestnpi0, fqtestpi0stage(fqtestmaxpi0)
     + , fqtestpi0pcflg(fqtestmaxpi0)
     + , fqtestpi0mom1(fqtestmaxpi0), fqtestpi0mom2(fqtestmaxpi0)
     + , fqtestpi0momtot(fqtestmaxpi0), fqtestpi0dconv1(fqtestmaxpi0), fqtestpi0dconv2(fqtestmaxpi0)
     + , fqtestpi0t0(fqtestmaxpi0), fqtestpi0totmu(fqtestmaxpi0), fqtestpi0nll(fqtestmaxpi0)
     + , fqtestpi0mass(fqtestmaxpi0), fqtestpi0photangle(fqtestmaxpi0), fqtestpi0pos(3,fqtestmaxpi0)
     + , fqtestpi0dir1(3,fqtestmaxpi0), fqtestpi0dir2(3,fqtestmaxpi0), fqtestpi0dirtot(3,fqtestmaxpi0)

c  proton decay to mu gamma fit information

      integer fqpmgpcflg
      real fqpmgmom1
      real fqpmgmom2
      real fqpmgt01
      real fqpmgt02
      real fqpmgtotmu
      real fqpmgnll
      real fqpmgpos1
      real fqpmgpos2
      real fqpmgdir1
      real fqpmgdir2

      common /fitqunpmg/ fqpmgpcflg(2)
     + , fqpmgmom1(2), fqpmgmom2(2), fqpmgt01(2), fqpmgt02(2)
     + , fqpmgtotmu(2), fqpmgnll(2)
     + , fqpmgpos1(3,2), fqpmgpos2(3,2), fqpmgdir1(3,2), fqpmgdir2(3,2)