RPCGeometryServTest Class Reference

#include <RPCGeometryServTest.h>

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Author

M. Maggi - INFN Bari

Definition at line 30 of file RPCGeometryServTest.h.

Constructor & Destructor Documentation

RPCGeometryServTest::RPCGeometryServTest

(

const edm::ParameterSet &

pset

)

Definition at line 35 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 43 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 48 of file RPCGeometryServTest.cc.

References a, bacRranges, barzranges, gather_cfg::cout, dashedLine_, forRranges, g, edm::EventSetup::get(), myName(), RPCChamber::nrolls(), alignCSCRings::r, RPCChamber::rolls(), alignCSCRings::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 ); 
      
      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 39 of file RPCGeometryServTest.h.

References myName_.

Referenced by analyze().

{ return myName_;}

Member Data Documentation

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

private

Definition at line 48 of file RPCGeometryServTest.h.

Referenced by analyze().

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

private

Definition at line 46 of file RPCGeometryServTest.h.

Referenced by analyze().

const std::string RPCGeometryServTest::dashedLine_

private

Definition at line 44 of file RPCGeometryServTest.h.

Referenced by analyze().

const int RPCGeometryServTest::dashedLineWidth_

private

Definition at line 43 of file RPCGeometryServTest.h.

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

private

Definition at line 47 of file RPCGeometryServTest.h.

Referenced by analyze().

const std::string RPCGeometryServTest::myName_

private

Definition at line 45 of file RPCGeometryServTest.h.

Referenced by myName().