TestGeometry.cxx

Go to the documentation of this file.

// **************************************************************************
// authors: L. Lavezzi (univ. of Torino & INFN, Italy)
//
 
//include system lib
#include <iostream>
#include <iomanip>
#include <string>
#include <cmath>
 
// Include files
#include "GaudiKernel/AlgFactory.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/IEventProcessor.h"
 
#include "GaudiKernel/Incident.h"
#include "GaudiKernel/IIncidentSvc.h"
#include "GaudiKernel/Memory.h"
 
#include <csignal>
 
// for save digit & cluster
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/RegistryEntry.h"
#include "GaudiKernel/MsgStream.h"
 
#include "CgemRawEvent/CgemDigi.h"
#include "CgemRecEvent/RecCgemCluster.h"
 
#include "Identifier/CgemID.h"
#include "RawEvent/RawDataUtil.h"
#include "RawEvent/DigiEvent.h"
#include "ReconEvent/ReconEvent.h"
#include "EventModel/EventHeader.h"
#include "GaudiKernel/SmartDataPtr.h"
 
#include "ReadCosmicRayData/TestGeometry.h"
#include "TMath.h"
//using namespace std;
                                                                                                                         
 
TestGeometry::TestGeometry(const std::string& name, ISvcLocator* pSvcLocator):
  Algorithm(name,pSvcLocator){
 
  declareProperty("LUTfile", lutfile = "/bes3fs/cgemCosmic/data/CGEM_cosmic_look_up_table_from_17_to_17.root");
 
}
 
TestGeometry::~TestGeometry(){}
 
StatusCode TestGeometry::initialize(){
  MsgStream log(msgSvc(), name());
 
  cout << "TestGeometry initialize" << endl;
  StatusCode sc = service("CgemGeomSvc", m_geomSvc);
  if(sc != StatusCode::SUCCESS) {
    log << MSG::ERROR << "can not use CgemGeomSvc" << endreq;
    return StatusCode::FAILURE;
  }
 
  // LUT 
  lutreader = new CgemLUTReader(lutfile);
  lutreader->ReadLUT();
 
  const int nlayer = 3; 
  const int nsheet[nlayer] = {1, 2, 2};
  const int nview = 2;
 
  for(int ilayer = 0; ilayer < nlayer; ilayer++) {
    for(int isheet = 0; isheet < nsheet[ilayer]; isheet++) {
      for(int iview=0; iview < nview; iview++) {
 
    CgemGeoReadoutPlane* anode = m_geomSvc->getReadoutPlane(ilayer, isheet);
    int layerid = anode->getLayerId();
    int sheetid = anode->getSheetId();
    int nxstrip = anode->getNXstrips();
    int nvstrip = anode->getNVstrips();
    
    int nstrip = nxstrip;
    if(iview==1) nstrip = nvstrip;
 
    for(int istrip=0; istrip < nstrip; istrip++) {
      const Identifier ident = CgemID::strip_id(layerid, sheetid, iview, istrip);
      bool isxstrip = CgemID::is_xstrip(ident);
      int  stripid = CgemID::strip(ident);
      double length = anode->getLength();
      if(isxstrip==false) length = anode->getVStripLength(stripid);
 
      cout << layerid << " " << sheetid << " " << iview << " " << stripid << " " << length << endl;
 
 
    }
      }
    }
  }
 
 
  return StatusCode::SUCCESS;
}
 
 
StatusCode TestGeometry::execute(){
    MsgStream log(msgSvc(), name());
    cout << "->TestGeometry::execute" << endl;
    cout << "geomtry service pointer " << m_geomSvc << endl;
    return StatusCode::SUCCESS;
}
 
 
 
StatusCode TestGeometry::finalize(){
  MsgStream log(msgSvc(),name());
  log << MSG::INFO << "TestGeometry finalize()" << endreq;
  
  return StatusCode::SUCCESS;
}