RPCGeometryServTest Class Reference

Inheritance diagram for RPCGeometryServTest:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
const std::string &	myName ()
	RPCGeometryServTest (const edm::ParameterSet &pset)
	~RPCGeometryServTest ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
std::map< int, std::pair < double, double > >	bacRranges
std::map< int, std::pair < double, double > >	barzranges
const std::string	dashedLine_
const int	dashedLineWidth_
std::map< int, std::pair < double, double > >	forRranges
const std::string	myName_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 31 of file RPCGeometryServTest.cc.

Constructor & Destructor Documentation

RPCGeometryServTest::RPCGeometryServTest

(

const edm::ParameterSet &

pset

)

Definition at line 52 of file RPCGeometryServTest.cc.

References gather_cfg::cout.

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

RPCGeometryServTest::~RPCGeometryServTest

(

)

Definition at line 60 of file RPCGeometryServTest.cc.

{
}

Member Function Documentation

void RPCGeometryServTest::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 65 of file RPCGeometryServTest.cc.

References a, bacRranges, barzranges, gather_cfg::cout, dashedLine_, forRranges, g, edm::EventSetup::get(), RPCChamber::id(), myName(), RPCChamber::nrolls(), csvReporter::r, RPCChamber::rolls(), asciidump::s, and edm::second().

{
  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 ); 
      RPCDetId detId=ch->id();
      
      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;
}

const std::string& RPCGeometryServTest::myName ( void  )

inline

Definition at line 40 of file RPCGeometryServTest.cc.

Referenced by analyze().

{ return myName_;}

Member Data Documentation

std::map<int, std::pair<double, double> > RPCGeometryServTest::bacRranges

private

Definition at line 49 of file RPCGeometryServTest.cc.

Referenced by analyze().

std::map<int, std::pair<double, double> > RPCGeometryServTest::barzranges

private

Definition at line 47 of file RPCGeometryServTest.cc.

Referenced by analyze().

const std::string RPCGeometryServTest::dashedLine_

private

Definition at line 45 of file RPCGeometryServTest.cc.

Referenced by analyze().

const int RPCGeometryServTest::dashedLineWidth_

private

Definition at line 44 of file RPCGeometryServTest.cc.

std::map<int, std::pair<double, double> > RPCGeometryServTest::forRranges

private

Definition at line 48 of file RPCGeometryServTest.cc.

Referenced by analyze().

const std::string RPCGeometryServTest::myName_

private

Definition at line 46 of file RPCGeometryServTest.cc.