BesMdcGeoParameter Class Reference

#include <BesMdcGeoParameter.hh>

Public Member Functions
	BesMdcGeoParameter (void)
	~BesMdcGeoParameter (void)
void	InitFromFile (void)
void	InitFromSvc (void)
void	Dump (void)
double	SignalWireR ()
double	FieldWireR ()
int	Signal2Global (int x)
BesMdcWire	Wire (int)
BesMdcWire	SignalWire (int, int)
int	LayerNo (void)
int	SignalLayerNo (void)
const BesMdcLayer &	Layer (int) const
const BesMdcLayer &	SignalLayer (int) const
int	SegmentNo (void)
BesMdcMember	Segment (int x)

Static Public Member Functions
static BesMdcGeoParameter *	GetGeo (void)

Detailed Description

Definition at line 113 of file BesMdcGeoParameter.hh.

Constructor & Destructor Documentation

BesMdcGeoParameter()

BesMdcGeoParameter::BesMdcGeoParameter

(

void

)

Definition at line 78 of file BesMdcGeoParameter.cc.

                                          {
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  IG4Svc* tmpSvc;
  G4Svc* m_G4Svc;
  StatusCode sc=svcLocator->service("G4Svc", tmpSvc);
  m_G4Svc=dynamic_cast<G4Svc *>(tmpSvc); 
  if (!sc.isSuccess())  
   std::cout<<"BesMdcGeoParameter::Could not open G4 Service"<<std::endl;
  if(m_G4Svc->GetMdcDataInput()== 0){ 
    cout<<"------- get BesMdcGeoParameter from file --------"<<endl;
    InitFromFile();
  }
  if(m_G4Svc->GetMdcDataInput()== 1) {
    cout<<"=======get BesMdcGeoParameter from MdcGeomSvc======="<<endl;
    InitFromSvc(); 
  }
  
  //Dump();
  if(fPointer)
  { G4Exception("BesMdcGeoParameter constructed twice."); }
  fPointer=this;
}

Referenced by GetGeo().

~BesMdcGeoParameter()

BesMdcGeoParameter::~BesMdcGeoParameter ( void )

inline

Definition at line 116 of file BesMdcGeoParameter.hh.

{};

Member Function Documentation

Dump()

void BesMdcGeoParameter::Dump

(

void

)

Definition at line 309 of file BesMdcGeoParameter.cc.

                              {
   //cout<<"------BesMdcGeoParameter info :--------"<<endl;
    cout<<" fLayerNo: "<<fLayerNo<<endl;
    cout<<" fWireNo: "<<fWireNo<<endl;
    cout<<" fSignalLayerNo: "<<fSignalLayerNo<<endl;
    cout<<" fSignalWireR: "<<fSignalWireR<<endl;
    cout<<" fFieldWireR: "<<fFieldWireR<<endl;
    
    cout<<"fSingalLayer:"<<endl;
    for(int i=0; i<fSignalLayerNo; i++){
      cout<<fSignalLayer[i]+1<<' '; }
    cout<<endl;
   
    for(int i=0;i<fLayerNo;i++){
      cout<<"Layer["<<i<<"]: "
          <<" name:"<<fLayer[i].Name() <<"  wireNo:"<<fLayer[i].WireNo()
          <<"  length: "<<fLayer[i].Length() <<"  r: "<<fLayer[i].R();
      if (i<75)  cout<<"  phi:"<<fLayer[i].Phi()*180/pi;
      else cout<<"  phi:"<<(fLayer[i].Phi()-fLayer[i].ShiftPhi())*180/pi;
      cout<<"  firstWire: "<<fLayer[i].FirstWire()
          <<"  rotateCell: "<<fLayer[i].RotateCell()<<endl;
    }    
 
    cout<<"fSegmentNo:"<<fSegmentNo<<endl;
    for(int j=0;j<fSegmentNo;j++){
      cout<<"length:"<<fMdcSegment[j].Length()
          <<"  innerR:"<<fMdcSegment[j].InnerR()
          <<"  outR:"<<fMdcSegment[j].OutR()
          <<"  z:"<<fMdcSegment[j].Z()
          <<"  name:"<<fMdcSegment[j].Name()<<endl;
   }
 
}

