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 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


#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 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 nclsy;
    int planey[NSINGMAX];
    float ys[NSINGMAX][2];
    float signlys[NSINGMAX];
    int clsy[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];
//      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;
            planey[i]  = 0;
            ys[i][0]   = 0;
            ys[i][1]   = 0;
            signlys[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;
            }
        }
    }
    
};
// ==================================================================
/**
 * \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];
    double xm_b[NPLANE],ym_b[NPLANE];
    double resx[NPLANE],resy[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;

/*     cMini2track(){ */
/* //  -------------------------- */
/* //  fitting routine parameters */
/* //  -------------------------- */
/*      zini = 23.5; */
/*      trackmode = 0; */
/*      istepmin = 3; */
/* //  -------------------------- */
/*      pfixed = 0.; */
/*      chi2  = 0; */
/*      nstep = 0; */
/*      for(int it1=0;it1<5;it1++){ */
/*          al[it1] = 0; */
/*          for(int it2=0;it2<5;it2++)cov[it1][it2] = 0; */
/*      }; */
/*      for(int ip=0;ip<6;ip++){ */
/*          xgood[ip]  = 0; */
/*          ygood[ip]  = 0; */
/*          xm[ip]     = 0; */
/*          ym[ip]     = 0; */
/*          xm_a[ip]     = 0; */
/*          ym_a[ip]     = 0; */
/*          xm_b[ip]     = 0; */
/*          ym_b[ip]     = 0; */
/*          zm[ip]     = 0; */
/*          resx[ip]   = 0; */
/*          resy[ip]   = 0; */
/*          xv[ip]     = 0; */
/*          yv[ip]     = 0; */
/*          zv[ip]     = 0; */
/*          axv[ip]    = 0; */
/*          ayv[ip]    = 0; */
/*          dedxtrk_x[ip] = 0; */
/*          dedxtrk_y[ip] = 0; */
/*      }; */
        
/*     } */


};
//
/**
 * \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;
    }
  }
  
};

/**
 * \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 );

};

/**
 * \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];

};
/**
 * \brief Struct to configure data reduction
 */
struct cTrkSW{
    int pfaid;
};

extern "C" {

    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 cMini2track track_; 
    extern struct cTrkETA     pfa_; 
    extern struct cTrkSW      sw_;

    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*);

}

