KalFitGemHit.cxx

Go to the documentation of this file.

#include "KalFitAlg/KalFitGemHit.h"
#include <cmath>
 
const double a_stero[3] = {(45.94*3.1415926/180),(31.10*3.1415926/180),(32.99*3.1415926/180)};
const double w_sheet[3] = {549.78,417.097,550.614};
const double r_layer[3] = {87.5026,132.7686,175.2686};
const double dmr_layer[3] = {79.838,125.104,167.604};
const int  l_layer[3] = {532,690,847};
const int  n_sheet[3] = {1,2,2};
const double x_reso[3]={0.1372,0.1476,0.1412};
const double v_reso[3]={0.1273,0.1326,0.1378};
KalFitGemHit::KalFitGemHit(double phi, double v, int layer,int sheet)
{
    /*
    myPhi=phi;
    myZ=z;
    myR=r;
    myErrPhi=0.023/r; //error in r*phi/r (rad)
    myErrZ=0.05; //error in z (cm) 
    myErrR=1.0; //error in R (cm)
    */
    myPhi=(phi+sheet*w_sheet[layer])/r_layer[layer];
    //cout << "pos phi=" << myPhi*dmr_layer[layer] << endl;
    //cout << "(v-(phi)*cos(a_stero[layer]))" << (v-(phi)*cos(a_stero[layer])) << endl;
    //cout << "pos v="   << myPhi*dmr_layer[layer]*cos(a_stero[layer])+(v-(phi)*cos(a_stero[layer])) << endl;
    //cout << "(phi+sheet*w_sheet[layer])/r_layer[layer]=" << (phi+sheet*w_sheet[layer])/r_layer[layer] << endl;
    //cout <<"phi=" << myPhi << endl;
    if(myPhi>3.1415926) myPhi=myPhi-6.2831852;
    //cout << "new phi=" << myPhi << endl;
    //cout << "first z=" << (v-(phi)*cos(a_stero[layer]))/sin(a_stero[layer]) << endl;
    myZ=(((v-(phi)*cos(a_stero[layer]))/sin(a_stero[layer]))-0.5*l_layer[layer])/10.;
    //cout << "second z="<< (v-(phi)*cos(a_stero[layer]))/sin(a_stero[layer])-0.5*l_layer[layer] << endl;
    //cout << "third z=" <<  myZ << endl;
    //cout << "v=" << v <<  endl;
    //cout << "l_layer=" << l_layer[layer] << endl;
    //cout << "0.5*l_layer=" << 0.5*l_layer[layer] << endl;
    
    myR=dmr_layer[layer]/10.;
    //double Dz,Cov,a,b;
    //a=1./(sin(a_stero[layer])*sin(a_stero[layer]));
    //b=1./(tan(a_stero[layer])*tan(a_stero[layer]));
    //Cov=-1./(sin(a_stero[layer])*sin(a_stero[layer]));
    //Dz=(a+2*Cov)*0.01467*0.01467+b*0.01349*0.01349;
    //myErrZ=sqrt(Dz);
    myErrZ=sqrt(0.01*v_reso[layer]*v_reso[layer]/(sin(a_stero[layer])*sin(a_stero[layer]))+0.01*x_reso[layer]*x_reso[layer]/(tan(a_stero[layer])*tan(a_stero[layer])));
    myErrPhi=x_reso[layer]/r_layer[layer];
}
 
KalFitGemHit::~KalFitGemHit() {;}