FieldWireR()

double BesMdcGeoParameter::FieldWireR ( )

inline

Definition at line 124 of file BesMdcGeoParameter.hh.

{return fFieldWireR;}

Referenced by BesMdcConstruction::Construct().

GetGeo()

BesMdcGeoParameter * BesMdcGeoParameter::GetGeo ( void )

static

Definition at line 34 of file BesMdcGeoParameter.cc.

                                                   {
      if (! fPointer) fPointer = new BesMdcGeoParameter();
     return fPointer;
}

Referenced by BesMdcConstruction::BesMdcConstruction(), BesMdcDigitizer::BesMdcDigitizer(), BesMdcSD::BesMdcSD(), and BesRawDataWriter::BesRawDataWriter().

InitFromFile()

void BesMdcGeoParameter::InitFromFile

(

void

)

Definition at line 148 of file BesMdcGeoParameter.cc.

                                         {
  int wireNo, firstWire;
  double length, phi, r, rotateCell,rotateAngle;
  double innerR, outR, z;
  string name, line;
 
 G4String geoPath = ReadBoostRoot::GetBoostRoot(); 
  if(!geoPath){ 
    G4Exception("BOOST environment not set!");
  }
  if(ReadBoostRoot::GetCgem()!=0){
       geoPath += "/dat/Mdc_cgem.txt";
       G4cout << "use CGEM" << endl;
  } else{
       geoPath += "/dat/Mdc.txt";
  }
  G4cout << "Read Mdc geometry file: " << geoPath << endl;
  
  ifstream inFile(geoPath);
  if(!inFile.good()){
    cout<<"Error, mdc parameters file not exist"<<endl; 
    return;
  }
 
  getline(inFile, line);
  inFile>>fLayerNo>>fWireNo>>fSignalLayerNo>>fSignalWireR>>fFieldWireR;
 
  inFile.seekg(1,ios::cur);
  getline(inFile, line);
  int i,signalLayer;
  for(i=0; i<fSignalLayerNo; i++){
    inFile>>signalLayer;
    fSignalLayer[i]=signalLayer-1;
 }
 
  inFile.seekg(1,ios::cur);
  getline(inFile, line);
  getline(inFile, line);
  for( i=0; i<fLayerNo; i++){
    inFile>>name>>wireNo>>length>>r>>phi>>firstWire>>rotateCell;
    getline(inFile, line);
    
    rotateAngle=2*pi*rotateCell/wireNo;
     
    fLayer[i].SetName(name);fLayer[i].SetRadius(r);    
    fLayer[i].SetLength(length); fLayer[i].SetRotateCell(rotateCell);
    fLayer[i].SetRotateAngle(rotateAngle); fLayer[i].SetWireNo(wireNo); 
    fLayer[i].SetShiftPhi(twopi/wireNo); fLayer[i].SetFirstWire(firstWire);    
    
    phi*=(pi/180); 
    if(phi<0)phi += fLayer[i].ShiftPhi();
    fLayer[i].SetPhi(phi); 
    
    if(i==0){
      fLayer[i].SetSumWireNo(wireNo); fLayer[i].SetBeginWireNo(0);
    }else{
      fLayer[i].SetBeginWireNo(fLayer[i-1].SumWireNo());
      fLayer[i].SetSumWireNo(fLayer[i-1].SumWireNo()+wireNo); 
    }
    
    for(int j=0; j<wireNo; j++){
      fWirePhi[fLayer[i].BeginWireNo()+j]=j*fLayer[i].ShiftPhi()+phi; 
    }
  }  
 
  if(fLayer[fLayerNo-1].SumWireNo()!= fWireNo){
    cout<<"Total wire number is not consistant!"<<endl;
  }    
  
  getline(inFile, line);
  inFile>>fSegmentNo;  
  inFile.seekg(1,ios::cur);
  getline(inFile, line);
  getline(inFile, line);
 
  for(i=0; i<fSegmentNo; i++){
    inFile>>length>>innerR>>outR>>z>>name;
    getline(inFile,line);
 
    fMdcSegment[i].SetLength(length); fMdcSegment[i].SetInnerR(innerR);
    fMdcSegment[i].SetOutR(outR); fMdcSegment[i].SetZ(z);
    fMdcSegment[i].SetName(name);
    }
 
}

