TrackerLevel2/src/ExtTrack.cpp

/**
 * \file ExtTrack.cpp
 * \author Elena Vannuccini
 */
#include <ExtTrack.h>
#include <iostream>
#include <math.h>
using namespace std;
//......................................
// F77 routines
//......................................
extern "C" {    
//     void dotrack_(int*, double*, double*, double*, double*, int*);
//     void dotrack2_(int*, double*, double*, double*, double*,double*, double*, double*,int*);
//     void dotrack3_(int*, double*, double*, double*, double*,double*, double*, double*,double*,int*);
//     void mini2_(int*,int*,int*); 
//     void guess_();
//     void gufld_(float*, float*);
//     float risxeta2_(float *);
//     float risxeta3_(float *);
//     float risxeta4_(float *);
//     float risyeta2_(float *);
}

//--------------------------------------
//
//
//--------------------------------------
ExtTrack::ExtTrack(){
    chi2  = 0;
    nstep = 0;
    for(int it1=0;it1<5;it1++){
        al[it1] = 0;
        for(int it2=0;it2<5;it2++)coval[it1][it2] = 0;
    };
    nplanes=0;

    xgood   = 0;
    ygood   = 0;
    xm      = 0;
    ym      = 0;
    zm      = 0;
    resx    = 0;
    resy    = 0;
    xv      = 0;
    yv      = 0;
    zv      = 0;
    axv     = 0;
    ayv     = 0;
    dedx_x  = 0;
    dedx_y  = 0;

    TrkParams::SetMiniDefault();

};
//--------------------------------------
//
//
//--------------------------------------
ExtTrack::ExtTrack(Int_t dim){
    chi2  = 0;
    nstep = 0;
    for(int it1=0;it1<5;it1++){
        al[it1] = 0;
        for(int it2=0;it2<5;it2++)coval[it1][it2] = 0;
    };
    
    SetDimension(dim);

    TrkParams::SetMiniDefault();
};
//--------------------------------------
//
//
//--------------------------------------
ExtTrack::ExtTrack(const ExtTrack& t){
    chi2  = t.chi2;
    nstep = t.nstep;
    for(int it1=0;it1<5;it1++){
        al[it1] = t.al[it1];
        for(int it2=0;it2<5;it2++)coval[it1][it2] = t.coval[it1][it2];
    };
    SetDimension(t.nplanes);
    for(int ip=0;ip<nplanes;ip++){
        xgood[ip]  = t.xgood[ip];
        ygood[ip]  = t.ygood[ip];
        xm[ip]     = t.xm[ip];
        ym[ip]     = t.ym[ip];
        zm[ip]     = t.zm[ip];
        resx[ip]   = t.resx[ip];
        resy[ip]   = t.resy[ip];
        xv[ip]     = t.xv[ip];
        yv[ip]     = t.yv[ip];
        zv[ip]     = t.zv[ip];
        axv[ip]    = t.axv[ip];
        ayv[ip]    = t.ayv[ip];
        dedx_x[ip] = t.dedx_x[ip];
        dedx_y[ip] = t.dedx_y[ip];
    };

    TrkParams::SetMiniDefault();
};
//--------------------------------------
//
//
//--------------------------------------
void ExtTrack::SetDimension(Int_t dim){

    nplanes = dim;

    xgood   = new Int_t[nplanes];
    ygood   = new Int_t[nplanes];
    xm      = new Float_t[nplanes];
    ym      = new Float_t[nplanes];
    zm      = new Float_t[nplanes];
    resx    = new Float_t[nplanes];
    resy    = new Float_t[nplanes];
    xv      = new Float_t[nplanes];
    yv      = new Float_t[nplanes];
    zv      = new Float_t[nplanes];
    axv     = new Float_t[nplanes];
    ayv     = new Float_t[nplanes];
    dedx_x  = new Float_t[nplanes];
    dedx_y  = new Float_t[nplanes];
    
    SetMeasure(0,0,0);
    SetResolution(0,0);
    SetGood(0,0);
    FitReset();
}
//--------------------------------------
//
//
//--------------------------------------
void ExtTrack::Clear(Option_t* option){

    FitReset();
    if(nplanes){
        delete []       xgood;
        delete []       ygood;
        delete []       xm;
        delete []       ym;
        delete []       zm;
        delete []       resx;
        delete []       resy;
        delete []       xv;
        delete []       yv;
        delete []       zv;
        delete []       axv;
        delete []       ayv;
        delete []       dedx_x;
        delete []       dedx_y;
    }

    nplanes = 0;
    xgood   = 0;
    ygood   = 0;
    xm      = 0;
    ym      = 0;
    zm      = 0;
    resx    = 0;
    resy    = 0;
    xv      = 0;
    yv      = 0;
    zv      = 0;
    axv     = 0;
    ayv     = 0;
    dedx_x  = 0;
    dedx_y  = 0;

};
//--------------------------------------
//
//
//--------------------------------------
void ExtTrack::Delete(){
    Clear();
};

