EMBAccordionDetails::Clockwork Class Reference

Collaboration diagram for EMBAccordionDetails::Clockwork:

Public Member Functions
int	phiGap (double radius, double xhit, double yhit) const
	compute number (0 to 1023) of closest electrode according to nominal accordion geometry More...
void	getRPhi ()
	initialization routine More...
double	phi0 (double radius) const
	phi of first absorber as function of radius for nominal accordion geometry (before sagging). More...
double	Distance_Ele (const double &x, const double &y, const int &PhiC, int &Num_Straight, const int &Num_Coude, double &xl)
	Compute the distance to the electrode. More...
double	Distance_Abs (const double &x, const double &y, const int &nabs, const int &Num_Straight, const int &Num_Coude)
	Compute the distance to the absorber. More...

Public Attributes
double	gam0
double	Rmin
double	Rmax
double	dR
double	Rphi [5000]
int	NRphi
int	Nbrt
int	Nbrt1
double	rint_eleFib
double	delta [15]
double	rc [15]
double	phic [15]
double	xc [15]
double	yc [15]
double	rMinAccordion
double	rMaxAccordion
double	zMinBarrel
double	zMaxBarrel
double	etaMaxBarrel
int	NCellTot
int	NCellMax
const GeoStraightAccSection *	electrodeStraightSection
const GeoStraightAccSection *	absorberStraightSection

Detailed Description

Definition at line 17 of file EMBAccordionDetails.cxx.

Member Function Documentation

Distance_Abs()

double EMBAccordionDetails::Clockwork::Distance_Abs	(	const double &	x,
		const double &	y,
		const int &	nabs,
		const int &	Num_Straight,
		const int &	Num_Coude
	)

Compute the distance to the absorber.

Distance_Ele()

double EMBAccordionDetails::Clockwork::Distance_Ele	(	const double &	x,
		const double &	y,
		const int &	PhiC,
		int &	Num_Straight,
		const int &	Num_Coude,
		double &	xl
	)

Compute the distance to the electrode.

getRPhi()

void EMBAccordionDetails::Clockwork::getRPhi

(

)

initialization routine

Definition at line 74 of file EMBAccordionDetails.cxx.

{
  const double dl=0.001;
  const double inv_dl = 1 / dl;
  double cenx[15],ceny[15];
  //double xl,xl2;
  double sum1[5000],sumx[5000];
  //xl=0;
  //xl2=0.;
  NRphi=5000;
  Rmin=1500.;
  dR=0.10;
  Rmax=0.;
 
  const double rint= rint_eleFib;
  const double inv_rint = 1. / rint;
  const double dt=dl * inv_rint;
  const double inv_dt = 1. / dt;
 
  for (int i=0;i<NRphi;i++) {
     sum1[i]=0.;
     sumx[i]=0.;
   }
  for (int i=0;i<15;i++) {
     cenx[i]=rc[i]*cos(phic[i]);
     ceny[i]=rc[i]*sin(phic[i]);
  }
 
  for (int i=0; i<15; i++) {
     
// fold
     double phi0,phi1;
     if (i==0) {
        phi0=-M_PI/2.;
        phi1=-M_PI/2.+delta[0];
     }
     else if (i==14) {
        phi0=-M_PI+delta[13];
        phi1=-M_PI/2.;
     }
     else {
        if (i%2==1) {
           phi0=delta[i];
           phi1=M_PI-delta[i-1];
        }
        else {
           phi0=-M_PI+delta[i-1];
           phi1=-delta[i];
        }
     }
     //xl2+=rint*fabs(phi1-phi0);
     int nstep=int((phi1-phi0) * inv_dt)+1;
//     std::cout << "fold " << phi0 << " " << phi1 << " " << nstep << std::endl;
     for (int ii=0;ii<nstep;ii++) {
        //xl+=dl;
        double phi=phi0+dt*((double)ii);
        double x=cenx[i]+rint*cos(phi);
        double y=ceny[i]+rint*sin(phi);
        double radius=sqrt(x*x+y*y);
        if (radius>Rmax) Rmax=radius;
        double phid=atan(y/x);
        int ir=((int) ((radius-Rmin)/dR) );
        if (ir>=0 && ir < NRphi) {
           sum1[ir]+=1.;
           sumx[ir]+=phid;
        }
     }
 
// straight section
     if (i<14) {
        double dx=cenx[i+1]-cenx[i];
        double dy=ceny[i+1]-ceny[i];
        double along=dx*dx+dy*dy-4.*rint*rint;
        along=sqrt(along);
        double x0=0.5*(cenx[i+1]+cenx[i]);
        double y0=0.5*(ceny[i+1]+ceny[i]);
        double phi;
        if (i%2==0) phi=M_PI/2-delta[i];
        else        phi=-M_PI/2.+delta[i];
        double x1=x0-0.5*along*cos(phi);
        double y1=y0-0.5*along*sin(phi);
        //xl2+=along;
        int nstep=int(along * inv_dl)+1;
//        std::cout << "straight" << x0 << " " << y0 << along << " " << nstep << std::endl;
        for (int ii=0;ii<nstep;ii++) {
           //xl+=dl;
           double x=x1+dl*((double)ii)*cos(phi);
           double y=y1+dl*((double)ii)*sin(phi);
           double radius=sqrt(x*x+y*y);
           if (radius>Rmax) Rmax=radius;
           double phid=atan(y/x);
           int ir=((int) ((radius-Rmin)/dR) );
           if (ir>=0 && ir < NRphi) {
              sum1[ir]+=1.;
              sumx[ir]+=phid;
           }
        }
     }
  }
//  std::cout << "total electrode length " << xl << " " << xl2 << std::endl;
//  std::cout << "rmax in accordion " << Rmax << std::endl;
  for (int i=0; i<NRphi; i++) {
    if (sum1[i]>0) {
     Rphi[i]=sumx[i]/sum1[i];
     // Not used: double radius = Rmin + ((double(i))+0.5)*dR;
//     std::cout << " r,phi0 " << radius << " " << Rphi[i] << std::endl;
    } 
    else Rphi[i]=0.;
  }
}