Referenced by BesMdcGeoParameter().

InitFromSvc()

void BesMdcGeoParameter::InitFromSvc

(

void

)

Definition at line 235 of file BesMdcGeoParameter.cc.

                                     {
    ISvcLocator* svcLocator = Gaudi::svcLocator();
    IMdcGeomSvc* ISvc;
    MdcGeomSvc* mdcGeomSvc;
    StatusCode sc=svcLocator->service("MdcGeomSvc", ISvc);
    mdcGeomSvc=dynamic_cast<MdcGeomSvc *>(ISvc);
    if (!sc.isSuccess()) 
    std::cout<<"BesMdcGeoParameter::Could not open Geometry Service"<<std::endl;
 
    fLayerNo= mdcGeomSvc->Misc()->LayerNo();
    fWireNo=mdcGeomSvc->Misc()->WireNo();
    fSignalLayerNo=mdcGeomSvc->Misc()->SLayerNo();
    fSignalWireR=mdcGeomSvc->Misc()->SWireR();
    fFieldWireR=mdcGeomSvc->Misc()->FWireR();
   
    int i,signalLayer;
    for(i=0; i<fSignalLayerNo; i++){
      signalLayer=mdcGeomSvc->Layer(i)->SLayer();  
      fSignalLayer[i]=signalLayer-1;
     }
     
    string name;
    int wireNo,firstWire;
    double length, r, phi,rotateCell,rotateAngle;
    for(i=0;i<fLayerNo;i++){    
      name=mdcGeomSvc->GeneralLayer(i)->LayerName();   
      wireNo=mdcGeomSvc->GeneralLayer(i)->NCell()*2;
      length= mdcGeomSvc->GeneralLayer(i)->Length();   
      r= mdcGeomSvc->GeneralLayer(i)->Radius();
      phi=mdcGeomSvc->GeneralLayer(i)->nomPhi();
      firstWire=mdcGeomSvc->GeneralLayer(i)->First();
      rotateCell= mdcGeomSvc->GeneralLayer(i)->nomShift();
      
      rotateAngle=2*pi*rotateCell/wireNo;
 
      fLayer[i].SetName(name);fLayer[i].SetRadius(r);             
      fLayer[i].SetLength(length); fLayer[i].SetRotateCell(rotateCell);
      fLayer[i].SetRotateAngle(rotateAngle); fLayer[i].SetWireNo(wireNo); 
      fLayer[i].SetShiftPhi(twopi/wireNo); fLayer[i].SetFirstWire(firstWire);
      fLayer[i].SetPhi(phi);
 
     if(i==0){
       fLayer[i].SetSumWireNo(wireNo); fLayer[i].SetBeginWireNo(0);
     }else{
       fLayer[i].SetBeginWireNo(fLayer[i-1].SumWireNo());
       fLayer[i].SetSumWireNo(fLayer[i-1].SumWireNo()+wireNo);
     }
 
    for(int j=0; j<wireNo; j++){
      fWirePhi[fLayer[i].BeginWireNo()+j]=j*fLayer[i].ShiftPhi()+phi;
      }
  }
 
  if(fLayer[fLayerNo-1].SumWireNo()!= fWireNo){
    cout<<"Total wire number is not consistant!"<<endl;
  }
  
 
  double innerR,outR,z;
  fSegmentNo=mdcGeomSvc->getSegmentNo();
  for(i=0;i<fSegmentNo;i++){    
    length=mdcGeomSvc->End(i)->Length();
    innerR=mdcGeomSvc->End(i)->InnerR();
    outR=mdcGeomSvc->End(i)->OutR();
    z=mdcGeomSvc->End(i)->Z();
    name=mdcGeomSvc->End(i)->Name();
 
    fMdcSegment[i].SetLength(length); fMdcSegment[i].SetInnerR(innerR);
    fMdcSegment[i].SetOutR(outR); fMdcSegment[i].SetZ(z);
    fMdcSegment[i].SetName(name);
    }
 
}