/**
 * Set the  position measurements 
 */
void ExtTrack::SetMeasure(Double_t *xmeas, Double_t *ymeas, Double_t *zmeas){
    for(Int_t i=0; i<nplanes; i++) xm[i]= (xmeas? *xmeas++ : 0.);
    for(Int_t i=0; i<nplanes; i++) ym[i]= (ymeas? *ymeas++ : 0.);
    for(Int_t i=0; i<nplanes; i++) zm[i]= (ymeas? *zmeas++ : 0.);
}
/**
 * Set the TrkTrack position resolution 
 */
void ExtTrack::SetResolution(Double_t *rx, Double_t *ry){
    for(Int_t i=0; i<nplanes; i++) resx[i]= (rx? *rx++ : 0.);
    for(Int_t i=0; i<nplanes; i++) resy[i]= (ry? *ry++ : 0.);
}
/**
 * Set the TrkTrack good measurement 
 */
void ExtTrack::SetGood(Int_t *xg, Int_t *yg){

    for(Int_t i=0; i<nplanes; i++) xgood[i]=(xg? *xg++ : 0 );
    for(Int_t i=0; i<nplanes; i++) ygood[i]=(yg? *yg++ : 0 );
}

//--------------------------------------
//
//
//--------------------------------------
void ExtTrack::Dump(){
    cout << endl << "========== Track " ;
    cout << endl << "al       : "; for(int i=0; i<5; i++)cout << al[i] << " ";
    cout << endl << "chi^2    : "<< chi2;
    cout << endl << "n.step   : "<< nstep;
    cout << endl << "xgood    : "; for(int i=0; i<nplanes; i++)cout << XGood(i) ;
    cout << endl << "ygood    : "; for(int i=0; i<nplanes; i++)cout << YGood(i) ;
    cout << endl << "xm       : "; for(int i=0; i<nplanes; i++)cout << xm[i] << " ";
    cout << endl << "ym       : "; for(int i=0; i<nplanes; i++)cout << ym[i] << " ";
    cout << endl << "zm       : "; for(int i=0; i<nplanes; i++)cout << zm[i] << " ";
    cout << endl << "xv       : "; for(int i=0; i<nplanes; i++)cout << xv[i] << " ";
    cout << endl << "yv       : "; for(int i=0; i<nplanes; i++)cout << yv[i] << " ";
    cout << endl << "zv       : "; for(int i=0; i<nplanes; i++)cout << zv[i] << " ";
    cout << endl << "resx     : "; for(int i=0; i<nplanes; i++)cout << resx[i] << " ";
    cout << endl << "resy     : "; for(int i=0; i<nplanes; i++)cout << resy[i] << " ";
    cout << endl << "coval    : "; for(int i=0; i<5; i++)cout << coval[0][i]<<" ";
    cout << endl << "           "; for(int i=0; i<5; i++)cout << coval[1][i]<<" ";
    cout << endl << "           "; for(int i=0; i<5; i++)cout << coval[2][i]<<" ";
    cout << endl << "           "; for(int i=0; i<5; i++)cout << coval[3][i]<<" ";
    cout << endl << "           "; for(int i=0; i<5; i++)cout << coval[4][i]<<" ";
    cout << endl << "dedx_x   : "; for(int i=0; i<nplanes; i++)cout << dedx_x[i] << " ";
    cout << endl << "dedx_y   : "; for(int i=0; i<nplanes; i++)cout << dedx_y[i] << " ";

    cout << endl;
}
/**
 * Reset the fit parameters 
 */
void ExtTrack::FitReset(){
    for(int i=0; i<5; i++) al[i]=-9999.;
    chi2=0.;
    nstep=0;
    for(int i=0; i<nplanes; i++) xv[i]=0.;
    for(int i=0; i<nplanes; i++) yv[i]=0.;
    for(int i=0; i<nplanes; i++) zv[i]=0.;
    for(int i=0; i<nplanes; i++) axv[i]=0.;
    for(int i=0; i<nplanes; i++) ayv[i]=0.;
    for(int i=0; i<5; i++) {
        for(int j=0; j<5; j++) coval[i][j]=0.;
    }
}
/**
 * Method to fill minimization-routine common
 * NB (zini=23.5 is implicitelly set)
 */
