CMS 3D CMS Logo

SiStripRegFEDSelector.cc
Go to the documentation of this file.
1 
7 
9 {
10  tok_seed_ = consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<edm::InputTag>("regSeedLabel"));
11  delta_=iConfig.getParameter<double>("delta");
12 
13  tok_raw_ = consumes<FEDRawDataCollection>(iConfig.getParameter<edm::InputTag>("rawInputLabel"));
14 
15  produces<FEDRawDataCollection>();
16 }
17 
18 
20 {
21 
22 }
23 
24 void
26 {
27  auto producedData = std::make_unique<FEDRawDataCollection>();
28 
30  iEvent.getByToken(tok_seed_,trigSeedTrks);
31 
32  std::vector< edm::Ref<reco::IsolatedPixelTrackCandidateCollection> > isoPixTrackRefs;
33  trigSeedTrks->getObjects(trigger::TriggerTrack, isoPixTrackRefs);
34 
36  iEvent.getByToken(tok_raw_,rawIn);
37 
39 
40  std::vector<int> stripFEDVec;
41 
42  //get vector of regions
44 
45  //size of region (eta,phi)
46  const std::pair<double,double> regDim=strip_cabling->regionDimensions();
47 
49 
50  bool fedSaved[1000];
51  for (int i=0; i<1000; i++) fedSaved[i]=false;
52 
53  //cycle on seeds
54  for (uint32_t p=0; p<isoPixTrackRefs.size(); p++)
55  {
56  double etaObj_=isoPixTrackRefs[p]->track()->eta();
57  double phiObj_=isoPixTrackRefs[p]->track()->phi();
58 
59  //cycle on regions
60  for (uint32_t i=0; i<ccab.size(); i++)
61  {
63  double dphi=fabs(pos.second-phiObj_);
64  if (dphi>acos(-1)) dphi=2*acos(-1)-dphi;
65  double R=sqrt(pow(pos.first-etaObj_,2)+dphi*dphi);
66  if (R-sqrt(pow(regDim.first/2,2)+pow(regDim.second/2,2))>delta_) continue;
67  //get vector of subdets within region
68  const SiStripRegionCabling::RegionCabling regSubdets= ccab[i];
69  //cycle on subdets
70  for (uint32_t idet=0; idet<SiStripRegionCabling::ALLSUBDETS; idet++)
71  {
72  //get vector of layers within subdet of region
73  const SiStripRegionCabling::WedgeCabling regSubdetLayers=regSubdets[idet];
74  for (uint32_t ilayer=0; ilayer<SiStripRegionCabling::ALLLAYERS; ilayer++)
75  {
76  //get map of vectors of feds withing the layer of subdet of region
77  const SiStripRegionCabling::ElementCabling fedVectorMap=regSubdetLayers[ilayer];
78  SiStripRegionCabling::ElementCabling::const_iterator it=fedVectorMap.begin();
79  for( ; it!=fedVectorMap.end(); it++)
80  {
81  for (uint32_t op=0; op<(it->second).size(); op++)
82  {
83  //get fed id
84  int fediid=(it->second)[op].fedId();
85  if (!fedSaved[fediid])
86  {
87  stripFEDVec.push_back(fediid);
88  }
89  fedSaved[fediid]=true;
90 
91  }
92  }
93 
94  }
95 
96  }
97  }
98  }
99 
100 
102 
103  const FEDRawDataCollection *rdc=rawIn.product();
104 
105  // if ( ( rawData[i].provenance()->processName() != e.processHistory().rbegin()->processName() ) )
106  // continue ; // skip all raw collections not produced by the current process
107 
108  for ( int j=0; j< FEDNumbering::MAXFEDID; ++j )
109  {
110  bool rightFED=false;
111  for (uint32_t k=0; k<stripFEDVec.size(); k++)
112  {
113  if (j==stripFEDVec[k])
114  {
115  rightFED=true;
116  }
117  }
118  if (!rightFED) continue;
119  const FEDRawData & fedData = rdc->FEDData(j);
120  size_t size=fedData.size();
121 
122  if ( size > 0 )
123  {
124  // this fed has data -- lets copy it
125  FEDRawData & fedDataProd = producedData->FEDData(j);
126  if ( fedDataProd.size() != 0 ) {
127 // std::cout << " More than one FEDRawDataCollection with data in FED ";
128 // std::cout << j << " Skipping the 2nd\n";
129  continue;
130  }
131  fedDataProd.resize(size);
132  unsigned char *dataProd=fedDataProd.data();
133  const unsigned char *data=fedData.data();
134  for ( unsigned int k=0; k<size; ++k ) {
135  dataProd[k]=data[k];
136  }
137  }
138  }
139 
140  iEvent.put(std::move(producedData));
141 
142 }
143 
144 void
146 }
147 
148 void
150 }
size
Write out results.
SiStripRegFEDSelector(const edm::ParameterSet &)
T getParameter(std::string const &) const
const std::pair< double, double > regionDimensions() const
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:137
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:579
void produce(edm::Event &, const edm::EventSetup &) override
std::vector< Element > ElementCabling
edm::ESHandle< SiStripRegionCabling > strip_cabling
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > tok_seed_
int iEvent
Definition: GenABIO.cc:230
const FEDRawData & FEDData(int fedid) const
retrieve data for fed
void resize(size_t newsize)
Definition: FEDRawData.cc:32
std::vector< ElementCabling > WedgeCabling
T sqrt(T t)
Definition: SSEVec.h:18
edm::EDGetTokenT< FEDRawDataCollection > tok_raw_
std::pair< double, double > Position
const Cabling & getRegionCabling() const
int k[5][pyjets_maxn]
const Position position(const uint32_t) const
T const * product() const
Definition: Handle.h:81
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
T get() const
Definition: EventSetup.h:63
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28
std::vector< WedgeCabling > RegionCabling
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
def move(src, dest)
Definition: eostools.py:510
std::vector< RegionCabling > Cabling