#include <SiStripRegFEDSelector.h>

Inheritance diagram for SiStripRegFEDSelector:

Public Member Functions
	SiStripRegFEDSelector (const edm::ParameterSet &)

	~SiStripRegFEDSelector ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Attributes
double	delta_

edm::InputTag	rawInLabel_

edm::InputTag	seedLabel_

edm::ESHandle < SiStripRegionCabling >	strip_cabling

Private Member Functions
virtual void	beginJob ()

virtual void	endJob ()

virtual void	produce (edm::Event &, const edm::EventSetup &)

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 23 of file SiStripRegFEDSelector.h.

Constructor & Destructor Documentation

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

Definition at line 9 of file SiStripRegFEDSelector.cc.

References delta_, edm::ParameterSet::getParameter(), rawInLabel_, and seedLabel_.

 {
   seedLabel_=iConfig.getParameter<edm::InputTag>("regSeedLabel");
   delta_=iConfig.getParameter<double>("delta");
   
   rawInLabel_=iConfig.getParameter<edm::InputTag>("rawInputLabel");
 
   produces<FEDRawDataCollection>();
 }

SiStripRegFEDSelector::~SiStripRegFEDSelector ( )

Definition at line 20 of file SiStripRegFEDSelector.cc.

 {
  
 }

Member Function Documentation

void SiStripRegFEDSelector::beginJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 146 of file SiStripRegFEDSelector.cc.

146 {

147 }

void SiStripRegFEDSelector::endJob ( void )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 150 of file SiStripRegFEDSelector.cc.

150 {

151 }

void SiStripRegFEDSelector::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDProducer.

Definition at line 26 of file SiStripRegFEDSelector.cc.

References SiStripRegionCabling::ALLLAYERS, SiStripRegionCabling::ALLSUBDETS, FEDRawData::data(), runTheMatrix::data, delta_, FEDRawDataCollection::FEDData(), edm::EventSetup::get(), edm::Event::getByLabel(), i, j, gen::k, FEDNumbering::MAXFEDID, L1TEmulatorMonitor_cff::p, pos, funct::pow(), edm::Handle< T >::product(), edm::Event::put(), dttmaxenums::R, rawInLabel_, FEDRawData::resize(), seedLabel_, FEDRawData::size(), findQualityFiles::size, mathSSE::sqrt(), strip_cabling, and trigger::TriggerTrack.

Referenced by python.JSONExport.JsonExport::export(), and python.HTMLExport.HTMLExport::export().

 {
   std::auto_ptr<FEDRawDataCollection> producedData(new FEDRawDataCollection);
 
   edm::Handle<trigger::TriggerFilterObjectWithRefs> trigSeedTrks;
   iEvent.getByLabel(seedLabel_,trigSeedTrks);
   
   std::vector< edm::Ref<reco::IsolatedPixelTrackCandidateCollection> > isoPixTrackRefs;
   trigSeedTrks->getObjects(trigger::TriggerTrack, isoPixTrackRefs);
  
   edm::Handle<FEDRawDataCollection> rawIn;
   iEvent.getByLabel(rawInLabel_,rawIn);
   
   iSetup.get<SiStripRegionCablingRcd>().get(strip_cabling);
   
   std::vector<int> stripFEDVec; 
   
   //get vector of regions
   const SiStripRegionCabling::Cabling ccab=strip_cabling->getRegionCabling();
   
   //size of region (eta,phi)
   const std::pair<double,double> regDim=strip_cabling->regionDimensions();
   
   const SiStripRegionCabling::ElementCabling elcabling;
   
   bool fedSaved[1000];
   for (int i=0; i<1000; i++) fedSaved[i]=false;
   
   //cycle on seeds
   for (uint32_t p=0; p<isoPixTrackRefs.size(); p++)
     {
       double etaObj_=isoPixTrackRefs[p]->track()->eta();
       double phiObj_=isoPixTrackRefs[p]->track()->phi();
       
       //cycle on regions
       for (uint32_t i=0; i<ccab.size(); i++)
     {
       SiStripRegionCabling::Position pos=strip_cabling->position(i);
       double dphi=fabs(pos.second-phiObj_);
       if (dphi>acos(-1)) dphi=2*acos(-1)-dphi;
       double R=sqrt(pow(pos.first-etaObj_,2)+dphi*dphi);
       if (R-sqrt(pow(regDim.first/2,2)+pow(regDim.second/2,2))>delta_) continue;
       //get vector of subdets within region
       const SiStripRegionCabling::RegionCabling regSubdets= ccab[i];
       //cycle on subdets
       for (uint32_t idet=0; idet<SiStripRegionCabling::ALLSUBDETS; idet++)
         {
           //get vector of layers within subdet of region
           const SiStripRegionCabling::WedgeCabling regSubdetLayers=regSubdets[idet];    
           for (uint32_t ilayer=0; ilayer<SiStripRegionCabling::ALLLAYERS; ilayer++)
         {
           //get map of vectors of feds withing the layer of subdet of region
           const SiStripRegionCabling::ElementCabling fedVectorMap=regSubdetLayers[ilayer];
           SiStripRegionCabling::ElementCabling::const_iterator it=fedVectorMap.begin();
           for( ; it!=fedVectorMap.end(); it++)
             {
               for (uint32_t op=0; op<(it->second).size(); op++)
             {
               //get fed id 
               int fediid=(it->second)[op].fedId(); 
               if (!fedSaved[fediid]) 
                 {
                   stripFEDVec.push_back(fediid);
                 }
               fedSaved[fediid]=true;
               
             }
             }   
           
         }
           
         }
     }
     }
   
   
   
   const FEDRawDataCollection *rdc=rawIn.product();
   
   //   if ( ( rawData[i].provenance()->processName() != e.processHistory().rbegin()->processName() ) )
   //       continue ; // skip all raw collections not produced by the current process
   
   for ( int j=0; j< FEDNumbering::MAXFEDID; ++j ) 
     {
       bool rightFED=false;
       for (uint32_t k=0; k<stripFEDVec.size(); k++) 
     {
       if (j==stripFEDVec[k])
         {
           rightFED=true;
         }
     }
       if (!rightFED) continue;
       const FEDRawData & fedData = rdc->FEDData(j);
       size_t size=fedData.size();
       
       if ( size > 0 ) 
     {
       // this fed has data -- lets copy it
       FEDRawData & fedDataProd = producedData->FEDData(j);
       if ( fedDataProd.size() != 0 ) {
 //      std::cout << " More than one FEDRawDataCollection with data in FED ";
 //      std::cout << j << " Skipping the 2nd\n";
         continue;
       }
       fedDataProd.resize(size);
       unsigned char *dataProd=fedDataProd.data();
       const unsigned char *data=fedData.data();
       for ( unsigned int k=0; k<size; ++k ) {
         dataProd[k]=data[k];
       }
     }
     }
   
   iEvent.put(producedData);  
   
 }