TrackerLevel2/inc/TrkStruct.h

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

#define NTRKMAX 10
#define NSINGMAX 100

#define NCLMAX_VIEW 50
#define NCLSTR 20
#define NCLMAX (NCLMAX_VIEW*12)
#define NCLBUFF (NCLMAX*NCLSTR)

#define NPLANE 6

#include <CalibTrk1Event.h>
#include <CalibTrk2Event.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];
    int nstep;
    int idcand;
};
//


extern "C" {

    extern struct cTrkCalib  pedsigbad_;
    extern struct cTrkLevel0 level0event_;
    extern struct cTrkLevel1 level1event_;
    extern struct cTrkLevel2 level2event_;
    extern struct cPath      path_;
    extern struct cBPath     bpath_;
    extern struct cDbg       dbg_;
    extern struct cMini2track track_; 

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


#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			#define NTRKMAX 10
10			#define NSINGMAX 100
11
12	pam-fi	1.8	#define NCLMAX_VIEW 50
13			#define NCLSTR 20
14	pam-fi	1.5	#define NCLMAX (NCLMAX_VIEW*12)
15			#define NCLBUFF (NCLMAX*NCLSTR)
16
17	pam-fi	1.8	#define NPLANE 6
18
19	pam-fi	1.2	#include <CalibTrk1Event.h>
20			#include <CalibTrk2Event.h>
21
22	mocchiut	1.1	#include <TString.h>
23			#include <TFile.h>
24	pam-fi	1.2	#include <TTree.h>
25	mocchiut	1.1
26	pam-fi	1.2	#include <fstream>
27			#include <iostream>
28	mocchiut	1.1	/**
29	pam-fi	1.2	* \brief Struct to pass calibration parameters to F77 routines
30	mocchiut	1.1	*/
31	pam-fi	1.2	struct cTrkCalib {
32
33	mocchiut	1.1	float pedestal[128][24][12];
34			float pedestal_t[128][24][12];
35			float sigma[128][24][12];
36			float sigma_t[128][24][12];
37			int bad[128][24][12];
38	pam-fi	1.2
39			void FillACalibFrom(TFile* , Int_t , Int_t );
40			void FillFCalibFrom(TFile* , Int_t , Int_t );
41			void FillTCalibFrom(TFile* , Int_t , Int_t );
42			void FillTCalibFrom(TString);
43
44	mocchiut	1.1	};
45			// ==================================================================
46			/**
47			* \brief Struct to pass tracker LEVEL0 data to F77 routines
48			*/
49			struct cTrkLevel0 {
50			int DAQmode[12];
51			int DSPnumber[12];
52			int DATAlength[12];
53			int eventn[12];
54			int nclust[12];
55			int cutc[12];
56			int cutcl[12];
57			int addrcluster[3][12];
58			int signcluster[3][12];
59			int fc[12];
60			int compressiontime[12];
61			int fl5[12];
62			int fl4[12];
63			int fl3[12];
64			int fl2[12];
65			int fl1[12];
66			int fl6[12];
67			int checksum[12];
68			int TOTDATAlength;
69			int datatracker[49152];
70			int pnum[12];
71			int cmdnum[12];
72			int bid[12];
73			int alarm[12];
74			int aswr[12];
75			int good0;
76			int crc[12];
77			};
78			// ==================================================================
79			/**
80			* \brief Struct to pass tracker LEVEL1 data to F77 routines
81			*/
82			struct cTrkLevel1 {
83	pam-fi	1.6	// int good1;
84	pam-fi	1.8	int good[12];
85	mocchiut	1.1	int nclstr1;
86	pam-fi	1.8	int view[NCLMAX];
87			int ladder[NCLMAX];
88			int maxs[NCLMAX];
89			int mult[NCLMAX];
90			float dedx[NCLMAX];
91			int indstart[NCLMAX];
92			int indmax[NCLMAX];
93	mocchiut	1.1	int totCLlength;
94	pam-fi	1.8	float clsignal[NCLBUFF];
95			float clsigma[NCLBUFF];
96			int cladc[NCLBUFF];
97			int clbad[NCLBUFF];
98			float cnev[24][12];
99			int cnnev[24][12];
100			float cnrmsev[24][12];
101			int fshower[12];
102			int whichtrack[NCLMAX];
103	mocchiut	1.1	};
104			// ==================================================================
105			/**
106			* \brief Struct to pass tracker LEVEL2 data to F77 routines
107			*/
108			struct cTrkLevel2 {
109	pam-fi	1.6	// int good2;
110			// int crc[12];
111			int good[12];
112	mocchiut	1.1	// ----------------------------
113			int nclsx;
114			int planex[NSINGMAX];
115			float xs[NSINGMAX][2];
116			float signlxs[NSINGMAX];
117	pam-fi	1.8	int clsx[NSINGMAX];
118	mocchiut	1.1	// ----------------------------
119			int nclsy;
120			int planey[NSINGMAX];
121			float ys[NSINGMAX][2];
122			float signlys[NSINGMAX];
123	pam-fi	1.8	int clsy[NSINGMAX];
124	mocchiut	1.1	// ----------------------------
125			int ntrk;
126			int image[NTRKMAX];
127			float xm_nt[NTRKMAX][6];
128			float ym_nt[NTRKMAX][6];
129			float zm_nt[NTRKMAX][6];
130			float resx_nt[NTRKMAX][6];
131			float resy_nt[NTRKMAX][6];
132			float al_nt[NTRKMAX][5];
133			float coval[NTRKMAX][5][5];
134			float chi2_nt[NTRKMAX];
135	pam-fi	1.8	int nstep_nt[NTRKMAX];
136			int xgood_nt[NTRKMAX][6];
137	mocchiut	1.1	int ygood_nt[NTRKMAX][6];
138			float xv_nt[NTRKMAX][6];
139			float yv_nt[NTRKMAX][6];
140			float zv_nt[NTRKMAX][6];
141			float axv_nt[NTRKMAX][6];
142			float ayv_nt[NTRKMAX][6];
143			float dedx_x[NTRKMAX][6];
144			float dedx_y[NTRKMAX][6];
145	pam-fi	1.8	int cltrx[NTRKMAX][6];
146			int cltry[NTRKMAX][6];
147	pam-fi	1.5	// float bdl[NTRKMAX];
148	pam-fi	1.3
149			void InitcTrkLevel2(){
150
151	pam-fi	1.6	// good2 = 0;
152	pam-fi	1.3	for(Int_t i=0; i<12 ; i++)
153	pam-fi	1.6	// crc[i] = 0;
154			good[i] = 0;
155	pam-fi	1.3	// ----------------------------
156			nclsx = 0;
157			nclsy = 0;
158			for(Int_t i=0; i<NSINGMAX ; i++){
159			planex[i] = 0;
160			xs[i][0] = 0;
161			xs[i][1] = 0;
162			signlxs[i] = 0;
163			planey[i] = 0;
164			ys[i][0] = 0;
165			ys[i][1] = 0;
166			signlys[i] = 0;
167			}
168			// ----------------------------
169			ntrk =0 ;
170
171			for(Int_t i=0; i<NTRKMAX ; i++){
172			image[i] = 0;
173			chi2_nt[i] = 0;
174	pam-fi	1.5	// bdl[i] = 0;
175	pam-fi	1.3
176			for(Int_t ii=0; ii<5 ; ii++){
177			al_nt[i][ii] = 0;
178			for(Int_t iii=0; iii<5 ; iii++)
179			coval[i][ii][iii] = 0;
180			}
181
182			for(Int_t ii=0; ii<6 ; ii++){
183			xm_nt[i][ii] = 0;
184			ym_nt[i][ii] = 0;
185			zm_nt[i][ii] = 0;
186			resx_nt[i][ii] = 0;
187			resy_nt[i][ii] = 0;
188			xgood_nt[i][ii] = 0;
189			ygood_nt[i][ii] = 0;
190			xv_nt[i][ii] = 0;
191			yv_nt[i][ii] = 0;
192			zv_nt[i][ii] = 0;
193			axv_nt[i][ii] = 0;
194			ayv_nt[i][ii] = 0;
195			dedx_x[i][ii] = 0;
196			dedx_y[i][ii] = 0;
197			}
198			}
199			}
200
201	mocchiut	1.1	};
202			// ==================================================================
203			/**
204			* \brief Struct to pass calibration/parameter file paths to F77 routines
205			*/
206			struct cPath {
207	pam-fi	1.8	char path[256];
208	mocchiut	1.1	int pathlen;
209	pam-fi	1.2	int error;
210	mocchiut	1.1	/**
211			* Fill the struct variables from a TString object
212			*/
213			void FillWith(TString s){
214	pam-fi	1.2	pathlen = s.Length();
215			const char *pc = s.Data();
216			for(Int_t i=0; i<=pathlen; i++) path[i] = *pc++;
217	mocchiut	1.1	};
218	pam-fi	1.2
219
220	mocchiut	1.1	};
221
222			// ==================================================================
223			/**
224			* \brief Struct to pass magnetic-field file paths to F77 routines
225			*/
226			struct cBPath {
227	pam-fi	1.8	char b_path[256];
228	mocchiut	1.1	int b_pathlen;
229			int b_error;
230			int b_loaded;
231			/**
232			* Fill the struct variables from a TString object and set
233			* the load flag to FALSE.
234			*/
235			void FillWith(TString s){
236			b_loaded = 0;
237			b_pathlen = s.Length();
238			const char *pc = s.Data();
239			for(Int_t i=0; i<=b_pathlen; i++) b_path[i] = *pc++;
240			};
241			/**
242			* Fill the struct variables from a TString object
243			*/
244			int BIsLoaded(){ return b_loaded; };
245
246			};
247
248	pam-fi	1.2	// ==================================================================
249			/**
250			* \brief Struct to set debug mode in F77 routines
251			*/
252			struct cDbg {
253			int debug;
254			int verbose;
255			int warning;
256			// bool debug;
257			// bool verbose;
258			void SetNone() {debug=0; verbose=0; warning=0;};
259			void SetWarning(){debug=0; verbose=0; warning=1;};
260			void SetVerbose(){debug=0; verbose=1; warning=1;};
261			void SetDebug() {debug=1; verbose=1; warning=1;};
262			};
263
264	pam-fi	1.8
265			// ==================================================================
266			/**
267			* \brief Struct to pass mini2 track parameters to F77 routines
268			*/
269			struct cMini2track {
270			double al[5];
271			double xm[NPLANE],ym[NPLANE],zm[NPLANE];
272			double xm_a[NPLANE],ym_a[NPLANE];
273			double xm_b[NPLANE],ym_b[NPLANE];
274			double resx[NPLANE],resy[NPLANE];
275			double xgood[NPLANE],ygood[NPLANE];
276			double dedxtrk_x[NPLANE];
277			double dedxtrk_y[NPLANE];
278			double zini;
279			double pfixed;
280			double chi2;
281			double xv[NPLANE],yv[NPLANE],zv[NPLANE];
282			double axv[NPLANE],ayv[NPLANE];
283			double cov[5][5];
284			int nstep;
285			int idcand;
286			};
287			//
288
289
290	pam-fi	1.2	extern "C" {
291	pam-fi	1.8
292	pam-fi	1.2	extern struct cTrkCalib pedsigbad_;
293			extern struct cTrkLevel0 level0event_;
294			extern struct cTrkLevel1 level1event_;
295			extern struct cTrkLevel2 level2event_;
296			extern struct cPath path_;
297			extern struct cBPath bpath_;
298			extern struct cDbg dbg_;
299	pam-fi	1.8	extern struct cMini2track track_;
300	pam-fi	1.7
301	pam-fi	1.2	void fillpedsigfromdefault_();
302			int readmipparam_();
303			int readchargeparam_();
304			int readvkmask_();
305			int readalignparam_();
306			int readetaparam_();
307			void reductionflight_(int*);
308			int analysisflight_();
309			}
310
311	mocchiut	1.1
312			#endif
313