|
|
|
#include <RPCGeometryServTest.h>
Public Member Functions | |
| virtual void | analyze (const edm::Event &, const edm::EventSetup &) |
| const std::string & | myName () |
| RPCGeometryServTest (const edm::ParameterSet &pset) | |
| ~RPCGeometryServTest () | |
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_ |
Definition at line 30 of file RPCGeometryServTest.h.
| 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.
{
}
| 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_;}
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().
1.7.3