#include <CgemGeoReadoutPlane.h>

Public Member Functions
	CgemGeoReadoutPlane (int iLayer, int iSheet, double rx, double rv, double phi_min, double dX_strip, double dV, double w, double zmin, double L, int NXStrip, int NVStrip, double Xpitch, double XstripWidth, double Vpitch, double VstripWidth, double stereoAngle, double midROfGap, double outROfGap)

	~CgemGeoReadoutPlane ()

void	print ()

int	getLayerId ()

int	getSheetId ()

double	getRX ()

double	getRV ()

double	getPhimin ()

double	getXmin ()

double	getWidth ()

double	getZmin ()

double	getLength ()

double	getXPitch ()

double	getVPitch ()

double	getStereoAngle ()

double	getMidRAtGap ()

double	getVmax ()

int	getNXstrips () const

int	getNVstrips () const

bool	OnThePlane (double phi, double z) const

double	getX (double phi) const

int	getXStripID (double phi) const

int	getXStripID (double x, double y) const

int	getXStripID (G4ThreeVector pos) const

int	getClosestXStripID (double phi, double &dist)

int	getClosestVStripID (G4ThreeVector pos, double &dist) const

double	getDist2ClosestXStripCenter (double phi, int &id)

double	getDist2ClosestVStripCenter (G4ThreeVector pos, int &id)

void	getStripID (G4ThreeVector pos, int &X_ID, int &V_ID) const

void	getFiredStripID (G4ThreeVector pos1, G4ThreeVector pos2, vector< int > &vecXID, vector< int > &vecVID) const

double	getVFromLocalXZ (double x, double zeta, bool checkRange=true) const

double	getVFromPhiZ (double phi, double z, bool checkRange=true) const

double	getVFromPhiZ_nearPhiMin (double phi, double z, bool checkRange=true) const

double	getVInNextSheetFromV (double v, double phiminNext) const

int	getVIDFromV (double V) const

int	getVIDInNextSheetFromVID (int vID, double phimin_next) const

double	getPhiMin_strip () const

double	getCentralXFromXID (int X_ID) const

double	getPhiFromXID (int X_ID) const

double	getCentralVFromVID (int V_ID) const

double	getVStripLength (int V_ID) const

double	getZFromXV (double X, double V, int checkXRange=1, int checkVRange=1) const

double	getZFromPhiV (double phi, double V, int checkXRange=1) const

double	incidentAngleX (Hep3Vector momentum, double phi) const

double	incidentAngleV (Hep3Vector momentum, double phi) const

double	dAngleAcute (double phi1, double phi2) const

Detailed Description

Definition at line 14 of file CgemGeoReadoutPlane.h.

Constructor & Destructor Documentation

◆ CgemGeoReadoutPlane()

CgemGeoReadoutPlane::CgemGeoReadoutPlane	(	int	iLayer,
		int	iSheet,
		double	rx,
		double	rv,
		double	phi_min,
		double	dX_strip,
		double	dV,
		double	w,
		double	zmin,
		double	L,
		int	NXStrip,
		int	NVStrip,
		double	Xpitch,
		double	XstripWidth,
		double	Vpitch,
		double	VstripWidth,
		double	stereoAngle,
		double	midROfGap,
		double	outROfGap
	)

Definition at line 6 of file CgemGeoReadoutPlane.cxx.

                                           {
 
    m_iLayer = iLayer;
    m_iSheet = iSheet;
 
    m_RX = rx;
    m_RV = rv;
    m_Rmid = 0.5*(rx+rv);
    
    m_W = w;
    m_Zmin = zmin;
    m_Zmax = zmin+L;
    m_L = L;
    
    m_XPitch = Xpitch;            
    m_XstripWidth = XstripWidth;  
    m_VPitch = Vpitch;            
    m_VstripWidth = VstripWidth;  
    //m_XPitch = Xpitch*m_RX/m_Rmid;            // scaling from plane to cylinder
    //m_XstripWidth = XstripWidth*m_RX/m_Rmid;  // scaling from plane to cylinder
    //m_VPitch = Vpitch*m_RV/m_Rmid;            // scaling from plane to cylinder
    //m_VstripWidth = VstripWidth*m_RV/m_Rmid;  // scaling from plane to cylinder
    m_StereoAngle = stereoAngle / 180.*CLHEP::pi; //degree -> radian
    m_k = stereoAngle/fabs(stereoAngle); // sign of the stereo angle
    m_MidRGap = midROfGap;
    m_OutRGap = outROfGap;
    m_NXstrips = NXStrip;
    m_NVstrips = NVStrip;
 
    m_Phimin = phiMin;
    while(m_Phimin<-M_PI) m_Phimin+=2*M_PI;
    while(m_Phimin>M_PI)  m_Phimin-=2*M_PI;
    //m_Xmin = m_Phimin*rx;
    m_Xmin = m_Phimin*m_Rmid;
 
    //double dPhi_sensitive = NXStrip*m_XPitch/m_Rmid;
    //m_Phimin_strip = dX_strip/m_Rmid+m_Phimin;
    m_Phimin_strip = dX_strip/m_Rmid;
    //while(m_Phimin_strip<0)      m_Phimin_strip+=2*M_PI;
    //while(m_Phimin_strip>2*M_PI) m_Phimin_strip-=2*M_PI;
    m_Xmin_strip = m_Phimin_strip*m_Rmid;
    m_Xmax_strip = m_Xmin_strip+NXStrip*m_XPitch;
    m_Phimax_strip = m_Xmax_strip/m_Rmid;
    //if(NXStrip*m_XPitch>2*M_PI*rx/) cout<<"CgemGeoReadoutPlane : "<<endl;
 
    //m_Vmax = m_L*fabs(sin(m_StereoAngle))+m_W*m_RV/m_RX*cos(m_StereoAngle);
    //m_Vmax = m_L*fabs(sin(m_StereoAngle))+m_W*cos(m_StereoAngle);
    //m_NXstrips = floor(m_W / m_XPitch);
    //m_NVstrips = floor(m_Vmax / m_VPitch);
    //m_Vmin = 0.;
    
    m_Vtot = NXStrip*m_XPitch*cos(m_StereoAngle)+m_L*fabs(sin(m_StereoAngle));
    dV=(m_Vtot-NVStrip*m_VPitch)/2.0;// FIXME?
    m_Vmin_strip = dV;
    m_Vmax = m_Vmin_strip+NVStrip*m_VPitch;
 
 
    cout<<"    dX = "<<dX_strip<<endl;
    print();
}

