TrackerLevel2/inc/TrkStruct.h

/**
 * \file TrkStruct.h
 * \author Elena Vannuccini
 * \date
 */
#ifndef trkstruct_h
#define trkstruct_h

#define NPLANE 6
#define NVIEW 12
#define NVK 24
#define NLADDER 3
#define NSENSOR 2
#define NSTRIPVK 128
#define NSTRIP NSTRIPVK*NVK

#define NCLMAX_VIEW 200         //max n.clusters per view
#define NCLSTR 20               //max n.strip per cluster
#define NCLMAX (NCLMAX_VIEW*12) //max n.clusters total
#define NCLBUFF (NCLMAX*NCLSTR)

#define NDBLT_MAX_NT 1000   // 0.2*ndblt_max 
#define NTRPT_MAX_NT 10000   // 0.2*ntrpt_max
#define NCLOYZ_MAX 200 
#define NCLOXZ_MAX 200    

#define NTRKMAX 10
#define NSINGMAX NCLMAX //100

#define NANGMAX 21
#define NETAVALMAX 500

#define NTRACKSMAX (NCLOYZ_MAX*NCLOXZ_MAX)

#include <CalibTrk1Event.h>
#include <CalibTrk2Event.h>
using namespace pamela;
//#include <TrkCalib.h>

#include <TString.h>
#include <TFile.h>
#include <TTree.h>


#include <fstream>
#include <iostream>
/**
 * \brief Struct to pass calibration parameters to F77 routines
 */
struct cTrkCalib {

    float pedestal[128][24][12];
    float pedestal_t[128][24][12];
    float sigma[128][24][12];
    float sigma_t[128][24][12];
    int bad[128][24][12];

    void Reset(){
        for(int is=0; is<128; is++){
            for(int ivk=0; ivk<24; ivk++){
                for(int iv=0; iv<12; iv++){
                    pedestal[is][ivk][iv]=0.;
                    pedestal_t[is][ivk][iv]=0.;
                    sigma[is][ivk][iv]=0.;
                    sigma_t[is][ivk][iv]=0.;
                    bad[is][ivk][iv]=1;
                }
            }
        }               
    }

/*     void FillACalibFrom(TFile* , Int_t , Int_t ); */
/*     void FillFCalibFrom(TFile* , Int_t , Int_t ); */
/*     void FillTCalibFrom(TFile* , Int_t , Int_t ); */
/*     void FillTCalibFrom(TString); */

};
// ==================================================================
/**
 * \brief Struct to pass tracker LEVEL0 data to F77 routines
 */
struct cTrkLevel0 {
    int DAQmode[12];
    int DSPnumber[12];
    int DATAlength[12];
    int eventn[12];
    int nclust[12];
    int cutc[12];
    int cutcl[12];
    int addrcluster[3][12];
    int signcluster[3][12];
    int fc[12];
    int compressiontime[12];
    int fl5[12]; 
    int fl4[12];
    int fl3[12];
    int fl2[12];
    int fl1[12];
    int fl6[12];
    int checksum[12];
    int TOTDATAlength;
    int datatracker[49152];
    int pnum[12];    
    int cmdnum[12];
    int bid[12];
    int alarm[12];
    int aswr[12]; 
    int good0;
    int crc[12]; 
};
// ==================================================================
/**
 * \brief Struct to pass tracker LEVEL1 data to F77 routines
 */
struct cTrkLevel1 {
//    int good1;
    int good[12];
    int nclstr1;
    int view[NCLMAX];
    int ladder[NCLMAX];
    int maxs[NCLMAX];
    int mult[NCLMAX];
    float dedx[NCLMAX];
    int indstart[NCLMAX];
    int indmax[NCLMAX];
    int totCLlength;
    float clsignal[NCLBUFF];
    float clsigma[NCLBUFF];
    int cladc[NCLBUFF];
    int clbad[NCLBUFF];
    float cnev[24][12];
    int cnnev[24][12];
    float cnrmsev[24][12];
    int fshower[12];
    int whichtrack[NCLMAX];

};
// ==================================================================
/**
 * \brief Struct to pass tracker LEVEL2 data to F77 routines
 */
struct cTrkLevel2 {
//    int good2;
//    int crc[12];
    int good[12];
    int vkflag[24][12];
//  ----------------------------
    int nclsx;
    int planex[NSINGMAX];
    float xs[NSINGMAX][2];
    float signlxs[NSINGMAX];
    int clsx[NSINGMAX];
    int sxbad[NSINGMAX];
    int multmaxsx[NSINGMAX];
//  ----------------------------
    int nclsy;
    int planey[NSINGMAX];
    float ys[NSINGMAX][2];
    float signlys[NSINGMAX];
    int clsy[NSINGMAX];
    int sybad[NSINGMAX];
    int multmaxsy[NSINGMAX];
//  ----------------------------
    int ntrk;
    int image[NTRKMAX];
    float xm_nt[NTRKMAX][6];
    float ym_nt[NTRKMAX][6];
    float zm_nt[NTRKMAX][6];
    float resx_nt[NTRKMAX][6];
    float resy_nt[NTRKMAX][6];
    float tailx[NTRKMAX][6];
    float taily[NTRKMAX][6];
    float al_nt[NTRKMAX][5];
    float coval[NTRKMAX][5][5];
    float chi2_nt[NTRKMAX];
    int nstep_nt[NTRKMAX];
    int xgood_nt[NTRKMAX][6];
    int ygood_nt[NTRKMAX][6];
    int ls[NTRKMAX][6];
    int xbad[NTRKMAX][6];
    int ybad[NTRKMAX][6];
    float xv_nt[NTRKMAX][6];
    float yv_nt[NTRKMAX][6];
    float zv_nt[NTRKMAX][6];
    float axv_nt[NTRKMAX][6];
    float ayv_nt[NTRKMAX][6];
    float dedx_x[NTRKMAX][6];
    float dedx_y[NTRKMAX][6];
    int cltrx[NTRKMAX][6];
    int cltry[NTRKMAX][6];
    int   multmaxx[NTRKMAX][6];    
    int   multmaxy[NTRKMAX][6];    
    float seedx[NTRKMAX][6];       
    float seedy[NTRKMAX][6];    
    float xpu[NTRKMAX][6];      
    float ypu[NTRKMAX][6];      
//      float bdl[NTRKMAX];

