#include <iostream>
#include <fstream>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#include "spliTChanOutC.h"
#include "fitqunoutC.h"

using namespace std;

extern "C"
{
  void clrfqcmns_(int *);
}

void clrfqcmns_(int *iClrMod) {
//   clear all fiTQun and spliTChan common block variables
  
  cout << "Clearing fiTQun common block variables..." << endl;
  
  if (*iClrMod>0) {
    cout << "Preserving subevent variables..." << endl;
  }
  else {
    
    for (int ise=0; ise<MAXSE; ise++) {
      fitqunse_.cluster_tstart[ise] = 0;
      fitqunse_.cluster_tend[ise] = 0;
      fitqunse_.cluster_nhits[ise] = 0;
      fitqunse_.cluster_totq[ise] = 0;
      fitqunse_.cluster_goodflag[ise] = 0;
      
      for (int imethod=0; imethod<NPEAKFINDS; imethod++) {
        for (int ipeak=0; ipeak<MAXSE; ipeak++) {
          fitqunse_.cluster_ipeak[ise][imethod][ipeak] = 0;
          fitqunse_.cluster_timeofpeak[ise][imethod][ipeak] = 0;
        }
        fitqunse_.cluster_npeaks[ise][imethod] = 0;
      }
    }
    fitqunse_.cluster_ncand = 0;
    
    for (int imethod=0; imethod<NPEAKFINDS; imethod++) {
      
      for (int ipeak=0; ipeak<MAXSE; ipeak++) {
        fitqunse_.muechk_tpeak[imethod][ipeak]    = 0;
        fitqunse_.muechk_bg[imethod][ipeak]       = 0;
        fitqunse_.muechk_mean[imethod][ipeak]     = 0;
        fitqunse_.muechk_excess[imethod][ipeak]   = 0;
        fitqunse_.muechk_signif[imethod][ipeak]   = 0;
        fitqunse_.muechk_icluster[imethod][ipeak] = 0;
      }
      fitqunse_.muechk_ncand[imethod] = 0;
      fitqunse_.muechk_toya[imethod] = 0;
    }
    
    for (int ise=0; ise<MAXSE; ise++) {
      fitquntwnd_.fqtwnd_iclstr[ise] = 0;
      fitquntwnd_.fqtwnd_npeak[ise] = 0;
      fitquntwnd_.fqtwnd_prftt0[ise] = 0;
      for (int ix=0; ix<3; ix++) {
        fitquntwnd_.fqtwnd_prftpos[ise][ix] = 0;
      }
      fitquntwnd_.fqtwnd[ise][0] = 0;
      fitquntwnd_.fqtwnd[ise][1] = 0;
      for (int ipeak=0; ipeak<MAXNPEAK; ipeak++) {
        fitquntwnd_.fqtwnd_peakt0[ise][ipeak] = 0;
        fitquntwnd_.fqtwnd_peakiness[ise][ipeak] = 0;
      }
    }
    fitquntwnd_.trgoff = 0;
    fitquntwnd_.fqntwnd = 0;
  }
  
  for (int ise=0; ise<MAXNPEAK; ise++) {
    for (int ipx=0; ipx<MAXPX; ipx++) {
      fitqun1r_.fq1rmom[ise][ipx] = 0;
      fitqun1r_.fq1rt0[ise][ipx] = 0;
      fitqun1r_.fq1rtotmu[ise][ipx] = 0;
      fitqun1r_.fq1rnll[ise][ipx] = 0;
      for (int ix=0; ix<3; ix++) {
        fitqun1r_.fq1rpos[ise][ipx][ix] = 0;
        fitqun1r_.fq1rdir[ise][ipx][ix] = 0;
      }
      fitqun1r_.fq1rdconv[ise][ipx] = 0;
      fitqun1r_.fq1reloss[ise][ipx] = 0;
      fitqun1r_.fq1rpcflg[ise][ipx] = 0;
    }
    fitqun1r_.fqq50[ise] = 0;
    fitqun1r_.fqn50[ise] = 0;
    fitqun1r_.fqitwnd[ise] = 0;
    fitqun1r_.fqipeak[ise] = 0;
    fitqun1r_.fqnhitpmt[ise] = 0;
    fitqun1r_.fqtotq[ise] = 0;
    fitqun1r_.fq0rtotmu[ise] = 0;
    fitqun1r_.fq0rnll[ise] = 0;
  }
  fitqun1r_.fqnse = 0;
  
  for (int itype=0; itype<2; itype++) {
    fitqunpi0_.fqpi0pcflg[itype] = 0;
    fitqunpi0_.fqpi0mom1[itype] = 0;
    fitqunpi0_.fqpi0mom2[itype] = 0;
    fitqunpi0_.fqpi0momtot[itype] = 0;
    fitqunpi0_.fqpi0dconv1[itype] = 0;
    fitqunpi0_.fqpi0dconv2[itype] = 0;
    fitqunpi0_.fqpi0t0[itype] = 0;
    fitqunpi0_.fqpi0totmu[itype] = 0;
    fitqunpi0_.fqpi0nll[itype] = 0;
    fitqunpi0_.fqpi0mass[itype] = 0;
    fitqunpi0_.fqpi0photangle[itype] = 0;
    for (int ix=0; ix<3; ix++) {
      fitqunpi0_.fqpi0pos[itype][ix] = 0;
      fitqunpi0_.fqpi0dir1[itype][ix] = 0;
      fitqunpi0_.fqpi0dir2[itype][ix] = 0;
      fitqunpi0_.fqpi0dirtot[itype][ix] = 0;
    }
  }
  
  for (int itype=0; itype<2; itype++) {
    fitqunpmg_.fqpmgpcflg[itype] = 0;
    fitqunpmg_.fqpmgmom1[itype] = 0;
    fitqunpmg_.fqpmgmom2[itype] = 0;
    fitqunpmg_.fqpmgt01[itype] = 0;
    fitqunpmg_.fqpmgt02[itype] = 0;
    fitqunpmg_.fqpmgtotmu[itype] = 0;
    fitqunpmg_.fqpmgnll[itype] = 0;
    for (int ix=0; ix<3; ix++) {
      fitqunpmg_.fqpmgpos1[itype][ix] = 0;
      fitqunpmg_.fqpmgpos2[itype][ix] = 0;
      fitqunpmg_.fqpmgdir1[itype][ix] = 0;
      fitqunpmg_.fqpmgdir2[itype][ix] = 0;
    }
  }

  for (int idx1=0; idx1<FQMAXNMRFIT; idx1++) {
    fitqunmr_.fqmrifit[idx1] = 0;
    fitqunmr_.fqmrnring[idx1] = 0;
    fitqunmr_.fqmrpcflg[idx1] = 0;
    fitqunmr_.fqmrnll[idx1] = 0;
    fitqunmr_.fqmrtotmu[idx1] = 0;
    for (int idx2=0; idx2<FQMAXNRING; idx2++) {
      
      fitqunmr_.fqmrpid[idx1][idx2] = 0;
      
      fitqunmr_.fqmrmom[idx1][idx2] = 0;
      fitqunmr_.fqmrdconv[idx1][idx2] = 0;
      fitqunmr_.fqmreloss[idx1][idx2] = 0;
      fitqunmr_.fqmrt0[idx1][idx2] = 0;
      
      for (int ix=0; ix<3; ix++) {
        fitqunmr_.fqmrpos[idx1][idx2][ix] = 0;
        fitqunmr_.fqmrdir[idx1][idx2][ix] = 0;
      }
    }
  }
  fitqunmr_.fqnmrfit = 0;
  
  for (int idx1=0; idx1<FQMSMAXFIT; idx1++) {
    fitqunms_.fqmsnseg[idx1] = 0;
    fitqunms_.fqmspid[idx1] = 0;
    fitqunms_.fqmsifit[idx1] = 0;
    fitqunms_.fqmsimer[idx1] = 0;
    fitqunms_.fqmspcflg[idx1] = 0;
    fitqunms_.fqmstotmu[idx1] = 0;
    fitqunms_.fqmsnll[idx1] = 0;
    for (int idx2=0; idx2<FQMSMAXRING; idx2++) {
      fitqunms_.fqmsmom[idx1][idx2] = 0;
      fitqunms_.fqmseloss[idx1][idx2] = 0;
      fitqunms_.fqmst0[idx1][idx2] = 0;
      for (int ix=0; ix<3; ix++) {
        fitqunms_.fqmspos[idx1][idx2][ix] = 0;
        fitqunms_.fqmsdir[idx1][idx2][ix] = 0;
      }
    }
  }
  fitqunms_.fqmsnfit = 0;
  
  for (int i=0; i<FQTESTMAX1R; i++) {
    fitqun1rtest_.fqtest1rstage[i] = 0;
    fitqun1rtest_.fqtest1rse[i] = 0;
    fitqun1rtest_.fqtest1rpid[i] = 0;
    fitqun1rtest_.fqtest1rpcflg[i] = 0;
    fitqun1rtest_.fqtest1rmom[i] = 0;
    fitqun1rtest_.fqtest1rt0[i] = 0;
    fitqun1rtest_.fqtest1rtotmu[i] = 0;
    fitqun1rtest_.fqtest1rnll[i] = 0;
    for (int ix=0; ix<3; ix++) {
      fitqun1rtest_.fqtest1rpos[i][ix] = 0;
      fitqun1rtest_.fqtest1rdir[i][ix] = 0;
    }
    fitqun1rtest_.fqtest1rdconv[i] = 0;
    fitqun1rtest_.fqtest1reloss[i] = 0;
  }
  fitqun1rtest_.fqtestn1r = 0;
  
  for (int i=0; i<FQTESTMAXPI0; i++) {
    fitqunpi0test_.fqtestpi0stage[i] = 0;
    fitqunpi0test_.fqtestpi0pcflg[i] = 0;
    fitqunpi0test_.fqtestpi0mom1[i] = 0;
    fitqunpi0test_.fqtestpi0mom2[i] = 0;
    fitqunpi0test_.fqtestpi0momtot[i] = 0;
    fitqunpi0test_.fqtestpi0dconv1[i] = 0;
    fitqunpi0test_.fqtestpi0dconv2[i] = 0;
    fitqunpi0test_.fqtestpi0t0[i] = 0;
    fitqunpi0test_.fqtestpi0totmu[i] = 0;
    fitqunpi0test_.fqtestpi0nll[i] = 0;
    fitqunpi0test_.fqtestpi0mass[i] = 0;
    fitqunpi0test_.fqtestpi0photangle[i] = 0;
    for (int ix=0; ix<3; ix++) {
      fitqunpi0test_.fqtestpi0pos[i][ix] = 0;
      fitqunpi0test_.fqtestpi0dir1[i][ix] = 0;
      fitqunpi0test_.fqtestpi0dir2[i][ix] = 0;
      fitqunpi0test_.fqtestpi0dirtot[i][ix] = 0;
    }
  }
  fitqunpi0test_.fqtestnpi0 = 0;
  
}