GeometricDetLoader Class Reference

#include <GeometricDetLoader.h>

Inheritance diagram for GeometricDetLoader:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob (edm::EventSetup const &)
virtual void	endJob ()
	GeometricDetLoader (const edm::ParameterSet &iConfig)
void	putOne (const GeometricDet *gd, PGeometricDet *pgd, int lev)
	~GeometricDetLoader ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

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
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
CurrentProcessingContext const *	currentContext () const
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

Definition at line 14 of file GeometricDetLoader.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 29 of file GeometricDetLoader.cc.

References gather_cfg::cout.

{
  std::cout<<"GeometricDetLoader::GeometricDetLoader"<<std::endl;
}

GeometricDetLoader::~GeometricDetLoader

(

)

Definition at line 34 of file GeometricDetLoader.cc.

References gather_cfg::cout.

{
  std::cout<<"GeometricDetLoader::~GeometricDetLoader"<<std::endl;
}

Member Function Documentation

virtual void GeometricDetLoader::analyze ( const edm::Event &  , const edm::EventSetup &    )

inlinevirtual

Implements edm::EDAnalyzer.

Definition at line 20 of file GeometricDetLoader.h.

{}

void GeometricDetLoader::beginJob ( edm::EventSetup const &  es )

virtual

Definition at line 40 of file GeometricDetLoader.cc.

References cond::service::PoolDBOutputService::beginOfTime(), GeometricDet::components(), prof2calltree::count, gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), GeometricDet::deepComponents(), cond::service::PoolDBOutputService::endOfTime(), edm::EventSetup::get(), edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), putOne(), and patCandidatesForDimuonsSequences_cff::tracker.