    void InitcTrkLevel2(){
        
//    good2 = 0;
        for(Int_t i=0; i<12 ; i++){
            good[i] = 0;
            for(Int_t ii=0; ii<24 ; ii++)vkflag[ii][i]=0;
        }
        //  ----------------------------
        nclsx = 0;
        nclsy = 0;
        for(Int_t i=0; i<NSINGMAX ; i++){
            planex[i]  = 0;
            xs[i][0]   = 0;
            xs[i][1]   = 0;
            signlxs[i] = 0;
            clsx[i]    = 0;
            sxbad[i]   = 0;
            planey[i]  = 0;
            ys[i][0]   = 0;
            ys[i][1]   = 0;
            signlys[i] = 0;
            clsy[i]    = 0;
            sybad[i]   = 0;
        }
        //  ----------------------------
        ntrk =0 ;
        
        for(Int_t i=0; i<NTRKMAX ; i++){
            image[i]   = 0; 
            chi2_nt[i] = 0;
//      bdl[i]     = 0;
            
            for(Int_t ii=0; ii<5 ; ii++){
                al_nt[i][ii] = 0;
                for(Int_t iii=0; iii<5 ; iii++)
                    coval[i][ii][iii] = 0;
            }
            
            for(Int_t ii=0; ii<6 ; ii++){
                xm_nt[i][ii]    = 0;
                ym_nt[i][ii]    = 0;
                zm_nt[i][ii]    = 0;
                resx_nt[i][ii]  = 0;
                resy_nt[i][ii]  = 0;
                xgood_nt[i][ii] = 0;
                ygood_nt[i][ii] = 0;
                xv_nt[i][ii]    = 0;
                yv_nt[i][ii]    = 0;
                zv_nt[i][ii]    = 0;
                axv_nt[i][ii]   = 0;
                ayv_nt[i][ii]   = 0;
                dedx_x[i][ii]   = 0;
                dedx_y[i][ii]   = 0;
                multmaxx[i][ii] = 0; 
                multmaxy[i][ii] = 0;
                seedx[i][ii]  = 0;  
                seedy[i][ii]  = 0;
                xpu[i][ii]    = 0;   
                ypu[i][ii]    = 0;   
            }
        }
    }
    
};
// ==================================================================
/**
 * \brief Struct to pass calibration/parameter file paths to F77 routines
 */
struct cPath {
    char path[256];
    int  pathlen;
    int  error;
/**
 * Fill the struct variables from a TString object
 */
    void FillWith(TString s){
        pathlen = s.Length();
        const char *pc = s.Data();
        for(Int_t i=0; i<=pathlen; i++) path[i] = *pc++;
    };
        
};

// ==================================================================
/* /\** */
/*  * \brief Struct to pass magnetic-field file paths to F77 routines */
/*  *\/ */
/* struct cBPath { */
/*     char b_path[256]; */
/*     int  b_pathlen; */
/*     int  b_error; */
/*     int  b_loaded; */
/* /\** */
/*  * Fill the struct variables from a TString object and set  */
/*  * the load flag to FALSE. */
/*  *\/ */
/*     void FillWith(TString s){ */
/*      b_loaded  = 0; */
/*      b_pathlen = s.Length(); */
/*      const char *pc = s.Data(); */
/*      for(Int_t i=0; i<=b_pathlen; i++) b_path[i] = *pc++; */
/*     }; */
/* /\** */
/*  * Fill the struct variables from a TString object */
/*  *\/ */
/*     int BIsLoaded(){ return b_loaded; }; */
        
/* }; */

// ==================================================================
/**
 * \brief Struct to set debug mode in F77 routines
 */
struct cDbg {
    int debug;
    int verbose;
    int warning;
//    bool debug;
//    bool verbose;
    void SetNone()   {debug=0; verbose=0; warning=0;};
    void SetWarning(){debug=0; verbose=0; warning=1;};
    void SetVerbose(){debug=0; verbose=1; warning=1;};
    void SetDebug()  {debug=1; verbose=1; warning=1;};
};


// ==================================================================
/**
 * \brief Struct to pass mini2 track parameters to F77 routines
 */
struct cMini2track {
    double al[5];
    double xm[NPLANE],ym[NPLANE],zm[NPLANE];
    double xm_a[NPLANE],ym_a[NPLANE],zm_a[NPLANE];
    double xm_b[NPLANE],ym_b[NPLANE],zm_b[NPLANE];
    double resx[NPLANE],resy[NPLANE];
    double tailx[NPLANE],taily[NPLANE];
    double xgood[NPLANE],ygood[NPLANE];
    double dedxtrk_x[NPLANE];
    double dedxtrk_y[NPLANE];
    double zini;
    double pfixed;
    double chi2;
    double xv[NPLANE],yv[NPLANE],zv[NPLANE];
    double axv[NPLANE],ayv[NPLANE];
    double cov[5][5];
    double fact;
    int nstep;
    int idcand;
    int trackmode;
    int istepmin;

};

