PhysicalPartsTree Class Reference

Inheritance diagram for PhysicalPartsTree:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginRun (const edm::Run &, const edm::EventSetup &)
	PhysicalPartsTree (const edm::ParameterSet &)
	~PhysicalPartsTree ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

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)
static void	prevalidate (ConfigurationDescriptions &)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Definition at line 26 of file PhysicalPartsTree.cc.

Constructor & Destructor Documentation

PhysicalPartsTree::PhysicalPartsTree ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 41 of file PhysicalPartsTree.cc.

{

}

PhysicalPartsTree::~PhysicalPartsTree

(

)

Definition at line 47 of file PhysicalPartsTree.cc.

{
}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 51 of file PhysicalPartsTree.cc.

References gather_cfg::cout.

                                                                                     { 
  std::cout << "analyze does nothing" << std::endl;
}

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

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 55 of file PhysicalPartsTree.cc.

References DDLogicalPart::category(), python.rootplot.argparse::category, DDEnums::categoryName(), DDExpandedView::copyno(), gather_cfg::cout, spr::find(), DDExpandedView::geoHistory(), edm::EventSetup::get(), estimatePileup::hist, DDExpandedView::logicalPart(), DDLogicalPart::material(), DDBase< N, C >::name(), DDExpandedView::next(), dbtoconf::parent, DDExpandedView::rotation(), DDLogicalPart::solid(), DDBase< N, C >::toString(), DDExpandedView::translation(), x, detailsBasic3DVector::y, and detailsBasic3DVector::z.

                                                                             {
  //set a tolerance for "near zero"
  double tolerance = 1.0e-3;

  std::string physicalPartsTreeFileName("PHYSICALPARTSTREE.dat");
  std::string logicalPartTypesFileName("LOGICALPARTTYPES.dat");
  std::string nominalPlacementsFileName("NOMINALPLACEMENTS.dat");
  std::string detectorPartsFileName("DETECTORPARTS.dat");
  std::string pospartsGraphFileName("POSPARTSGRAPH.dat");

  std::ofstream physicalPartsTreeOS(physicalPartsTreeFileName.c_str());
  std::ofstream logicalPartTypesOS(logicalPartTypesFileName.c_str());
  std::ofstream nominalPlacementsOS(nominalPlacementsFileName.c_str());
  std::ofstream detectorPartsOS(detectorPartsFileName.c_str());
  std::ofstream pospartsGraphOS(pospartsGraphFileName.c_str());


  std::string slpname;
  std::string scategory;
  std::vector<std::string> lpname_vec;


  std::cout << "PhysicalPartsTree Analyzer..." << std::endl;
  edm::ESTransientHandle<DDCompactView> pDD;

  iSetup.get<IdealGeometryRecord>().get( "", pDD );
  
  const DDCompactView & cpv = *pDD;
  DDExpandedView epv(cpv);
  int pospartid=0;
  size_t lastfound;
  std::string lgname;
   while(epv.next()){
 
    //for table physicalpartstree
    std::ostringstream parent, child,logicalpartid,parentid;
    child << epv.logicalPart().name();
    const DDGeoHistory & hist = epv.geoHistory();
    parent << hist[hist.size()-2].logicalPart().name();
    logicalpartid << epv.geoHistory();
    parentid<<hist[hist.size()-2];
    lgname=logicalpartid.str();
    lastfound=lgname.find_last_of("/\\");
    lgname=lgname.substr(lastfound+1,lgname.size());
    physicalPartsTreeOS<<lgname<<","<< parentid.str() << "," << child.str()<<std::endl;
//    std::cout << "Logical Part Name= "<<logicalpartid.str() << std::endl;
     
    //for table nominalPlacements
    bool reflection = false;

    DD3Vector x, y, z;
    epv.rotation().GetComponents(x, y, z);
    //Hep3Vector xv = epv.rotation().colX();
    //Hep3Vector yv = epv.rotation().colY();
    //Hep3Vector zv = epv.rotation().colZ();
    //if ( xv.cross(yv) * zv  < 0) {
    //            reflection = true;
    //          }
    if ( (1.0 + (x.Cross(y)).Dot(z)) <= tolerance ) {
      reflection = true;
    }
    std::vector<double> comps(9);
    epv.rotation().GetComponents(comps.begin(), comps.end());
    nominalPlacementsOS<< logicalpartid.str()<<","
               << epv.translation().X() << "," 
               << epv.translation().Y() << ","  
               << epv.translation().Z()<< "," 
//             << epv.rotation().xx()<< "," 
//             << epv.rotation().xy()<< "," 
//             << epv.rotation().xz()<< "," 
//             << epv.rotation().yx()<< "," 
//             << epv.rotation().yy()<< "," 
//             << epv.rotation().yz()<< "," 
//             << epv.rotation().zx()<< "," 
//             << epv.rotation().zy()<< "," 
//             << epv.rotation().zz()<< "," 
               << comps[0]<<","
               << comps[1]<<","
               << comps[2]<<","
               << comps[3]<<","
               << comps[4]<<","
               << comps[5]<<","
               << comps[6]<<","
               << comps[7]<<","
               << comps[8]<<","
               << (int)reflection
               <<std::endl;


    //for table DetectorParts
    detectorPartsOS <<epv.logicalPart().solid().toString() << ","
            <<epv.logicalPart().material().toString() << ","
            <<logicalpartid.str()
            <<std::endl;
   
    //for table PosPartsGraph 

    pospartsGraphOS<<++pospartid<< ","<<epv.copyno()<< ","
    <<parent.str() << "," << child.str()<<std::endl;


    //for table logicalPartTypes
    //this table is special because the DPNAMEs are subset of Physicalpart_id
    //while we're looping all the Physicalpart_ids in epv, so we just 
    //keep the unique DPNAMEs into a vector and save them in looping epv.
    // if there's no lpname in lpname_vec
    //>> the OS file.
 
    slpname=epv.logicalPart().toString();

    std::vector<std::string>::iterator it = find(lpname_vec.begin(), lpname_vec.end(), slpname);    
   
    if (it==lpname_vec.end()){
      if(DDEnums::categoryName(epv.logicalPart().category())){     
    std::ostringstream category;
    category<<DDEnums::categoryName(epv.logicalPart().category());
    scategory=category.str();
    std::string unspec="unspecified";
    if(scategory==unspec){
      lpname_vec.push_back(slpname);
      logicalPartTypesOS <<slpname<< ","
                 <<"DETECTORPARTS"<<std::endl;
    }
    else{
      lpname_vec.push_back(slpname);
      logicalPartTypesOS <<slpname<< ","<<scategory<<std::endl;

    }
      }
      else{
    lpname_vec.push_back(slpname);
    logicalPartTypesOS <<slpname<< ","
               <<"DETECTORPARTS"<<std::endl;

      }
    }



  } //finish looping through expanded view.
 

  physicalPartsTreeOS.close();
  logicalPartTypesOS.close();
  nominalPlacementsOS.close();
  detectorPartsOS.close();
  pospartsGraphOS.close();
}