{
  std::cout<<"GeometricDetLoader::beginJob"<<std::endl;
  PGeometricDet* pgd = new PGeometricDet;
  edm::Service<cond::service::PoolDBOutputService> mydbservice;
  if( !mydbservice.isAvailable() ){
    std::cout<<"PoolDBOutputService unavailable"<<std::endl;
    return;
  }
  edm::ESHandle<DDCompactView> pDD;
  edm::ESHandle<GeometricDet> rDD;
  es.get<IdealGeometryRecord>().get( pDD );
  es.get<IdealGeometryRecord>().get( rDD );
  const GeometricDet* tracker = &(*rDD);

  // so now I have the tracker itself. loop over all its components to store them.
  putOne(tracker, pgd, 0);
  std::vector<const GeometricDet*> tc = tracker->components();
  std::cout <<"Tracker has " << tc.size() << " components." << std::endl; //, lets go through them." << std::endl;
  std::vector<const GeometricDet*>::const_iterator git = tc.begin();
  std::vector<const GeometricDet*>::const_iterator egit = tc.end();
  int count=0;
  int lev = 1;
  for (; git!= egit; ++git) {  // one level below "tracker"
    putOne(*git, pgd, lev);
    std::vector<const GeometricDet*> inone = (*git)->components();
    //    << ctste.name((*git)->type())
    //    std::cout << lev << " type " << (*git)->type() << " " << int((*git)->geographicalId()) << std::endl; // << " has " << inone.size() << " components." << std::endl;
    if ( inone.size() == 0 )  ++count;
    std::vector<const GeometricDet*>::const_iterator git2 = inone.begin();
    std::vector<const GeometricDet*>::const_iterator egit2 = inone.end();
    ++lev;
    for (; git2 != egit2; ++git2) { // level 2
      putOne(*git2, pgd, lev);
      std::vector<const GeometricDet*> intwo= (*git2)->components();
      //      std::cout << lev << "\ttype " << (*git2)->type() << " " << int((*git2)->geographicalId()) << std::endl; // << " has " << intwo.size() << " components." << std::endl;
      if ( intwo.size() == 0 )  ++count;
      std::vector<const GeometricDet*>::const_iterator git3 = intwo.begin();
      std::vector<const GeometricDet*>::const_iterator egit3 = intwo.end();
      ++lev;
      for (; git3 != egit3; ++git3) { // level 3
    putOne(*git3, pgd, lev);
    std::vector<const GeometricDet*> inthree= (*git3)->components();
    //  std::cout << lev << "\t\ttype " << (*git3)->type() << " " << int((*git3)->geographicalId()) << std::endl; // << " has " << inthree.size() << " components." << std::endl;
    if ( inthree.size() == 0 )  ++count;
    std::vector<const GeometricDet*>::const_iterator git4 = inthree.begin();
    std::vector<const GeometricDet*>::const_iterator egit4 = inthree.end();
    ++lev;
    for (; git4 != egit4; ++git4) { //level 4
      putOne(*git4, pgd, lev);
      std::vector<const GeometricDet*> infour= (*git4)->components();
      //      std::cout << lev << "\t\t\ttype " << (*git4)->type() << " " << int((*git4)->geographicalId()) << std::endl; // << " has " << infour.size() << " components." << std::endl;
      if ( infour.size() == 0 )  ++count;
      std::vector<const GeometricDet*>::const_iterator git5 = infour.begin();
      std::vector<const GeometricDet*>::const_iterator egit5 = infour.end();
      ++lev;
      for (; git5 != egit5; ++git5) { // level 5
        putOne(*git5, pgd, lev);
        std::vector<const GeometricDet*> infive= (*git5)->components();
        //      std::cout << lev << "\t\t\t\ttype " << (*git5)->type() << " " << int((*git5)->geographicalId()) << std::endl; // << " has " << infive.size() << " components." << std::endl;
        if ( infive.size() == 0 )  ++count;
        std::vector<const GeometricDet*>::const_iterator git6 = infive.begin();
        std::vector<const GeometricDet*>::const_iterator egit6 = infive.end();
        ++lev;
        for (; git6 != egit6; ++git6) { //level 6
          putOne(*git6, pgd, lev);
          std::vector<const GeometricDet*> insix= (*git6)->components();
          //          std::cout << lev << "\t\t\t\t\ttype " << (*git6)->type() << " " << int((*git6)->geographicalId()) << std::endl; // << " has " << insix.size() << " components." << std::endl;
          if ( insix.size() == 0 )  ++count;
        } // level 6
        --lev;
      } // level 5
      --lev;
    } // level 4
    --lev;
      } //level 3
      --lev;
    } // level 2
    --lev;
  }
  std::vector<const GeometricDet*> modules =  tracker->deepComponents();
  std::cout << " No. of Tracker components \"deepComponents\" = " << modules.size() << std::endl;
  std::cout << " Counted # of lowest \"leaves\" = " << count << std::endl; 
  if ( mydbservice->isNewTagRequest("PGeometricDetRcd") ) {
    mydbservice->createNewIOV<PGeometricDet>( pgd
                          , mydbservice->beginOfTime()
                          , mydbservice->endOfTime()
                          , "PGeometricDetRcd");
  } else {
    std::cout << "PGeometricDetRcd Tag is already present." << std::endl;
  }
}

virtual void GeometricDetLoader::endJob ( void  )

inlinevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 21 of file GeometricDetLoader.h.

{};

void GeometricDetLoader::putOne	(	const GeometricDet *	gd,
		PGeometricDet *	pgd,
		int	lev
	)

Definition at line 133 of file GeometricDetLoader.cc.