◆ ~CgemGeoReadoutPlane()

CgemGeoReadoutPlane::~CgemGeoReadoutPlane ( )

Definition at line 73 of file CgemGeoReadoutPlane.cxx.

74{}

Member Function Documentation

◆ dAngleAcute()

double CgemGeoReadoutPlane::dAngleAcute	(	double	phi1,
		double	phi2
	)		const

Definition at line 469 of file CgemGeoReadoutPlane.cxx.

{
    double dPhi = phi1-phi2;
    while(dPhi<-CLHEP::pi) dPhi+=CLHEP::twopi;
    while(dPhi> CLHEP::pi) dPhi-=CLHEP::twopi;
    if(dPhi<-0.5*CLHEP::pi) dPhi+=CLHEP::pi;
    else if(dPhi>0.5*CLHEP::pi) dPhi-=CLHEP::pi;
    return dPhi;
}

Referenced by incidentAngleV(), and incidentAngleX().

◆ getCentralVFromVID()

double CgemGeoReadoutPlane::getCentralVFromVID ( int V_ID ) const

inline

Definition at line 78 of file CgemGeoReadoutPlane.h.

        { 
            //V_ID=V_ID-m_NVstrips*m_iSheet; 
            return (V_ID+0.5)*m_VPitch+m_Vmin_strip; 
        }

Referenced by getClosestVStripID(), getDist2ClosestVStripCenter(), getVIDInNextSheetFromVID(), getVStripLength(), CgemROOTGeo::InitGeometry(), and TestMapping::PositionCgemBoss().

◆ getCentralXFromXID()

double CgemGeoReadoutPlane::getCentralXFromXID ( int X_ID ) const

inline

Definition at line 73 of file CgemGeoReadoutPlane.h.

        {   //X_ID=X_ID-m_NXstrips*m_iSheet; 
            return (X_ID+0.5)*m_XPitch+m_Xmin_strip; 
        }// central X with m_Xmin as origin

Referenced by getClosestXStripID(), getDist2ClosestXStripCenter(), getPhiFromXID(), and TestMapping::PositionCgemBoss().

◆ getClosestVStripID()

int CgemGeoReadoutPlane::getClosestVStripID	(	G4ThreeVector	pos,
		double &	dist
	)		const

Definition at line 193 of file CgemGeoReadoutPlane.cxx.

{
    double z = pos.z();
    dist = -9999.;
    int V_ID=-99;
    if(z < m_Zmin){
        V_ID = -3;
    } else if(z > m_Zmin+m_L){
        V_ID = -4;
    } else {
        double phi = pos.phi(); 
        double v   = getVFromPhiZ(phi, z);
        //V_ID = floor((v-m_Vmin_strip) / m_VPitch);
        V_ID = getVIDFromV(v);
        int vid2 = V_ID;
        if(V_ID==-1) vid2=0;
        if(V_ID==-2) vid2=m_NVstrips-1;
        if(v>=0.) dist       = fabs(v-getCentralVFromVID(vid2))-0.5*m_VstripWidth;
    }
 
    return V_ID;
}

Referenced by getStripID().

◆ getClosestXStripID()

int CgemGeoReadoutPlane::getClosestXStripID	(	double	phi,
		double &	dist
	)

Definition at line 126 of file CgemGeoReadoutPlane.cxx.

{
    //cout<<"getClosestXStripID: phi0="<<phi<<endl;
    //cout<<"getClosestXStripID: phi="<<phi<<endl;
    //cout<<"getXStripID: "<<xID<<endl;
    int xID = getXStripID(phi);
    int xID2 = xID;
    if(xID2==-1) xID2=0;
    if(xID2==-2) xID2=m_NXstrips-1;
    double X = getX(phi);
    dist = fabs(X-getCentralXFromXID(xID2))-0.5*m_XstripWidth;
    //cout<<"getClosestXStripID: m_Xmin="<<m_Xmin<<", m_RX="<<m_RX<<", X="<<X
    //  <<", getCentralXFromXID(xID)="<<getCentralXFromXID(xID)
    //  <<", m_XstripWidth="<<m_XstripWidth<<endl;
    return xID;
}

◆ getDist2ClosestVStripCenter()

double CgemGeoReadoutPlane::getDist2ClosestVStripCenter	(	G4ThreeVector	pos,
		int &	id
	)

Definition at line 217 of file CgemGeoReadoutPlane.cxx.

                                                                                  {
    double z = pos.z();
    double dist=-9999.;
    int V_ID=-99;
    if(z < m_Zmin){
        V_ID = -3;
    } else if(z > m_Zmin+m_L){
        V_ID = -4;
    } else {
        double phi = pos.phi(); 
        double v   = getVFromPhiZ(phi, z);
        V_ID = getVIDFromV(v);
        int vid2 = V_ID;
        if(V_ID==-1) vid2=0;
        if(V_ID==-2) vid2=m_NVstrips-1;
        dist       = v-getCentralVFromVID(vid2);
    }
    id=V_ID;
 
    return dist;
}

