49 PSimHitContainer::const_iterator simIt;
52 std::map<RPCDetId, std::vector<double> > allsims;
54 for (simIt = simHit->begin(); simIt != simHit->end(); simIt++)
57 const RPCRoll* soll =
dynamic_cast<const RPCRoll*
> (rpcGeom->roll(Rsid));
58 int ptype = simIt->particleType();
60 if (ptype == 13 || ptype == -13)
63 std::vector<double> buff;
64 if (allsims.find(Rsid) != allsims.end())
69 buff.push_back(simIt->localPosition().x());
87 else if (Rsid.
region() == (-1))
98 for (detUnitIt = rpcDigis->begin(); detUnitIt != rpcDigis->end(); ++detUnitIt)
100 const RPCDetId Rsid = (*detUnitIt).first;
101 const RPCRoll* roll =
dynamic_cast<const RPCRoll*
> (rpcGeom->roll(Rsid));
104 std::vector<double> sims;
105 if (allsims.find(Rsid) != allsims.end())
107 sims = allsims[Rsid];
116 if (Rsid.
region() == (+1))
121 else if (Rsid.
region() == (-1))
129 if (sims.size() == 0)
149 if (Rsid.
ring() == 2)
153 if (Rsid.
roll() == 1)
155 if (Rsid.
roll() == 2)
157 if (Rsid.
roll() == 3)
162 if (Rsid.
roll() == 1)
164 if (Rsid.
roll() == 2)
166 if (Rsid.
roll() == 3)
170 if (Rsid.
ring() == 3)
172 if (Rsid.
roll() == 1)
174 if (Rsid.
roll() == 2)
176 if (Rsid.
roll() == 3)
183 if (sims.size() == 1 && ndigi == 1)
185 double dis = roll->
centreOfStrip(range.first->strip()).
x() - sims[0];
190 if (Rsid.
ring() == -2)
192 else if (Rsid.
ring() == -1)
194 else if (Rsid.
ring() == 0)
196 else if (Rsid.
ring() == 1)
198 else if (Rsid.
ring() == 2)
217 if (Rsid.
ring() == 2)
221 if (Rsid.
roll() == 1)
223 if (Rsid.
roll() == 2)
225 if (Rsid.
roll() == 3)
230 if (Rsid.
roll() == 1)
232 if (Rsid.
roll() == 2)
234 if (Rsid.
roll() == 3)
238 if (Rsid.
ring() == 3)
240 if (Rsid.
roll() == 1)
242 if (Rsid.
roll() == 2)
244 if (Rsid.
roll() == 3)
249 if (Rsid.
region() == (+1))
261 if (Rsid.
region() == (-1))
MonitorElement * CLS_Endcap_1_Ring2_A
MonitorElement * CLS_Endcap_123_Ring3_A
MonitorElement * CLS_Endcap_1_Ring2_B
LocalPoint centreOfStrip(int strip) const
MonitorElement * CLS_Endcap_23_Ring2_A
MonitorElement * noiseCLSEndcaps
MonitorElement * ResLayer6_barrel
MonitorElement * Res_Endcap1_Ring2_A
MonitorElement * Res_Endcap123_Ring3_A
MonitorElement * ResLayer3_barrel
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
MonitorElement * CLS_Endcap_4
MonitorElement * ResWplu2
MonitorElement * BxDisc_4Min
MonitorElement * StripProf
MonitorElement * ResWzer0
MonitorElement * ResDmin3
MonitorElement * ResLayer5_barrel
MonitorElement * ResWmin1
MonitorElement * ResLayer2_barrel
MonitorElement * ResDplu4
MonitorElement * ResDplu3
MonitorElement * ResDmin1
MonitorElement * ResWplu1
MonitorElement * Res_Endcap123_Ring3_C
MonitorElement * BxDisc_4Plus
Abs< T >::type abs(const T &t)
MonitorElement * Res_Endcap123_Ring3_B
MonitorElement * CLS_Endcap_123_Ring3_B
MonitorElement * ResWmin2
MonitorElement * xyvDmin4
MonitorElement * CLS_Endcap_1_Ring2_C
MonitorElement * ResDplu2
MonitorElement * Res_Endcap23_Ring2_C
MonitorElement * ResLayer1_barrel
MonitorElement * Res_Endcap1_Ring2_B
MonitorElement * noiseCLS
MonitorElement * ResDmin2
std::vector< RPCDigi >::const_iterator const_iterator
edm::EDGetTokenT< RPCDigiCollection > rpcDigiToken
MonitorElement * Res_Endcap23_Ring2_A
MonitorElement * clsBarrel
MonitorElement * ResDmin4
MonitorElement * ResLayer4_barrel
MonitorElement * CLS_Endcap_123_Ring3_C
std::pair< const_iterator, const_iterator > Range
MonitorElement * ResDplu1
edm::EDGetTokenT< edm::PSimHitContainer > simHitToken
MonitorElement * Res_Endcap23_Ring2_B
MonitorElement * xyvDplu4
MonitorElement * CLS_Endcap_23_Ring2_C
MonitorElement * Res_Endcap1_Ring2_C
MonitorElement * CLS_Endcap_23_Ring2_B
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.
MonitorElement * noiseCLSBarrel