CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
CSCSegtoRPC Class Reference

#include <CSCSegtoRPC.h>

Public Member Functions

 CSCSegtoRPC (edm::ConsumesCollector iC, const edm::ParameterSet &)
 
std::unique_ptr< RPCRecHitCollection > thePoints (CSCSegmentCollection const *allCSCSegments, edm::EventSetup const &iSetup, bool debug, double eyr)
 

Private Attributes

edm::ESGetToken< CSCGeometry, MuonGeometryRecord > cscGeoToken_
 
edm::ESGetToken< CSCObjectMap, MuonGeometryRecord > cscMapToken_
 
int maxBX
 
int minBX
 
edm::ESGetToken< RPCGeometry, MuonGeometryRecord > rpcGeoToken_
 

Detailed Description

Definition at line 18 of file CSCSegtoRPC.h.

Constructor & Destructor Documentation

◆ CSCSegtoRPC()

CSCSegtoRPC::CSCSegtoRPC ( edm::ConsumesCollector  iC,
const edm::ParameterSet &  iConfig 
)
explicit

Definition at line 16 of file CSCSegtoRPC.cc.

References edm::ParameterSet::getParameter(), maxBX, and minBX.

17  : rpcGeoToken_(iC.esConsumes()), cscGeoToken_(iC.esConsumes()), cscMapToken_(iC.esConsumes()) {
18  minBX = iConfig.getParameter<int>("minBX");
19  maxBX = iConfig.getParameter<int>("maxBX");
20 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
CSCSegtoRPC::cscGeoToken_
edm::ESGetToken< CSCGeometry, MuonGeometryRecord > cscGeoToken_
Definition: CSCSegtoRPC.h:28
CSCSegtoRPC::cscMapToken_
edm::ESGetToken< CSCObjectMap, MuonGeometryRecord > cscMapToken_
Definition: CSCSegtoRPC.h:29
CSCSegtoRPC::rpcGeoToken_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord > rpcGeoToken_
Definition: CSCSegtoRPC.h:27

Member Function Documentation

◆ thePoints()

std::unique_ptr< RPCRecHitCollection > CSCSegtoRPC::thePoints ( CSCSegmentCollection const *  allCSCSegments,
edm::EventSetup const &  iSetup,
bool  debug,
double  eyr 
)

Definition at line 22 of file CSCSegtoRPC.cc.

References cms::cuda::assert(), PV3DBase< T, PVType, FrameType >::barePhi(), CSCDetId::chamber(), CSCGeometry::chamber(), cscGeoToken_, cscMapToken_, debug, isotrackApplyRegressor::df, l1ctLayer1_cff::dr, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, CSCDetId::endcap(), edm::EventSetup::getHandle(), CSCObjectMap::getRolls(), RPCRoll::id(), CSCChamber::id(), TrapezoidalStripTopology::localPosition(), LogDebug, RPCResidualsHLT_cfi::MaxD, minBX, RPCRoll::nstrips(), PV3DBase< T, PVType, FrameType >::perp(), TrapezoidalStripTopology::pitch(), RPCDetId::ring(), CSCDetId::ring(), RPCGeometry::roll(), rpcGeoToken_, RPCGeomServ::segment(), RPCRecHit::setTimeAndError(), mathSSE::sqrt(), RPCDetId::station(), CSCDetId::station(), TrapezoidalStripTopology::stripLength(), GeomDet::surface(), GeomDet::toGlobal(), GeomDet::toLocal(), RPCRoll::topology(), X, PV3DBase< T, PVType, FrameType >::x(), TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, beamSpotPI::Y, PV3DBase< T, PVType, FrameType >::y(), beamSpotPI::Z, and PV3DBase< T, PVType, FrameType >::z().

25  {
26  edm::ESHandle<RPCGeometry> rpcGeo = iSetup.getHandle(rpcGeoToken_);
27  edm::ESHandle<CSCGeometry> cscGeo = iSetup.getHandle(cscGeoToken_);
28  edm::ESHandle<CSCObjectMap> cscMap = iSetup.getHandle(cscMapToken_);
29 
30  double MaxD = 80.;
31 
32  if (debug)
33  LogDebug("CSCSegtoRPC") << "CSC \t Number of CSC Segments in this event = " << allCSCSegments->size() << std::endl;
34 
35  auto _ThePoints = std::make_unique<RPCRecHitCollection>();
36  edm::OwnVector<RPCRecHit> RPCPointVector;
37 
38  if (allCSCSegments->size() == 0) {
39  if (debug)
40  LogDebug("CSCSegtoRPC") << "CSC 0 segments skiping event" << std::endl;
41  } else {
42  std::map<CSCDetId, int> CSCSegmentsCounter;
43  CSCSegmentCollection::const_iterator segment;
44 
45  for (segment = allCSCSegments->begin(); segment != allCSCSegments->end(); ++segment) {
46  CSCSegmentsCounter[segment->cscDetId()]++;
47  }
48 
49  float myTime = -9999.;
50  float myTimeErr = -9999.;
51  int myBx = -99;
52  const float bunchCrossTimeDiff = 25.; // time between bunch crossings
53 
54  if (debug)
55  LogDebug("CSCSegtoRPC") << "CSC \t loop over all the CSCSegments " << std::endl;
56  for (segment = allCSCSegments->begin(); segment != allCSCSegments->end(); ++segment) {
57  CSCDetId CSCId = segment->cscDetId();
58 
59  myTime = segment->time();
60  myBx = round(myTime / bunchCrossTimeDiff);
61  if (!(myBx <= maxBX && myBx >= minBX)) // rpc read in bx in [-2, 2]
62  continue;
63 
64  if (debug)
65  LogDebug("CSCSegtoRPC") << "CSC \t \t This Segment is in Chamber id: " << CSCId << std::endl;
66  if (debug)
67  LogDebug("CSCSegtoRPC") << "CSC \t \t Number of segments in this CSC = " << CSCSegmentsCounter[CSCId]
68  << std::endl;
69  if (debug)
70  LogDebug("CSCSegtoRPC")
71  << "CSC \t \t Is the only one in this CSC? is not ind the ring 1? Are there more than 2 "
72  "segments in the event?"
73  << std::endl;
74 
75  if (CSCSegmentsCounter[CSCId] == 1 && CSCId.ring() != 1 && allCSCSegments->size() >= 2) {
76  if (debug)
77  LogDebug("CSCSegtoRPC") << "CSC \t \t yes" << std::endl;
78  int cscEndCap = CSCId.endcap();
79  int cscStation = CSCId.station();
80  int cscRing = CSCId.ring();
81  int rpcRegion = 1;
82  if (cscEndCap == 2)
83  rpcRegion = -1; //Relacion entre las endcaps
84  int rpcRing = cscRing;
85  if (cscRing == 4)
86  rpcRing = 1;
87  int rpcStation = cscStation;
88  int rpcSegment = CSCId.chamber();
89 
90  LocalPoint segmentPosition = segment->localPosition();
91  LocalVector segmentDirection = segment->localDirection();
92  float dz = segmentDirection.z();
93 
94  if (debug)
95  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Information about the segment"
96  << "RecHits =" << segment->nRecHits() << "Angle =" << acos(dz) * 180 / 3.1415926
97  << std::endl;
98 
99  if (debug)
100  LogDebug("CSCSegtoRPC")
101  << "CSC \t \t Is a good Segment? dim = 4, 4 <= nRecHits <= 10 Incident angle int range 45 < "
102  << acos(dz) * 180 / 3.1415926 << " < 135? " << std::endl;
103 
104  if ((segment->dimension() == 4) && (segment->nRecHits() <= 10 && segment->nRecHits() >= 4)) {
105  if (debug)
106  LogDebug("CSCSegtoRPC") << "CSC \t \t yes" << std::endl;
107  if (debug)
108  LogDebug("CSCSegtoRPC") << "CSC \t \t CSC Segment Dimension " << segment->dimension() << std::endl;
109 
110  float Xo = segmentPosition.x();
111  float Yo = segmentPosition.y();
112  float Zo = segmentPosition.z();
113  float dx = segmentDirection.x();
114  float dy = segmentDirection.y();
115  float dz = segmentDirection.z();
116 
117  if (debug)
118  LogDebug("CSCSegtoRPC") << "Creating the CSCIndex" << std::endl;
119  CSCStationIndex theindex(rpcRegion, rpcStation, rpcRing, rpcSegment);
120  if (debug)
121  LogDebug("CSCSegtoRPC") << "Getting the Rolls for the given index" << std::endl;
122  std::set<RPCDetId> rollsForThisCSC = cscMap->getRolls(theindex);
123  if (debug)
124  LogDebug("CSCSegtoRPC") << "CSC \t \t Getting chamber from Geometry" << std::endl;
125  const CSCChamber* TheChamber = cscGeo->chamber(CSCId);
126  if (debug)
127  LogDebug("CSCSegtoRPC") << "CSC \t \t Getting ID from Chamber" << std::endl;
128  const CSCDetId TheId = TheChamber->id();
129 
130  if (debug)
131  LogDebug("CSCSegtoRPC") << "CSC \t \t Number of rolls for this CSC = " << rollsForThisCSC.size()
132  << std::endl;
133 
134  if (debug)
135  LogDebug("CSCSegtoRPC") << "CSC \t \t Printing The Id" << TheId << std::endl;
136 
137  if (rpcRing != 1) { //They don't exist in Run3!
138 
139  assert(!rollsForThisCSC.empty());
140 
141  if (debug)
142  LogDebug("CSCSegtoRPC") << "CSC \t \t Loop over all the rolls asociated to this CSC" << std::endl;
143  for (std::set<RPCDetId>::iterator iteraRoll = rollsForThisCSC.begin(); iteraRoll != rollsForThisCSC.end();
144  iteraRoll++) {
145  const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll);
146  RPCDetId rpcId = rollasociated->id();
147 
148  if (debug)
149  LogDebug("CSCSegtoRPC") << "CSC \t \t \t We are in the roll getting the surface" << rpcId << std::endl;
150  const BoundPlane& RPCSurface = rollasociated->surface();
151 
152  if (debug)
153  LogDebug("CSCSegtoRPC") << "CSC \t \t \t RollID: " << rpcId << std::endl;
154 
155  if (debug)
156  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Doing the extrapolation to this roll" << std::endl;
157  if (debug)
158  LogDebug("CSCSegtoRPC") << "CSC \t \t \t CSC Segment Direction in CSCLocal " << segmentDirection
159  << std::endl;
160  if (debug)
161  LogDebug("CSCSegtoRPC") << "CSC \t \t \t CSC Segment Point in CSCLocal " << segmentPosition
162  << std::endl;
163 
164  GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0, 0, 0));
165  if (debug)
166  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Center (0,0,0) of the Roll in Global" << CenterPointRollGlobal
167  << std::endl;
168  GlobalPoint CenterPointCSCGlobal = TheChamber->toGlobal(LocalPoint(0, 0, 0));
169  if (debug)
170  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Center (0,0,0) of the CSC in Global" << CenterPointCSCGlobal
171  << std::endl;
172  GlobalPoint segmentPositionInGlobal =
173  TheChamber->toGlobal(segmentPosition); //new way to convert to global
174  if (debug)
175  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Segment Position in Global" << segmentPositionInGlobal
176  << std::endl;
177  LocalPoint CenterRollinCSCFrame = TheChamber->toLocal(CenterPointRollGlobal);
178 
179  if (debug) { //to check CSC RPC phi relation!
180  float rpcphi = 0;
181  float cscphi = 0;
182 
183  (CenterPointRollGlobal.barePhi() < 0) ? rpcphi = 2 * 3.141592 + CenterPointRollGlobal.barePhi()
184  : rpcphi = CenterPointRollGlobal.barePhi();
185 
186  (CenterPointCSCGlobal.barePhi() < 0) ? cscphi = 2 * 3.1415926536 + CenterPointCSCGlobal.barePhi()
187  : cscphi = CenterPointCSCGlobal.barePhi();
188 
189  float df = fabs(cscphi - rpcphi);
190  float dr = fabs(CenterPointRollGlobal.perp() - CenterPointCSCGlobal.perp());
191  float diffz = CenterPointRollGlobal.z() - CenterPointCSCGlobal.z();
192  float dfg = df * 180. / 3.14159265;
193 
194  if (debug)
195  LogDebug("CSCSegtoRPC") << "CSC \t \t \t z of RPC=" << CenterPointRollGlobal.z() << "z of CSC"
196  << CenterPointCSCGlobal.z() << " dfg=" << dfg << std::endl;
197 
198  RPCGeomServ rpcsrv(rpcId);
199 
200  if (dr > 200. || fabs(dz) > 55. || dfg > 1.) {
201  if (debug)
202  LogDebug("CSCSegtoRPC")
203  << "\t \t \t CSC Station= " << CSCId.station() << " Ring= " << CSCId.ring()
204  << " Chamber= " << CSCId.chamber() << " cscphi=" << cscphi * 180 / 3.14159265
205  << "\t RPC Station= " << rpcId.station() << " ring= " << rpcId.ring() << " segment =-> "
206  << rpcsrv.segment() << " rollphi=" << rpcphi * 180 / 3.14159265 << "\t dfg=" << dfg
207  << " dz=" << diffz << " dr=" << dr << std::endl;
208  }
209  }
210 
211  float D = CenterRollinCSCFrame.z();
212 
213  float X = Xo + dx * D / dz;
214  float Y = Yo + dy * D / dz;
215  float Z = D;
216 
217  const TrapezoidalStripTopology* top_ =
218  dynamic_cast<const TrapezoidalStripTopology*>(&(rollasociated->topology()));
219  LocalPoint xmin = top_->localPosition(0.);
220  if (debug)
221  LogDebug("CSCSegtoRPC") << "CSC \t \t \t xmin of this Roll " << xmin << "cm" << std::endl;
222  LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
223  if (debug)
224  LogDebug("CSCSegtoRPC") << "CSC \t \t \t xmax of this Roll " << xmax << "cm" << std::endl;
225  float rsize = fabs(xmax.x() - xmin.x());
226  if (debug)
227  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Roll Size " << rsize << "cm" << std::endl;
228  float stripl = top_->stripLength();
229  float stripw = top_->pitch();
230 
231  if (debug)
232  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Strip Lenght " << stripl << "cm" << std::endl;
233  if (debug)
234  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Strip Width " << stripw << "cm" << std::endl;
235 
236  if (debug)
237  LogDebug("CSCSegtoRPC") << "CSC \t \t \t X Predicted in CSCLocal= " << X << "cm" << std::endl;
238  if (debug)
239  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Y Predicted in CSCLocal= " << Y << "cm" << std::endl;
240  if (debug)
241  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Z Predicted in CSCLocal= " << Z << "cm" << std::endl;
242 
243  float extrapolatedDistance = sqrt((X - Xo) * (X - Xo) + (Y - Yo) * (Y - Yo) + (Z - Zo) * (Z - Zo));
244 
245  if (debug)
246  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Is the distance of extrapolation less than MaxD? ="
247  << extrapolatedDistance << "cm"
248  << " MaxD=" << MaxD << "cm" << std::endl;
249 
250  if (extrapolatedDistance <= MaxD) {
251  if (debug)
252  LogDebug("CSCSegtoRPC") << "CSC \t \t \t yes" << std::endl;
253 
254  GlobalPoint GlobalPointExtrapolated = TheChamber->toGlobal(LocalPoint(X, Y, Z));
255  if (debug)
256  LogDebug("CSCSegtoRPC")
257  << "CSC \t \t \t Point ExtraPolated in Global" << GlobalPointExtrapolated << std::endl;
258 
259  LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
260 
261  if (debug)
262  LogDebug("CSCSegtoRPC")
263  << "CSC \t \t \t Point Extrapolated in RPCLocal" << PointExtrapolatedRPCFrame << std::endl;
264  if (debug)
265  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Corner of the Roll = (" << rsize * eyr << "," << stripl * eyr
266  << ")" << std::endl;
267  if (debug)
268  LogDebug("CSCSegtoRPC")
269  << "CSC \t \t \t Info About the Point Extrapolated in X Abs ("
270  << fabs(PointExtrapolatedRPCFrame.x()) << "," << fabs(PointExtrapolatedRPCFrame.y()) << ","
271  << fabs(PointExtrapolatedRPCFrame.z()) << ")" << std::endl;
272  if (debug)
273  LogDebug("CSCSegtoRPC") << "CSC \t \t \t dz=" << fabs(PointExtrapolatedRPCFrame.z())
274  << " dx=" << fabs(PointExtrapolatedRPCFrame.x())
275  << " dy=" << fabs(PointExtrapolatedRPCFrame.y()) << std::endl;
276 
277  if (debug)
278  LogDebug("CSCSegtoRPC") << "CSC \t \t \t Does the extrapolation go inside this roll????" << std::endl;
279 
280  if (fabs(PointExtrapolatedRPCFrame.z()) < 1. && fabs(PointExtrapolatedRPCFrame.x()) < rsize * eyr &&
281  fabs(PointExtrapolatedRPCFrame.y()) < stripl * eyr) {
282  if (debug)
283  LogDebug("CSCSegtoRPC") << "CSC \t \t \t \t yes" << std::endl;
284  if (debug)
285  LogDebug("CSCSegtoRPC") << "CSC \t \t \t \t Creating the RecHit" << std::endl;
286 
287  RPCRecHit RPCPoint(rpcId, myBx, PointExtrapolatedRPCFrame);
288  RPCPoint.setTimeAndError(myTime, myTimeErr);
289 
290  if (debug)
291  LogDebug("CSCSegtoRPC") << "CSC \t \t \t \t Clearing the vector" << std::endl;
292  RPCPointVector.clear();
293  if (debug)
294  LogDebug("CSCSegtoRPC") << "CSC \t \t \t \t Pushing back" << std::endl;
295  RPCPointVector.push_back(RPCPoint);
296  if (debug)
297  LogDebug("CSCSegtoRPC") << "CSC \t \t \t \t Putting the vector" << std::endl;
298  _ThePoints->put(rpcId, RPCPointVector.begin(), RPCPointVector.end());
299  }
300  }
301  }
302  }
303  }
304  }
305  }
306  }
307  return _ThePoints;
308 }
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:69
align::LocalPoint
Point3DBase< Scalar, LocalTag > LocalPoint
Definition: Definitions.h:30
CSCGeometry::chamber
const CSCChamber * chamber(CSCDetId id) const
Return the chamber corresponding to given DetId.
Definition: CSCGeometry.cc:100
GeomDet::toLocal
LocalPoint toLocal(const GlobalPoint &gp) const
Conversion to the R.F. of the GeomDet.
Definition: GeomDet.h:58
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
TrackerOfflineValidation_Dqm_cff.xmin
xmin
Definition: TrackerOfflineValidation_Dqm_cff.py:10
X
#define X(str)
Definition: MuonsGrabber.cc:38
RPCDetId::ring
int ring() const
Definition: RPCDetId.h:59
CSCSegtoRPC::cscGeoToken_
edm::ESGetToken< CSCGeometry, MuonGeometryRecord > cscGeoToken_
Definition: CSCSegtoRPC.h:28
RPCGeometry::roll
const RPCRoll * roll(RPCDetId id) const
Return a roll given its id.
Definition: RPCGeometry.cc:50
RPCRoll::nstrips
int nstrips() const
Definition: RPCRoll.cc:24
cms::cuda::assert
assert(be >=bs)
PV3DBase::barePhi
T barePhi() const
Definition: PV3DBase.h:65
edm::RangeMap::const_iterator
C::const_iterator const_iterator
constant access iterator type
Definition: RangeMap.h:43
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:23
CSCDetId::chamber
int chamber() const
Definition: CSCDetId.h:62
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:49
RPCResidualsHLT_cfi.MaxD
MaxD
Definition: RPCResidualsHLT_cfi.py:18
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
debug
#define debug
Definition: HDRShower.cc:19
TrackerOfflineValidation_Dqm_cff.xmax
xmax
Definition: TrackerOfflineValidation_Dqm_cff.py:11
funct::D
DecomposeProduct< arg, typename Div::arg > D
Definition: Factorize.h:141
CSCDetId::station
int station() const
Definition: CSCDetId.h:79
CSCDetId::endcap
int endcap() const
Definition: CSCDetId.h:85
RPCDetId::station
int station() const
Definition: RPCDetId.h:78
l1ctLayer1_cff.dr
dr
Definition: l1ctLayer1_cff.py:70
RPCRoll::topology
const Topology & topology() const override
Definition: RPCRoll.cc:18
RPCRoll::id
RPCDetId id() const
Definition: RPCRoll.cc:16
CSCDetId::ring
int ring() const
Definition: CSCDetId.h:68
TrapezoidalStripTopology::stripLength
float stripLength() const override
det heigth (strip length in the middle)
Definition: TrapezoidalStripTopology.h:63
CSCSegtoRPC::cscMapToken_
edm::ESGetToken< CSCObjectMap, MuonGeometryRecord > cscMapToken_
Definition: CSCSegtoRPC.h:29
CSCSegtoRPC::rpcGeoToken_
edm::ESGetToken< RPCGeometry, MuonGeometryRecord > rpcGeoToken_
Definition: CSCSegtoRPC.h:27
CSCObjectMap::getRolls
std::set< RPCDetId > const & getRolls(CSCStationIndex index) const
Definition: CSCObjectMap.cc:42
TrapezoidalStripTopology::localPosition
LocalPoint localPosition(float strip) const override
Definition: TrapezoidalStripTopology.cc:30
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:241
CSCChamber::id
CSCDetId id() const
Get the (concrete) DetId.
Definition: CSCChamber.h:34

Member Data Documentation

◆ cscGeoToken_

edm::ESGetToken<CSCGeometry, MuonGeometryRecord> CSCSegtoRPC::cscGeoToken_
private

Definition at line 28 of file CSCSegtoRPC.h.

Referenced by thePoints().

◆ cscMapToken_

edm::ESGetToken<CSCObjectMap, MuonGeometryRecord> CSCSegtoRPC::cscMapToken_
private

Definition at line 29 of file CSCSegtoRPC.h.

Referenced by thePoints().

◆ maxBX

int CSCSegtoRPC::maxBX
private

Definition at line 31 of file CSCSegtoRPC.h.

Referenced by CSCSegtoRPC().

◆ minBX

int CSCSegtoRPC::minBX
private

Definition at line 30 of file CSCSegtoRPC.h.

Referenced by CSCSegtoRPC(), and thePoints().

◆ rpcGeoToken_

edm::ESGetToken<RPCGeometry, MuonGeometryRecord> CSCSegtoRPC::rpcGeoToken_
private

Definition at line 27 of file CSCSegtoRPC.h.

Referenced by thePoints().

Generated for CMSSW Reference Manual by doxygen 1.8.14