Referenced by InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ getDist2ClosestXStripCenter()

double CgemGeoReadoutPlane::getDist2ClosestXStripCenter	(	double	phi,
		int &	id
	)

Definition at line 143 of file CgemGeoReadoutPlane.cxx.

{
    //cout<<"getClosestXStripID: phi0="<<phi<<endl;
    //phi = phi - m_Xmin/m_RX;
    //while(phi<0) phi += CLHEP::twopi;
    //while(phi>CLHEP::twopi)  phi-=CLHEP::twopi;
    //cout<<"getClosestXStripID: phi="<<phi<<endl;
    //double X = m_RX * phi;// in X-zeta plane (X-V plane)
    //if(X == m_W) X -= 0.001 * m_XPitch;
    //int xID = floor((X-(m_Xmin_strip-m_Xmin))/m_XPitch);
    //cout<<"getXStripID: "<<xID<<endl;
    //if(xID==m_NXstrips) xID=m_NXstrips-1;
    //if(xID<0) xID=0;
    //if(xID>=m_NXstrips) xID=m_NXstrips-1;
    double dist=-9999.;
    int xID = getXStripID(phi);
    id=xID;
    int xID2 = xID;
    if(xID2==-1) xID2=0;
    if(xID2==-2) xID2=m_NXstrips-1;
    if(xID2>=0&&xID2<m_NXstrips) {
        double X = getX(phi);
        dist = X-getCentralXFromXID(xID2);
        //cout<<"getClosestXStripID: m_Xmin="<<m_Xmin<<", m_RX="<<m_RX<<", X="<<X
        //  <<", getCentralXFromXID(xID)="<<getCentralXFromXID(xID)
        //  <<", m_XstripWidth="<<m_XstripWidth<<endl;
    }
    return dist;
}

Referenced by InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ getFiredStripID()

void CgemGeoReadoutPlane::getFiredStripID	(	G4ThreeVector	pos1,
		G4ThreeVector	pos2,
		vector< int > &	vecXID,
		vector< int > &	vecVID
	)		const