#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.9
14			#define NCLMAX_VIEW 200 //max n.clusters per view
15			#define NCLSTR 20 //max n.strip per cluster
16			#define NCLMAX (NCLMAX_VIEW*12) //max n.clusters total
17			#define NCLBUFF (NCLMAX*NCLSTR)
18
19			#define NDBLT_MAX_NT 1000 // 0.2*ndblt_max
20			#define NTRPT_MAX_NT 10000 // 0.2*ntrpt_max
21			#define NCLOYZ_MAX 200
22			#define NCLOXZ_MAX 200
23
24	mocchiut	1.1	#define NTRKMAX 10
25	pam-fi	1.9	#define NSINGMAX NCLMAX //100
26	mocchiut	1.1
27	pam-fi	1.15	#define NANGMAX 21
28	pam-fi	1.16	#define NETAVALMAX 500
29	pam-fi	1.5
30	pam-fi	1.8
31	pam-fi	1.16	#include <CalibTrk1Event.h>
32			#include <CalibTrk2Event.h>
33			using namespace pamela;
34			//#include <TrkCalib.h>
35	pam-fi	1.2
36	mocchiut	1.1	#include <TString.h>
37			#include <TFile.h>
38	pam-fi	1.2	#include <TTree.h>
39	mocchiut	1.1
40	pam-fi	1.16
41	pam-fi	1.2	#include <fstream>
42			#include <iostream>
43	mocchiut	1.1	/**
44	pam-fi	1.2	* \brief Struct to pass calibration parameters to F77 routines
45	mocchiut	1.1	*/
46	pam-fi	1.2	struct cTrkCalib {
47
48	mocchiut	1.1	float pedestal[128][24][12];
49			float pedestal_t[128][24][12];
50			float sigma[128][24][12];
51			float sigma_t[128][24][12];
52			int bad[128][24][12];
53	pam-fi	1.2
54	pam-fi	1.16	/* void FillACalibFrom(TFile* , Int_t , Int_t ); */
55			/* void FillFCalibFrom(TFile* , Int_t , Int_t ); */
56			/* void FillTCalibFrom(TFile* , Int_t , Int_t ); */
57			/* void FillTCalibFrom(TString); */
58	pam-fi	1.2
59	mocchiut	1.1	};
60			// ==================================================================
61			/**
62			* \brief Struct to pass tracker LEVEL0 data to F77 routines
63			*/
64			struct cTrkLevel0 {
65			int DAQmode[12];
66			int DSPnumber[12];
67			int DATAlength[12];
68			int eventn[12];
69			int nclust[12];
70			int cutc[12];
71			int cutcl[12];
72			int addrcluster[3][12];
73			int signcluster[3][12];
74			int fc[12];
75			int compressiontime[12];
76			int fl5[12];
77			int fl4[12];
78			int fl3[12];
79			int fl2[12];
80			int fl1[12];
81			int fl6[12];
82			int checksum[12];
83			int TOTDATAlength;
84			int datatracker[49152];
85			int pnum[12];
86			int cmdnum[12];
87			int bid[12];
88			int alarm[12];
89			int aswr[12];
90			int good0;
91			int crc[12];
92			};
93			// ==================================================================
94			/**
95			* \brief Struct to pass tracker LEVEL1 data to F77 routines
96			*/
97			struct cTrkLevel1 {
98	pam-fi	1.6	// int good1;
99	pam-fi	1.8	int good[12];
100	mocchiut	1.1	int nclstr1;
101	pam-fi	1.8	int view[NCLMAX];
102			int ladder[NCLMAX];
103			int maxs[NCLMAX];
104			int mult[NCLMAX];
105			float dedx[NCLMAX];
106			int indstart[NCLMAX];
107			int indmax[NCLMAX];
108	mocchiut	1.1	int totCLlength;
109	pam-fi	1.8	float clsignal[NCLBUFF];
110			float clsigma[NCLBUFF];
111			int cladc[NCLBUFF];
112			int clbad[NCLBUFF];
113			float cnev[24][12];
114			int cnnev[24][12];
115			float cnrmsev[24][12];
116			int fshower[12];
117			int whichtrack[NCLMAX];
118	pam-fi	1.16
119	mocchiut	1.1	};
120			// ==================================================================
121			/**
122			* \brief Struct to pass tracker LEVEL2 data to F77 routines
123			*/
124			struct cTrkLevel2 {
125	pam-fi	1.6	// int good2;
126			// int crc[12];
127	pam-fi	1.19	int good[12];
128			int vkflag[24][12];
129	mocchiut	1.1	// ----------------------------
130			int nclsx;
131			int planex[NSINGMAX];
132			float xs[NSINGMAX][2];
133			float signlxs[NSINGMAX];
134	pam-fi	1.8	int clsx[NSINGMAX];
135	mocchiut	1.1	// ----------------------------
136			int nclsy;
137			int planey[NSINGMAX];
138			float ys[NSINGMAX][2];
139			float signlys[NSINGMAX];
140	pam-fi	1.8	int clsy[NSINGMAX];
141	mocchiut	1.1	// ----------------------------
142			int ntrk;
143			int image[NTRKMAX];
144			float xm_nt[NTRKMAX][6];
145			float ym_nt[NTRKMAX][6];
146			float zm_nt[NTRKMAX][6];
147			float resx_nt[NTRKMAX][6];
148			float resy_nt[NTRKMAX][6];
149	pam-fi	1.19	float tailx[NTRKMAX][6];
150			float taily[NTRKMAX][6];
151	mocchiut	1.1	float al_nt[NTRKMAX][5];
152			float coval[NTRKMAX][5][5];
153			float chi2_nt[NTRKMAX];
154	pam-fi	1.8	int nstep_nt[NTRKMAX];
155			int xgood_nt[NTRKMAX][6];
156	mocchiut	1.1	int ygood_nt[NTRKMAX][6];
157	pam-fi	1.19	int ls[NTRKMAX][6];
158			int xbad[NTRKMAX][6];
159			int ybad[NTRKMAX][6];
160	mocchiut	1.1	float xv_nt[NTRKMAX][6];
161			float yv_nt[NTRKMAX][6];
162			float zv_nt[NTRKMAX][6];
163			float axv_nt[NTRKMAX][6];
164			float ayv_nt[NTRKMAX][6];
165			float dedx_x[NTRKMAX][6];
166			float dedx_y[NTRKMAX][6];
167	pam-fi	1.8	int cltrx[NTRKMAX][6];
168			int cltry[NTRKMAX][6];
169	pam-fi	1.5	// float bdl[NTRKMAX];
170	pam-fi	1.3
171	pam-fi	1.19	void InitcTrkLevel2(){
172
173	pam-fi	1.6	// good2 = 0;
174	pam-fi	1.19	for(Int_t i=0; i<12 ; i++){
175			good[i] = 0;
176			for(Int_t ii=0; ii<24 ; ii++)vkflag[ii][i]=0;
177			}
178			// ----------------------------
179			nclsx = 0;
180			nclsy = 0;
181			for(Int_t i=0; i<NSINGMAX ; i++){
182			planex[i] = 0;
183			xs[i][0] = 0;
184			xs[i][1] = 0;
185			signlxs[i] = 0;
186			planey[i] = 0;
187			ys[i][0] = 0;
188			ys[i][1] = 0;
189			signlys[i] = 0;
190			}
191			// ----------------------------
192			ntrk =0 ;
193
194			for(Int_t i=0; i<NTRKMAX ; i++){
195			image[i] = 0;
196			chi2_nt[i] = 0;
197			// bdl[i] = 0;
198
199			for(Int_t ii=0; ii<5 ; ii++){
200			al_nt[i][ii] = 0;
201			for(Int_t iii=0; iii<5 ; iii++)
202			coval[i][ii][iii] = 0;
203			}
204
205			for(Int_t ii=0; ii<6 ; ii++){
206			xm_nt[i][ii] = 0;
207			ym_nt[i][ii] = 0;
208			zm_nt[i][ii] = 0;
209			resx_nt[i][ii] = 0;
210			resy_nt[i][ii] = 0;
211			xgood_nt[i][ii] = 0;
212			ygood_nt[i][ii] = 0;
213			xv_nt[i][ii] = 0;
214			yv_nt[i][ii] = 0;
215			zv_nt[i][ii] = 0;
216			axv_nt[i][ii] = 0;
217			ayv_nt[i][ii] = 0;
218			dedx_x[i][ii] = 0;
219			dedx_y[i][ii] = 0;
220			}
221			}
222	pam-fi	1.3	}
223
224	mocchiut	1.1	};
225			// ==================================================================
226			/**
227			* \brief Struct to pass calibration/parameter file paths to F77 routines
228			*/
229			struct cPath {
230	pam-fi	1.8	char path[256];
231	mocchiut	1.1	int pathlen;
232	pam-fi	1.2	int error;
233	mocchiut	1.1	/**
234			* Fill the struct variables from a TString object
235			*/
236			void FillWith(TString s){
237	pam-fi	1.2	pathlen = s.Length();
238			const char *pc = s.Data();
239			for(Int_t i=0; i<=pathlen; i++) path[i] = *pc++;
240	mocchiut	1.1	};
241	pam-fi	1.11
242	mocchiut	1.1	};
243
244			// ==================================================================
245	pam-fi	1.11	/* /\** */
246			/* * \brief Struct to pass magnetic-field file paths to F77 routines */
247			/* \/ /
248			/* struct cBPath { */
249			/* char b_path[256]; */
250			/* int b_pathlen; */
251			/* int b_error; */
252			/* int b_loaded; */
253			/* /\** */
254			/* * Fill the struct variables from a TString object and set */
255			/* * the load flag to FALSE. */
256			/* \/ /
257			/* void FillWith(TString s){ */
258			/* b_loaded = 0; */
259			/* b_pathlen = s.Length(); */
260			/* const char pc = s.Data(); /
261			/* for(Int_t i=0; i<=b_pathlen; i++) b_path[i] = pc++; /
262			/* }; */
263			/* /\** */
264			/* * Fill the struct variables from a TString object */
265			/* \/ /
266			/* int BIsLoaded(){ return b_loaded; }; */
267	mocchiut	1.1
268	pam-fi	1.11	/* }; */
269	mocchiut	1.1
270	pam-fi	1.2	// ==================================================================
271			/**
272			* \brief Struct to set debug mode in F77 routines
273			*/
274			struct cDbg {
275			int debug;
276			int verbose;
277			int warning;
278			// bool debug;
279			// bool verbose;
280			void SetNone() {debug=0; verbose=0; warning=0;};
281			void SetWarning(){debug=0; verbose=0; warning=1;};
282			void SetVerbose(){debug=0; verbose=1; warning=1;};
283			void SetDebug() {debug=1; verbose=1; warning=1;};
284			};
285
286	pam-fi	1.8
287			// ==================================================================
288			/**
289			* \brief Struct to pass mini2 track parameters to F77 routines
290			*/
291			struct cMini2track {
292			double al[5];
293			double xm[NPLANE],ym[NPLANE],zm[NPLANE];
294			double xm_a[NPLANE],ym_a[NPLANE];
295			double xm_b[NPLANE],ym_b[NPLANE];
296			double resx[NPLANE],resy[NPLANE];
297			double xgood[NPLANE],ygood[NPLANE];
298			double dedxtrk_x[NPLANE];
299			double dedxtrk_y[NPLANE];
300			double zini;
301			double pfixed;
302			double chi2;
303			double xv[NPLANE],yv[NPLANE],zv[NPLANE];
304			double axv[NPLANE],ayv[NPLANE];
305			double cov[5][5];
306	pam-fi	1.18	double fact;
307	pam-fi	1.8	int nstep;
308			int idcand;
309	pam-fi	1.17	int trackmode;
310	pam-fi	1.18	int istepmin;
311	pam-fi	1.20
312			/* cMini2track(){ */
313			/* // -------------------------- */
314			/* // fitting routine parameters */
315			/* // -------------------------- */
316			/* zini = 23.5; */
317			/* trackmode = 0; */
318			/* istepmin = 3; */
319			/* // -------------------------- */
320			/* pfixed = 0.; */
321			/* chi2 = 0; */
322			/* nstep = 0; */
323			/* for(int it1=0;it1<5;it1++){ */
324			/* al[it1] = 0; */
325			/* for(int it2=0;it2<5;it2++)cov[it1][it2] = 0; */
326			/* }; */
327			/* for(int ip=0;ip<6;ip++){ */
328			/* xgood[ip] = 0; */
329			/* ygood[ip] = 0; */
330			/* xm[ip] = 0; */
331			/* ym[ip] = 0; */
332			/* xm_a[ip] = 0; */
333			/* ym_a[ip] = 0; */
334			/* xm_b[ip] = 0; */
335			/* ym_b[ip] = 0; */
336			/* zm[ip] = 0; */
337			/* resx[ip] = 0; */
338			/* resy[ip] = 0; */
339			/* xv[ip] = 0; */
340			/* yv[ip] = 0; */
341			/* zv[ip] = 0; */
342			/* axv[ip] = 0; */
343			/* ayv[ip] = 0; */
344			/* dedxtrk_x[ip] = 0; */
345			/* dedxtrk_y[ip] = 0; */
346			/* }; */
347
348			/* } */
349
350
351	pam-fi	1.8	};
352			//
353	pam-fi	1.9	/**
354			* \Struct for the hough transform variables
355			*/
356
357			struct cTrkHough {
358
359	pam-fi	1.12	int ndblt_nt;
360			float alfayz1_nt[NDBLT_MAX_NT];
361			float alfayz2_nt[NDBLT_MAX_NT];
362			int db_cloud_nt[NDBLT_MAX_NT];
363			int ntrpt_nt;
364			float alfaxz1_nt[NTRPT_MAX_NT];
365			float alfaxz2_nt[NTRPT_MAX_NT];
366			float alfaxz3_nt[NTRPT_MAX_NT];
367			int tr_cloud_nt[NTRPT_MAX_NT];
368			int nclouds_yz_nt;
369			float alfayz1_av_nt[NCLOYZ_MAX];
370			float alfayz2_av_nt[NCLOYZ_MAX];
371			int ptcloud_yz_nt[NCLOYZ_MAX];
372			int nclouds_xz_nt;
373			float alfaxz1_av_nt[NCLOXZ_MAX];
374			float alfaxz2_av_nt[NCLOXZ_MAX];
375			float alfaxz3_av_nt[NCLOXZ_MAX];
376			int ptcloud_xz_nt[NCLOXZ_MAX];
377	pam-fi	1.9	int nclstr;
378			float totaltime;
379			float houghtime;
380			float fittime;
381
382
383			void InitcTrkHough(){
384
385	pam-fi	1.12	ndblt_nt = 0;
386			ntrpt_nt = 0;
387			nclouds_yz_nt = 0;
388			nclouds_xz_nt = 0;
389	pam-fi	1.9	nclstr = 0;
390			totaltime = 0;
391			houghtime = 0;
392			fittime = 0;
393			for(int i=0;i<NDBLT_MAX_NT ;i++){
394	pam-fi	1.12	alfayz1_nt[i] = 0;
395			alfayz2_nt[i] = 0;
396			db_cloud_nt[i] = 0;
397	pam-fi	1.9	}
398			for(int i=0;i<NTRPT_MAX_NT ;i++){
399	pam-fi	1.12	alfaxz1_nt[i] = 0;
400			alfaxz2_nt[i] = 0;
401			alfaxz3_nt[i] = 0;
402			tr_cloud_nt[i] = 0;
403	pam-fi	1.9	}
404			for(int i=0;i<NCLOYZ_MAX ;i++){
405	pam-fi	1.12	alfayz1_av_nt[i] = 0;
406			alfayz2_av_nt[i] = 0;
407			ptcloud_yz_nt[i] = 0;
408	pam-fi	1.9	}
409			for(int i=0;i<NCLOXZ_MAX ;i++){
410	pam-fi	1.12	ptcloud_xz_nt[i] = 0;
411			alfaxz1_av_nt[i] = 0;
412			alfaxz2_av_nt[i] = 0;
413			alfaxz3_av_nt[i] = 0;
414	pam-fi	1.9	}
415			}
416
417			};
418	pam-fi	1.8
419	pam-fi	1.14	/**
420			* \brief Struct to pass VA1-mask to F77 routines
421			*/
422			struct cTrkMask {
423
424			int mask_vk[NVK][NVIEW];
425			int mask_vk_run[NVK][NVIEW];
426
427	pam-fi	1.16	// void Set(TFile* , Int_t , Int_t );
428	pam-fi	1.14
429			};
430	pam-fi	1.8
431	pam-fi	1.15	/**
432			* \brief Struct for pfa parameters
433			*/
434			struct cTrkETA {
435
436			int nangbin;
437			float angL[NANGMAX],angR[NANGMAX];
438			int netaval;
439			float eta2[NANGMAX][NETAVALMAX];
440			float feta2[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
441			float eta3[NANGMAX][NETAVALMAX];
442			float feta3[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
443			float eta4[NANGMAX][NETAVALMAX];
444			float feta4[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
445
446			};
447	pam-fi	1.20	/**
448			* \brief Struct to configure data reduction
449			*/
450			struct cTrkSW{
451			int pfaid;
452			};
453	pam-fi	1.15
454	pam-fi	1.2	extern "C" {
455	pam-fi	1.8
456	pam-fi	1.9	extern struct cTrkCalib pedsigbad_;
457	pam-fi	1.16	extern struct cTrkMask mask_;
458	pam-fi	1.9	extern struct cTrkLevel0 level0event_;
459			extern struct cTrkLevel1 level1event_;
460			extern struct cTrkLevel2 level2event_;
461			extern struct cPath path_;
462			extern struct cDbg dbg_;
463	pam-fi	1.12	extern struct cTrkHough houghevent_;
464	pam-fi	1.8	extern struct cMini2track track_;
465	pam-fi	1.15	extern struct cTrkETA pfa_;
466	pam-fi	1.20	extern struct cTrkSW sw_;
467	pam-fi	1.7
468	pam-fi	1.2	void fillpedsigfromdefault_();
469			int readmipparam_();
470			int readchargeparam_();
471			int readvkmask_();
472			int readalignparam_();
473			int readetaparam_();
474			void reductionflight_(int*);
475	pam-fi	1.20	int analysisflight_();
476	pam-fi	1.16	int readb_();
477	pam-fi	1.20	void gufld_(float, float);
478			void xyzpam_(int,int,int,int,int,float,float,float,float*);
479	pam-fi	1.21	float riseta_(int,float);
480	pam-fi	1.16
481	pam-fi	1.2	}
482
483	mocchiut	1.1	#endif