GeometricTimingDetLoader Class Reference

#include <GeometricTimingDetLoader.h>

Inheritance diagram for GeometricTimingDetLoader:

Public Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &) override
void	beginJob () override
void	beginRun (edm::Run const &iEvent, edm::EventSetup const &) override
void	endJob () override
void	endRun (edm::Run const &iEvent, edm::EventSetup const &) override
	GeometricTimingDetLoader (const edm::ParameterSet &iConfig)
	~GeometricTimingDetLoader () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns >
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	putOne (const GeometricTimingDet *gd, PGeometricTimingDet *pgd, int lev)

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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 GeometricTimingDetLoader.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 26 of file GeometricTimingDetLoader.cc.

References LogDebug.

{
  LogDebug("GeometricTimingDetLoader") 
    <<"GeometricTimingDetLoader::GeometricTimingDetLoader"<<std::endl;
}

GeometricTimingDetLoader::~GeometricTimingDetLoader ( )

override

Definition at line 32 of file GeometricTimingDetLoader.cc.

References LogDebug.

{
  LogDebug("GeometricTimingDetLoader") 
    <<"GeometricTimingDetLoader::~GeometricTimingDetLoader"<<std::endl;
}

Member Function Documentation

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

inlineoverridevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 22 of file GeometricTimingDetLoader.h.

{}

void GeometricTimingDetLoader::beginJob ( void  )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 20 of file GeometricTimingDetLoader.h.

References beginRun(), and iEvent.

{}

void GeometricTimingDetLoader::beginRun ( edm::Run const &  iEvent, edm::EventSetup const &  es  )

override

Definition at line 39 of file GeometricTimingDetLoader.cc.

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

Referenced by beginJob().

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

  // so now I have the tracker itself. loop over all its components to store them.
  putOne(tracker, pgd, 0);
  std::vector<const GeometricTimingDet*> tc = tracker->components();
  std::cout <<"Tracker has " << tc.size() << " components." << std::endl; //, lets go through them." << std::endl;
  std::vector<const GeometricTimingDet*>::const_iterator git = tc.begin();
  std::vector<const GeometricTimingDet*>::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 GeometricTimingDet*> 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.empty() )  ++count;
    std::vector<const GeometricTimingDet*>::const_iterator git2 = inone.begin();
    std::vector<const GeometricTimingDet*>::const_iterator egit2 = inone.end();
    ++lev;
    for (; git2 != egit2; ++git2) { // level 2
      putOne(*git2, pgd, lev);
      std::vector<const GeometricTimingDet*> intwo= (*git2)->components();
      //      std::cout << lev << "\ttype " << (*git2)->type() << " " << int((*git2)->geographicalId()) << std::endl; // << " has " << intwo.size() << " components." << std::endl;
      if ( intwo.empty() )  ++count;
      std::vector<const GeometricTimingDet*>::const_iterator git3 = intwo.begin();
      std::vector<const GeometricTimingDet*>::const_iterator egit3 = intwo.end();
      ++lev;
      for (; git3 != egit3; ++git3) { // level 3
    putOne(*git3, pgd, lev);
    std::vector<const GeometricTimingDet*> inthree= (*git3)->components();
    //  std::cout << lev << "\t\ttype " << (*git3)->type() << " " << int((*git3)->geographicalId()) << std::endl; // << " has " << inthree.size() << " components." << std::endl;
    if ( inthree.empty() )  ++count;
    std::vector<const GeometricTimingDet*>::const_iterator git4 = inthree.begin();
    std::vector<const GeometricTimingDet*>::const_iterator egit4 = inthree.end();
    ++lev;
    for (; git4 != egit4; ++git4) { //level 4
      putOne(*git4, pgd, lev);
      std::vector<const GeometricTimingDet*> infour= (*git4)->components();
      //      std::cout << lev << "\t\t\ttype " << (*git4)->type() << " " << int((*git4)->geographicalId()) << std::endl; // << " has " << infour.size() << " components." << std::endl;
      if ( infour.empty() )  ++count;
      std::vector<const GeometricTimingDet*>::const_iterator git5 = infour.begin();
      std::vector<const GeometricTimingDet*>::const_iterator egit5 = infour.end();
      ++lev;
      for (; git5 != egit5; ++git5) { // level 5
        putOne(*git5, pgd, lev);
        std::vector<const GeometricTimingDet*> 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.empty() )  ++count;
        std::vector<const GeometricTimingDet*>::const_iterator git6 = infive.begin();
        std::vector<const GeometricTimingDet*>::const_iterator egit6 = infive.end();
        ++lev;
        for (; git6 != egit6; ++git6) { //level 6
          putOne(*git6, pgd, lev);
          std::vector<const GeometricTimingDet*> 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.empty() )  ++count;
        } // level 6
        --lev;
      } // level 5
      --lev;
    } // level 4
    --lev;
      } //level 3
      --lev;
    } // level 2
    --lev;
  }
  std::vector<const GeometricTimingDet*> 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("PGeometricTimingDetRcd") ) {
    mydbservice->createNewIOV<PGeometricTimingDet>( pgd
                          , mydbservice->beginOfTime()
                          , mydbservice->endOfTime()
                          , "PGeometricTimingDetRcd");
  } else {
    std::cout << "PGeometricTimingDetRcd Tag is already present." << std::endl;
  }
}