// ==================================================================
/**
 * \brief Struct to pass DELTAB parameters to F77 routines
 */
struct cDeltaB {
    double delta0;
    double delta1;
    double dlt;
};

//
/**
 * \Struct for the hough transform variables
 */

struct cTrkHough {

  int ndblt_nt;
  float alfayz1_nt[NDBLT_MAX_NT];
  float alfayz2_nt[NDBLT_MAX_NT];
  int db_cloud_nt[NDBLT_MAX_NT];
  int ntrpt_nt;
  float alfaxz1_nt[NTRPT_MAX_NT];
  float alfaxz2_nt[NTRPT_MAX_NT];
  float alfaxz3_nt[NTRPT_MAX_NT];
  int tr_cloud_nt[NTRPT_MAX_NT];
  int nclouds_yz_nt;
  float alfayz1_av_nt[NCLOYZ_MAX];
  float alfayz2_av_nt[NCLOYZ_MAX];
  int ptcloud_yz_nt[NCLOYZ_MAX];
  int nclouds_xz_nt;
  float alfaxz1_av_nt[NCLOXZ_MAX];
  float alfaxz2_av_nt[NCLOXZ_MAX];
  float alfaxz3_av_nt[NCLOXZ_MAX];
  int ptcloud_xz_nt[NCLOXZ_MAX];
  int nclstr;
  float totaltime;
  float houghtime;
  float fittime;


  void InitcTrkHough(){

    ndblt_nt = 0;
    ntrpt_nt = 0;
    nclouds_yz_nt = 0;
    nclouds_xz_nt = 0;
    nclstr = 0;
    totaltime = 0;
    houghtime = 0;
    fittime = 0;
    for(int i=0;i<NDBLT_MAX_NT ;i++){
      alfayz1_nt[i] = 0;
      alfayz2_nt[i] = 0;
      db_cloud_nt[i] = 0;
    }
    for(int i=0;i<NTRPT_MAX_NT ;i++){
      alfaxz1_nt[i] = 0;
      alfaxz2_nt[i] = 0;
      alfaxz3_nt[i] = 0;
      tr_cloud_nt[i] = 0;
    }
    for(int i=0;i<NCLOYZ_MAX ;i++){
      alfayz1_av_nt[i] = 0;
      alfayz2_av_nt[i] = 0;
      ptcloud_yz_nt[i] = 0;
    }
    for(int i=0;i<NCLOXZ_MAX ;i++){
      ptcloud_xz_nt[i] = 0;
      alfaxz1_av_nt[i] = 0;
      alfaxz2_av_nt[i] = 0;
      alfaxz3_av_nt[i] = 0;
    }
  }
  
};

/**
 * \Struct for track candidates after hough transform
 */
struct cTrkCandidates {
    int ntracks;
    float al[NTRACKSMAX][5];
    float xm[NTRACKSMAX][NPLANE];
    float ym[NTRACKSMAX][NPLANE];
    float zm[NTRACKSMAX][NPLANE];
    float resx[NTRACKSMAX][NPLANE];
    float resy[NTRACKSMAX][NPLANE];
    float xv[NTRACKSMAX][NPLANE];
    float yv[NTRACKSMAX][NPLANE];
    float zv[NTRACKSMAX][NPLANE];
    float axv[NTRACKSMAX][NPLANE];
    float ayv[NTRACKSMAX][NPLANE];
    float xgood[NTRACKSMAX][NPLANE];
    float ygood[NTRACKSMAX][NPLANE];
    int cp[NTRACKSMAX][NPLANE];
    int cls[NTRACKSMAX][NPLANE];
    int sensor[NTRACKSMAX][NPLANE];
    int ladder[NTRACKSMAX][NPLANE];
    float bx[NTRACKSMAX][NPLANE];
    float by[NTRACKSMAX][NPLANE];
    float chi2[NTRACKSMAX];
    
    void Init(){
        ntracks=0;
        for(int i=0; i<NTRACKSMAX; i++){
            for(int ii=0; ii<5; ii++)al[i][ii]=0;
            chi2[i]=0;
            for(int ip=0; ip<NPLANE; ip++){
                xm[i][ip]=0;
                ym[i][ip]=0;
                zm[i][ip]=0;
                resx[i][ip]=0;
                resy[i][ip]=0;
                xv[i][ip]=0;
                yv[i][ip]=0;
                zv[i][ip]=0;
                axv[i][ip]=0;
                ayv[i][ip]=0;
                xgood[i][ip]=0;
                ygood[i][ip]=0;
                cp[i][ip]=0;
                cls[i][ip]=0;
                sensor[i][ip]=0;
                ladder[i][ip]=0;
                bx[i][ip]=0;
                by[i][ip]=0;
            }
        }

    };


};

/**
 * \brief Struct to pass VA1-mask to F77 routines
 */
struct cTrkMask {

    int mask_vk[NVK][NVIEW];
    int mask_vk_run[NVK][NVIEW];

//    void Set(TFile* , Int_t , Int_t );

    void Reset(){
        for(int ivk=0; ivk<NVK; ivk++){
            for(int iv=0; iv<NVIEW; iv++){
//              mask_vk[ivk][iv]=0;
                mask_vk_run[ivk][iv]=0;
            }
        }
    }

    void Dump(){
    }

};

/**
 * \brief Struct for pfa parameters
 */
struct cTrkETA {

