MdcGeomSvc Class Reference

#include <MdcGeomSvc.h>

Inheritance diagram for MdcGeomSvc:

List of all members.

Public Member Functions
	MdcGeomSvc (const std::string &name, ISvcLocator *svcloc)
	~MdcGeomSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	handle (const Incident &inc)
	this handle function is prepared for special use
const MdcGeoWire *const	Wire (unsigned id)
const MdcGeoWire *const	Wire (unsigned lyrid, unsigned wirid)
const MdcGeoLayer *const	Layer (unsigned id)
const MdcGeoSuper *const	SuperLayer (unsigned id)
const MdcGeoGeneral *const	GeneralLayer (unsigned id)
const MdcGeoMisc *const	Misc (void)
const MdcGeoEnd *const	End (unsigned id)
const int	getWireSize ()
const int	getLayerSize ()
const int	getSuperLayerSize ()
const int	getGeneralLayerSize ()
const int	getSegmentNo ()
void	Dump ()
Static Public Member Functions
static bool	getSagFlag (void)
static const InterfaceID &	interfaceID ()
Static Public Attributes
static bool	m_doSag = true
static bool	m_readAlignParDataBase = true
static bool	m_nomcalignment = true
Private Member Functions
void	clean ()
void	ReadFilePar ()
void	ReadTensionDataBase (std::vector< double > &wireTensionVec)
void	ReadWirePosDataBase (std::vector< vector< double > > &wirePosVec)
void	ReadAliParDataBase (vector< double > &Sx, vector< double > &Sy, vector< double > &Sz, vector< double > &Rx, vector< double > &Ry, vector< double > &Rz)
void	Fill ()
const int	getAlignParIndexEast (int lyr) const
const int	getAlignParIndexWest (int lyr) const
Private Attributes
double	m_wholeShiftX
double	m_wholeShiftY
double	m_wholeShiftZ
double	m_wholeRotatX
double	m_wholeRotatY
double	m_wholeRotatZ
int	m_mindex
std::string	m_alignFilePath
std::string	m_wirePosFilePath
std::string	m_wireTensionFilePath
bool	m_updataalign
vector< MdcGeoWire * >	fWires
vector< MdcGeoGeneral >	fGenerals
vector< MdcGeoLayer * >	fLayers
vector< MdcGeoSuper * >	fSupers
MdcGeoMisc	fMisc
vector< MdcGeoEnd * >	fEnd
IDataProviderSvc *	m_pCalibDataSvc

---

Detailed Description

Definition at line 25 of file MdcGeomSvc.h.

---

Constructor & Destructor Documentation

MdcGeomSvc::MdcGeomSvc	(	const std::string &	name,
		ISvcLocator *	svcloc
	)

Definition at line 30 of file MdcGeomSvc.cxx.

References m_alignFilePath, m_doSag, m_nomcalignment, m_readAlignParDataBase, m_wholeRotatX, m_wholeRotatY, m_wholeRotatZ, m_wholeShiftX, m_wholeShiftY, m_wholeShiftZ, m_wirePosFilePath, m_wireTensionFilePath, and deljobs::string.

                                                                   : Service(name, svcloc) {
    if(getenv("MDCGEOMSVCROOT")){
        m_alignFilePath = std::string(getenv("MDCGEOMSVCROOT"))+std::string("/share/MdcAlignPar.dat");
        //std::cout<<" the MDC alignment file: "<<m_alignFilePath<<std::endl;

        m_wirePosFilePath = std::string(getenv("MDCGEOMSVCROOT"))+std::string("/share/WirePosCalib.dat");
        //std::cout<<" the MDC wire position  file: "<<m_wirePosFilePath<<std::endl;

        m_wireTensionFilePath = std::string(getenv("MDCGEOMSVCROOT"))+std::string("/share/mdcWireTension.dat");
        //std::cout<<" the MDC wire tension file: "<<m_wireTensionFilePath<<std::endl;

    }
    else {
        std::cout<<"A fatal error, contact wangjk..."<<std::endl;
    }

    declareProperty("doSag",      m_doSag = true);
    declareProperty("readAlignParDataBase", m_readAlignParDataBase = true);
    declareProperty("mcnoalignment",m_nomcalignment=true);
    declareProperty("wholeShiftX",m_wholeShiftX = 0.);
    declareProperty("wholeShiftY",m_wholeShiftY = 0.);
    declareProperty("wholeShiftZ",m_wholeShiftZ = 0.);
    declareProperty("wholeRotatX",m_wholeRotatX = 0.);
    declareProperty("wholeRotatY",m_wholeRotatY = 0.);
    declareProperty("wholeRotatZ",m_wholeRotatZ = 0.);
    declareProperty("alignFilePath",m_alignFilePath);
    declareProperty("wirePosFilePath",m_wirePosFilePath);
    declareProperty("wireTensionFilePath",m_wireTensionFilePath);
    

}

MdcGeomSvc::~MdcGeomSvc

(

)

Definition at line 112 of file MdcGeomSvc.cxx.

References fEnd, fGenerals, fLayers, fSupers, and fWires.

                       {
    for(vector<MdcGeoLayer*>::iterator it1 = fLayers.begin(); it1 != fLayers.end(); it1++) delete *it1;
    for(vector<MdcGeoSuper*>::iterator it2 = fSupers.begin(); it2 != fSupers.end(); it2++) delete *it2;
    for(vector<MdcGeoWire*>::iterator it3 = fWires.begin(); it3 != fWires.end(); it3++) delete *it3;
    for(vector<MdcGeoEnd*>::iterator it4 = fEnd.begin(); it4 != fEnd.end(); it4++) delete *it4;  
    fGenerals.clear();
    fWires.clear();
    fLayers.clear();
    fSupers.clear();
    fEnd.clear();  
}

---

Member Function Documentation

void MdcGeomSvc::clean

(

)

[private]

Definition at line 124 of file MdcGeomSvc.cxx.

References fEnd, fGenerals, fLayers, fSupers, and fWires.