Definition at line 290 of file CgemGeoReadoutPlane.cxx.

                                                        {
    //cout << "get fired strip id" << endl;
    bool printDebug = false;
    double z[2];
    z[0] = pos1.z();
    z[1] = pos2.z();
    if( (z[0]<m_Zmin && z[1]<m_Zmin) || (z[0]>(m_Zmin+m_L) && z[1]>(m_Zmin+m_L)) ) return;
 
    double phi[2];
    double phimin=-CLHEP::pi;
    phi[0]=pos1.phi();
    phi[1]=pos2.phi();
    if(fabs(phi[1]-phi[0])>CLHEP::pi){
        phimin=0.0;
        for(int i=0; i<2; i++) if(phi[i]<0.0) phi[i]+=CLHEP::twopi;
    }
    if(phi[0]>phi[1]){
        double tmp=phi[0];
        phi[0]=phi[1];
        phi[1]=tmp;
        tmp=z[1];
        z[1]=z[0];
        z[0]=tmp;
    }
    if(printDebug) cout<<"phimin="<<phimin<<", phi[0,1]="<<phi[0]<<","<<phi[1]<<", z[0,1]="<<z[0]<<","<<z[1]<<endl;
    for(int i=0; i<2; i++){
        if(z[i]<m_Zmin){
            double phiTmp=(phi[1]-phi[0])/(z[1]-z[0])*(m_Zmin-z[0])+phi[0];
            phi[i]=phiTmp;
            z[i]=m_Zmin;
        }
        if(z[i]>(m_Zmin+m_L)){
            double phiTmp=(phi[1]-phi[0])/(z[1]-z[0])*(m_Zmin+m_L-z[0])+phi[0];
            phi[i]=phiTmp;
            z[i]=m_Zmin+m_L;
        }   
    }
    double phiActiveMin=m_Phimin+m_Phimin_strip+10e-10;// 2019-05-05 wangll, +10e-10 2019-10-30
    double phiActiveMax=m_Phimin+m_Phimax_strip-10e-10;// 2019-05-05 wangll, -10e-10 2019-10-30
    if(printDebug) cout<<"phimin="<<phimin<<", phiActiveMin="<<phiActiveMin<<", phiActiveMax="<<phiActiveMax<<", phi[0,1]="<<phi[0]<<","<<phi[1]<<", z[0,1]="<<z[0]<<","<<z[1]<<endl;
    while(phiActiveMax>=CLHEP::twopi) phiActiveMax-=CLHEP::twopi;// [0,2pi]
    while(phiActiveMin>=CLHEP::twopi) phiActiveMin-=CLHEP::twopi;// [0,2pi]
    while(phiActiveMax<0            ) phiActiveMax+=CLHEP::twopi;// [0,2pi]
    while(phiActiveMin<0            ) phiActiveMin+=CLHEP::twopi;// [0,2pi]
    if(phimin==-CLHEP::pi) 
    {
        /*if(phiActiveMin>CLHEP::pi||phiActiveMax>CLHEP::pi) {
          phiActiveMin-=CLHEP::twopi;
          phiActiveMax-=CLHEP::twopi;
          }*/
        while(phiActiveMin>CLHEP::pi) phiActiveMin-=CLHEP::twopi;
        while(phiActiveMax>CLHEP::pi) phiActiveMax-=CLHEP::twopi;
    }// [-pi,pi]
    if(printDebug) cout<<"phimin="<<phimin<<", phiActiveMin="<<phiActiveMin<<", phiActiveMax="<<phiActiveMax<<", phi[0,1]="<<phi[0]<<","<<phi[1]<<", z[0,1]="<<z[0]<<","<<z[1]<<endl;
 
    int iCase=-1;
    if(phiActiveMin<=phi[0])
    {
        if(phiActiveMax<=phiActiveMin) iCase=1;
        else if(phiActiveMax<=phi[0])  iCase=2;
        else if(phiActiveMax<=phi[1])  iCase=3;
        else                           iCase=4;
    }
    else if(phiActiveMin<=phi[1])
    {
        if(phiActiveMax<=phi[0])             iCase=5;
        else if(phiActiveMax<=phiActiveMin)  iCase=6;
        else if(phiActiveMax<=phi[1])        iCase=7;
        else                                 iCase=8;
    }
    else
    {
        if(phiActiveMax<=phi[0])             iCase=9;
        else if(phiActiveMax<=phi[1])        iCase=10;
        else if(phiActiveMax<=phiActiveMin)  iCase=11;
        else                                 iCase=12;
    }
    if(printDebug) cout<<"iCase = "<<iCase<<endl;
 
    double zBoundary[2];
    zBoundary[0]=z[0]+(z[1]-z[0])/(phi[1]-phi[0])*(phiActiveMin-phi[0]);
    zBoundary[1]=z[0]+(z[1]-z[0])/(phi[1]-phi[0])*(phiActiveMax-phi[0]);
    switch(iCase)
    {
        case 2:
        case 9:
        case 12:
            return;// no cluster
        case 1:
        case 4:
        case 6:
        case 11:
            break;// no change in points
        case 3:
        case 10:
            phi[1]=phiActiveMax;
            z[1]=zBoundary[1];
            break;
        case 5:
        case 8:
            phi[0]=phiActiveMin;
            z[0]=zBoundary[0];
            break;
        case 7:
            phi[0]=phiActiveMin;
            phi[1]=phiActiveMax;
            z[0]=zBoundary[0];
            z[1]=zBoundary[1];
            break;
        default:
            cout<<"Error in CgemGeoReadoutPlane::getFiredStripID: iCase = "<<iCase<<endl;
            return;
    }
    if(printDebug) cout<<"phimin="<<phimin<<", phiActiveMin="<<phiActiveMin<<", phiActiveMax="<<phiActiveMax<<", phi[0,1]="<<phi[0]<<","<<phi[1]<<", z[0,1]="<<z[0]<<","<<z[1]<<endl;
 
    G4ThreeVector pos_1,pos_2;
    pos_1.setRhoPhiZ(m_Rmid,phi[0],z[0]);
    pos_2.setRhoPhiZ(m_Rmid,phi[1],z[1]);
    int XID_min(-9),XID_max(-9);
    int VID_min(-9),VID_max(-9);
    getStripID(pos_1,XID_min,VID_min);
    getStripID(pos_2,XID_max,VID_max);
    if( iCase!=6 && VID_max<VID_min ) {
        int tmp=VID_max;
        VID_max=VID_min;
        VID_min=tmp;
    }
    if(printDebug) {
        cout<<"on the sheet: phi "<<phi[0]<<" to "<<phi[1]<<", z "<<z[0]<<" to "<<z[1]<<endl;
        cout<<"              XID "<<XID_min<<" to "<<XID_max<<", VID "<<VID_min<<" to "<<VID_max<<endl;
    }
 
    if(iCase!=6)
    {
        for(int i=0; i<=abs(XID_min-XID_max);i++) vecXID.push_back(XID_min+i);
        for(int i=0; i<=abs(VID_min-VID_max);i++) vecVID.push_back(VID_min+i);
    }
    else
    {
        pos_1.setRhoPhiZ(m_Rmid,phiActiveMax,zBoundary[1]);
        pos_2.setRhoPhiZ(m_Rmid,phiActiveMin,zBoundary[0]);
        int XID_1(-9),XID_2(-9);
        int VID_1(-9),VID_2(-9);
        getStripID(pos_1,XID_1,VID_1);
        getStripID(pos_2,XID_2,VID_2);
        for(int i=0; i<=abs(XID_min-XID_1);i++) vecXID.push_back(XID_min+i);
        for(int i=0; i<=abs(XID_max-XID_2);i++) vecXID.push_back(XID_2+i);
        int VID_0=min(VID_1,VID_min);
        for(int i=0; i<=abs(VID_min-VID_1);i++) vecVID.push_back(VID_0+i);
        VID_0=min(VID_2,VID_max);
        for(int i=0; i<=abs(VID_max-VID_2);i++) vecVID.push_back(VID_0+i);
    }
}

Referenced by BesCgemDigitizer::Digitize_v2().

◆ getLayerId()

int CgemGeoReadoutPlane::getLayerId ( )

inline

Definition at line 28 of file CgemGeoReadoutPlane.h.

28{return m_iLayer;};

Referenced by TestMapping::CheckCrossingCgemBoss(), TestMapping::CheckFEBPositionCgemBoss(), TestGeometry::initialize(), and TestMapping::PositionCgemBoss().

◆ getLength()

double CgemGeoReadoutPlane::getLength ( )

inline

Definition at line 36 of file CgemGeoReadoutPlane.h.

36{return m_L;};

Referenced by TestMapping::CheckCrossingCgemBoss(), TestHit::execute(), TestTrack::execute(), CgemROOTGeo::InitGeometry(), and TestGeometry::initialize().

◆ getMidRAtGap()

double CgemGeoReadoutPlane::getMidRAtGap ( )

inline

Definition at line 40 of file CgemGeoReadoutPlane.h.

40{return m_MidRGap;};

Referenced by TestCluster::execute(), TestClusterWithHit::execute(), TestTrack::execute(), CgemHitOnTrack::getFitStuff(), CgemROOTGeo::InitGeometry(), CgemCosmicRayQA::initialize(), and TestTrack::initialize().

