CMS 3D CMS Logo
RPCSynchronizer.cc
Go to the documentation of this file.
1 #include "Geometry/RPCGeometry/interface/RPCRoll.h"
2 #include "Geometry/RPCGeometry/interface/RPCRollSpecs.h"
3 #include "SimMuon/RPCDigitizer/src/RPCSynchronizer.h"
4 #include "Geometry/CommonTopologies/interface/RectangularStripTopology.h"
5 #include "Geometry/CommonTopologies/interface/TrapezoidalStripTopology.h"
6 
7 #include "FWCore/Framework/interface/Frameworkfwd.h"
8 #include "FWCore/Framework/interface/EventSetup.h"
9 #include "FWCore/Framework/interface/EDAnalyzer.h"
10 #include "FWCore/Framework/interface/Event.h"
11 #include "FWCore/ParameterSet/interface/ParameterSet.h"
12 
13 #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
14 #include "SimDataFormats/TrackingHit/interface/PSimHit.h"
15 #include "Geometry/RPCGeometry/interface/RPCGeometry.h"
16 #include "Geometry/Records/interface/MuonGeometryRecord.h"
17 #include "DataFormats/MuonDetId/interface/RPCDetId.h"
18 #include "SimMuon/RPCDigitizer/src/RPCSimSetUp.h"
19 #include "SimDataFormats/TrackingHit/interface/PSimHit.h"
20 
21 #include "CLHEP/Random/RandGaussQ.h"
22 
23 #include <cstring>
24 #include <iostream>
25 #include <fstream>
26 #include <string>
27 #include <vector>
28 #include <cstdlib>
29 #include <cmath>
30 
31 using namespace std;
32 
33 RPCSynchronizer::RPCSynchronizer(const edm::ParameterSet& config) {
34  resRPC = config.getParameter<double>("timeResolution");
35  timOff = config.getParameter<double>("timingRPCOffset");
36  dtimCs = config.getParameter<double>("deltatimeAdjacentStrip");
37  resEle = config.getParameter<double>("timeJitter");
38  sspeed = config.getParameter<double>("signalPropagationSpeed");
39  lbGate = config.getParameter<double>("linkGateWidth");
40  LHCGate = config.getParameter<double>("Gate");
41  cosmics = config.getParameter<bool>("cosmics");
42  irpc_timing_res = config.getParameter<double>("IRPC_time_resolution");
43  irpc_electronics_jitter = config.getParameter<double>("IRPC_electronics_jitter");
44  N_BX = config.getParameter<int>("BX_range");
45  //"magic" parameter for cosmics
46  cosmicPar = 37.62;
47 
48  double c = 299792458; // [m/s]
49  //light speed in [cm/ns]
50  cspeed = c * 1e+2 * 1e-9;
51  //signal propagation speed [cm/ns]
52  sspeed = sspeed * cspeed;
53 }
54 
55 RPCSynchronizer::~RPCSynchronizer() {}
56 
57 int RPCSynchronizer::getSimHitBx(const PSimHit* simhit, CLHEP::HepRandomEngine* engine) {
58  RPCSimSetUp* simsetup = this->getRPCSimSetUp();
59  const RPCGeometry* geometry = simsetup->getGeometry();
60  float timeref = simsetup->getTime(simhit->detUnitId());
61 
62  int bx = -999;
63  LocalPoint simHitPos = simhit->localPosition();
64  float tof = simhit->timeOfFlight();
65 
66  //automatic variable to prevent memory leak
67 
68  float rr_el = CLHEP::RandGaussQ::shoot(engine, 0., resEle);
69 
70  RPCDetId SimDetId(simhit->detUnitId());
71 
72  const RPCRoll* SimRoll = nullptr;
73 
74  for (TrackingGeometry::DetContainer::const_iterator it = geometry->dets().begin(); it != geometry->dets().end();
75  it++) {
76  if (dynamic_cast<const RPCChamber*>(*it) != nullptr) {
77  auto ch = dynamic_cast<const RPCChamber*>(*it);
78 
79  std::vector<const RPCRoll*> rollsRaf = (ch->rolls());
80  for (std::vector<const RPCRoll*>::iterator r = rollsRaf.begin(); r != rollsRaf.end(); ++r) {
81  if ((*r)->id() == SimDetId) {
82  SimRoll = &(*(*r));
83  break;
84  }
85  }
86  }
87  }
88 
89  if (SimRoll != nullptr) {
90  float distanceFromEdge = 0;
91  float half_stripL = 0.;
92 
93  if (SimRoll->id().region() == 0) {
94  const RectangularStripTopology* top_ = dynamic_cast<const RectangularStripTopology*>(&(SimRoll->topology()));
95  half_stripL = top_->stripLength() / 2;
96  distanceFromEdge = half_stripL + simHitPos.y();
97  } else {
98  const TrapezoidalStripTopology* top_ = dynamic_cast<const TrapezoidalStripTopology*>(&(SimRoll->topology()));
99  half_stripL = top_->stripLength() / 2;
100  distanceFromEdge = half_stripL - simHitPos.y();
101  }
102 
103  float prop_time = distanceFromEdge / sspeed;
104 
105  double rr_tim1 = CLHEP::RandGaussQ::shoot(engine, 0., resRPC);
106  double total_time = tof + prop_time + timOff + rr_tim1 + rr_el;
107 
108  // Bunch crossing assignment
109  double time_differ = 0.;
110 
111  if (cosmics) {
112  time_differ = (total_time - (timeref + ((half_stripL / sspeed) + timOff))) / cosmicPar;
113  } else if (!cosmics) {
114  time_differ = total_time - (timeref + (half_stripL / sspeed) + timOff);
115  }
116 
117  double inf_time = 0;
118  double sup_time = 0;
119 
120  for (int n = -N_BX; n <= N_BX; ++n) {
121  if (cosmics) {
122  inf_time = (-lbGate / 2 + n * LHCGate) / cosmicPar;
123  sup_time = (lbGate / 2 + n * LHCGate) / cosmicPar;
124  } else if (!cosmics) {
125  inf_time = -lbGate / 2 + n * LHCGate;
126  sup_time = lbGate / 2 + n * LHCGate;
127  }
128 
129  if (inf_time < time_differ && time_differ < sup_time) {
130  bx = n;
131  break;
132  }
133  }
134  }
135 
136  return bx;
137 }
138 
139 int RPCSynchronizer::getSimHitBxAndTimingForIRPC(const PSimHit* simhit, CLHEP::HepRandomEngine* engine) {
140  RPCSimSetUp* simsetup = this->getRPCSimSetUp();
141  const RPCGeometry* geometry = simsetup->getGeometry();
142  float timeref = simsetup->getTime(simhit->detUnitId());
143 
144  int bx = -999;
145  LocalPoint simHitPos = simhit->localPosition();
146  float tof = simhit->timeOfFlight();
147 
148  //automatic variable to prevent memory leak
149 
150  // float rr_el = CLHEP::RandGaussQ::shoot(engine, 0.,resEle);
151  float rr_el = CLHEP::RandGaussQ::shoot(engine, 0., irpc_electronics_jitter);
152 
153  RPCDetId SimDetId(simhit->detUnitId());
154 
155  const RPCRoll* SimRoll = nullptr;
156 
157  for (TrackingGeometry::DetContainer::const_iterator it = geometry->dets().begin(); it != geometry->dets().end();
158  it++) {
159  if (dynamic_cast<const RPCChamber*>(*it) != nullptr) {
160  auto ch = dynamic_cast<const RPCChamber*>(*it);
161 
162  std::vector<const RPCRoll*> rollsRaf = (ch->rolls());
163  for (std::vector<const RPCRoll*>::iterator r = rollsRaf.begin(); r != rollsRaf.end(); ++r) {
164  if ((*r)->id() == SimDetId) {
165  SimRoll = &(*(*r));
166  break;
167  }
168  }
169  }
170  }
171 
172  if (SimRoll != nullptr) {
173  float distanceFromEdge = 0;
174  float half_stripL = 0.;
175 
176  if (SimRoll->id().region() == 0) {
177  const RectangularStripTopology* top_ = dynamic_cast<const RectangularStripTopology*>(&(SimRoll->topology()));
178  half_stripL = top_->stripLength() / 2;
179  distanceFromEdge = half_stripL + simHitPos.y();
180  } else {
181  const TrapezoidalStripTopology* top_ = dynamic_cast<const TrapezoidalStripTopology*>(&(SimRoll->topology()));
182  half_stripL = top_->stripLength() / 2;
183  distanceFromEdge = half_stripL - simHitPos.y();
184  }
185 
186  float prop_time = distanceFromEdge / sspeed;
187 
188  // double rr_tim1 = CLHEP::RandGaussQ::shoot(engine, 0.,resRPC);
189  double rr_tim1 = CLHEP::RandGaussQ::shoot(engine, 0., irpc_timing_res);
190 
191  double total_time = tof + prop_time + timOff + rr_tim1 + rr_el;
192 
193  // Bunch crossing assignment
194  double time_differ = 0.;
195 
196  if (cosmics) {
197  time_differ = (total_time - (timeref + ((half_stripL / sspeed) + timOff))) / cosmicPar;
198  } else if (!cosmics) {
199  time_differ = total_time - (timeref + (half_stripL / sspeed) + timOff);
200  }
201 
202  double exact_total_time = tof + prop_time + timOff;
203  double exact_time_differ = 0.;
204 
205  if (cosmics) {
206  exact_time_differ = (exact_total_time - (timeref + ((half_stripL / sspeed) + timOff))) / cosmicPar;
207  } else if (!cosmics) {
208  exact_time_differ = exact_total_time - (timeref + (half_stripL / sspeed) + timOff);
209  }
210 
211  double inf_time = 0;
212  double sup_time = 0;
213 
214  for (int n = -N_BX; n <= N_BX; ++n) {
215  if (cosmics) {
216  inf_time = (-lbGate / 2 + n * LHCGate) / cosmicPar;
217  sup_time = (lbGate / 2 + n * LHCGate) / cosmicPar;
218  } else if (!cosmics) {
219  inf_time = -lbGate / 2 + n * LHCGate;
220  sup_time = lbGate / 2 + n * LHCGate;
221  }
222 
223  if (inf_time < time_differ && time_differ < sup_time) {
224  bx = n;
225  the_exact_time = exact_time_differ;
226  the_smeared_time = time_differ;
227 
228  //cout<<"Debug\t"<<inf_time<<'\t'<<sup_time<<endl;
230  // cout<<"Bingo\t"<<time_differ<<'\t'<<bx<<'\t'<<exact_time_differ<<'\t'<<exact_time_differ-time_differ<<'\t'<<exact_time_differ-bx*25.<<endl;
231  break;
232  }
233  }
234  }
235  return bx;
236 }
RPCSimSetUp::getTime
float getTime(uint32_t id)
Definition: RPCSimSetUp.cc:440
RPCRoll
Definition: RPCRoll.h:12
PSimHit::timeOfFlight
float timeOfFlight() const
Definition: PSimHit.h:73
TrapezoidalStripTopology::stripLength
float stripLength() const override
det heigth (strip length in the middle)
Definition: TrapezoidalStripTopology.h:63
dqmiodumpmetadata.n
n
Definition: dqmiodumpmetadata.py:28
RPCDetId
Definition: RPCDetId.h:16
geometry
Definition: geometry.py:1
PSimHitContainer.h
l1GtPatternGenerator_cfi.bx
bx
Definition: l1GtPatternGenerator_cfi.py:18
RPCSynchronizer::getSimHitBx
int getSimHitBx(const PSimHit *, CLHEP::HepRandomEngine *)
Definition: RPCSynchronizer.cc:57
RPCRoll.h
EDAnalyzer.h
RPCSynchronizer::~RPCSynchronizer
~RPCSynchronizer()
Definition: RPCSynchronizer.cc:55
RPCSimSetUp::getGeometry
const RPCGeometry * getGeometry()
Definition: RPCSimSetUp.h:49
RPCSynchronizer::RPCSynchronizer
RPCSynchronizer(const edm::ParameterSet &config)
Definition: RPCSynchronizer.cc:33
RectangularStripTopology.h
PSimHit::detUnitId
unsigned int detUnitId() const
Definition: PSimHit.h:97
config
Definition: config.py:1
PSimHit.h
RectangularStripTopology
Definition: RectangularStripTopology.h:11
TrapezoidalStripTopology
Definition: TrapezoidalStripTopology.h:21
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:52
Point3DBase< float, LocalTag >
RPCDetId.h
RectangularStripTopology::stripLength
float stripLength() const override
Definition: RectangularStripTopology.h:44
ecal_dqm_sourceclient-live_cfg.cosmics
cosmics
Definition: ecal_dqm_sourceclient-live_cfg.py:55
RPCRollSpecs.h
edm::ParameterSet
Definition: ParameterSet.h:47
Event.h
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
alignCSCRings.r
r
Definition: alignCSCRings.py:93
std
Definition: JetResolutionObject.h:76
TrapezoidalStripTopology.h
Frameworkfwd.h
EventSetup.h
RPCSimSetUp.h
RPCGeometry
Definition: RPCGeometry.h:20
RPCSynchronizer::getSimHitBxAndTimingForIRPC
int getSimHitBxAndTimingForIRPC(const PSimHit *, CLHEP::HepRandomEngine *)
Definition: RPCSynchronizer.cc:139
RPCSimSetUp
Definition: RPCSimSetUp.h:28
ParameterSet.h
PSimHit
Definition: PSimHit.h:15
MuonGeometryRecord.h
c
auto & c
Definition: CAHitNtupletGeneratorKernelsImpl.h:56
RPCGeometry.h
RPCSynchronizer.h
lumi::N_BX
const unsigned int N_BX
Definition: ConstantDef.h:5
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
Generated for CMSSW Reference Manual by doxygen 1.8.16