SiStripDetInfoFileWriter.cc

Go to the documentation of this file.

// -*- C++ -*-
// Package:    SiStripCommon
// Class:      SiStripDetInfoFileWriter
// Original Author:  G. Bruno
//         Created:  Mon May 20 10:04:31 CET 2007

#include "CalibTracker/SiStripCommon/plugins/SiStripDetInfoFileWriter.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h" 
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"


using namespace cms;
using namespace std;


SiStripDetInfoFileWriter::SiStripDetInfoFileWriter(const edm::ParameterSet& iConfig) {

  
  edm::LogInfo("SiStripDetInfoFileWriter::SiStripDetInfoFileWriter");

  filePath_ = iConfig.getUntrackedParameter<std::string>("FilePath",std::string("SiStripDetInfo.dat"));

}


SiStripDetInfoFileWriter::~SiStripDetInfoFileWriter(){

   edm::LogInfo("SiStripDetInfoFileWriter::~SiStripDetInfoFileWriter");
}



void SiStripDetInfoFileWriter::beginRun(const edm::Run&, const edm::EventSetup& iSetup){


  outputFile_.open(filePath_.c_str());

  if (outputFile_.is_open()){


    edm::ESHandle<TrackerGeometry> pDD;

    iSetup.get<TrackerDigiGeometryRecord>().get( pDD );

    edm::LogInfo("SiStripDetInfoFileWriter::beginRun - got geometry  ")<<std::endl;

    
    edm::LogInfo("SiStripDetInfoFileWriter") <<" There are "<<pDD->detUnits().size() <<" detectors"<<std::endl;
    
    for( const auto& it : pDD->detUnits()) {
  
      const StripGeomDetUnit* mit = dynamic_cast<StripGeomDetUnit const *>(it);

      if(mit!=nullptr){

    uint32_t detid=(mit->geographicalId()).rawId();
    double stripLength = mit->specificTopology().stripLength();
    unsigned short numberOfAPVs= mit->specificTopology().nstrips()/128;
    float thickness=mit->specificSurface().bounds().thickness();


    if(numberOfAPVs<1 || numberOfAPVs>6 ) {
      edm::LogError("SiStripDetInfoFileWriter")<<" Problem with Number of strips in detector.. "<< mit->specificTopology().nstrips() <<  "Will not write this entry to file"<< endl;
      continue;
    }

    outputFile_ << detid << " "<< numberOfAPVs << " " << stripLength << " "<< thickness << "\n";

      }

    }
    
    outputFile_.close();

  }
  
  else {

    edm::LogError("SiStripDetInfoFileWriter::beginRun - Unable to open file")<<endl;
    return;
  
  }

}