void GeometricTimingDetLoader::endJob ( void  )

inlineoverridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 24 of file GeometricTimingDetLoader.h.

References putOne().

{}

void GeometricTimingDetLoader::endRun ( edm::Run const &  iEvent, edm::EventSetup const &    )

inlineoverride

Definition at line 23 of file GeometricTimingDetLoader.h.

{}

void GeometricTimingDetLoader::putOne ( const GeometricTimingDet *  gd, PGeometricTimingDet *  pgd, int  lev  )

private

Definition at line 132 of file GeometricTimingDetLoader.cc.

References PGeometricTimingDet::Item::a11_, PGeometricTimingDet::Item::a12_, PGeometricTimingDet::Item::a13_, PGeometricTimingDet::Item::a21_, PGeometricTimingDet::Item::a22_, PGeometricTimingDet::Item::a23_, PGeometricTimingDet::Item::a31_, PGeometricTimingDet::Item::a32_, PGeometricTimingDet::Item::a33_, ddbox, ddtrap, DEFINE_FWK_MODULE, GeometricTimingDet::geographicalID(), PGeometricTimingDet::Item::geographicalID_, PGeometricTimingDet::Item::level_, DDName::name(), GeometricTimingDet::name(), PGeometricTimingDet::Item::name_, GeometricTimingDet::params(), PGeometricTimingDet::Item::params_0, PGeometricTimingDet::Item::params_1, PGeometricTimingDet::Item::params_10, PGeometricTimingDet::Item::params_2, PGeometricTimingDet::Item::params_3, PGeometricTimingDet::Item::params_4, PGeometricTimingDet::Item::params_5, PGeometricTimingDet::Item::params_6, PGeometricTimingDet::Item::params_7, PGeometricTimingDet::Item::params_8, PGeometricTimingDet::Item::params_9, PGeometricTimingDet::pgeomdets_, GeometricTimingDet::phi(), PGeometricTimingDet::Item::phi_, GeometricTimingDet::pixROCCols(), PGeometricTimingDet::Item::pixROCCols_, GeometricTimingDet::pixROCRows(), PGeometricTimingDet::Item::pixROCRows_, GeometricTimingDet::pixROCx(), PGeometricTimingDet::Item::pixROCx_, GeometricTimingDet::pixROCy(), PGeometricTimingDet::Item::pixROCy_, GeometricTimingDet::radLength(), PGeometricTimingDet::Item::radLength_, GeometricTimingDet::rho(), PGeometricTimingDet::Item::rho_, makeMuonMisalignmentScenario::rot, GeometricTimingDet::rotation(), GeometricTimingDet::shape(), PGeometricTimingDet::Item::shape_, GeometricTimingDet::siliconAPVNum(), PGeometricTimingDet::Item::siliconAPVNum_, GeometricTimingDet::stereo(), PGeometricTimingDet::Item::stereo_, GeometricTimingDet::translation(), GeometricTimingDet::type(), PGeometricTimingDet::Item::type_, x, PGeometricTimingDet::Item::x_, GeometricTimingDet::xi(), PGeometricTimingDet::Item::xi_, y, PGeometricTimingDet::Item::y_, z, and PGeometricTimingDet::Item::z_.

Referenced by beginRun(), and endJob().

                                                                                                        {

//   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;
  PGeometricTimingDet::Item item;
  const DDTranslation& tran = gd->translation();
  const 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_          = static_cast<int>(gd->shape());
  item.type_           = gd->type();
  if(gd->shape()==DDSolidShape::ddbox){
    item.params_0=gd->params()[0];
    item.params_1=gd->params()[1];
    item.params_2=gd->params()[2];
    item.params_3=0;
    item.params_4=0;
    item.params_5=0;
    item.params_6=0;
    item.params_7=0;
    item.params_8=0;
    item.params_9=0;
    item.params_10=0;
  }else if(gd->shape()==DDSolidShape::ddtrap){
    item.params_0=gd->params()[0];
    item.params_1=gd->params()[1];
    item.params_2=gd->params()[2];
    item.params_3=gd->params()[3];
    item.params_4=gd->params()[4];
    item.params_5=gd->params()[5];
    item.params_6=gd->params()[6];
    item.params_7=gd->params()[7];
    item.params_8=gd->params()[8];
    item.params_9=gd->params()[9];
    item.params_10=gd->params()[10];
  }else{
    item.params_0=0;
    item.params_1=0;
    item.params_2=0;
    item.params_3=0;
    item.params_4=0;
    item.params_5=0;
    item.params_6=0;
    item.params_7=0;
    item.params_8=0;
    item.params_9=0;
    item.params_10=0;
  } 
  item.geographicalID_ = gd->geographicalID();
  // FIXME: These are moved to PGeometricTimingDetExtra:
  //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 );
}