void ExtTrack::FillMiniStruct(cMini2track& track){

    for(int i=0; i<nplanes; i++){ 

        track.xgood[i]=XGood(i);
        track.ygood[i]=YGood(i);
        
        track.xm[i]=xm[i];
        track.ym[i]=ym[i];
        track.zm[i]=zm[i];
        
        float segment = 100.; //temporaneo
        track.xm_a[i]=xm[i];
        track.xm_b[i]=xm[i];
        track.ym_a[i]=ym[i];
        track.ym_b[i]=ym[i];
        if(       XGood(i) && !YGood(i) ){
            track.ym_a[i] = track.ym_a[i]+segment;
            track.ym_b[i] = track.ym_b[i]-segment;
        }else if( !XGood(i) && YGood(i)){
            track.xm_a[i] = track.xm_a[i]+segment;
            track.xm_b[i] = track.xm_b[i]-segment;
        }

        
        track.resx[i]=resx[i];
        track.resy[i]=resy[i];
    }

    for(int i=0; i<5; i++) track.al[i]=al[i]; 

    track.zini = 55.;

//    track.zini = 23.5; 
// ZINI = 23.5 !!! it should be the same parameter in all codes
    
}
/**
 * Method to set values from  minimization-routine common
 */
void ExtTrack::SetFromMiniStruct(cMini2track *track){

    for(int i=0; i<5; i++) {
        al[i]=track->al[i];
        for(int j=0; j<5; j++) coval[i][j]=track->cov[i][j];
    }
    chi2  = track->chi2;
    nstep = track->nstep;
    for(int i=0; i<6; i++){ 
        xv[i]  = track->xv[i];
        yv[i]  = track->yv[i];
        zv[i]  = track->zv[i];
        xm[i]  = track->xm[i];
        ym[i]  = track->ym[i];
        zm[i]  = track->zm[i];
        axv[i] = track->axv[i];
        ayv[i] = track->ayv[i]; 
    }
    
}