    int nangbin;
    float angL[NANGMAX],angR[NANGMAX];
    int netaval;
    float eta2[NANGMAX][NETAVALMAX];
    float feta2[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
    float eta3[NANGMAX][NETAVALMAX];
    float feta3[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
    float eta4[NANGMAX][NETAVALMAX];
    float feta4[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
    float fcorr[NANGMAX][NLADDER][NVIEW];
    float e234ax[6];//F77: e2fax,e2tax,e3fax,e3tax,e4fax,e4tax
    float e234ay[6];//F77: e2fay,e2tay,e3fay,e3tay,e4fay,e4tay

};
/**
 * \brief Struct to configure data reduction
 */
struct cTrkSW{
    int pfaid;
};
/**
 * \brief Struct to read alignment parameters
 */
struct cTrkAlign {
    
    double omega[NSENSOR][NLADDER][NPLANE];
    double beta[NSENSOR][NLADDER][NPLANE];
    double gamma[NSENSOR][NLADDER][NPLANE];
    double dx[NSENSOR][NLADDER][NPLANE];
    double dy[NSENSOR][NLADDER][NPLANE];
    double dz[NSENSOR][NLADDER][NPLANE];

};

/**
 * \brief Struct to check reduction procedure: CN evaluation
 */
struct cTrkCN {

    float cn[NVK][NVIEW];
    float cnrms[NVK][NVIEW];
    int cnn[NVK][NVIEW];
    int clstr[NSTRIPVK][NVK][NVIEW];
    int strange[NSTRIPVK][NVK][NVIEW];
 
};
/**
 * \brief Struct to check reduction procedure: pedestal subtraction
 */
struct cTrkEvent { 
    float value[NSTRIP];    
};

/**
 * \brief Struct to check reduction procedure: pedestal subtraction
 */
struct cTrkMIP {
    float mip[NLADDER][NVIEW];
};

extern "C" {

    extern struct cTrkMIP     adc2mip_;
    extern struct cTrkCalib   pedsigbad_;
    extern struct cTrkMask    mask_; 
    extern struct cTrkLevel0  level0event_;
    extern struct cTrkLevel1  level1event_;
    extern struct cTrkLevel2  level2event_;
    extern struct cPath       path_;
    extern struct cDbg        dbg_;
    extern struct cTrkHough   houghevent_;
    extern struct cTrkCandidates trackcandidates_;
    extern struct cMini2track track_; 
    extern struct cTrkETA     pfa_; 
    extern struct cTrkSW      sw_;
    extern struct cDeltaB     deltab_;
    extern struct cTrkAlign   alignparameters_;
    extern struct cTrkEvent   calibratedsignal_;

    void fillpedsigfromdefault_();
    int readmipparam_(); 
    int readchargeparam_(); 
    int readvkmask_(); 
    int readalignparam_(); 
    int readetaparam_(); 
    void reductionflight_(int*);
    int analysisflight_();
    int  readb_();
    void gufld_(float*, float*);
    void xyzpam_(int*,int*,int*,int*,int*,float*,float*,float*,float*);
    float riseta_(int*,float*);
    void chisq_(int*, int*);

};


#endif
1	mocchiut	1.1	/**
2			* \file TrkStruct.h
3			* \author Elena Vannuccini
4			* \date
5			*/
6			#ifndef trkstruct_h
7			#define trkstruct_h
8
9	pam-fi	1.9	#define NPLANE 6
10	pam-fi	1.14	#define NVIEW 12
11			#define NVK 24
12	pam-fi	1.15	#define NLADDER 3
13	pam-fi	1.32	#define NSENSOR 2
14	pam-fi	1.33	#define NSTRIPVK 128
15			#define NSTRIP NSTRIPVK*NVK
16	pam-fi	1.9
17			#define NCLMAX_VIEW 200 //max n.clusters per view
18			#define NCLSTR 20 //max n.strip per cluster
19			#define NCLMAX (NCLMAX_VIEW*12) //max n.clusters total
20			#define NCLBUFF (NCLMAX*NCLSTR)
21
22			#define NDBLT_MAX_NT 1000 // 0.2*ndblt_max
23			#define NTRPT_MAX_NT 10000 // 0.2*ntrpt_max
24			#define NCLOYZ_MAX 200
25			#define NCLOXZ_MAX 200
26
27	mocchiut	1.1	#define NTRKMAX 10
28	pam-fi	1.9	#define NSINGMAX NCLMAX //100
29	mocchiut	1.1
30	pam-fi	1.15	#define NANGMAX 21
31	pam-fi	1.16	#define NETAVALMAX 500
32	pam-fi	1.5
33	pam-fi	1.28	#define NTRACKSMAX (NCLOYZ_MAX*NCLOXZ_MAX)
34	pam-fi	1.8
35	pam-fi	1.16	#include <CalibTrk1Event.h>
36			#include <CalibTrk2Event.h>
37			using namespace pamela;
38			//#include <TrkCalib.h>
39	pam-fi	1.2
40	mocchiut	1.1	#include <TString.h>
41			#include <TFile.h>
42	pam-fi	1.2	#include <TTree.h>
43	mocchiut	1.1
44	pam-fi	1.16
45	pam-fi	1.28
46	pam-fi	1.2	#include <fstream>
47			#include <iostream>
48	mocchiut	1.1	/**
49	pam-fi	1.2	* \brief Struct to pass calibration parameters to F77 routines
50	mocchiut	1.1	*/
51	pam-fi	1.2	struct cTrkCalib {
52
53	mocchiut	1.1	float pedestal[128][24][12];
54			float pedestal_t[128][24][12];
55			float sigma[128][24][12];
56			float sigma_t[128][24][12];
57			int bad[128][24][12];
58	pam-fi	1.2
59	pam-fi	1.23	void Reset(){
60			for(int is=0; is<128; is++){
61			for(int ivk=0; ivk<24; ivk++){
62			for(int iv=0; iv<12; iv++){
63			pedestal[is][ivk][iv]=0.;
64			pedestal_t[is][ivk][iv]=0.;
65			sigma[is][ivk][iv]=0.;
66			sigma_t[is][ivk][iv]=0.;
67			bad[is][ivk][iv]=1;
68			}
69			}
70			}
71			}
72
73	pam-fi	1.16	/* void FillACalibFrom(TFile* , Int_t , Int_t ); */
74			/* void FillFCalibFrom(TFile* , Int_t , Int_t ); */
75			/* void FillTCalibFrom(TFile* , Int_t , Int_t ); */
76			/* void FillTCalibFrom(TString); */
77	pam-fi	1.2
78	mocchiut	1.1	};
79			// ==================================================================
80			/**
81			* \brief Struct to pass tracker LEVEL0 data to F77 routines
82			*/
83			struct cTrkLevel0 {
84			int DAQmode[12];
85			int DSPnumber[12];
86			int DATAlength[12];
87			int eventn[12];
88			int nclust[12];
89			int cutc[12];
90			int cutcl[12];
91			int addrcluster[3][12];
92			int signcluster[3][12];
93			int fc[12];
94			int compressiontime[12];
95			int fl5[12];
96			int fl4[12];
97			int fl3[12];
98			int fl2[12];
99			int fl1[12];
100			int fl6[12];
101			int checksum[12];
102			int TOTDATAlength;
103			int datatracker[49152];
104			int pnum[12];
105			int cmdnum[12];
106			int bid[12];
107			int alarm[12];
108			int aswr[12];
109			int good0;
110			int crc[12];
111			};
112			// ==================================================================
113			/**
114			* \brief Struct to pass tracker LEVEL1 data to F77 routines
115			*/
116			struct cTrkLevel1 {
117	pam-fi	1.6	// int good1;
118	pam-fi	1.8	int good[12];
119	mocchiut	1.1	int nclstr1;
120	pam-fi	1.8	int view[NCLMAX];
121			int ladder[NCLMAX];
122			int maxs[NCLMAX];
123			int mult[NCLMAX];
124			float dedx[NCLMAX];
125			int indstart[NCLMAX];
126			int indmax[NCLMAX];
127	mocchiut	1.1	int totCLlength;
128	pam-fi	1.8	float clsignal[NCLBUFF];
129			float clsigma[NCLBUFF];
130			int cladc[NCLBUFF];
131			int clbad[NCLBUFF];
132			float cnev[24][12];
133			int cnnev[24][12];
134			float cnrmsev[24][12];
135			int fshower[12];
136			int whichtrack[NCLMAX];
137	pam-fi	1.16
138	mocchiut	1.1	};
139			// ==================================================================
140			/**
141			* \brief Struct to pass tracker LEVEL2 data to F77 routines
142			*/
143			struct cTrkLevel2 {
144	pam-fi	1.6	// int good2;
145			// int crc[12];
146	pam-fi	1.19	int good[12];
147			int vkflag[24][12];
148	mocchiut	1.1	// ----------------------------
149			int nclsx;
150			int planex[NSINGMAX];
151			float xs[NSINGMAX][2];
152			float signlxs[NSINGMAX];
153	pam-fi	1.8	int clsx[NSINGMAX];
154	pam-fi	1.31	int sxbad[NSINGMAX];
155			int multmaxsx[NSINGMAX];
156	mocchiut	1.1	// ----------------------------
157			int nclsy;
158			int planey[NSINGMAX];
159			float ys[NSINGMAX][2];
160			float signlys[NSINGMAX];
161	pam-fi	1.8	int clsy[NSINGMAX];
162	pam-fi	1.31	int sybad[NSINGMAX];
163			int multmaxsy[NSINGMAX];
164	mocchiut	1.1	// ----------------------------
165			int ntrk;
166			int image[NTRKMAX];
167			float xm_nt[NTRKMAX][6];
168			float ym_nt[NTRKMAX][6];
169			float zm_nt[NTRKMAX][6];
170			float resx_nt[NTRKMAX][6];
171			float resy_nt[NTRKMAX][6];
172	pam-fi	1.19	float tailx[NTRKMAX][6];
173			float taily[NTRKMAX][6];
174	mocchiut	1.1	float al_nt[NTRKMAX][5];
175			float coval[NTRKMAX][5][5];
176			float chi2_nt[NTRKMAX];
177	pam-fi	1.8	int nstep_nt[NTRKMAX];
178			int xgood_nt[NTRKMAX][6];
179	mocchiut	1.1	int ygood_nt[NTRKMAX][6];
180	pam-fi	1.19	int ls[NTRKMAX][6];
181			int xbad[NTRKMAX][6];
182			int ybad[NTRKMAX][6];
183	mocchiut	1.1	float xv_nt[NTRKMAX][6];
184			float yv_nt[NTRKMAX][6];
185			float zv_nt[NTRKMAX][6];
186			float axv_nt[NTRKMAX][6];
187			float ayv_nt[NTRKMAX][6];
188			float dedx_x[NTRKMAX][6];
189			float dedx_y[NTRKMAX][6];
190	pam-fi	1.8	int cltrx[NTRKMAX][6];
191			int cltry[NTRKMAX][6];
192	pam-fi	1.27	int multmaxx[NTRKMAX][6];
193			int multmaxy[NTRKMAX][6];
194			float seedx[NTRKMAX][6];
195			float seedy[NTRKMAX][6];
196			float xpu[NTRKMAX][6];
197			float ypu[NTRKMAX][6];
198	pam-fi	1.5	// float bdl[NTRKMAX];
199	pam-fi	1.3
200	pam-fi	1.19	void InitcTrkLevel2(){
201
202	pam-fi	1.6	// good2 = 0;
203	pam-fi	1.19	for(Int_t i=0; i<12 ; i++){
204			good[i] = 0;
205			for(Int_t ii=0; ii<24 ; ii++)vkflag[ii][i]=0;
206			}
207			// ----------------------------
208			nclsx = 0;
209			nclsy = 0;
210			for(Int_t i=0; i<NSINGMAX ; i++){
211			planex[i] = 0;
212			xs[i][0] = 0;
213			xs[i][1] = 0;
214			signlxs[i] = 0;
215	pam-fi	1.31	clsx[i] = 0;
216			sxbad[i] = 0;
217	pam-fi	1.19	planey[i] = 0;
218			ys[i][0] = 0;
219			ys[i][1] = 0;
220			signlys[i] = 0;
221	pam-fi	1.31	clsy[i] = 0;
222			sybad[i] = 0;
223	pam-fi	1.19	}
224			// ----------------------------
225			ntrk =0 ;
226
227			for(Int_t i=0; i<NTRKMAX ; i++){
228			image[i] = 0;
229			chi2_nt[i] = 0;
230			// bdl[i] = 0;
231
232			for(Int_t ii=0; ii<5 ; ii++){
233			al_nt[i][ii] = 0;
234			for(Int_t iii=0; iii<5 ; iii++)
235			coval[i][ii][iii] = 0;
236			}
237
238			for(Int_t ii=0; ii<6 ; ii++){
239			xm_nt[i][ii] = 0;
240			ym_nt[i][ii] = 0;
241			zm_nt[i][ii] = 0;
242			resx_nt[i][ii] = 0;
243			resy_nt[i][ii] = 0;
244			xgood_nt[i][ii] = 0;
245			ygood_nt[i][ii] = 0;
246			xv_nt[i][ii] = 0;
247			yv_nt[i][ii] = 0;
248			zv_nt[i][ii] = 0;
249			axv_nt[i][ii] = 0;
250			ayv_nt[i][ii] = 0;
251			dedx_x[i][ii] = 0;
252			dedx_y[i][ii] = 0;
253	pam-fi	1.27	multmaxx[i][ii] = 0;
254			multmaxy[i][ii] = 0;
255			seedx[i][ii] = 0;
256			seedy[i][ii] = 0;
257			xpu[i][ii] = 0;
258			ypu[i][ii] = 0;
259	pam-fi	1.19	}
260			}
261	pam-fi	1.3	}
262
263	mocchiut	1.1	};
264			// ==================================================================
265			/**
266			* \brief Struct to pass calibration/parameter file paths to F77 routines
267			*/
268			struct cPath {
269	pam-fi	1.8	char path[256];
270	mocchiut	1.1	int pathlen;
271	pam-fi	1.2	int error;
272	mocchiut	1.1	/**
273			* Fill the struct variables from a TString object
274			*/
275			void FillWith(TString s){
276	pam-fi	1.2	pathlen = s.Length();
277			const char *pc = s.Data();
278			for(Int_t i=0; i<=pathlen; i++) path[i] = *pc++;
279	mocchiut	1.1	};
280	pam-fi	1.11
281	mocchiut	1.1	};
282
283			// ==================================================================
284	pam-fi	1.11	/* /\** */
285			/* * \brief Struct to pass magnetic-field file paths to F77 routines */
286			/* \/ /
287			/* struct cBPath { */
288			/* char b_path[256]; */
289			/* int b_pathlen; */
290			/* int b_error; */
291			/* int b_loaded; */
292			/* /\** */
293			/* * Fill the struct variables from a TString object and set */
294			/* * the load flag to FALSE. */
295			/* \/ /
296			/* void FillWith(TString s){ */
297			/* b_loaded = 0; */
298			/* b_pathlen = s.Length(); */
299			/* const char pc = s.Data(); /
300			/* for(Int_t i=0; i<=b_pathlen; i++) b_path[i] = pc++; /
301			/* }; */
302			/* /\** */
303			/* * Fill the struct variables from a TString object */
304			/* \/ /
305			/* int BIsLoaded(){ return b_loaded; }; */
306	mocchiut	1.1
307	pam-fi	1.11	/* }; */
308	mocchiut	1.1
309	pam-fi	1.2	// ==================================================================
310			/**
311			* \brief Struct to set debug mode in F77 routines
312			*/
313			struct cDbg {
314			int debug;
315			int verbose;
316			int warning;
317			// bool debug;
318			// bool verbose;
319			void SetNone() {debug=0; verbose=0; warning=0;};
320			void SetWarning(){debug=0; verbose=0; warning=1;};
321			void SetVerbose(){debug=0; verbose=1; warning=1;};
322			void SetDebug() {debug=1; verbose=1; warning=1;};
323			};
324
325	pam-fi	1.8
326			// ==================================================================
327			/**
328			* \brief Struct to pass mini2 track parameters to F77 routines
329			*/
330			struct cMini2track {
331			double al[5];
332			double xm[NPLANE],ym[NPLANE],zm[NPLANE];
333	pam-fi	1.32	double xm_a[NPLANE],ym_a[NPLANE],zm_a[NPLANE];
334			double xm_b[NPLANE],ym_b[NPLANE],zm_b[NPLANE];
335	pam-fi	1.8	double resx[NPLANE],resy[NPLANE];
336	pam-fi	1.22	double tailx[NPLANE],taily[NPLANE];
337	pam-fi	1.8	double xgood[NPLANE],ygood[NPLANE];
338			double dedxtrk_x[NPLANE];
339			double dedxtrk_y[NPLANE];
340			double zini;
341			double pfixed;
342			double chi2;
343			double xv[NPLANE],yv[NPLANE],zv[NPLANE];
344			double axv[NPLANE],ayv[NPLANE];
345			double cov[5][5];
346	pam-fi	1.18	double fact;
347	pam-fi	1.8	int nstep;
348			int idcand;
349	pam-fi	1.17	int trackmode;
350	pam-fi	1.18	int istepmin;
351	pam-fi	1.20
352	pam-fi	1.8	};
353	pam-fi	1.29
354			// ==================================================================
355			/**
356			* \brief Struct to pass DELTAB parameters to F77 routines
357			*/
358			struct cDeltaB {
359			double delta0;
360			double delta1;
361	pam-fi	1.30	double dlt;
362	pam-fi	1.29	};
363
364	pam-fi	1.8	//
365	pam-fi	1.9	/**
366			* \Struct for the hough transform variables
367			*/
368
369			struct cTrkHough {
370
371	pam-fi	1.12	int ndblt_nt;
372			float alfayz1_nt[NDBLT_MAX_NT];
373			float alfayz2_nt[NDBLT_MAX_NT];
374			int db_cloud_nt[NDBLT_MAX_NT];
375			int ntrpt_nt;
376			float alfaxz1_nt[NTRPT_MAX_NT];
377			float alfaxz2_nt[NTRPT_MAX_NT];
378			float alfaxz3_nt[NTRPT_MAX_NT];
379			int tr_cloud_nt[NTRPT_MAX_NT];
380			int nclouds_yz_nt;
381			float alfayz1_av_nt[NCLOYZ_MAX];
382			float alfayz2_av_nt[NCLOYZ_MAX];
383			int ptcloud_yz_nt[NCLOYZ_MAX];
384			int nclouds_xz_nt;
385			float alfaxz1_av_nt[NCLOXZ_MAX];
386			float alfaxz2_av_nt[NCLOXZ_MAX];
387			float alfaxz3_av_nt[NCLOXZ_MAX];
388			int ptcloud_xz_nt[NCLOXZ_MAX];
389	pam-fi	1.9	int nclstr;
390			float totaltime;
391			float houghtime;
392			float fittime;
393
394
395			void InitcTrkHough(){
396
397	pam-fi	1.12	ndblt_nt = 0;
398			ntrpt_nt = 0;
399			nclouds_yz_nt = 0;
400			nclouds_xz_nt = 0;
401	pam-fi	1.9	nclstr = 0;
402			totaltime = 0;
403			houghtime = 0;
404			fittime = 0;
405			for(int i=0;i<NDBLT_MAX_NT ;i++){
406	pam-fi	1.12	alfayz1_nt[i] = 0;
407			alfayz2_nt[i] = 0;
408			db_cloud_nt[i] = 0;
409	pam-fi	1.9	}
410			for(int i=0;i<NTRPT_MAX_NT ;i++){
411	pam-fi	1.12	alfaxz1_nt[i] = 0;
412			alfaxz2_nt[i] = 0;
413			alfaxz3_nt[i] = 0;
414			tr_cloud_nt[i] = 0;
415	pam-fi	1.9	}
416			for(int i=0;i<NCLOYZ_MAX ;i++){
417	pam-fi	1.12	alfayz1_av_nt[i] = 0;
418			alfayz2_av_nt[i] = 0;
419			ptcloud_yz_nt[i] = 0;
420	pam-fi	1.9	}
421			for(int i=0;i<NCLOXZ_MAX ;i++){
422	pam-fi	1.12	ptcloud_xz_nt[i] = 0;
423			alfaxz1_av_nt[i] = 0;
424			alfaxz2_av_nt[i] = 0;
425			alfaxz3_av_nt[i] = 0;
426	pam-fi	1.9	}
427			}
428
429			};
430	pam-fi	1.8
431	pam-fi	1.14	/**
432	pam-fi	1.28	* \Struct for track candidates after hough transform
433			*/
434			struct cTrkCandidates {
435			int ntracks;
436			float al[NTRACKSMAX][5];
437			float xm[NTRACKSMAX][NPLANE];
438			float ym[NTRACKSMAX][NPLANE];
439			float zm[NTRACKSMAX][NPLANE];
440			float resx[NTRACKSMAX][NPLANE];
441			float resy[NTRACKSMAX][NPLANE];
442			float xv[NTRACKSMAX][NPLANE];
443			float yv[NTRACKSMAX][NPLANE];
444			float zv[NTRACKSMAX][NPLANE];
445			float axv[NTRACKSMAX][NPLANE];
446			float ayv[NTRACKSMAX][NPLANE];
447			float xgood[NTRACKSMAX][NPLANE];
448			float ygood[NTRACKSMAX][NPLANE];
449			int cp[NTRACKSMAX][NPLANE];
450			int cls[NTRACKSMAX][NPLANE];
451			int sensor[NTRACKSMAX][NPLANE];
452			int ladder[NTRACKSMAX][NPLANE];
453			float bx[NTRACKSMAX][NPLANE];
454			float by[NTRACKSMAX][NPLANE];
455			float chi2[NTRACKSMAX];
456
457			void Init(){
458			ntracks=0;
459			for(int i=0; i<NTRACKSMAX; i++){
460			for(int ii=0; ii<5; ii++)al[i][ii]=0;
461			chi2[i]=0;
462			for(int ip=0; ip<NPLANE; ip++){
463			xm[i][ip]=0;
464			ym[i][ip]=0;
465			zm[i][ip]=0;
466			resx[i][ip]=0;
467			resy[i][ip]=0;
468			xv[i][ip]=0;
469			yv[i][ip]=0;
470			zv[i][ip]=0;
471			axv[i][ip]=0;
472			ayv[i][ip]=0;
473			xgood[i][ip]=0;
474			ygood[i][ip]=0;
475			cp[i][ip]=0;
476			cls[i][ip]=0;
477			sensor[i][ip]=0;
478			ladder[i][ip]=0;
479			bx[i][ip]=0;
480			by[i][ip]=0;
481			}
482			}
483
484			};
485
486
487			};
488
489			/**
490	pam-fi	1.14	* \brief Struct to pass VA1-mask to F77 routines
491			*/
492			struct cTrkMask {
493
494			int mask_vk[NVK][NVIEW];
495			int mask_vk_run[NVK][NVIEW];
496
497	pam-fi	1.16	// void Set(TFile* , Int_t , Int_t );
498	pam-fi	1.14
499	pam-fi	1.23	void Reset(){
500			for(int ivk=0; ivk<NVK; ivk++){
501			for(int iv=0; iv<NVIEW; iv++){
502			// mask_vk[ivk][iv]=0;
503			mask_vk_run[ivk][iv]=0;
504			}
505			}
506			}
507
508	pam-fi	1.24	void Dump(){
509			}
510
511	pam-fi	1.14	};
512	pam-fi	1.8
513	pam-fi	1.15	/**
514			* \brief Struct for pfa parameters
515			*/
516			struct cTrkETA {
517
518			int nangbin;
519			float angL[NANGMAX],angR[NANGMAX];
520			int netaval;
521			float eta2[NANGMAX][NETAVALMAX];
522			float feta2[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
523			float eta3[NANGMAX][NETAVALMAX];
524			float feta3[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
525			float eta4[NANGMAX][NETAVALMAX];
526			float feta4[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
527	pam-fi	1.25	float fcorr[NANGMAX][NLADDER][NVIEW];
528	pam-fi	1.26	float e234ax[6];//F77: e2fax,e2tax,e3fax,e3tax,e4fax,e4tax
529			float e234ay[6];//F77: e2fay,e2tay,e3fay,e3tay,e4fay,e4tay
530	pam-fi	1.15
531			};
532	pam-fi	1.20	/**
533			* \brief Struct to configure data reduction
534			*/
535			struct cTrkSW{
536			int pfaid;
537			};
538	pam-fi	1.32	/**
539			* \brief Struct to read alignment parameters
540			*/
541			struct cTrkAlign {
542
543			double omega[NSENSOR][NLADDER][NPLANE];
544			double beta[NSENSOR][NLADDER][NPLANE];
545			double gamma[NSENSOR][NLADDER][NPLANE];
546			double dx[NSENSOR][NLADDER][NPLANE];
547			double dy[NSENSOR][NLADDER][NPLANE];
548			double dz[NSENSOR][NLADDER][NPLANE];
549
550			};
551	pam-fi	1.15
552	pam-fi	1.33	/**
553			* \brief Struct to check reduction procedure: CN evaluation
554			*/
555			struct cTrkCN {
556
557			float cn[NVK][NVIEW];
558			float cnrms[NVK][NVIEW];
559			int cnn[NVK][NVIEW];
560			int clstr[NSTRIPVK][NVK][NVIEW];
561			int strange[NSTRIPVK][NVK][NVIEW];
562
563			};
564			/**
565			* \brief Struct to check reduction procedure: pedestal subtraction
566			*/
567			struct cTrkEvent {
568			float value[NSTRIP];
569			};
570
571	pam-fi	1.36	/**
572			* \brief Struct to check reduction procedure: pedestal subtraction
573			*/
574			struct cTrkMIP {
575			float mip[NLADDER][NVIEW];
576			};
577
578	pam-fi	1.2	extern "C" {
579	pam-fi	1.8
580	pam-fi	1.36	extern struct cTrkMIP adc2mip_;
581	pam-fi	1.9	extern struct cTrkCalib pedsigbad_;
582	pam-fi	1.16	extern struct cTrkMask mask_;
583	pam-fi	1.9	extern struct cTrkLevel0 level0event_;
584			extern struct cTrkLevel1 level1event_;
585			extern struct cTrkLevel2 level2event_;
586			extern struct cPath path_;
587			extern struct cDbg dbg_;
588	pam-fi	1.12	extern struct cTrkHough houghevent_;
589	pam-fi	1.28	extern struct cTrkCandidates trackcandidates_;
590	pam-fi	1.8	extern struct cMini2track track_;
591	pam-fi	1.15	extern struct cTrkETA pfa_;
592	pam-fi	1.20	extern struct cTrkSW sw_;
593	pam-fi	1.29	extern struct cDeltaB deltab_;
594	pam-fi	1.32	extern struct cTrkAlign alignparameters_;
595	pam-fi	1.33	extern struct cTrkEvent calibratedsignal_;
596	pam-fi	1.7
597	pam-fi	1.2	void fillpedsigfromdefault_();
598			int readmipparam_();
599			int readchargeparam_();
600			int readvkmask_();
601			int readalignparam_();
602			int readetaparam_();
603			void reductionflight_(int*);
604	pam-fi	1.20	int analysisflight_();
605	pam-fi	1.16	int readb_();
606	pam-fi	1.20	void gufld_(float, float);
607			void xyzpam_(int,int,int,int,int,float,float,float,float*);
608	pam-fi	1.21	float riseta_(int,float);
609	pam-fi	1.29	void chisq_(int, int);
610	pam-fi	1.16
611	pam-fi	1.32	};
612	pam-fi	1.2
613	pam-fi	1.33
614	mocchiut	1.1	#endif