◆ getNVstrips()

int CgemGeoReadoutPlane::getNVstrips ( ) const

inline

Definition at line 43 of file CgemGeoReadoutPlane.h.

43{ return m_NVstrips; };

Referenced by TestMapping::CheckCrossingCgemBoss(), TestMapping::CheckFEBPositionCgemBoss(), CgemROOTGeo::InitGeometry(), TestGeometry::initialize(), TestMapping::PositionCgemBoss(), InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ getNXstrips()

int CgemGeoReadoutPlane::getNXstrips ( ) const

inline

Definition at line 42 of file CgemGeoReadoutPlane.h.

42{ return m_NXstrips; };

Referenced by TestMapping::CheckCrossingCgemBoss(), TestMapping::CheckFEBPositionCgemBoss(), CgemROOTGeo::InitGeometry(), TestGeometry::initialize(), TestMapping::PositionCgemBoss(), InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ getPhiFromXID()

double CgemGeoReadoutPlane::getPhiFromXID ( int X_ID ) const

Definition at line 267 of file CgemGeoReadoutPlane.cxx.

                                                        { // from -pi to pi
    double X=getCentralXFromXID(X_ID);// central X with m_Xmin as origin
    double phi=(X+m_Xmin)/m_Rmid;
    while(phi<-M_PI) phi+=CLHEP::twopi;
    while(phi>M_PI)  phi-=CLHEP::twopi;
    return phi;
}

Referenced by TestMapping::CheckFEBPositionCgemBoss(), CgemROOTGeo::InitGeometry(), and TestMapping::PositionCgemBoss().

◆ getPhimin()

double CgemGeoReadoutPlane::getPhimin ( )

inline

Definition at line 32 of file CgemGeoReadoutPlane.h.

32{return m_Phimin;};

◆ getPhiMin_strip()

double CgemGeoReadoutPlane::getPhiMin_strip ( ) const

inline

Definition at line 72 of file CgemGeoReadoutPlane.h.

72{return m_Phimin_strip+m_Phimin;}

◆ getRV()

double CgemGeoReadoutPlane::getRV ( )

inline

Definition at line 31 of file CgemGeoReadoutPlane.h.

31{return m_RV;};

Referenced by TestMapping::CheckCrossingCgemBoss(), TestCluster::execute(), TestClusterWithHit::execute(), CgemHitOnTrack::getFitStuff(), CgemROOTGeo::InitGeometry(), CgemSegmentFitAlg::initialize(), and TestMapping::PositionCgemBoss().

◆ getRX()

double CgemGeoReadoutPlane::getRX ( )

inline

Definition at line 30 of file CgemGeoReadoutPlane.h.

30{return m_RX;};

Referenced by TestMapping::CheckCrossingCgemBoss(), TestCluster::execute(), TestClusterWithHit::execute(), CgemHitOnTrack::getFitStuff(), CgemROOTGeo::InitGeometry(), CgemSegmentFitAlg::initialize(), and TestMapping::PositionCgemBoss().

◆ getSheetId()

int CgemGeoReadoutPlane::getSheetId ( )

inline

Definition at line 29 of file CgemGeoReadoutPlane.h.

29{return m_iSheet;};

Referenced by TestMapping::CheckCrossingCgemBoss(), TestMapping::CheckFEBPositionCgemBoss(), TestGeometry::initialize(), and TestMapping::PositionCgemBoss().

◆ getStereoAngle()

double CgemGeoReadoutPlane::getStereoAngle ( )

inline

Definition at line 39 of file CgemGeoReadoutPlane.h.

39{return m_StereoAngle;};

Referenced by TestMapping::CheckCrossingCgemBoss(), CgemHitOnTrack::getFitStuff(), CgemLineFit::GetIntersection(), and CgemROOTGeo::InitGeometry().

◆ getStripID()

void CgemGeoReadoutPlane::getStripID	(	G4ThreeVector	pos,
		int &	X_ID,
		int &	V_ID
	)		const

Definition at line 240 of file CgemGeoReadoutPlane.cxx.

                                                                                  {
    /*
    double z = pos.z();
    if(z < m_Zmin){
        V_ID = -2;
    } else if(z > m_Zmin+m_L){
        V_ID = -1;
    } else if(X_ID >= 0){
        double x = pos.x();
        double y = pos.y();
        double phi = atan2(y,x);
        phi = phi - m_Xmin/m_RX;
        while(phi < 0)             phi+=CLHEP::twopi;
        while(phi > CLHEP::twopi)  phi-=CLHEP::twopi;
        double X = m_RX * phi;// in X-zeta plane (X-V plane) @ m_RX
        double zeta = z - m_Zmin;
        V_ID = floor((getVFromLocalXZ(X, zeta)-m_Vmin_strip) / m_VPitch);
    } else{
        V_ID = -3;
    }
    */
    
    X_ID = getXStripID(pos);
    double dist=-9999.;
    V_ID = getClosestVStripID(pos, dist);
}

Referenced by TestMapping::CheckCrossingCgemBoss(), getFiredStripID(), and InductionGTS::setMultiElectrons().

◆ getVFromLocalXZ()

double CgemGeoReadoutPlane::getVFromLocalXZ	(	double	x,
		double	zeta,
		bool	checkRange = `true`
	)		const

inline

Definition at line 163 of file CgemGeoReadoutPlane.h.

