RegSelSiLUT.cxx File Reference

#include "RegSelLUT/RegSelSiLUT.h"
#include "RegSelLUT/RegSelModule.h"
#include "RegSelLUT/RegSelTimer.h"
#include "RegSelLUT/RegSelROB.h"
#include <iostream>
#include <fstream>
#include <set>

Include dependency graph for RegSelSiLUT.cxx:

Go to the source code of this file.

Functions
double	phiprime (double x, double y, double phi, double r, bool &status)
bool	philimits (double x, double y, double rmin, double rmax, double &phimin, double &phimax)

Function Documentation

philimits()

bool philimits	(	double	x,
		double	y,
		double	rmin,
		double	rmax,
		double &	phimin,
		double &	phimax
	)

Definition at line 403 of file RegSelSiLUT.cxx.

                                                                                                {
 
  double phi[4];
 
  bool status = true;
 
  //  std::cout << "phimin=" << phimin << "\tphimax=" << phimax << "\tlimits" << std::endl;
 
  phi[0] = phiprime(x, y, phimin, rmin, status);
  phi[1] = phiprime(x, y, phimin, rmax, status);
 
  phi[2] = phiprime(x, y, phimax, rmin, status);
  phi[3] = phiprime(x, y, phimax, rmax, status);
 
  //  for ( int i=0 ; i<4 ; i++  ) std::cout << "phi[" << i << "]=" << phi[i] << std::endl;
 
  if ( status ) { 
 
    // this is needed to take account of the phi=pi boundary, it's not rigorous, but it works
 
    if ( phi[0]*phi[1]<0 && (phi[0]>0.5*M_PI||phi[1]>0.5*M_PI) ) phimin = ( phi[0]<phi[1] ? phi[1] : phi[0] ); 
    else                                                         phimin = ( phi[0]<phi[1] ? phi[0] : phi[1] );  
  
    if ( phi[2]*phi[3]<0 && (phi[2]>0.5*M_PI||phi[3]>0.5*M_PI) ) phimax = ( phi[2]>phi[3] ? phi[3] : phi[2] ); 
    else                                                         phimax = ( phi[2]>phi[3] ? phi[2] : phi[3] );  
 
  }
  
  //  std::cout << "phimin=" << phimin << "\tphimax=" << phimax << "\tlimits (out)" << std::endl;
 
  return status;
 
}

phiprime()

double phiprime	(	double	x,
		double	y,
		double	phi,
		double	r,
		bool &	status
	)

Definition at line 371 of file RegSelSiLUT.cxx.

                                                                        {
 
  if ( x*x+y*y>r*r ) {
    //    std::cerr << "duffer x=" << x << "\ty=" << y << "\tr=" << r << std::endl; 
    status = false; 
    return 0;
  } 
 
  double a = (y*cos(phi)-x*sin(phi))/r;
 
  // this is bad, getting the return value to issue an error code.
  // ouch! should maybe define an exception and catch it?
  // this should never happen
  if ( fabs(a)>1 ) { status = false; return 0; } 
 
  status &= true;
 
  double phip = asin(a)+phi;
 
  //  std::cout << "phiprime()\tx=" << x << "\ty=" << y << "\tphi=" << phi << "\tr=" << r << "\tphiprime=" << phip << " before" << std::endl;
 
  // make sure wrapped between -pi to pi
  if ( phip>M_PI )  phip = phip-2*M_PI;
  if ( phip<-M_PI ) phip = phip+2*M_PI;
 
  //  std::cout << "phiprime()\tx=" << x << "\ty=" << y << "\tphi=" << phi << "\tr=" << r << "\tphiprime=" << phip << std::endl;
 
  return phip;
 
}