phi0()

double EMBAccordionDetails::Clockwork::phi0

(

double

radius

)

const

phi of first absorber as function of radius for nominal accordion geometry (before sagging).

Definition at line 286 of file EMBAccordionDetails.cxx.

{
 int ir;
 if (radius < Rmin) ir=0;
 else {
   if (radius > Rmax) radius=Rmax-0.0001;
   ir=((int) ((radius-Rmin)/dR) );
 }
 return Rphi[ir];
}

phiGap()

int EMBAccordionDetails::Clockwork::phiGap	(	double	radius,
		double	xhit,
		double	yhit
	)		const

compute number (0 to 1023) of closest electrode according to nominal accordion geometry

Definition at line 189 of file EMBAccordionDetails.cxx.

{
  const double m2pi = 2.0*Gaudi::Units::pi;
  double phi_0=phi0(radius)+gam0;   // from -M_PI to M_PI
  double phi_hit=atan2(yhit,xhit);  // from -M_PI to M_PI
  double dphi=phi_hit-phi_0;
// bring back to 0-2pi
  if (dphi<0) dphi=dphi+m2pi;
  if (dphi>=m2pi) dphi=dphi-m2pi;
  dphi=dphi*(1024/m2pi);
  int ngap=((int) dphi);
  return ngap;
}

Member Data Documentation

absorberStraightSection

const GeoStraightAccSection* EMBAccordionDetails::Clockwork::absorberStraightSection

Definition at line 70 of file EMBAccordionDetails.cxx.

delta

double EMBAccordionDetails::Clockwork::delta[15]

Definition at line 56 of file EMBAccordionDetails.cxx.

dR

double EMBAccordionDetails::Clockwork::dR

Definition at line 50 of file EMBAccordionDetails.cxx.

electrodeStraightSection

const GeoStraightAccSection* EMBAccordionDetails::Clockwork::electrodeStraightSection

Definition at line 69 of file EMBAccordionDetails.cxx.

etaMaxBarrel

double EMBAccordionDetails::Clockwork::etaMaxBarrel

Definition at line 65 of file EMBAccordionDetails.cxx.

gam0

double EMBAccordionDetails::Clockwork::gam0

Definition at line 47 of file EMBAccordionDetails.cxx.

Nbrt

int EMBAccordionDetails::Clockwork::Nbrt

Definition at line 53 of file EMBAccordionDetails.cxx.

Nbrt1

int EMBAccordionDetails::Clockwork::Nbrt1

Definition at line 54 of file EMBAccordionDetails.cxx.

NCellMax

int EMBAccordionDetails::Clockwork::NCellMax

Definition at line 67 of file EMBAccordionDetails.cxx.

NCellTot

int EMBAccordionDetails::Clockwork::NCellTot

Definition at line 66 of file EMBAccordionDetails.cxx.

NRphi

int EMBAccordionDetails::Clockwork::NRphi

Definition at line 52 of file EMBAccordionDetails.cxx.

phic

double EMBAccordionDetails::Clockwork::phic[15]

Definition at line 58 of file EMBAccordionDetails.cxx.

rc

double EMBAccordionDetails::Clockwork::rc[15]

Definition at line 57 of file EMBAccordionDetails.cxx.

rint_eleFib

double EMBAccordionDetails::Clockwork::rint_eleFib

Definition at line 55 of file EMBAccordionDetails.cxx.

Rmax

double EMBAccordionDetails::Clockwork::Rmax

Definition at line 49 of file EMBAccordionDetails.cxx.

rMaxAccordion

double EMBAccordionDetails::Clockwork::rMaxAccordion

Definition at line 62 of file EMBAccordionDetails.cxx.

Rmin

double EMBAccordionDetails::Clockwork::Rmin

Definition at line 48 of file EMBAccordionDetails.cxx.

rMinAccordion

double EMBAccordionDetails::Clockwork::rMinAccordion

Definition at line 61 of file EMBAccordionDetails.cxx.

Rphi

double EMBAccordionDetails::Clockwork::Rphi[5000]

Definition at line 51 of file EMBAccordionDetails.cxx.

xc

double EMBAccordionDetails::Clockwork::xc[15]

Definition at line 59 of file EMBAccordionDetails.cxx.

yc

double EMBAccordionDetails::Clockwork::yc[15]

Definition at line 60 of file EMBAccordionDetails.cxx.

zMaxBarrel

double EMBAccordionDetails::Clockwork::zMaxBarrel

Definition at line 64 of file EMBAccordionDetails.cxx.

zMinBarrel

double EMBAccordionDetails::Clockwork::zMinBarrel

Definition at line 63 of file EMBAccordionDetails.cxx.

---

The documentation for this class was generated from the following file:

• EMBAccordionDetails.cxx