ClassImp(ExtTrack);
1	pam-fi	1.1	/**
2			* \file ExtTrack.cpp
3			* \author Elena Vannuccini
4			*/
5			#include <ExtTrack.h>
6			#include <iostream>
7			#include <math.h>
8			using namespace std;
9			//......................................
10			// F77 routines
11			//......................................
12			extern "C" {
13			// void dotrack_(int, double, double, double, double, int);
14			// void dotrack2_(int, double, double, double, double,double, double, double,int*);
15			// void dotrack3_(int, double, double, double, double,double, double, double,double,int);
16			// void mini2_(int,int,int*);
17			// void guess_();
18			// void gufld_(float, float);
19			// float risxeta2_(float *);
20			// float risxeta3_(float *);
21			// float risxeta4_(float *);
22			// float risyeta2_(float *);
23			}
24
25			//--------------------------------------
26			//
27			//
28			//--------------------------------------
29			ExtTrack::ExtTrack(){
30			chi2 = 0;
31			nstep = 0;
32			for(int it1=0;it1<5;it1++){
33			al[it1] = 0;
34			for(int it2=0;it2<5;it2++)coval[it1][it2] = 0;
35			};
36			nplanes=0;
37
38			xgood = 0;
39			ygood = 0;
40			xm = 0;
41			ym = 0;
42			zm = 0;
43			resx = 0;
44			resy = 0;
45			xv = 0;
46			yv = 0;
47			zv = 0;
48			axv = 0;
49			ayv = 0;
50			dedx_x = 0;
51			dedx_y = 0;
52
53			TrkParams::SetMiniDefault();
54
55			};
56			//--------------------------------------
57			//
58			//
59			//--------------------------------------
60			ExtTrack::ExtTrack(Int_t dim){
61			chi2 = 0;
62			nstep = 0;
63			for(int it1=0;it1<5;it1++){
64			al[it1] = 0;
65			for(int it2=0;it2<5;it2++)coval[it1][it2] = 0;
66			};
67
68			SetDimension(dim);
69
70			TrkParams::SetMiniDefault();
71			};
72			//--------------------------------------
73			//
74			//
75			//--------------------------------------
76			ExtTrack::ExtTrack(const ExtTrack& t){
77			chi2 = t.chi2;
78			nstep = t.nstep;
79			for(int it1=0;it1<5;it1++){
80			al[it1] = t.al[it1];
81			for(int it2=0;it2<5;it2++)coval[it1][it2] = t.coval[it1][it2];
82			};
83			SetDimension(t.nplanes);
84			for(int ip=0;ip<nplanes;ip++){
85			xgood[ip] = t.xgood[ip];
86			ygood[ip] = t.ygood[ip];
87			xm[ip] = t.xm[ip];
88			ym[ip] = t.ym[ip];
89			zm[ip] = t.zm[ip];
90			resx[ip] = t.resx[ip];
91			resy[ip] = t.resy[ip];
92			xv[ip] = t.xv[ip];
93			yv[ip] = t.yv[ip];
94			zv[ip] = t.zv[ip];
95			axv[ip] = t.axv[ip];
96			ayv[ip] = t.ayv[ip];
97			dedx_x[ip] = t.dedx_x[ip];
98			dedx_y[ip] = t.dedx_y[ip];
99			};
100
101			TrkParams::SetMiniDefault();
102			};
103			//--------------------------------------
104			//
105			//
106			//--------------------------------------
107			void ExtTrack::SetDimension(Int_t dim){
108
109			nplanes = dim;
110
111			xgood = new Int_t[nplanes];
112			ygood = new Int_t[nplanes];
113			xm = new Float_t[nplanes];
114			ym = new Float_t[nplanes];
115			zm = new Float_t[nplanes];
116			resx = new Float_t[nplanes];
117			resy = new Float_t[nplanes];
118			xv = new Float_t[nplanes];
119			yv = new Float_t[nplanes];
120			zv = new Float_t[nplanes];
121			axv = new Float_t[nplanes];
122			ayv = new Float_t[nplanes];
123			dedx_x = new Float_t[nplanes];
124			dedx_y = new Float_t[nplanes];
125
126			SetMeasure(0,0,0);
127			SetResolution(0,0);
128			SetGood(0,0);
129			FitReset();
130			}
131			//--------------------------------------
132			//
133			//
134			//--------------------------------------
135			void ExtTrack::Clear(Option_t* option){
136
137			FitReset();
138			if(nplanes){
139			delete [] xgood;
140			delete [] ygood;
141			delete [] xm;
142			delete [] ym;
143			delete [] zm;
144			delete [] resx;
145			delete [] resy;
146			delete [] xv;
147			delete [] yv;
148			delete [] zv;
149			delete [] axv;
150			delete [] ayv;
151			delete [] dedx_x;
152			delete [] dedx_y;
153			}
154
155			nplanes = 0;
156			xgood = 0;
157			ygood = 0;
158			xm = 0;
159			ym = 0;
160			zm = 0;
161			resx = 0;
162			resy = 0;
163			xv = 0;
164			yv = 0;
165			zv = 0;
166			axv = 0;
167			ayv = 0;
168			dedx_x = 0;
169			dedx_y = 0;
170
171			};
172			//--------------------------------------
173			//
174			//
175			//--------------------------------------
176			void ExtTrack::Delete(){
177			Clear();
178			};
179
180			/**
181			* Set the position measurements
182			*/
183			void ExtTrack::SetMeasure(Double_t xmeas, Double_t ymeas, Double_t *zmeas){
184			for(Int_t i=0; i<nplanes; i++) xm[i]= (xmeas? *xmeas++ : 0.);
185			for(Int_t i=0; i<nplanes; i++) ym[i]= (ymeas? *ymeas++ : 0.);
186			for(Int_t i=0; i<nplanes; i++) zm[i]= (ymeas? *zmeas++ : 0.);
187			}
188			/**
189			* Set the TrkTrack position resolution
190			*/
191			void ExtTrack::SetResolution(Double_t rx, Double_t ry){
192			for(Int_t i=0; i<nplanes; i++) resx[i]= (rx? *rx++ : 0.);
193			for(Int_t i=0; i<nplanes; i++) resy[i]= (ry? *ry++ : 0.);
194			}
195			/**
196			* Set the TrkTrack good measurement
197			*/
198			void ExtTrack::SetGood(Int_t xg, Int_t yg){
199
200			for(Int_t i=0; i<nplanes; i++) xgood[i]=(xg? *xg++ : 0 );
201			for(Int_t i=0; i<nplanes; i++) ygood[i]=(yg? *yg++ : 0 );
202			}
203
204			//--------------------------------------
205			//
206			//
207			//--------------------------------------
208			void ExtTrack::Dump(){
209			cout << endl << "========== Track " ;
210			cout << endl << "al : "; for(int i=0; i<5; i++)cout << al[i] << " ";
211			cout << endl << "chi^2 : "<< chi2;
212			cout << endl << "n.step : "<< nstep;
213			cout << endl << "xgood : "; for(int i=0; i<nplanes; i++)cout << XGood(i) ;
214			cout << endl << "ygood : "; for(int i=0; i<nplanes; i++)cout << YGood(i) ;
215			cout << endl << "xm : "; for(int i=0; i<nplanes; i++)cout << xm[i] << " ";
216			cout << endl << "ym : "; for(int i=0; i<nplanes; i++)cout << ym[i] << " ";
217			cout << endl << "zm : "; for(int i=0; i<nplanes; i++)cout << zm[i] << " ";
218			cout << endl << "xv : "; for(int i=0; i<nplanes; i++)cout << xv[i] << " ";
219			cout << endl << "yv : "; for(int i=0; i<nplanes; i++)cout << yv[i] << " ";
220			cout << endl << "zv : "; for(int i=0; i<nplanes; i++)cout << zv[i] << " ";
221			cout << endl << "resx : "; for(int i=0; i<nplanes; i++)cout << resx[i] << " ";
222			cout << endl << "resy : "; for(int i=0; i<nplanes; i++)cout << resy[i] << " ";
223			cout << endl << "coval : "; for(int i=0; i<5; i++)cout << coval[0][i]<<" ";
224			cout << endl << " "; for(int i=0; i<5; i++)cout << coval[1][i]<<" ";
225			cout << endl << " "; for(int i=0; i<5; i++)cout << coval[2][i]<<" ";
226			cout << endl << " "; for(int i=0; i<5; i++)cout << coval[3][i]<<" ";
227			cout << endl << " "; for(int i=0; i<5; i++)cout << coval[4][i]<<" ";
228			cout << endl << "dedx_x : "; for(int i=0; i<nplanes; i++)cout << dedx_x[i] << " ";
229			cout << endl << "dedx_y : "; for(int i=0; i<nplanes; i++)cout << dedx_y[i] << " ";
230
231			cout << endl;
232			}
233			/**
234			* Reset the fit parameters
235			*/
236			void ExtTrack::FitReset(){
237			for(int i=0; i<5; i++) al[i]=-9999.;
238			chi2=0.;
239			nstep=0;
240			for(int i=0; i<nplanes; i++) xv[i]=0.;
241			for(int i=0; i<nplanes; i++) yv[i]=0.;
242			for(int i=0; i<nplanes; i++) zv[i]=0.;
243			for(int i=0; i<nplanes; i++) axv[i]=0.;
244			for(int i=0; i<nplanes; i++) ayv[i]=0.;
245			for(int i=0; i<5; i++) {
246			for(int j=0; j<5; j++) coval[i][j]=0.;
247			}
248			}
249			/**
250			* Method to fill minimization-routine common
251			* NB (zini=23.5 is implicitelly set)
252			*/
253			void ExtTrack::FillMiniStruct(cMini2track& track){
254
255			for(int i=0; i<nplanes; i++){
256
257			track.xgood[i]=XGood(i);
258			track.ygood[i]=YGood(i);
259
260			track.xm[i]=xm[i];
261			track.ym[i]=ym[i];
262			track.zm[i]=zm[i];
263
264			float segment = 100.; //temporaneo
265			track.xm_a[i]=xm[i];
266			track.xm_b[i]=xm[i];
267			track.ym_a[i]=ym[i];
268			track.ym_b[i]=ym[i];
269			if( XGood(i) && !YGood(i) ){
270			track.ym_a[i] = track.ym_a[i]+segment;
271			track.ym_b[i] = track.ym_b[i]-segment;
272			}else if( !XGood(i) && YGood(i)){
273			track.xm_a[i] = track.xm_a[i]+segment;
274			track.xm_b[i] = track.xm_b[i]-segment;
275			}
276
277
278			track.resx[i]=resx[i];
279			track.resy[i]=resy[i];
280			}
281
282			for(int i=0; i<5; i++) track.al[i]=al[i];
283
284			track.zini = 55.;
285
286			// track.zini = 23.5;
287			// ZINI = 23.5 !!! it should be the same parameter in all codes
288
289			}
290			/**
291			* Method to set values from minimization-routine common
292			*/
293			void ExtTrack::SetFromMiniStruct(cMini2track *track){
294
295			for(int i=0; i<5; i++) {
296			al[i]=track->al[i];
297			for(int j=0; j<5; j++) coval[i][j]=track->cov[i][j];
298			}
299			chi2 = track->chi2;
300			nstep = track->nstep;
301			for(int i=0; i<6; i++){
302			xv[i] = track->xv[i];
303			yv[i] = track->yv[i];
304			zv[i] = track->zv[i];
305			xm[i] = track->xm[i];
306			ym[i] = track->ym[i];
307			zm[i] = track->zm[i];
308			axv[i] = track->axv[i];
309			ayv[i] = track->ayv[i];
310			}
311
312			}
313
314			ClassImp(ExtTrack);