Referenced by handle().

                      {
    for(vector<MdcGeoLayer*>::iterator it1 = fLayers.begin(); it1 != fLayers.end(); it1++) delete *it1;
    for(vector<MdcGeoSuper*>::iterator it2 = fSupers.begin(); it2 != fSupers.end(); it2++) delete *it2;
    for(vector<MdcGeoWire*>::iterator it3 = fWires.begin(); it3 != fWires.end(); it3++) delete *it3;
    for(vector<MdcGeoEnd*>::iterator it4 = fEnd.begin(); it4 != fEnd.end(); it4++) delete *it4;
    fGenerals.clear();
    fWires.clear();
    fLayers.clear();
    fSupers.clear();
    fEnd.clear();
}

void MdcGeomSvc::Dump

(

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 694 of file MdcGeomSvc.cxx.

Referenced by main().

{}

const MdcGeoEnd *const MdcGeomSvc::End

(

unsigned

id

)

Definition at line 813 of file MdcGeomSvc.cxx.

References fEnd.

Referenced by BesMdcGeoParameter::InitFromSvc().

                              {
        if (id < fEnd.size())     
            return  fEnd[id];

        return 0;
    }  

void MdcGeomSvc::Fill

(

)

[private]

StatusCode MdcGeomSvc::finalize

(

)

[virtual]

Definition at line 106 of file MdcGeomSvc.cxx.

References Bes_Common::INFO.

Referenced by main().

                                  {
    MsgStream log(messageService(), name());
    log << MSG::INFO << name() << ": End of Run" << endreq;
    return StatusCode::SUCCESS;
}

const MdcGeoGeneral *const MdcGeomSvc::GeneralLayer

(

unsigned

id

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 800 of file MdcGeomSvc.cxx.

References fGenerals.

Referenced by BesMdcGeoParameter::InitFromSvc().

                                       {
        if (id < fGenerals.size()) 
            return & fGenerals[id];

        return 0;
    }

const int MdcGeomSvc::getAlignParIndexEast

(

int

lyr

)

const [private]

Definition at line 696 of file MdcGeomSvc.cxx.

Referenced by ReadFilePar().

                                                       {
    if((lyr>=0) && (lyr<=16))  return 0;  
    else if((lyr>=17) && (lyr<=20)) return 1;
    else if((lyr>=21) && (lyr<=24)) return 2;
    else if((lyr>=25) && (lyr<=28)) return 3;
    else if((lyr>=29) && (lyr<=32)) return 4;
    else if((lyr>=33) && (lyr<=36)) return 5;
    else if((lyr>=37) && (lyr<=41)) return 6;
    else if(lyr>=42) return 7; 
    else std::cout<<" Hi, ERROR OCCUR !!!"<<std::endl;
    return -1;
}

const int MdcGeomSvc::getAlignParIndexWest

(

int

lyr

)

const [private]

Definition at line 718 of file MdcGeomSvc.cxx.

Referenced by ReadFilePar().

                                                       {
    if((lyr>=0) && (lyr<=16))  return 8;
    else if((lyr>=17) && (lyr<=20)) return 9;
    else if((lyr>=21) && (lyr<=24)) return 10;
    else if((lyr>=25) && (lyr<=28)) return 11;
    else if((lyr>=29) && (lyr<=32)) return 12;
    else if((lyr>=33) && (lyr<=36)) return 13;
    else if((lyr>=37) && (lyr<=41)) return 14;
    else if(lyr>=42) return 15;
    else std::cout<<" Hi, ERROR OCCUR !!!"<<std::endl;
    return -1;
}

const int MdcGeomSvc::getGeneralLayerSize

(

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 684 of file MdcGeomSvc.cxx.

References fGenerals.

{
    return fGenerals.size();
}

const int MdcGeomSvc::getLayerSize

(

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 674 of file MdcGeomSvc.cxx.

References fLayers.

Referenced by MdcDetector::MdcDetector(), and KalFitAlg::set_Mdc().

{
    return fLayers.size();
}

bool MdcGeomSvc::getSagFlag

(

void

)

[static]

Definition at line 821 of file MdcGeomSvc.cxx.

References m_doSag.

Referenced by ReadFilePar(), and MdcGeoWire::Sag().

                                {

    return m_doSag;
}  

const int MdcGeomSvc::getSegmentNo

(

)

Definition at line 689 of file MdcGeomSvc.cxx.

References fEnd.

Referenced by BesMdcGeoParameter::InitFromSvc().

{
    return fEnd.size();
}   

const int MdcGeomSvc::getSuperLayerSize

(

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 679 of file MdcGeomSvc.cxx.

References fSupers.

Referenced by MdcDetector::MdcDetector(), and KalFitAlg::set_Mdc().

{
    return fSupers.size();
}

const int MdcGeomSvc::getWireSize

(

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 669 of file MdcGeomSvc.cxx.

References fWires.

Referenced by MdcDetector::MdcDetector(), and KalFitAlg::set_Mdc().

{
    return fWires.size();
}

void MdcGeomSvc::handle

(

const Incident &

inc

)

this handle function is prepared for special use

Definition at line 741 of file MdcGeomSvc.cxx.

References clean(), Bes_Common::DEBUG, false, Bes_Common::FATAL, m_mindex, m_nomcalignment, m_readAlignParDataBase, m_updataalign, and ReadFilePar().

                                          {
    MsgStream log( messageService(), name() );
    log << MSG::DEBUG << "handle: " << inc.type() << endreq;
    IDataProviderSvc* m_eventSvc;
    Gaudi::svcLocator()->service("EventDataSvc", m_eventSvc, true);
    SmartDataPtr<Event::EventHeader> eventHeader(m_eventSvc,"/Event/EventHeader");
    if (!eventHeader) { 
    log << MSG::FATAL << "Could not find Event Header" << endreq;
      }
    if (m_updataalign) return;
    if (inc.type() == "NewRun" ){
        log << MSG::DEBUG << "Begin Event" << endreq;
        clean();
        m_updataalign = true;
    if(m_nomcalignment&&m_mindex==0) {  
        int   RunNo=eventHeader->runNumber();
        if(RunNo<0)  m_readAlignParDataBase=false ;
           else   m_readAlignParDataBase=true;
        m_mindex+=1;
    cout<<"m__RunNo="<<RunNo<<"m_mindex="<<m_mindex<<endl;
        }
   //std::cout<<"############"<<m_readAlignParDataBase<<std::endl;
        ReadFilePar();          
    }
}

StatusCode MdcGeomSvc::initialize

(

)

[virtual]

Definition at line 72 of file MdcGeomSvc.cxx.

References Bes_Common::DEBUG, calibUtil::ERROR, Bes_Common::INFO, m_mindex, m_pCalibDataSvc, m_updataalign, and ReadFilePar().

Referenced by main().

                                    {
    MsgStream log(messageService(), name());
    log << MSG::INFO << name() << ": Start of run initialisation" << endreq;

    StatusCode sc = Service::initialize();
    if ( sc.isFailure() ) return sc;
    m_mindex=0; 
    m_updataalign = false;
    IIncidentSvc* incsvc;
    sc = service("IncidentSvc", incsvc);
    int priority = 100;
    if( sc.isSuccess() ){
        incsvc -> addListener(this, "NewRun", priority);
    }

    //    ReadFilePar(); //  get geometry data from file SimUtil/dat/Mdc.txt
    // Fill();  //  get geometry data from Database

    sc = service("CalibDataSvc", m_pCalibDataSvc, true);

    if ( !sc.isSuccess() ) {
        log << MSG::ERROR 
            << "Could not get IDataProviderSvc interface of CalibXmlCnvSvc" 
            << endreq;
        return sc;
    } else {
        log << MSG::DEBUG 
            << "Retrieved IDataProviderSvc interface of CalibXmlCnvSvc" 
            << endreq;
    }
     ReadFilePar(); 
    return StatusCode::SUCCESS;
}

static const InterfaceID& IMdcGeomSvc::interfaceID

(

)

[inline, static, inherited]

Definition at line 17 of file IMdcGeomSvc.h.

References IID_IMdcGeomSvc().

{ return IID_IMdcGeomSvc; } 

const MdcGeoLayer *const MdcGeomSvc::Layer

(

unsigned

id

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 784 of file MdcGeomSvc.cxx.

References fLayers.

Referenced by BesMdcSD::Distance(), KalFitAlg::execute(), BesMdcGeoParameter::InitFromSvc(), MdcDetector::MdcDetector(), MdcTrackUtil::nLayerTrackPassed(), BesMdcSD::ProcessHits(), KalFitAlg::set_Mdc(), and Wire().

                                {
        if (id < fLayers.size())
            return  fLayers[id];

        return 0;
    }

const MdcGeoMisc *const MdcGeomSvc::Misc

(

void

)

Definition at line 808 of file MdcGeomSvc.cxx.

References fMisc.

Referenced by BesMdcGeoParameter::InitFromSvc().

                    {
    return & fMisc;
}

StatusCode MdcGeomSvc::queryInterface	(	const InterfaceID &	riid,
		void **	ppvUnknown
	)			[virtual]

Definition at line 62 of file MdcGeomSvc.cxx.

References IID_IMdcGeomSvc().

                                                                                   {

    if ( IID_IMdcGeomSvc.versionMatch(riid) ) { 
        *ppvInterface = static_cast<IMdcGeomSvc*> (this); 
    } else { 
        return Service::queryInterface(riid, ppvInterface) ; 
    }
    return StatusCode::SUCCESS;
}

void MdcGeomSvc::ReadAliParDataBase	(	vector< double > &	Sx,
		vector< double > &	Sy,
		vector< double > &	Sz,
		vector< double > &	Rx,
		vector< double > &	Ry,
		vector< double > &	Rz
	)			[private]

Definition at line 638 of file MdcGeomSvc.cxx.

References calibUtil::ERROR, Bes_Common::INFO, m_pCalibDataSvc, m_wholeRotatX, m_wholeRotatY, m_wholeRotatZ, m_wholeShiftX, m_wholeShiftY, m_wholeShiftZ, and deljobs::string.

Referenced by ReadFilePar().

                                                                   {
    MsgStream log(messageService(), name());
    std::string fullPath = "/Calib/MdcAlign";
    log << MSG::INFO<<"ReadAliParDataBase() alignParPath = "<<fullPath<< endreq;
    cout << "Read alignment parameters from Calibration Database!" << endl;

    SmartDataPtr<CalibData::MdcAlignData> alignpar(m_pCalibDataSvc, fullPath);
    if( ! alignpar){
        log << MSG::ERROR << "can not get MdcAlignConst via SmartPtr" 
            << endreq;
        return;
    }
    double tmp1,tmp2,tmp3,tmp4,tmp5,tmp6;
    for(int k=0;k<16;k++){
        tmp1 = alignpar->getdxEP(k);
        tmp2 = alignpar->getdyEP(k);
        tmp3 = alignpar->getdzEP(k);
        tmp4 = alignpar->getrxEP(k);
        tmp5 = alignpar->getryEP(k);
        tmp6 = alignpar->getrzEP(k);

        Sx.push_back(m_wholeShiftX+tmp1);
        Sy.push_back(m_wholeShiftY+tmp2);
        Sz.push_back(m_wholeShiftZ+tmp3);
        Rx.push_back(m_wholeRotatX+tmp4);
        Ry.push_back(m_wholeRotatY+tmp5);
        Rz.push_back(m_wholeRotatZ+tmp6);
    }

}

void MdcGeomSvc::ReadFilePar

(

)

[private]

Definition at line 137 of file MdcGeomSvc.cxx.

References abs, cos(), DBL_MAX, deljobs::end, fEnd, fGenerals, MdcGeoGeneral::First(), dchain::children::First, fLayers, fMisc, fSupers, MdcGeoMisc::FWireR(), fWires, getAlignParIndexEast(), getAlignParIndexWest(), getSagFlag(), genRecEmupikp::i, MdcGeoSuper::Id(), MdcGeoGeneral::Id(), MdcGeoMisc::InnerR(), MdcGeoMisc::InnerTk(), ganga-rec::j, MdcGeoGeneral::LayerName(), MdcGeoMisc::LayerNo(), MdcGeoGeneral::Length(), genRecEmupikp::line, MdcGeoSuper::LowerR(), m_alignFilePath, m_readAlignParDataBase, m_updataalign, m_wholeRotatX, m_wholeRotatY, m_wholeRotatZ, m_wholeShiftX, m_wholeShiftY, m_wholeShiftZ, m_wirePosFilePath, m_wireTensionFilePath, MdcGeoGeneral::NCell(), MdcGeoMisc::NFWire(), MdcGeoGeneral::nomOffset(), MdcGeoGeneral::nomPhi(), MdcGeoGeneral::nomShift(), MdcGeoMisc::NSWire(), MdcGeoGeneral::OffB(), MdcGeoGeneral::OffF(), MdcGeoGeneral::Offset(), MdcGeoMisc::OuterR(), MdcGeoMisc::OuterTk(), MdcGeoGeneral::Phi(), Phi(), phi0, pi, rb::R(), MdcGeoGeneral::Radius(), ReadAliParDataBase(), ReadTensionDataBase(), ReadWirePosDataBase(), MdcGeoGeneral::RxEast(), MdcGeoGeneral::RxWest(), MdcGeoGeneral::RyEast(), MdcGeoGeneral::RyWest(), MdcGeoGeneral::RzEast(), MdcGeoGeneral::RzWest(), MdcGeoGeneral::Shift(), sin(), MdcGeoMisc::SLayerNo(), deljobs::string, MdcGeoMisc::SWireR(), MdcGeoGeneral::SxEast(), MdcGeoGeneral::SxWest(), MdcGeoGeneral::SyEast(), MdcGeoGeneral::SyWest(), MdcGeoGeneral::SzEast(), MdcGeoGeneral::SzWest(), title, MdcGeoGeneral::TwistB(), MdcGeoGeneral::TwistF(), MdcGeoSuper::Type(), MdcGeoSuper::UpperR(), MdcGeoMisc::WireNo(), and wireNo.

Referenced by handle(), and initialize().

                            {
    std::string geometryFilePath = getenv("MDCSIMROOT");
    geometryFilePath += "/dat/Mdc.txt";
    std::ifstream inFile(geometryFilePath.c_str() );

    if(!inFile.good()){
        std::cout<<"Error, mdc parameters file not exist"<<std::endl; 
        return;
    }
    std::string alignFilePath;
    std::string wirePosFilePath;
    std::string wireTensionFilePath;
    if (!m_updataalign) {
        alignFilePath = std::string(getenv("MDCGEOMSVCROOT"))+std::string("/share/MdcAlignPar.dat");
        wirePosFilePath = std::string(getenv("MDCGEOMSVCROOT"))+std::string("/share/WirePosCalib.dat");
        wireTensionFilePath = std::string(getenv("MDCGEOMSVCROOT"))+std::string("/share/mdcWireTension.dat");
    } else {
        alignFilePath = m_alignFilePath;
        wirePosFilePath = m_wirePosFilePath;
        wireTensionFilePath = m_wireTensionFilePath;
    }
    std::ifstream alignFile(alignFilePath.c_str()); 
    std::ifstream wirePosFile(wirePosFilePath.c_str());
    std::ifstream wireTensionFile(wireTensionFilePath.c_str());

    if(!wirePosFile){
        std::cout<<"Error, mdc wire position file not exist..."<<std::endl;
    }

    if(!wireTensionFile){
        std::cout<<"Error, mdc wire tension file not exist..."<<std::endl;
    }

      std::vector<double> wireTensionVec;
    //cout << __FILE__ << "    " << __LINE__ << "  readdb = " << m_readAlignParDataBase
    //    << "   update align = " << m_updataalign << endl;
    if (m_readAlignParDataBase && m_updataalign) {
        ReadTensionDataBase(wireTensionVec);
    }else {
        cout << "Read wireTension from file:" << wireTensionFilePath.c_str() << endl;
        int idd;
        double tension;
        for(int ii=0; ii<6796; ii++){
            wireTensionFile>>idd>>tension;

            if(getSagFlag()){
                wireTensionVec.push_back(tension);
            }
            else { 
                tension = DBL_MAX;
                wireTensionVec.push_back(tension);
            }
        }
    }

    std::vector<vector<double> > wirePosVec(6796);
    if (m_readAlignParDataBase && m_updataalign) {
        ReadWirePosDataBase(wirePosVec);
    } else {
        cout << "Read wirePosition from file:" << wirePosFilePath.c_str() << endl;
        double dxe,dye,dze,dxw,dyw,dzw;
        int wid;
        std::string title; 
        getline(wirePosFile, title);
        wirePosFile.seekg(1,ios::cur);
        for(int j=0; j<6796; j++){
            wirePosFile>>wid>>dxe>>dye>>dze>>dxw>>dyw>>dzw; 
            wirePosVec[j].push_back(dxe);
            wirePosVec[j].push_back(dye);
            wirePosVec[j].push_back(dze);
            wirePosVec[j].push_back(dxw);
            wirePosVec[j].push_back(dyw);
            wirePosVec[j].push_back(dzw);
            /*
               for(int jj=0;jj<6;jj++){
               std::cout<<" wid: "<<wid<<" dxe: "<<wirePosVec[j][0]<<" dye: "<<wirePosVec[j][1]<<" dze: "<<wirePosVec[j][2]<<" dxw: "<<wirePosVec[j][3]<<" dyw: "<<wirePosVec[j][4]<<" dzw: "<<wirePosVec[j][5]<<std::endl;
               }
               */
        }
    }

    double signalWireR, fieldWireR;
    int TLayerNo, TWireNo,TSignalLayerNo,fSignalLayer[50];  
    int i,wireNo, firstWire,segmentNo,signalLayer;
    double length, nomphi, phi, r,nomaadiv2, aadiv2,
           nomaa, aa, nomrotateCell, rotateCell, wlength, innerR, outR, z;
    std::string layer, name, line;

    vector<int> WireNo, fSegmentNo;
    vector<double> wLength, R, nomPhi, Phi,nomadiv2, adiv2, First, 
        noma, a, nomebusing, ebusing, nomsigma, sigma, nomdeltaz, deltaz,
        nomRotate, Rotate, fLength, fInnerR, fOutR, fZ;
    vector<string> LayerName, fName;  
    vector<double> Sx,Sy,Sz,Rx,Ry,Rz;
    double tmp1,tmp2,tmp3,tmp4,tmp5,tmp6;


    if (m_readAlignParDataBase && m_updataalign) {
        ReadAliParDataBase(Sx, Sy, Sz, Rx, Ry, Rz);
    } else {
        cout << "Read alignment parameters from file:" << alignFilePath.c_str() << endl;
        getline(alignFile, line);
        alignFile.seekg(1,ios::cur);
        for(int k=0;k<16;k++){
            alignFile>>line>>tmp1>>tmp2>>tmp3>>tmp4>>tmp5>>tmp6;   
            //changed by luot 10.03.12
            //if (!m_updataalign) {
            //    tmp1 = 0.0;
            //    tmp2 = 0.0;
            //    tmp3 = 0.0;
            //    tmp4 = 0.0;
            //    tmp5 = 0.0;
            //    tmp6 = 0.0;
            //}
            Sx.push_back(m_wholeShiftX+tmp1);
            Sy.push_back(m_wholeShiftY+tmp2);
            Sz.push_back(m_wholeShiftZ+tmp3);
            Rx.push_back(m_wholeRotatX+tmp4);
            Ry.push_back(m_wholeRotatY+tmp5);
            Rz.push_back(m_wholeRotatZ+tmp6);
        }
    }

    getline(inFile, line);
    inFile>>TLayerNo>>TWireNo>>TSignalLayerNo>>signalWireR>>fieldWireR;
    inFile.seekg(1,ios::cur);
    getline(inFile, line);
    for( i=0; i<TSignalLayerNo; i++){
        inFile>>signalLayer;
        fSignalLayer[i]=signalLayer-1;
    }

    inFile.seekg(1,ios::cur);
    getline(inFile, line);
    getline(inFile, line);

    int i_east,i_west;   
    double delta_rotateCell;
    for( i=0; i<TLayerNo; i++){
        i_east=getAlignParIndexEast(i);
        i_west=getAlignParIndexWest(i); 

        inFile>>layer>>wireNo>>wlength>>r>>nomphi>>firstWire>>nomrotateCell;
        getline(inFile, line);

        nomaadiv2=2*pi*nomrotateCell/wireNo; 
        //aadiv2=2*pi*rotateCell/wireNo+0.5*(Rz[i_west]-Rz[i_east]);
        nomaa=2*r*sin(nomaadiv2);
        delta_rotateCell=(Rz[i_west]-Rz[i_east])*wireNo/(4*pi);
        rotateCell=nomrotateCell+delta_rotateCell;
        aadiv2=2*pi*rotateCell/wireNo;
        aa=2*r*sin(aadiv2);

        double shiftPhi=twopi/wireNo;
        nomphi*=deg;
        if(nomphi<0) nomphi += shiftPhi; //By definition, phi of first wire must >=0
        phi=nomphi+Rz[i_east];
        LayerName.push_back(layer);
        WireNo.push_back(wireNo);
        wLength.push_back(wlength);
        R.push_back(r);
        nomPhi.push_back(nomphi);
        Phi.push_back(phi);
        //nomadiv2.push_back(nomaadiv2);
        //adiv2.push_back(aadiv2);
        First.push_back(firstWire);
        //noma.push_back(nomaa);
        a.push_back(aa);

        nomebusing.push_back(atan(nomaa/wlength));
        ebusing.push_back(atan(aa/wlength));

        //if(aa==0.0) sigma.push_back(0.0);
        //else  sigma.push_back(0.13*wlength/aa);
        //nomdeltaz.push_back(r*(1.0-cos(nomaadiv2)));
        //deltaz.push_back(r*(1.0-cos(aadiv2)));
        nomRotate.push_back(nomrotateCell);
        Rotate.push_back(rotateCell);
    }

    getline(inFile, line);
    inFile>>segmentNo;  
    inFile.seekg(1,ios::cur);
    getline(inFile, line);
    getline(inFile, line);

    for(i=0; i<segmentNo; i++){
        inFile>>length>>innerR>>outR>>z>>name;
        getline(inFile,line);
        fLength.push_back(length);
        fInnerR.push_back(innerR);
        fOutR.push_back(outR);
        fZ.push_back(z);
        fName.push_back(name);
    }


    //MdcGeoGeneral include signal wire and field wire
    MdcGeoGeneral hlh;
    double a1,a2,a3,nomphi0,phi0,cellphi;
    double noma1,noma2,noma3;
    for(i=0;i<TLayerNo;i++){
        i_east=getAlignParIndexEast(i);
        i_west=getAlignParIndexWest(i);
        hlh.Id(i);
        hlh.LayerName(LayerName[i]);    
        hlh.Radius(R[i]);
        hlh.Length(wLength[i]);
        hlh.NCell(WireNo[i]/2);
        //east endcap rotates away from designed position by Phi considering misalignment,
        //nomPhi without considering misalignment 
        hlh.Offset(Phi[i]);
        hlh.nomOffset(nomPhi[i]);
        hlh.Phi(Phi[i]);    
        hlh.nomPhi(nomPhi[i]);
        hlh.First((int)First[i]); 
        HepPoint3D offF(Sx[i_west],Sy[i_west],Sz[i_west]); 
        HepPoint3D offB(Sx[i_east],Sy[i_east],Sz[i_east]);
        hlh.OffF(offF);
        hlh.OffB(offB);
        hlh.TwistF(Rz[i_west]);
        hlh.TwistB(Rz[i_east]);
        //west endcap rotates away from designed position by Rotate considering misalignment,
        //nomRotate without considering misalignment 
        hlh.Shift(Rotate[i]);
        hlh.nomShift(nomRotate[i]);
        hlh.SxEast(Sx[i_east]);
        hlh.SxWest(Sx[i_west]);
        hlh.SyEast(Sy[i_east]);
        hlh.SyWest(Sy[i_west]);
        hlh.SzEast(Sz[i_east]);
        hlh.SzWest(Sz[i_west]);

        hlh.RxEast(Rx[i_east]);
        hlh.RxWest(Rx[i_west]);
        hlh.RyEast(Ry[i_east]);
        hlh.RyWest(Ry[i_west]);
        hlh.RzEast(Rz[i_east]);
        hlh.RzWest(Rz[i_west]);
        fGenerals.push_back(hlh);
    }

    MdcGeoSuper *suph;
    for(i=0;i<TSignalLayerNo;i++){
        //MdcGeoLayer only represent signal wires
        MdcGeoLayer *lh = new MdcGeoLayer();
        // lh is a signalwire layer 
        int lyr=fSignalLayer[i];
        // the value of lyr starts from 1,because the first layer of MDC is field
        // wire
        i_east=getAlignParIndexEast(lyr);
        i_west=getAlignParIndexWest(lyr);
        lh->SLayer(lyr+1);  
        if(i%4==0){
            suph = new MdcGeoSuper();
            suph->LowerR(R[lyr-1]);
            if(i==40){
                suph->UpperR(R[lyr+5]);
            }
            else {
                suph->UpperR(R[lyr+7]);
            }
            switch(i/4){
                case 0: suph->Type(-1); break;
                case 1: suph->Type( 1); break;
                case 2: suph->Type( 0); break;
                case 3: suph->Type( 0); break;
                case 4: suph->Type( 0); break;
                case 5: suph->Type(-1); break;
                case 6: suph->Type( 1); break;
                case 7: suph->Type(-1); break;
                case 8: suph->Type( 1); break;
                case 9: suph->Type( 0); break;
                case 10: suph->Type(0); break;
                default: break;
            }
            suph->Id(i/4);
            fSupers.push_back(suph);
        }
        cellphi=2*twopi/WireNo[lyr];
        phi0=Phi[lyr]+abs(First[lyr]-1)*(cellphi/2);
        nomphi0=nomPhi[lyr]+abs(First[lyr]-1)*(cellphi/2);
        lh->Id(i);
        int wircount=0;
        for(int wk=0; wk<i; wk++){
            int lyrwk=fSignalLayer[wk];
            wircount=wircount + WireNo[lyrwk]/2;
        };
        lh->Wirst(wircount);
        lh->Slant(ebusing[lyr]);
        lh->nomSlant(nomebusing[lyr]);
        lh->Radius(R[lyr]);
        lh->Length(wLength[lyr]);
        lh->RCSiz1(R[lyr]-R[lyr-1]);
        lh->RCSiz2(R[lyr+1]-R[lyr]);
        lh->PCSiz(R[lyr]*cellphi);
        lh->NCell(WireNo[lyr]/2);
        lh->Offset(phi0);
        lh->Shift(Rotate[lyr]);
        lh->nomOffset(nomphi0);
        lh->nomShift(nomRotate[lyr]);

        vector<MdcGeoSuper*>::iterator it1 = fSupers.end()-1;
        lh->Sup(* it1);
        fLayers.push_back(lh);
        int geo = -1;
        if (i < 8) {
            geo = i * 2 + 1;
        } else if (i < 20) {
            geo = i * 2 + 2;
        } else if (i < 36) {
            geo = i * 2 + 3;
        } else {
            geo = i * 2 + 4;
        } 
        lh->Gen(geo);

        for(int j=0;j<WireNo[lyr]/2;j++){
            MdcGeoWire *wh = new MdcGeoWire();
            const int wireId = wircount+j;
            wh->Id(wireId);

            a1 = R[lyr]*cos(phi0+j*cellphi)+Sx[i_east]+wirePosVec[wireId][0];
            a2 = R[lyr]*sin(phi0+j*cellphi)+Sy[i_east]+wirePosVec[wireId][1];
            a3 = wLength[lyr]/2+wirePosVec[wireId][2]+Sz[i_east];
            //a3 = wLength[lyr]/2+wirePosVec[wireId][2];
            //MDC rotate around x firstly;
            double xb = a1;
            double yb = a2 * cos(Rx[i_east]) + a3 * sin(Rx[i_east]);
            double zb = a3 * cos(Rx[i_east]) - a2 * sin(Rx[i_east]);

            //then MDC rotate around y;
            double xbnew = xb * cos(Ry[i_east]) - zb * sin(Ry[i_east]);
            double ybnew = yb;
            double zbnew = xb * sin(Ry[i_east]) + zb * cos(Ry[i_east]);
            noma1 = R[lyr]*cos(nomphi0+j*cellphi);
            noma2 = R[lyr]*sin(nomphi0+j*cellphi);
            noma3 = wLength[lyr]/2;
            //HepPoint3D pb(a1,a2,a3);
            HepPoint3D pb(xbnew,ybnew,zbnew);
            HepPoint3D nompb(noma1,noma2,noma3);
            HepPoint3D wireposb(wirePosVec[wireId][0],wirePosVec[wireId][1],wirePosVec[wireId][2]);
            wh->Backward(pb);
            wh->nomBackward(nompb);
            wh->BWirePos(wireposb);

            //double g1 = noma1*cos(Rz[i_east]) - noma2*sin(Rz[i_east]) + Sx[i_east] + wirePosVec[wireId][0];
            //double g2 = noma2*cos(Rz[i_east]) + noma1*sin(Rz[i_east]) + Sy[i_east] + wirePosVec[wireId][1];
            //double g3 = noma3 + wirePosVec[wireId][2];
            //HepPoint3D backward_test(g1,g2,g3);

            a1 = R[lyr]*cos(phi0+(j+lh->Shift())*cellphi)+Sx[i_west]+wirePosVec[wireId][3];
            a2 = R[lyr]*sin(phi0+(j+lh->Shift())*cellphi)+Sy[i_west]+wirePosVec[wireId][4];
            //a3 = -wLength[lyr]/2+wirePosVec[wireId][5];
            a3 = -wLength[lyr]/2+wirePosVec[wireId][5]+Sz[i_west];
            //MDC rotate around x firstly;
            double xf = a1;
            double yf = a2 * cos(Rx[i_west]) + a3 * sin(Rx[i_west]);
            double zf = a3 * cos(Rx[i_west]) - a2 * sin(Rx[i_west]);

            //then MDC rotate around y;
            double xfnew = xf * cos(Ry[i_west]) - zf * sin(Ry[i_west]);
            double yfnew = yf;
            double zfnew = xf * sin(Ry[i_west]) + zf * cos(Ry[i_west]);

            noma1 = R[lyr]*cos(nomphi0+(j+lh->nomShift())*cellphi);
            noma2 = R[lyr]*sin(nomphi0+(j+lh->nomShift())*cellphi);
            noma3 = -wLength[lyr]/2;
            //HepPoint3D pf(a1,a2,a3);
            HepPoint3D pf(xfnew, yfnew, zfnew);
            HepPoint3D nompf(noma1,noma2,noma3);
            HepPoint3D wireposf(wirePosVec[wireId][3],wirePosVec[wireId][4],wirePosVec[wireId][5]);
            wh->Forward(pf);
            wh->nomForward(nompf);
            wh->FWirePos(wireposf);

            //double f1 = noma1*cos(Rz[i_west]) - noma2*sin(Rz[i_west]) + Sx[i_west] + wirePosVec[wireId][3];
            //double f2 = noma2*cos(Rz[i_west]) + noma1*sin(Rz[i_west]) + Sy[i_west] + wirePosVec[wireId][4];
            //double f3 = noma3 + wirePosVec[wireId][5];
            //HepPoint3D forward_test(f1,f2,f3);
            wh->Layer(i);
            wh->Cell(j);
            wh->Stat(0);
            wh->Slant(ebusing[lyr]);
            wh->nomSlant(nomebusing[lyr]);

            wh->Tension(wireTensionVec[wireId]);
            //std::cout<<"  wh->Tension: "<< wh->Tension()<<std::endl;
            //std::cout<<" wire sag: "<<wh->Sag(wireId)<<std::endl;

            vector<MdcGeoLayer*>::iterator it2 = fLayers.end()-1;
            wh->Lyr(* it2);
            fWires.push_back(wh);
        }
    }

    fMisc.OuterR(fOutR[1]);
    fMisc.InnerR(fInnerR[2]);
    fMisc.OuterTk(fOutR[1]-fInnerR[1]);
    fMisc.InnerTk(fOutR[2]-fInnerR[2]);
    fMisc.NSWire(fWires.size());
    fMisc.NFWire(TWireNo-fWires.size());

    fMisc.LayerNo(TLayerNo);
    fMisc.WireNo(TWireNo);
    fMisc.SLayerNo(TSignalLayerNo);
    fMisc.SWireR(signalWireR);
    fMisc.FWireR(fieldWireR);

    for(i=0;i<segmentNo;i++){
        MdcGeoEnd * end=new MdcGeoEnd();
        end->Id(i);
        end->Length(fLength[i]);
        end->InnerR(fInnerR[i]);
        end->OutR(fOutR[i]);
        end->Z(fZ[i]);
        end->Name(fName[i]);
        fEnd.push_back(end);
    }
}

void MdcGeomSvc::ReadTensionDataBase

(

std::vector< double > &

wireTensionVec

)

[private]

Definition at line 590 of file MdcGeomSvc.cxx.

References calibUtil::ERROR, genRecEmupikp::i, Bes_Common::INFO, m_pCalibDataSvc, and deljobs::string.

Referenced by ReadFilePar().

                                                                      {
    std::string fullPath = "/Calib/MdcAlign";
    MsgStream log(messageService(), name());
    log << MSG::INFO<<"ReadTensionDataBase() wireTensionPath = "<<fullPath<< endreq;
    cout << "Read wireTension from Calibration Database!" << endl;

    SmartDataPtr<CalibData::MdcAlignData> tension(m_pCalibDataSvc, fullPath);
    if( ! tension){
        log << MSG::ERROR << "can not get MdcAlignConst via SmartPtr" 
            << endreq;
        return;
    }
    for(int i=0;i<6796;i++){
        double tens = tension->gettension(i);
        wireTensionVec.push_back(tens);
    }
}

void MdcGeomSvc::ReadWirePosDataBase

(

std::vector< vector< double > > &

wirePosVec

)

[private]

Definition at line 607 of file MdcGeomSvc.cxx.

References calibUtil::ERROR, Bes_Common::INFO, ganga-rec::j, m_pCalibDataSvc, and deljobs::string.

Referenced by ReadFilePar().

                                                                           {
    std::string fullPath = "/Calib/MdcAlign";
    MsgStream log(messageService(), name());
    log << MSG::INFO<<"ReadWirePosDataBase() wirePositionPath = "<<fullPath<< endreq;

    cout << "Read wirePosition from Calibration Database!" << endl;
    SmartDataPtr<CalibData::MdcAlignData> wirepos(m_pCalibDataSvc, fullPath);
    if( ! wirepos){
        log << MSG::ERROR << "can not get MdcAlignConst via SmartPtr" 
            << endreq;
        return;
    }
    double dxe,dye,dze,dxw,dyw,dzw;
    for(int j=0; j<6796; j++){

        dxe = wirepos->getdxWireEast(j);
        dye = wirepos->getdyWireEast(j);
        dze = wirepos->getdzWireEast(j);
        dxw = wirepos->getdxWireWest(j);
        dyw = wirepos->getdyWireWest(j);
        dzw = wirepos->getdzWireWest(j);

        wirePosVec[j].push_back(dxe);
        wirePosVec[j].push_back(dye);
        wirePosVec[j].push_back(dze);
        wirePosVec[j].push_back(dxw);
        wirePosVec[j].push_back(dyw);
        wirePosVec[j].push_back(dzw);
    }

}

const MdcGeoSuper *const MdcGeomSvc::SuperLayer

(

unsigned

id

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 792 of file MdcGeomSvc.cxx.

References fSupers.

                                     {
        if (id < fSupers.size())
            return  fSupers[id];

        return 0;
    }

const MdcGeoWire *const MdcGeomSvc::Wire	(	unsigned	lyrid,
		unsigned	wirid
	)			[virtual]

Implements IMdcGeomSvc.

Definition at line 776 of file MdcGeomSvc.cxx.

References fLayers, fWires, Layer(), and MdcGeoLayer::Wirst().

                                                  {
        if ((lyrid <fLayers.size()) && ((int) wirid < Layer(lyrid)->NCell()))
            return fWires[Layer(lyrid)->Wirst() + wirid];

        return 0;
    }

const MdcGeoWire *const MdcGeomSvc::Wire

(

unsigned

id

)

[virtual]

Implements IMdcGeomSvc.

Definition at line 768 of file MdcGeomSvc.cxx.

References fWires.

Referenced by BesMdcSD::Distance(), TrkReco::execute(), MdcHoughFinder::execute(), HoughValidUpdate::execute(), KalFitAlg::execute(), KalFitAlg::filter_fwd_anal(), KalFitAlg::filter_fwd_calib(), main(), MdcDetector::MdcDetector(), KalFitAlg::set_Mdc(), KalFitAlg::smoother_anal(), and KalFitAlg::smoother_calib().

                               {
        if (id < fWires.size())
            return  fWires[id];

        return 0;
    }

---

Member Data Documentation

vector<MdcGeoEnd*> MdcGeomSvc::fEnd [private]

Definition at line 91 of file MdcGeomSvc.h.

Referenced by clean(), End(), getSegmentNo(), ReadFilePar(), and ~MdcGeomSvc().

vector<MdcGeoGeneral> MdcGeomSvc::fGenerals [private]

Definition at line 87 of file MdcGeomSvc.h.

Referenced by clean(), GeneralLayer(), getGeneralLayerSize(), ReadFilePar(), and ~MdcGeomSvc().

vector<MdcGeoLayer*> MdcGeomSvc::fLayers [private]

Definition at line 88 of file MdcGeomSvc.h.

Referenced by clean(), getLayerSize(), Layer(), ReadFilePar(), Wire(), and ~MdcGeomSvc().

MdcGeoMisc MdcGeomSvc::fMisc [private]

Definition at line 90 of file MdcGeomSvc.h.

Referenced by Misc(), and ReadFilePar().

vector<MdcGeoSuper*> MdcGeomSvc::fSupers [private]

Definition at line 89 of file MdcGeomSvc.h.

Referenced by clean(), getSuperLayerSize(), ReadFilePar(), SuperLayer(), and ~MdcGeomSvc().

vector<MdcGeoWire*> MdcGeomSvc::fWires [private]

Definition at line 86 of file MdcGeomSvc.h.

Referenced by clean(), getWireSize(), ReadFilePar(), Wire(), and ~MdcGeomSvc().

std::string MdcGeomSvc::m_alignFilePath [private]

Definition at line 81 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), and ReadFilePar().

bool MdcGeomSvc::m_doSag = true [static]

Definition at line 54 of file MdcGeomSvc.h.

Referenced by getSagFlag(), and MdcGeomSvc().

int MdcGeomSvc::m_mindex [private]

Definition at line 80 of file MdcGeomSvc.h.

Referenced by handle(), and initialize().

bool MdcGeomSvc::m_nomcalignment = true [static]

Definition at line 56 of file MdcGeomSvc.h.

Referenced by handle(), and MdcGeomSvc().

IDataProviderSvc* MdcGeomSvc::m_pCalibDataSvc [private]

Definition at line 92 of file MdcGeomSvc.h.

Referenced by initialize(), ReadAliParDataBase(), ReadTensionDataBase(), and ReadWirePosDataBase().

bool MdcGeomSvc::m_readAlignParDataBase = true [static]

Definition at line 55 of file MdcGeomSvc.h.

Referenced by handle(), MdcGeomSvc(), and ReadFilePar().

bool MdcGeomSvc::m_updataalign [private]

Definition at line 84 of file MdcGeomSvc.h.

Referenced by handle(), initialize(), and ReadFilePar().

double MdcGeomSvc::m_wholeRotatX [private]

Definition at line 77 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), ReadAliParDataBase(), and ReadFilePar().

double MdcGeomSvc::m_wholeRotatY [private]

Definition at line 78 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), ReadAliParDataBase(), and ReadFilePar().

double MdcGeomSvc::m_wholeRotatZ [private]

Definition at line 79 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), ReadAliParDataBase(), and ReadFilePar().

double MdcGeomSvc::m_wholeShiftX [private]

Definition at line 74 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), ReadAliParDataBase(), and ReadFilePar().

double MdcGeomSvc::m_wholeShiftY [private]

Definition at line 75 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), ReadAliParDataBase(), and ReadFilePar().

double MdcGeomSvc::m_wholeShiftZ [private]

Definition at line 76 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), ReadAliParDataBase(), and ReadFilePar().

std::string MdcGeomSvc::m_wirePosFilePath [private]

Definition at line 82 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), and ReadFilePar().

std::string MdcGeomSvc::m_wireTensionFilePath [private]

Definition at line 83 of file MdcGeomSvc.h.

Referenced by MdcGeomSvc(), and ReadFilePar().

---