{
    double V = -9999.0;
    if( (!checkRange) || (x>=0. && x<=m_NXstrips*m_XPitch && zeta>=0 && zeta<=m_L) )
    {
        //double X = x/m_RX*m_RV;// X at V plane
        //if(m_StereoAngle>0) return zeta*sin(m_StereoAngle) + X*cos(m_StereoAngle);
        //else return (zeta-m_L)*sin(m_StereoAngle)+X*cos(m_StereoAngle);
        double z0=m_Zmax; if(m_k<0) z0=m_Zmin;
        //return X*cos(m_StereoAngle)+((m_k-1.0)*m_L/2.0-zeta)*sin(m_StereoAngle);
        V = x*cos(m_StereoAngle)+(z0-zeta-m_Zmin)*sin(m_StereoAngle);
    } 
    return V;
}

Referenced by getVFromPhiZ(), and getVFromPhiZ_nearPhiMin().

◆ getVFromPhiZ()

double CgemGeoReadoutPlane::getVFromPhiZ	(	double	phi,
		double	z,
		bool	checkRange = `true`
	)		const

inline

Definition at line 178 of file CgemGeoReadoutPlane.h.

{
    double dPhi=phi-m_Phimin;
    while(dPhi<0) dPhi+=CLHEP::twopi;
    while(dPhi>=CLHEP::twopi) dPhi-=CLHEP::twopi;
    double X=dPhi*m_Rmid;// X with m_Xmin as origin
    double x=X-m_Xmin_strip;// x with m_Xmin_strip as origin
    double zeta=z-m_Zmin;// local zeta [0,L]
    return getVFromLocalXZ(x,zeta, checkRange);
}

Referenced by DotsHelixFitter::calculateNewHelix(), getClosestVStripID(), getDist2ClosestVStripCenter(), HoughHit::residual(), and CgemMidDriftPlane::xyzToPhiVIdealGeom().

◆ getVFromPhiZ_nearPhiMin()

double CgemGeoReadoutPlane::getVFromPhiZ_nearPhiMin	(	double	phi,
		double	z,
		bool	checkRange = `true`
	)		const

inline

Definition at line 189 of file CgemGeoReadoutPlane.h.

{
    double dPhi=phi-m_Phimin;
    while(dPhi<-CLHEP::pi) dPhi+=CLHEP::twopi;
    while(dPhi>=CLHEP::pi) dPhi-=CLHEP::twopi;
    double X=dPhi*m_Rmid;// X with m_Xmin as origin
    double x=X-m_Xmin_strip;// x with m_Xmin_strip as origin
    double zeta=z-m_Zmin;// local zeta [0,L]
    return getVFromLocalXZ(x,zeta, checkRange);
}

Referenced by DotsHelixFitter::calculateNewHelix().

◆ getVIDFromV()

int CgemGeoReadoutPlane::getVIDFromV ( double V ) const

inline

Definition at line 220 of file CgemGeoReadoutPlane.h.

                                                          {
        double Vid = floor((V-m_Vmin_strip)/m_VPitch);
        if(Vid<0)           Vid=-1;
        if(Vid>=m_NVstrips) Vid=-2;
        //Vid=Vid+m_NVstrips*m_iSheet;
        return Vid;
}

Referenced by getClosestVStripID(), getDist2ClosestVStripCenter(), and getVIDInNextSheetFromVID().

◆ getVIDInNextSheetFromVID()

int CgemGeoReadoutPlane::getVIDInNextSheetFromVID	(	int	vID,
		double	phimin_next
	)		const

inline

Definition at line 200 of file CgemGeoReadoutPlane.h.

                                                                                          {
    int vId_next = -1;
    //vID=vID-m_NVstrips*m_iSheet;
    if( (vID>=0) && (vID<m_NVstrips) ){
        double v = getCentralVFromVID(vID);
        double v2 = getVInNextSheetFromV(v, phimin_next);
        //if(v2>=0.) return floor(v2/m_VPitch);
        vId_next = getVIDFromV(v2);// assume the next sheet with the same configuration
        //if(m_iLayer>0) { vId_next=vId_next+(1-2*m_iSheet)*m_NVstrips;// with assumption Nsheet=2 for layerID>0 FIXME}
    }
    return vId_next;
}

◆ getVInNextSheetFromV()

double CgemGeoReadoutPlane::getVInNextSheetFromV	(	double	v,
		double	phiminNext
	)		const

inline

Definition at line 61 of file CgemGeoReadoutPlane.h.

                                                                       { 
            double dPhi = phiminNext-m_Phimin;// assume the next m_Phimin_strip is the same
            while(dPhi<=0) dPhi+=CLHEP::twopi;
            while(dPhi>CLHEP::twopi) dPhi-=CLHEP::twopi;
            //cout<<"getVInNextSheetFromV: dphi = "<<dPhi<<endl;
            double dX = dPhi*m_Rmid;
            return v-dX*cos(m_StereoAngle); 
        }

Referenced by getVIDInNextSheetFromVID().

◆ getVmax()

double CgemGeoReadoutPlane::getVmax ( )

inline

Definition at line 41 of file CgemGeoReadoutPlane.h.

41{return m_Vmax;};

Referenced by CgemROOTGeo::InitGeometry().

◆ getVPitch()

double CgemGeoReadoutPlane::getVPitch ( )

inline

Definition at line 38 of file CgemGeoReadoutPlane.h.

38{return m_VPitch;};

Referenced by CgemROOTGeo::InitGeometry(), InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ getVStripLength()

double CgemGeoReadoutPlane::getVStripLength ( int V_ID ) const

inline

Definition at line 83 of file CgemGeoReadoutPlane.h.

        {
            double vStripLegth(0.);
            if(V_ID>=0&&V_ID<m_NVstrips) 
            {
                double V=getCentralVFromVID(V_ID);
                double Lz=V/fabs(sin(m_StereoAngle));
                if(Lz<=m_L) 
                {
                    vStripLegth=Lz/cos(m_StereoAngle);
                }
                else if(V/cos(m_StereoAngle)<=m_NXstrips*m_XPitch)
                {
                    vStripLegth=m_L/cos(m_StereoAngle);
                }
                else vStripLegth=(m_Vtot-V)/fabs(sin(m_StereoAngle)*cos(m_StereoAngle));
            }
            return vStripLegth;
        }

