BeamSpot Namespace Reference

Classes
struct	BeamSpotDB
	A simple struct to hold vertex info. More...
struct	Event
class	ID
struct	VrtHolder

Functions
void	myFCN_LLsolver (Int_t &, Double_t *, Double_t &, Double_t *, Int_t)
void	myFCN_LLsolverNorm (Int_t &, Double_t *, Double_t &, Double_t *, Int_t)
double	norm_xMin (-1e8)
double	norm_xMax (1e8)
double	norm_yMin (-1e8)
double	norm_yMax (1e8)
double	norm_zMin (-1e8)
double	norm_zMax (1e8)
double	pdfxy (double *x, double *p)

Variables
std::mutex	mutex
const std::vector< BeamSpot::VrtHolder > *vertexData	ATLAS_THREAD_SAFE

Function Documentation

myFCN_LLsolver()

void BeamSpot::myFCN_LLsolver	(	Int_t &	,
		Double_t *	,
		Double_t &	f,
		Double_t *	par,
		Int_t	)

Definition at line 1032 of file InDetBeamSpotVertex.cxx.

                                                                                     {
  constexpr double Pi = M_PI;
  //par[*]
  //0,  1,  2,   3,  4,  5,    6,       7,   8    9
  //X0, Y0, Ax.  Ay, sx, sy, rhoxy,   k     z0    sigma(z)
  
  f = 0;
  
  using Vertices = std::vector<BeamSpot::VrtHolder>;
  std::lock_guard<std::mutex> lock(BeamSpot::mutex);
  Vertices::const_iterator vit = BeamSpot::vertexData->begin();
  
  double temp =0;
  double x=0,y=0,z=0;
  double vxx,vyy, vxy;
  double covXX,covYY,covXY;
  double det,k2;
  
  // ln L = Sum[ ln(F) ]
  for ( ; vit != vertexData->end(); ++vit) {
    if (!vit->valid) continue; // don't use non-valid vertices
    temp =0;
    x = vit->x;
    y = vit->y;
    z = vit->z;
    vxx = vit->vxx;
    vxy = vit->vxy;
    vyy = vit->vyy;
    
    
    k2 = par[7]*par[7];
    
    covXX = k2 *vxx +  par[4]*par[4];
    covYY = k2 *vyy +  par[5]*par[5];
    covXY = k2 *vxy +  par[6] *par[4]* par[5];
    
    det = covXX * covYY - covXY*covXY;
    double recDet = 1./det;
 
    //temp =  TMath::Log(2*Pi * sqrt(std::abs(det)));
    temp  =  2*TMath::Log(2*Pi);
    temp  += TMath::Log(det);
    
    covXY = -covXY * recDet;
    double t = covXX *recDet;
    covXX = covYY *recDet;
    covYY = t;
    
    temp += ( 
              ( x - par[0] - par[2]*z) * covXX *  ( x - par[0] - par[2]*z) 
              +   ( y - par[1] - par[3]*z) * covYY *  ( y - par[1] - par[3]*z) 
              + 2*( x - par[0] - par[2]*z) * covXY *  ( y - par[1] - par[3]*z)
              );
    
    temp +=  TMath::Log( 2*Pi * par[9]*par[9] )  +  ( z - par[8]) * (z-par[8]) / (par[9] * par[9] );
    f+= 0.5*temp;
  }//for
  
  
}//myFCN

myFCN_LLsolverNorm()

void BeamSpot::myFCN_LLsolverNorm	(	Int_t &	,
		Double_t *	,
		Double_t &	,
		Double_t *	,
		Int_t	)

Definition at line 1100 of file InDetBeamSpotVertex.cxx.

                                                                                         {
 
}

norm_xMax()

double BeamSpot::norm_xMax

(

1e8

)

norm_xMin()

double BeamSpot::norm_xMin

(

-

1e8

)

norm_yMax()

double BeamSpot::norm_yMax

(

1e8

)

norm_yMin()

double BeamSpot::norm_yMin

(

-

1e8

)

norm_zMax()

double BeamSpot::norm_zMax

(

1e8

)

norm_zMin()

double BeamSpot::norm_zMin

(

-

1e8

)

pdfxy()

double BeamSpot::pdfxy	(	double *	x,
		double *	p )

Definition at line 1093 of file InDetBeamSpotVertex.cxx.

                                         {
  
  return 0; // TBD dlete
}

Variable Documentation

ATLAS_THREAD_SAFE

const std::vector< BeamSpot::VrtHolder >* vertexData BeamSpot::ATLAS_THREAD_SAFE

Definition at line 20 of file InDetBeamSpotVertex.cxx.

mutex

std::mutex BeamSpot::mutex

Definition at line 19 of file InDetBeamSpotVertex.cxx.