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 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 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];
    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++)
//    crc[i] = 0;       
        good[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;
};
//
/**
 * \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];

};

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

    void fillpedsigfromdefault_();
    int readmipparam_(); 
    int readchargeparam_(); 
    int readvkmask_(); 
    int readalignparam_(); 
    int readetaparam_(); 
    void reductionflight_(int*);
    int analysisflight_(int*);
    int  readb_();

}

#endif
1	/**
2	* \file TrkStruct.h
3	* \author Elena Vannuccini
4	* \date
5	*/
6	#ifndef trkstruct_h
7	#define trkstruct_h
8
9	#define NPLANE 6
10	#define NVIEW 12
11	#define NVK 24
12	#define NLADDER 3
13
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	#define NTRKMAX 10
25	#define NSINGMAX NCLMAX //100
26
27	#define NANGMAX 21
28	#define NETAVALMAX 500
29
30
31	#include <CalibTrk1Event.h>
32	#include <CalibTrk2Event.h>
33	using namespace pamela;
34	//#include <TrkCalib.h>
35
36	#include <TString.h>
37	#include <TFile.h>
38	#include <TTree.h>
39
40
41	#include <fstream>
42	#include <iostream>
43	/**
44	* \brief Struct to pass calibration parameters to F77 routines
45	*/
46	struct cTrkCalib {
47
48	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
54	/* 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
59	};
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	// int good1;
99	int good[12];
100	int nclstr1;
101	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	int totCLlength;
109	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
119	};
120	// ==================================================================
121	/**
122	* \brief Struct to pass tracker LEVEL2 data to F77 routines
123	*/
124	struct cTrkLevel2 {
125	// int good2;
126	// int crc[12];
127	int good[12];
128	// ----------------------------
129	int nclsx;
130	int planex[NSINGMAX];
131	float xs[NSINGMAX][2];
132	float signlxs[NSINGMAX];
133	int clsx[NSINGMAX];
134	// ----------------------------
135	int nclsy;
136	int planey[NSINGMAX];
137	float ys[NSINGMAX][2];
138	float signlys[NSINGMAX];
139	int clsy[NSINGMAX];
140	// ----------------------------
141	int ntrk;
142	int image[NTRKMAX];
143	float xm_nt[NTRKMAX][6];
144	float ym_nt[NTRKMAX][6];
145	float zm_nt[NTRKMAX][6];
146	float resx_nt[NTRKMAX][6];
147	float resy_nt[NTRKMAX][6];
148	float al_nt[NTRKMAX][5];
149	float coval[NTRKMAX][5][5];
150	float chi2_nt[NTRKMAX];
151	int nstep_nt[NTRKMAX];
152	int xgood_nt[NTRKMAX][6];
153	int ygood_nt[NTRKMAX][6];
154	float xv_nt[NTRKMAX][6];
155	float yv_nt[NTRKMAX][6];
156	float zv_nt[NTRKMAX][6];
157	float axv_nt[NTRKMAX][6];
158	float ayv_nt[NTRKMAX][6];
159	float dedx_x[NTRKMAX][6];
160	float dedx_y[NTRKMAX][6];
161	int cltrx[NTRKMAX][6];
162	int cltry[NTRKMAX][6];
163	// float bdl[NTRKMAX];
164
165	void InitcTrkLevel2(){
166
167	// good2 = 0;
168	for(Int_t i=0; i<12 ; i++)
169	// crc[i] = 0;
170	good[i] = 0;
171	// ----------------------------
172	nclsx = 0;
173	nclsy = 0;
174	for(Int_t i=0; i<NSINGMAX ; i++){
175	planex[i] = 0;
176	xs[i][0] = 0;
177	xs[i][1] = 0;
178	signlxs[i] = 0;
179	planey[i] = 0;
180	ys[i][0] = 0;
181	ys[i][1] = 0;
182	signlys[i] = 0;
183	}
184	// ----------------------------
185	ntrk =0 ;
186
187	for(Int_t i=0; i<NTRKMAX ; i++){
188	image[i] = 0;
189	chi2_nt[i] = 0;
190	// bdl[i] = 0;
191
192	for(Int_t ii=0; ii<5 ; ii++){
193	al_nt[i][ii] = 0;
194	for(Int_t iii=0; iii<5 ; iii++)
195	coval[i][ii][iii] = 0;
196	}
197
198	for(Int_t ii=0; ii<6 ; ii++){
199	xm_nt[i][ii] = 0;
200	ym_nt[i][ii] = 0;
201	zm_nt[i][ii] = 0;
202	resx_nt[i][ii] = 0;
203	resy_nt[i][ii] = 0;
204	xgood_nt[i][ii] = 0;
205	ygood_nt[i][ii] = 0;
206	xv_nt[i][ii] = 0;
207	yv_nt[i][ii] = 0;
208	zv_nt[i][ii] = 0;
209	axv_nt[i][ii] = 0;
210	ayv_nt[i][ii] = 0;
211	dedx_x[i][ii] = 0;
212	dedx_y[i][ii] = 0;
213	}
214	}
215	}
216
217	};
218	// ==================================================================
219	/**
220	* \brief Struct to pass calibration/parameter file paths to F77 routines
221	*/
222	struct cPath {
223	char path[256];
224	int pathlen;
225	int error;
226	/**
227	* Fill the struct variables from a TString object
228	*/
229	void FillWith(TString s){
230	pathlen = s.Length();
231	const char *pc = s.Data();
232	for(Int_t i=0; i<=pathlen; i++) path[i] = *pc++;
233	};
234
235	};
236
237	// ==================================================================
238	/* /\** */
239	/* * \brief Struct to pass magnetic-field file paths to F77 routines */
240	/* \/ /
241	/* struct cBPath { */
242	/* char b_path[256]; */
243	/* int b_pathlen; */
244	/* int b_error; */
245	/* int b_loaded; */
246	/* /\** */
247	/* * Fill the struct variables from a TString object and set */
248	/* * the load flag to FALSE. */
249	/* \/ /
250	/* void FillWith(TString s){ */
251	/* b_loaded = 0; */
252	/* b_pathlen = s.Length(); */
253	/* const char pc = s.Data(); /
254	/* for(Int_t i=0; i<=b_pathlen; i++) b_path[i] = pc++; /
255	/* }; */
256	/* /\** */
257	/* * Fill the struct variables from a TString object */
258	/* \/ /
259	/* int BIsLoaded(){ return b_loaded; }; */
260
261	/* }; */
262
263	// ==================================================================
264	/**
265	* \brief Struct to set debug mode in F77 routines
266	*/
267	struct cDbg {
268	int debug;
269	int verbose;
270	int warning;
271	// bool debug;
272	// bool verbose;
273	void SetNone() {debug=0; verbose=0; warning=0;};
274	void SetWarning(){debug=0; verbose=0; warning=1;};
275	void SetVerbose(){debug=0; verbose=1; warning=1;};
276	void SetDebug() {debug=1; verbose=1; warning=1;};
277	};
278
279
280	// ==================================================================
281	/**
282	* \brief Struct to pass mini2 track parameters to F77 routines
283	*/
284	struct cMini2track {
285	double al[5];
286	double xm[NPLANE],ym[NPLANE],zm[NPLANE];
287	double xm_a[NPLANE],ym_a[NPLANE];
288	double xm_b[NPLANE],ym_b[NPLANE];
289	double resx[NPLANE],resy[NPLANE];
290	double xgood[NPLANE],ygood[NPLANE];
291	double dedxtrk_x[NPLANE];
292	double dedxtrk_y[NPLANE];
293	double zini;
294	double pfixed;
295	double chi2;
296	double xv[NPLANE],yv[NPLANE],zv[NPLANE];
297	double axv[NPLANE],ayv[NPLANE];
298	double cov[5][5];
299	double fact;
300	int nstep;
301	int idcand;
302	int trackmode;
303	int istepmin;
304	};
305	//
306	/**
307	* \Struct for the hough transform variables
308	*/
309
310	struct cTrkHough {
311
312	int ndblt_nt;
313	float alfayz1_nt[NDBLT_MAX_NT];
314	float alfayz2_nt[NDBLT_MAX_NT];
315	int db_cloud_nt[NDBLT_MAX_NT];
316	int ntrpt_nt;
317	float alfaxz1_nt[NTRPT_MAX_NT];
318	float alfaxz2_nt[NTRPT_MAX_NT];
319	float alfaxz3_nt[NTRPT_MAX_NT];
320	int tr_cloud_nt[NTRPT_MAX_NT];
321	int nclouds_yz_nt;
322	float alfayz1_av_nt[NCLOYZ_MAX];
323	float alfayz2_av_nt[NCLOYZ_MAX];
324	int ptcloud_yz_nt[NCLOYZ_MAX];
325	int nclouds_xz_nt;
326	float alfaxz1_av_nt[NCLOXZ_MAX];
327	float alfaxz2_av_nt[NCLOXZ_MAX];
328	float alfaxz3_av_nt[NCLOXZ_MAX];
329	int ptcloud_xz_nt[NCLOXZ_MAX];
330	int nclstr;
331	float totaltime;
332	float houghtime;
333	float fittime;
334
335
336	void InitcTrkHough(){
337
338	ndblt_nt = 0;
339	ntrpt_nt = 0;
340	nclouds_yz_nt = 0;
341	nclouds_xz_nt = 0;
342	nclstr = 0;
343	totaltime = 0;
344	houghtime = 0;
345	fittime = 0;
346	for(int i=0;i<NDBLT_MAX_NT ;i++){
347	alfayz1_nt[i] = 0;
348	alfayz2_nt[i] = 0;
349	db_cloud_nt[i] = 0;
350	}
351	for(int i=0;i<NTRPT_MAX_NT ;i++){
352	alfaxz1_nt[i] = 0;
353	alfaxz2_nt[i] = 0;
354	alfaxz3_nt[i] = 0;
355	tr_cloud_nt[i] = 0;
356	}
357	for(int i=0;i<NCLOYZ_MAX ;i++){
358	alfayz1_av_nt[i] = 0;
359	alfayz2_av_nt[i] = 0;
360	ptcloud_yz_nt[i] = 0;
361	}
362	for(int i=0;i<NCLOXZ_MAX ;i++){
363	ptcloud_xz_nt[i] = 0;
364	alfaxz1_av_nt[i] = 0;
365	alfaxz2_av_nt[i] = 0;
366	alfaxz3_av_nt[i] = 0;
367	}
368	}
369
370	};
371
372	/**
373	* \brief Struct to pass VA1-mask to F77 routines
374	*/
375	struct cTrkMask {
376
377	int mask_vk[NVK][NVIEW];
378	int mask_vk_run[NVK][NVIEW];
379
380	// void Set(TFile* , Int_t , Int_t );
381
382	};
383
384	/**
385	* \brief Struct for pfa parameters
386	*/
387	struct cTrkETA {
388
389	int nangbin;
390	float angL[NANGMAX],angR[NANGMAX];
391	int netaval;
392	float eta2[NANGMAX][NETAVALMAX];
393	float feta2[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
394	float eta3[NANGMAX][NETAVALMAX];
395	float feta3[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
396	float eta4[NANGMAX][NETAVALMAX];
397	float feta4[NANGMAX][NLADDER][NVIEW][NETAVALMAX];
398
399	};
400
401	extern "C" {
402
403	extern struct cTrkCalib pedsigbad_;
404	extern struct cTrkMask mask_;
405	extern struct cTrkLevel0 level0event_;
406	extern struct cTrkLevel1 level1event_;
407	extern struct cTrkLevel2 level2event_;
408	extern struct cPath path_;
409	extern struct cDbg dbg_;
410	extern struct cTrkHough houghevent_;
411	extern struct cMini2track track_;
412	extern struct cTrkETA pfa_;
413
414	void fillpedsigfromdefault_();
415	int readmipparam_();
416	int readchargeparam_();
417	int readvkmask_();
418	int readalignparam_();
419	int readetaparam_();
420	void reductionflight_(int*);
421	int analysisflight_(int*);
422	int readb_();
423
424	}
425
426	#endif