Referenced by TestHit::execute(), TestTrack::execute(), and TestGeometry::initialize().

◆ getWidth()

double CgemGeoReadoutPlane::getWidth ( )

inline

Definition at line 34 of file CgemGeoReadoutPlane.h.

34{return m_W;};

Referenced by CgemROOTGeo::InitGeometry().

◆ getX()

double CgemGeoReadoutPlane::getX ( double phi ) const

Definition at line 91 of file CgemGeoReadoutPlane.cxx.

{
    double dPhi = phi - m_Phimin;
    while(dPhi<0)              dPhi += CLHEP::twopi;
    while(dPhi>=CLHEP::twopi)  dPhi -= CLHEP::twopi;
    if(dPhi>m_Phimax_strip) {
        double dphi_start = phi-(m_Phimin+m_Phimin_strip);
        while(dphi_start> 0)        dphi_start-=2.*M_PI;
        while(dphi_start<=-2.*M_PI) dphi_start+=2.*M_PI;
        //double dphi_end   = phi-(m_Phimin+m_Phimax_strip);
        //while(dphi_end< 0)      dphi_end+=2.*M_PI;
        //while(dphi_end>=2*M_PI) dphi_end-=2.*M_PI;
        double dphi_end = dPhi-m_Phimax_strip;
        if(dphi_start+dphi_end>0) dPhi=dphi_start+m_Phimin_strip;
    }
    double X = m_Rmid*dPhi;// in X-zeta plane (X-V plane)
    return X;
}

Referenced by getClosestXStripID(), getDist2ClosestXStripCenter(), and getXStripID().

◆ getXmin()

double CgemGeoReadoutPlane::getXmin ( )

inline

Definition at line 33 of file CgemGeoReadoutPlane.h.

33{return m_Xmin;};

Referenced by CgemROOTGeo::InitGeometry().

◆ getXPitch()

double CgemGeoReadoutPlane::getXPitch ( )

inline

Definition at line 37 of file CgemGeoReadoutPlane.h.

37{return m_XPitch;};

Referenced by CgemROOTGeo::InitGeometry(), InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ getXStripID() [1/3]

int CgemGeoReadoutPlane::getXStripID ( double phi ) const

Definition at line 111 of file CgemGeoReadoutPlane.cxx.

                                                     {
    double X = getX(phi);
    //cout<<"X="<<X<<", Xmin_strip="<<m_Xmin_strip<<endl;
    X -= m_Xmin_strip;
    //cout<<"X="<<X<<endl;
    int xID = floor(X/m_XPitch);
    //cout<<"getXStripID: "<<xID<<endl;
    if(xID<0) xID=-1;
    //if(xID>=m_NXstrips) xID=m_NXstrips-1;
    if(xID>=m_NXstrips) xID=-2;
    //if(xID>=0&&xID<m_NXstrips) return xID;
    //else return -1;
    return xID;
}

Referenced by getClosestXStripID(), getDist2ClosestXStripCenter(), getStripID(), and getXStripID().

◆ getXStripID() [2/3]

int CgemGeoReadoutPlane::getXStripID	(	double	x,
		double	y
	)		const

Definition at line 175 of file CgemGeoReadoutPlane.cxx.

                                                             {
    if(x!=0||y!=0){
        double phi=atan2(y,x);
        return getXStripID(phi);
    } else{
        return -3;
    }
}

◆ getXStripID() [3/3]

int CgemGeoReadoutPlane::getXStripID ( G4ThreeVector pos ) const

Definition at line 185 of file CgemGeoReadoutPlane.cxx.

                                                            {
    double x = pos.x();
    double y = pos.y();
    return getXStripID(x,y);
}

◆ getZFromPhiV()

double CgemGeoReadoutPlane::getZFromPhiV	(	double	phi,
		double	V,
		int	checkXRange = `1`
	)		const

inline

Definition at line 213 of file CgemGeoReadoutPlane.h.

                                                                                           {
    phi = phi-m_Phimin;
    while(phi<0)             phi+=CLHEP::twopi;
    while(phi>=CLHEP::twopi) phi-=CLHEP::twopi;
    return getZFromXV(phi*m_Rmid, V, checkXRange);
}

◆ getZFromXV()

double CgemGeoReadoutPlane::getZFromXV	(	double	X,
		double	V,
		int	checkXRange = `1`,
		int	checkVRange = `1`
	)		const

Definition at line 275 of file CgemGeoReadoutPlane.cxx.

                                                                                                 {
    double localX = X-m_Xmin_strip;// with m_Xmin_strip as origin
    double z = -9999.0;
    if( (checkVRange==0||((V>=m_Vmin_strip)&&(V<m_Vmax))) && ((checkXRange==0)||((localX>=0.0)&&(localX<=m_Xmax_strip-m_Xmin_strip))) ){
        //double localXOfVstrip = V / cos(m_StereoAngle);
        //double zeta = (localXOfVstrip - localX) / tan(m_StereoAngle);
        //if( m_StereoAngle < 0 ) zeta += m_L;
        //if( (zeta>=0) && (zeta<=m_L) ) 
        //z = zeta + m_Zmin;
        double z0=m_Zmax; if(m_k<0) z0=m_Zmin;
        z = localX/tan(m_StereoAngle) - V/sin(m_StereoAngle) + z0;
    }
    return z;
}

Referenced by getZFromPhiV().