Referenced by BesMdcGeoParameter().

Layer()

const BesMdcLayer & BesMdcGeoParameter::Layer

(

int

layerNumber

)

const

Definition at line 133 of file BesMdcGeoParameter.cc.

                                                                  {
    if(layerNumber<0 || layerNumber>89){
    cout<<"Error: Wrong layerNo: "<<layerNumber<<endl;
    }
    return fLayer[layerNumber];
}

Referenced by BesMdcConstruction::Construct().

LayerNo()

int BesMdcGeoParameter::LayerNo ( void )

inline

Definition at line 130 of file BesMdcGeoParameter.hh.

{return fLayerNo;}  

Segment()

BesMdcMember BesMdcGeoParameter::Segment ( int x )

inline

Definition at line 136 of file BesMdcGeoParameter.hh.

{return fMdcSegment[x];}

Referenced by BesMdcConstruction::Construct().

SegmentNo()

int BesMdcGeoParameter::SegmentNo ( void )

inline

Definition at line 135 of file BesMdcGeoParameter.hh.

{return fSegmentNo;}

Referenced by BesMdcConstruction::Construct().

Signal2Global()

int BesMdcGeoParameter::Signal2Global ( int x )

inline

Definition at line 125 of file BesMdcGeoParameter.hh.

{return  fSignalLayer[x];}

Referenced by BesMdcConstruction::Construct().

SignalLayer()

const BesMdcLayer & BesMdcGeoParameter::SignalLayer

(

int

layerNumber

)

const

Definition at line 140 of file BesMdcGeoParameter.cc.

                                                                        {
    if(layerNumber<0 || layerNumber>42){
    cout<<"Error: Wrong SignallayerNo: "<<layerNumber<<endl;
    }
    return fLayer[fSignalLayer[layerNumber]];
}

SignalLayerNo()

int BesMdcGeoParameter::SignalLayerNo ( void )

inline

Definition at line 131 of file BesMdcGeoParameter.hh.

{return fSignalLayerNo;}

Referenced by BesMdcDigitizer::Digitize().

SignalWire()

BesMdcWire BesMdcGeoParameter::SignalWire	(	int	signalLayerNo,
		int	wireNo )

Definition at line 118 of file BesMdcGeoParameter.cc.

{
 
    int i=fSignalLayer[signalLayerNo];
    int wireNoInLayer=2*wireNo+1-fLayer[i].FirstWire();//FirstWire():0,field;1,signal
    double phi=fLayer[i].Phi();
    double shiftPhi=fLayer[i].ShiftPhi();
    double wirePhi;
    wirePhi= wireNoInLayer*shiftPhi+phi;
 
    BesMdcWire temp(fLayer[i].Length(), fWirePhi[fLayer[i].BeginWireNo()+wireNoInLayer], fLayer[i].R(),fLayer[i].RotateAngle());
    return temp;
}

Referenced by BesMdcSD::ProcessHits().

SignalWireR()

double BesMdcGeoParameter::SignalWireR ( )

inline

Definition at line 123 of file BesMdcGeoParameter.hh.

{return fSignalWireR;}

Referenced by BesMdcConstruction::Construct().

Wire()

BesMdcWire BesMdcGeoParameter::Wire

(

int

wireNo

)

Definition at line 102 of file BesMdcGeoParameter.cc.

                                             {
 
    int i;
    for(i=0; i<fLayerNo; i++){
      if(fLayer[i].BeginWireNo()<=wireNo && wireNo<fLayer[i].SumWireNo()){
 
    break;
      }
    }
 
   BesMdcWire temp(fLayer[i].Length(), fWirePhi[wireNo], fLayer[i].R(), fLayer[i].RotateAngle());
    return temp;
}

---

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

• BesMdcGeoParameter.hh
• BesMdcGeoParameter.cc