References PGeometricDet::Item::_a11, PGeometricDet::Item::_a12, PGeometricDet::Item::_a13, PGeometricDet::Item::_a21, PGeometricDet::Item::_a22, PGeometricDet::Item::_a23, PGeometricDet::Item::_a31, PGeometricDet::Item::_a32, PGeometricDet::Item::_a33, PGeometricDet::Item::_geographicalID, PGeometricDet::Item::_level, PGeometricDet::Item::_name, PGeometricDet::Item::_params0, PGeometricDet::Item::_params1, PGeometricDet::Item::_params10, PGeometricDet::Item::_params2, PGeometricDet::Item::_params3, PGeometricDet::Item::_params4, PGeometricDet::Item::_params5, PGeometricDet::Item::_params6, PGeometricDet::Item::_params7, PGeometricDet::Item::_params8, PGeometricDet::Item::_params9, PGeometricDet::Item::_phi, PGeometricDet::Item::_pixROCCols, PGeometricDet::Item::_pixROCRows, PGeometricDet::Item::_pixROCx, PGeometricDet::Item::_pixROCy, PGeometricDet::Item::_radLength, PGeometricDet::Item::_rho, PGeometricDet::Item::_shape, PGeometricDet::Item::_siliconAPVNum, PGeometricDet::Item::_stereo, PGeometricDet::Item::_type, PGeometricDet::Item::_x, PGeometricDet::Item::_xi, PGeometricDet::Item::_y, PGeometricDet::Item::_z, GeometricDet::geographicalID(), DDName::name(), GeometricDet::name(), GeometricDet::params(), PGeometricDet::pgeomdets_, GeometricDet::phi(), GeometricDet::pixROCCols(), GeometricDet::pixROCRows(), GeometricDet::pixROCx(), GeometricDet::pixROCy(), GeometricDet::radLength(), GeometricDet::rho(), makeMuonMisalignmentScenario::rot, GeometricDet::rotation(), GeometricDet::shape(), GeometricDet::siliconAPVNum(), GeometricDet::stereo(), GeometricDet::translation(), GeometricDet::type(), x, GeometricDet::xi(), detailsBasic3DVector::y, and detailsBasic3DVector::z.

Referenced by beginJob().

                                                                                      {

//   std::cout << "putting name: " << gd->name().name();
//   std::cout << " gid: " << gd->geographicalID();
//   std::cout << " type: " << gd->type() << std::endl;
//  std::cout << "shape = " << gd->shape()<<"; name = "<<gd->name().name()<<"; parameter number = "<<gd->params().size()<<std::endl;
  PGeometricDet::Item item;
  DDTranslation tran = gd->translation();
  DDRotationMatrix rot = gd->rotation();
  DD3Vector x, y, z;
  rot.GetComponents(x, y, z);
  item._name           = gd->name().name();
  item._level          = lev;
  item._x              = tran.X();
  item._y              = tran.Y();
  item._z              = tran.Z();
  item._phi            = gd->phi();
  item._rho            = gd->rho();
  item._a11            = x.X();
  item._a12            = y.X();
  item._a13            = z.X();
  item._a21            = x.Y();
  item._a22            = y.Y();
  item._a23            = z.Y();
  item._a31            = x.Z();
  item._a32            = y.Z();
  item._a33            = z.Z();
  item._shape          = gd->shape();
  item._type           = gd->type();
  if(gd->shape()==1){
    item._params0=gd->params()[0];
    item._params1=gd->params()[1];
    item._params2=gd->params()[2];
    item._params3=0;
    item._params4=0;
    item._params5=0;
    item._params6=0;
    item._params7=0;
    item._params8=0;
    item._params9=0;
    item._params10=0;
  }else if(gd->shape()==3){
    item._params0=gd->params()[0];
    item._params1=gd->params()[1];
    item._params2=gd->params()[2];
    item._params3=gd->params()[3];
    item._params4=gd->params()[4];
    item._params5=gd->params()[5];
    item._params6=gd->params()[6];
    item._params7=gd->params()[7];
    item._params8=gd->params()[8];
    item._params9=gd->params()[9];
    item._params10=gd->params()[10];
  }else{
    item._params0=0;
    item._params1=0;
    item._params2=0;
    item._params3=0;
    item._params4=0;
    item._params5=0;
    item._params6=0;
    item._params7=0;
    item._params8=0;
    item._params9=0;
    item._params10=0;
  } 
  item._geographicalID = gd->geographicalID();
  // FIXME: These are moved to PGeometricDetExtra:
  //item._volume         = gd->volume();
  //item._density        = gd->density();
  //item._weight         = gd->weight();
  //item._copy           = gd->copyno();
  //item._material       = gd->material();
  item._radLength      = gd->radLength();
  item._xi             = gd->xi();
  item._pixROCRows     = gd->pixROCRows();
  item._pixROCCols     = gd->pixROCCols();
  item._pixROCx        = gd->pixROCx();
  item._pixROCy        = gd->pixROCy();
  item._stereo         = gd->stereo();
  item._siliconAPVNum = gd->siliconAPVNum();
  pgd->pgeomdets_.push_back ( item );
}