◆ getZmin()

double CgemGeoReadoutPlane::getZmin ( )

inline

Definition at line 35 of file CgemGeoReadoutPlane.h.

35{return m_Zmin;};

Referenced by TestMapping::CheckCrossingCgemBoss(), and CgemROOTGeo::InitGeometry().

◆ incidentAngleV()

double CgemGeoReadoutPlane::incidentAngleV	(	Hep3Vector	momentum,
		double	phi
	)		const

Definition at line 458 of file CgemGeoReadoutPlane.cxx.

{
    Hep3Vector zaxis(0,0,1);
    momentum.rotate(zaxis, -phi);
    Hep3Vector xaxis(1,0,0);
    momentum.rotate(xaxis, m_StereoAngle);
    double phiTrk = momentum.phi();
    double dPhi = dAngleAcute(phiTrk,0);
    return dPhi;
}

◆ incidentAngleX()

double CgemGeoReadoutPlane::incidentAngleX	(	Hep3Vector	momentum,
		double	phi
	)		const

Definition at line 445 of file CgemGeoReadoutPlane.cxx.

{
    double phiTrk = momentum.phi();
    /*double dPhi = phiTrk-phi;
    while(dPhi<-CLHEP::pi) dPhi+=CLHEP::twopi;
    while(dPhi> CLHEP::pi) dPhi-=CLHEP::twopi;
    if(dPhi<-0.5*CLHEP::pi) dPhi+=CLHEP::pi;
    else if(dPhi>0.5*CLHEP::pi) dPhi-=CLHEP::pi;
    */
    double dPhi = dAngleAcute(phiTrk,phi);
    return dPhi;
}

◆ OnThePlane()

bool CgemGeoReadoutPlane::OnThePlane	(	double	phi,
		double	z
	)		const

inline

Definition at line 154 of file CgemGeoReadoutPlane.h.

{
    double dPhi=phi-m_Phimin;// [0, 2pi) with m_Phimin as origin
    while(dPhi<0)       dPhi+=CLHEP::twopi;
    while(dPhi>=2*M_PI)  dPhi-=CLHEP::twopi;
    if(dPhi>=m_Phimin_strip && dPhi<=m_Xmax_strip/m_Rmid && z>=m_Zmin && z<=m_Zmax) return true;
    else return false;
}

Referenced by InductionGar2::setMultiElectrons(), and InductionGar::setMultiElectrons().

◆ print()

void CgemGeoReadoutPlane::print ( )

Definition at line 76 of file CgemGeoReadoutPlane.cxx.

                               {
    cout<<"construct one readout-plane with: "<<endl;
    cout<<"    iLayer "<<m_iLayer<<", iSheet "<<m_iSheet<<endl;
    cout<<"    Rx, Rv, Xmin, width = "<< m_RX<<", "<<m_RV<<", "<<m_Xmin<<", "<<m_W<<endl;
    cout<<"    Zmin, L = "<< m_Zmin <<", "<<m_L<<endl;
    cout<<"    XPitch, XstripWidth ="<<m_XPitch<<", "<<m_XstripWidth<<endl;
    cout<<"    sensitive X (Xmin as origin): "<<m_Xmin_strip<<" to "<<m_Xmax_strip<<endl;
    cout<<"    sensitive range in phi (phimin as origin): "<<m_Phimin_strip<<" to "<<m_Xmax_strip/m_Rmid<<endl;
    //cout<<"    VPitch, VstripWidth, StereoAngle, Vmax ="<<m_VPitch<<", "<<m_VstripWidth<<", "<<m_StereoAngle<<", "<<m_Vmax<<endl;
    cout<<"    VPitch, VstripWidth, StereoAngle, k ="<<m_VPitch<<", "<<m_VstripWidth<<", "<<m_StereoAngle<<", "<<m_k<<endl;
    cout<<"    Vtot, Vmin, Vmax = "<<m_Vtot<<", "<<m_Vmin_strip<<", "<<m_Vmax<<endl;
    cout<<"    N_Xstrips, N_Vstrips="<<m_NXstrips<<", "<<m_NVstrips<<endl;
    cout<<"    MidRGap="<<m_MidRGap<<", OutRGap="<<m_OutRGap<<endl;
}

Referenced by CgemGeoReadoutPlane().

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ CgemGeoReadoutPlane()

◆ ~CgemGeoReadoutPlane()

Member Function Documentation

◆ dAngleAcute()

◆ getCentralVFromVID()

◆ getCentralXFromXID()

◆ getClosestVStripID()

◆ getClosestXStripID()

◆ getDist2ClosestVStripCenter()

◆ getDist2ClosestXStripCenter()

◆ getFiredStripID()

◆ getLayerId()

◆ getLength()

◆ getMidRAtGap()

◆ getNVstrips()

◆ getNXstrips()

◆ getPhiFromXID()

◆ getPhimin()

◆ getPhiMin_strip()

◆ getRV()

◆ getRX()

◆ getSheetId()

◆ getStereoAngle()

◆ getStripID()

◆ getVFromLocalXZ()

◆ getVFromPhiZ()

◆ getVFromPhiZ_nearPhiMin()

◆ getVIDFromV()

◆ getVIDInNextSheetFromVID()

◆ getVInNextSheetFromV()

◆ getVmax()

◆ getVPitch()

◆ getVStripLength()

◆ getWidth()

◆ getX()

◆ getXmin()

◆ getXPitch()

◆ getXStripID() [1/3]

◆ getXStripID() [2/3]

◆ getXStripID() [3/3]

◆ getZFromPhiV()

◆ getZFromXV()

◆ getZmin()

◆ incidentAngleV()

◆ incidentAngleX()

◆ OnThePlane()

◆ print()