RPCGeometryServTest.cc

Go to the documentation of this file.

#include <memory>
#include <fstream>
#include <FWCore/Framework/interface/Frameworkfwd.h>

#include <FWCore/Framework/interface/EDAnalyzer.h>
#include <FWCore/Framework/interface/Event.h>
#include <FWCore/Framework/interface/EventSetup.h>
#include <FWCore/Framework/interface/ESHandle.h>
#include <FWCore/ParameterSet/interface/ParameterSet.h>

#include "Geometry/RPCGeometry/interface/RPCGeometry.h"
#include "Geometry/RPCGeometry/interface/RPCGeomServ.h"
#include <Geometry/Records/interface/MuonGeometryRecord.h>
#include <Geometry/CommonTopologies/interface/RectangularStripTopology.h>
#include <Geometry/CommonTopologies/interface/TrapezoidalStripTopology.h>


#include "Validation/MuonRPCGeometry/plugins/RPCGeometryServTest.h"

#include <string>
#include <cmath>
#include <vector>
#include <map>
#include <iomanip>
#include <set>

using namespace std;


RPCGeometryServTest::RPCGeometryServTest( const edm::ParameterSet& iConfig )
  : dashedLineWidth_(104), dashedLine_( std::string(dashedLineWidth_, '-') ), 
    myName_( "RPCGeometryServTest" ) 
{ 
  std::cout <<"======================== Opening output file"<< std::endl;
}


RPCGeometryServTest::~RPCGeometryServTest() 
{
}

void
RPCGeometryServTest::analyze( const edm::Event& iEvent, const edm::EventSetup& iSetup )
{
  edm::ESHandle<RPCGeometry> pDD;
  iSetup.get<MuonGeometryRecord>().get( pDD );     

  std::cout << myName() << ": Analyzer..." << std::endl;
  std::cout << "start " << dashedLine_ << std::endl;


  std::cout << " Geometry node for RPCGeom is  " << &(*pDD) << std::endl;   
  cout << " I have "<<pDD->detTypes().size()    << " detTypes" << endl;
  cout << " I have "<<pDD->detUnits().size()    << " detUnits" << endl;
  cout << " I have "<<pDD->dets().size()        << " dets" << endl;
  cout << " I have "<<pDD->rolls().size()       << " rolls" << endl;
  cout << " I have "<<pDD->chambers().size()    << " chambers" << endl;

  std::cout << myName() << ": Begin iteration over geometry..." << std::endl;
  std::cout << "iter " << dashedLine_ << std::endl;

  

  int iRPCCHcount = 0;
  LocalPoint a(0.,0.,0.);
  for(TrackingGeometry::DetContainer::const_iterator it = 
    pDD->dets().begin(); it != pDD->dets().end(); it++){
    
    //----------------------- RPCCHAMBER TEST ---------------------------

    if( dynamic_cast< RPCChamber* >( *it ) != 0 ){
      ++iRPCCHcount;
      RPCChamber* ch = dynamic_cast< RPCChamber* >( *it ); 
      
      std::vector< const RPCRoll*> rollsRaf = (ch->rolls());
      for(std::vector<const RPCRoll*>::iterator r = rollsRaf.begin();
      r != rollsRaf.end(); ++r){
    std::cout << dashedLine_ << " NEW ROLL" << std::endl;
    std::cout<<"Region = "<<(*r)->id().region()
         <<"  Ring = "<<(*r)->id().ring()
         <<"  Station = "<<(*r)->id().station()
         <<"  Sector = "<<(*r)->id().sector()
         <<"  Layer = "<<(*r)->id().layer()
         <<"  Subsector = "<<(*r)->id().subsector()
         <<"  Roll = "<<(*r)->id().roll()<<std::endl;
    RPCGeomServ s( (*r)->id());
    GlobalPoint g= (*r)->toGlobal(a);
    std::cout <<s.name()<<" eta partition "<<s.eta_partition()<< " nroll="<<ch->nrolls()
          <<" z="<<g.z()<<" phi="<<g.phi()<<" R="<<g.perp()
          <<std::endl;

    if((*r)->id().region() ==  0 ){
      if (barzranges.find(s.eta_partition()) != barzranges.end()){
        std::pair<double, double> cic=barzranges.find(s.eta_partition())->second;
        double cmin = cic.first;
        double cmax = cic.second;
        if (g.z()<cmin)
          cmin = g.z();
        if (g.z()>cmax)
          cmax = g.z();
        std::pair<double,double>  cic2(cmin,cmax);
        barzranges[s.eta_partition()]=cic2;
      }else{
        std::pair<double, double> cic(g.z(),g.z());
        barzranges[s.eta_partition()]=cic;
      }
    }
    else if((*r)->id().region() ==  +1 ){
      if (forRranges.find(s.eta_partition()) != forRranges.end()){
        std::pair<double, double> cic=forRranges.find(s.eta_partition())->second;
        double cmin = cic.first;
        double cmax = cic.second;
        if (g.perp()<cmin)
          cmin = g.perp();
        if (g.perp()>cmax)
          cmax = g.perp();
        std::pair<double,double>  cic2(cmin,cmax);
        forRranges[s.eta_partition()]=cic2;
      }else{
        std::pair<double, double> cic(g.perp(),g.perp());
        forRranges[s.eta_partition()]=cic;
      }
    }
    else if((*r)->id().region() ==  -1 ){
      if (bacRranges.find(s.eta_partition()) != bacRranges.end()){
        std::pair<double, double> cic=bacRranges.find(s.eta_partition())->second;
        double cmin = cic.first;
        double cmax = cic.second;
        if (g.perp()<cmin)
          cmin = g.perp();
        if (g.perp()>cmax)
          cmax = g.perp();
        std::pair<double,double>  cic2(cmin,cmax);
        bacRranges[s.eta_partition()]=cic2;
      }else{
        std::pair<double, double> cic(g.perp(),g.perp());
        bacRranges[s.eta_partition()]=cic;
      }
    }
      }
    }
  }

  std::cout <<std::endl;
  std::map<int, std::pair<double, double> >::iterator ieta;

  for (ieta=bacRranges.begin(); ieta != bacRranges.end(); ieta++){
    std::cout <<" Eta "<<ieta->first<<" Radii = ( "<<ieta->second.first
          <<", "<<ieta->second.second<<")"<<std::endl;
  }

  for (ieta=barzranges.begin(); ieta != barzranges.end(); ieta++){
    std::cout <<" Eta "<<ieta->first<<" Z = ( "<<ieta->second.first
          <<", "<<ieta->second.second<<")"<<std::endl;
  }

  for (ieta=forRranges.begin(); ieta != forRranges.end(); ieta++){
    std::cout <<" Eta "<<ieta->first<<" Radii = ( "<<ieta->second.first
          <<", "<<ieta->second.second<<")"<<std::endl;
  }


  std::cout << dashedLine_ << " end" << std::endl;
}