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

#include <RPCSynchronizer.h>

Public Member Functions

double getExactTime () const
 
RPCSimSetUp * getRPCSimSetUp ()
 
int getSimHitBx (const PSimHit *, CLHEP::HepRandomEngine *)
 
int getSimHitBxAndTimingForIRPC (const PSimHit *, CLHEP::HepRandomEngine *)
 
double getSmearedTime () const
 
 RPCSynchronizer (const edm::ParameterSet &config)
 
void setRPCSimSetUp (RPCSimSetUp *simsetup)
 
 ~RPCSynchronizer ()
 

Private Attributes

double cosmicPar
 
bool cosmics
 
double cspeed
 
double dtimCs
 
double irpc_electronics_jitter
 
double irpc_timing_res
 
double lbGate
 
double lbGateNew
 
double LHCGate
 
int N_BX
 
double resEle
 
double resRPC
 
double sspeed
 
double the_exact_time
 
double the_smeared_time
 
RPCSimSetUp * theSimSetUp
 
double timOff
 

Detailed Description

Class for the RPC strip response simulation based on a parametrized model (ORCA-based)

Author
Raffaello Trentadue – INFN Bari

Definition at line 35 of file RPCSynchronizer.h.

Constructor & Destructor Documentation

◆ RPCSynchronizer()

RPCSynchronizer::RPCSynchronizer ( const edm::ParameterSet &  config)

Definition at line 32 of file RPCSynchronizer.cc.

References DummyCfis::c, ecal_dqm_sourceclient-live_cfg::cosmics, MillePedeFileConverter_cfg::e, and lumi::N_BX.

32  {
33  resRPC = config.getParameter<double>("timeResolution");
34  timOff = config.getParameter<double>("timingRPCOffset");
35  dtimCs = config.getParameter<double>("deltatimeAdjacentStrip");
36  resEle = config.getParameter<double>("timeJitter");
37  sspeed = config.getParameter<double>("signalPropagationSpeed");
38  lbGate = config.getParameter<double>("linkGateWidth");
39  LHCGate = config.getParameter<double>("Gate");
40  cosmics = config.getParameter<bool>("cosmics");
41  irpc_timing_res = config.getParameter<double>("IRPC_time_resolution");
42  irpc_electronics_jitter = config.getParameter<double>("IRPC_electronics_jitter");
43  N_BX = config.getParameter<int>("BX_range");
44  //"magic" parameter for cosmics
45  cosmicPar = 37.62;
46 
47  double c = 299792458; // [m/s]
48  //light speed in [cm/ns]
49  cspeed = c * 1e+2 * 1e-9;
50  //signal propagation speed [cm/ns]
51  sspeed = sspeed * cspeed;
52 }
config
Definition: config.py:1
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
RPCSynchronizer::irpc_electronics_jitter
double irpc_electronics_jitter
Definition: RPCSynchronizer.h:60

◆ ~RPCSynchronizer()

RPCSynchronizer::~RPCSynchronizer ( )

Definition at line 54 of file RPCSynchronizer.cc.

54 {}

Member Function Documentation

◆ getExactTime()

double RPCSynchronizer::getExactTime ( ) const
inline

Definition at line 44 of file RPCSynchronizer.h.

References the_exact_time.

44 { return the_exact_time; }

◆ getRPCSimSetUp()

RPCSimSetUp* RPCSynchronizer::getRPCSimSetUp ( )
inline

Definition at line 43 of file RPCSynchronizer.h.

References theSimSetUp.

43 { return theSimSetUp; }
RPCSynchronizer::theSimSetUp
RPCSimSetUp * theSimSetUp
Definition: RPCSynchronizer.h:63

◆ getSimHitBx()

int RPCSynchronizer::getSimHitBx ( const PSimHit *  simhit,
CLHEP::HepRandomEngine *  engine 
)

Definition at line 56 of file RPCSynchronizer.cc.

References nano_mu_digi_cff::bx, ecal_dqm_sourceclient-live_cfg::cosmics, PSimHit::detUnitId(), RPCSimSetUp::getGeometry(), RPCSimSetUp::getTime(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, PSimHit::localPosition(), dqmiodumpmetadata::n, lumi::N_BX, RectangularStripTopology::stripLength(), TrapezoidalStripTopology::stripLength(), PSimHit::timeOfFlight(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by RPCSimSimple::simulate(), RPCSimParam::simulate(), RPCSimAverageNoise::simulate(), RPCSimAverage::simulate(), RPCSimAverageNoiseEff::simulate(), RPCSimAverageNoiseEffCls::simulate(), and RPCSimAsymmetricCls::simulate().

56  {
57  RPCSimSetUp* simsetup = this->getRPCSimSetUp();
58  const RPCGeometry* geometry = simsetup->getGeometry();
59  float timeref = simsetup->getTime(simhit->detUnitId());
60 
61  int bx = -999;
62  LocalPoint simHitPos = simhit->localPosition();
63  float tof = simhit->timeOfFlight();
64 
65  //automatic variable to prevent memory leak
66 
67  float rr_el = CLHEP::RandGaussQ::shoot(engine, 0., resEle);
68 
69  RPCDetId SimDetId(simhit->detUnitId());
70 
71  const RPCRoll* SimRoll = nullptr;
72 
73  for (TrackingGeometry::DetContainer::const_iterator it = geometry->dets().begin(); it != geometry->dets().end();
74  it++) {
75  if (dynamic_cast<const RPCChamber*>(*it) != nullptr) {
76  auto ch = dynamic_cast<const RPCChamber*>(*it);
77 
78  std::vector<const RPCRoll*> rollsRaf = (ch->rolls());
79  for (std::vector<const RPCRoll*>::iterator r = rollsRaf.begin(); r != rollsRaf.end(); ++r) {
80  if ((*r)->id() == SimDetId) {
81  SimRoll = &(*(*r));
82  break;
83  }
84  }
85  }
86  }
87 
88  if (SimRoll != nullptr) {
89  float distanceFromEdge = 0;
90  float half_stripL = 0.;
91 
92  if (SimRoll->id().region() == 0) {
93  const RectangularStripTopology* top_ = dynamic_cast<const RectangularStripTopology*>(&(SimRoll->topology()));
94  half_stripL = top_->stripLength() / 2;
95  distanceFromEdge = half_stripL + simHitPos.y();
96  } else {
97  const TrapezoidalStripTopology* top_ = dynamic_cast<const TrapezoidalStripTopology*>(&(SimRoll->topology()));
98  half_stripL = top_->stripLength() / 2;
99  distanceFromEdge = half_stripL - simHitPos.y();
100  }
101 
102  float prop_time = distanceFromEdge / sspeed;
103 
104  double rr_tim1 = CLHEP::RandGaussQ::shoot(engine, 0., resRPC);
105  double total_time = tof + prop_time + timOff + rr_tim1 + rr_el;
106 
107  // Bunch crossing assignment
108  double time_differ = 0.;
109 
110  if (cosmics) {
111  time_differ = (total_time - (timeref + ((half_stripL / sspeed) + timOff))) / cosmicPar;
112  } else if (!cosmics) {
113  time_differ = total_time - (timeref + (half_stripL / sspeed) + timOff);
114  }
115 
116  double inf_time = 0;
117  double sup_time = 0;
118 
119  for (int n = -N_BX; n <= N_BX; ++n) {
120  if (cosmics) {
121  inf_time = (-lbGate / 2 + n * LHCGate) / cosmicPar;
122  sup_time = (lbGate / 2 + n * LHCGate) / cosmicPar;
123  } else if (!cosmics) {
124  inf_time = -lbGate / 2 + n * LHCGate;
125  sup_time = lbGate / 2 + n * LHCGate;
126  }
127 
128  if (inf_time < time_differ && time_differ < sup_time) {
129  bx = n;
130  break;
131  }
132  }
133  }
134 
135  return bx;
136 }
PSimHit::detUnitId
unsigned int detUnitId() const
Definition: PSimHit.h:99
RPCSimSetUp::getTime
float getTime(uint32_t id)
Definition: RPCSimSetUp.cc:444
nano_mu_digi_cff.bx
bx
Definition: nano_mu_digi_cff.py:41
RectangularStripTopology::stripLength
float stripLength() const override
Definition: RectangularStripTopology.h:44
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
RPCSimSetUp::getGeometry
const RPCGeometry * getGeometry()
Definition: RPCSimSetUp.h:48
RPCSynchronizer::getRPCSimSetUp
RPCSimSetUp * getRPCSimSetUp()
Definition: RPCSynchronizer.h:43
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:54
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
PSimHit::timeOfFlight
float timeOfFlight() const
Definition: PSimHit.h:75
TrapezoidalStripTopology::stripLength
float stripLength() const override
det heigth (strip length in the middle)
Definition: TrapezoidalStripTopology.h:63

◆ getSimHitBxAndTimingForIRPC()

int RPCSynchronizer::getSimHitBxAndTimingForIRPC ( const PSimHit *  simhit,
CLHEP::HepRandomEngine *  engine 
)

Definition at line 138 of file RPCSynchronizer.cc.

References nano_mu_digi_cff::bx, ecal_dqm_sourceclient-live_cfg::cosmics, PSimHit::detUnitId(), RPCSimSetUp::getGeometry(), RPCSimSetUp::getTime(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, PSimHit::localPosition(), dqmiodumpmetadata::n, lumi::N_BX, RectangularStripTopology::stripLength(), TrapezoidalStripTopology::stripLength(), PSimHit::timeOfFlight(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by RPCSimModelTiming::simulate().

138  {
139  RPCSimSetUp* simsetup = this->getRPCSimSetUp();
140  const RPCGeometry* geometry = simsetup->getGeometry();
141  float timeref = simsetup->getTime(simhit->detUnitId());
142 
143  int bx = -999;
144  LocalPoint simHitPos = simhit->localPosition();
145  float tof = simhit->timeOfFlight();
146 
147  //automatic variable to prevent memory leak
148 
149  // float rr_el = CLHEP::RandGaussQ::shoot(engine, 0.,resEle);
150  float rr_el = CLHEP::RandGaussQ::shoot(engine, 0., irpc_electronics_jitter);
151 
152  RPCDetId SimDetId(simhit->detUnitId());
153 
154  const RPCRoll* SimRoll = nullptr;
155 
156  for (TrackingGeometry::DetContainer::const_iterator it = geometry->dets().begin(); it != geometry->dets().end();
157  it++) {
158  if (dynamic_cast<const RPCChamber*>(*it) != nullptr) {
159  auto ch = dynamic_cast<const RPCChamber*>(*it);
160 
161  std::vector<const RPCRoll*> rollsRaf = (ch->rolls());
162  for (std::vector<const RPCRoll*>::iterator r = rollsRaf.begin(); r != rollsRaf.end(); ++r) {
163  if ((*r)->id() == SimDetId) {
164  SimRoll = &(*(*r));
165  break;
166  }
167  }
168  }
169  }
170 
171  if (SimRoll != nullptr) {
172  float distanceFromEdge = 0;
173  float half_stripL = 0.;
174 
175  if (SimRoll->id().region() == 0) {
176  const RectangularStripTopology* top_ = dynamic_cast<const RectangularStripTopology*>(&(SimRoll->topology()));
177  half_stripL = top_->stripLength() / 2;
178  distanceFromEdge = half_stripL + simHitPos.y();
179  } else {
180  const TrapezoidalStripTopology* top_ = dynamic_cast<const TrapezoidalStripTopology*>(&(SimRoll->topology()));
181  half_stripL = top_->stripLength() / 2;
182  distanceFromEdge = half_stripL - simHitPos.y();
183  }
184 
185  float prop_time = distanceFromEdge / sspeed;
186 
187  // double rr_tim1 = CLHEP::RandGaussQ::shoot(engine, 0.,resRPC);
188  double rr_tim1 = CLHEP::RandGaussQ::shoot(engine, 0., irpc_timing_res);
189 
190  double total_time = tof + prop_time + timOff + rr_tim1 + rr_el;
191 
192  // Bunch crossing assignment
193  double time_differ = 0.;
194 
195  if (cosmics) {
196  time_differ = (total_time - (timeref + ((half_stripL / sspeed) + timOff))) / cosmicPar;
197  } else if (!cosmics) {
198  time_differ = total_time - (timeref + (half_stripL / sspeed) + timOff);
199  }
200 
201  double exact_total_time = tof + prop_time + timOff;
202  double exact_time_differ = 0.;
203 
204  if (cosmics) {
205  exact_time_differ = (exact_total_time - (timeref + ((half_stripL / sspeed) + timOff))) / cosmicPar;
206  } else if (!cosmics) {
207  exact_time_differ = exact_total_time - (timeref + (half_stripL / sspeed) + timOff);
208  }
209 
210  double inf_time = 0;
211  double sup_time = 0;
212 
213  for (int n = -N_BX; n <= N_BX; ++n) {
214  if (cosmics) {
215  inf_time = (-lbGate / 2 + n * LHCGate) / cosmicPar;
216  sup_time = (lbGate / 2 + n * LHCGate) / cosmicPar;
217  } else if (!cosmics) {
218  inf_time = -lbGate / 2 + n * LHCGate;
219  sup_time = lbGate / 2 + n * LHCGate;
220  }
221 
222  if (inf_time < time_differ && time_differ < sup_time) {
223  bx = n;
224  the_exact_time = exact_time_differ;
225  the_smeared_time = time_differ;
226 
227  //cout<<"Debug\t"<<inf_time<<'\t'<<sup_time<<endl;
229  // cout<<"Bingo\t"<<time_differ<<'\t'<<bx<<'\t'<<exact_time_differ<<'\t'<<exact_time_differ-time_differ<<'\t'<<exact_time_differ-bx*25.<<endl;
230  break;
231  }
232  }
233  }
234  return bx;
235 }
PSimHit::detUnitId
unsigned int detUnitId() const
Definition: PSimHit.h:99
RPCSimSetUp::getTime
float getTime(uint32_t id)
Definition: RPCSimSetUp.cc:444
nano_mu_digi_cff.bx
bx
Definition: nano_mu_digi_cff.py:41
RectangularStripTopology::stripLength
float stripLength() const override
Definition: RectangularStripTopology.h:44
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
RPCSimSetUp::getGeometry
const RPCGeometry * getGeometry()
Definition: RPCSimSetUp.h:48
RPCSynchronizer::getRPCSimSetUp
RPCSimSetUp * getRPCSimSetUp()
Definition: RPCSynchronizer.h:43
RPCSynchronizer::irpc_electronics_jitter
double irpc_electronics_jitter
Definition: RPCSynchronizer.h:60
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:54
ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it
auto & it
Definition: splitVertices.h:48
PSimHit::timeOfFlight
float timeOfFlight() const
Definition: PSimHit.h:75
TrapezoidalStripTopology::stripLength
float stripLength() const override
det heigth (strip length in the middle)
Definition: TrapezoidalStripTopology.h:63

◆ getSmearedTime()

double RPCSynchronizer::getSmearedTime ( ) const
inline

Definition at line 45 of file RPCSynchronizer.h.

References the_smeared_time.

Referenced by RPCSimModelTiming::simulate().

45 { return the_smeared_time; }

◆ setRPCSimSetUp()

void RPCSynchronizer::setRPCSimSetUp ( RPCSimSetUp *  simsetup)
inline

Definition at line 42 of file RPCSynchronizer.h.

References theSimSetUp.

Referenced by RPCSimSimple::simulate(), RPCSimParam::simulate(), RPCSimAverageNoise::simulate(), RPCSimAverage::simulate(), RPCSimAverageNoiseEffCls::simulate(), RPCSimAverageNoiseEff::simulate(), RPCSimModelTiming::simulate(), and RPCSimAsymmetricCls::simulate().

42 { theSimSetUp = simsetup; }
RPCSynchronizer::theSimSetUp
RPCSimSetUp * theSimSetUp
Definition: RPCSynchronizer.h:63

Member Data Documentation

◆ cosmicPar

double RPCSynchronizer::cosmicPar
private

Definition at line 56 of file RPCSynchronizer.h.

◆ cosmics

bool RPCSynchronizer::cosmics
private

Definition at line 58 of file RPCSynchronizer.h.

◆ cspeed

double RPCSynchronizer::cspeed
private

Definition at line 53 of file RPCSynchronizer.h.

◆ dtimCs

double RPCSynchronizer::dtimCs
private

Definition at line 50 of file RPCSynchronizer.h.

◆ irpc_electronics_jitter

double RPCSynchronizer::irpc_electronics_jitter
private

Definition at line 60 of file RPCSynchronizer.h.

◆ irpc_timing_res

double RPCSynchronizer::irpc_timing_res
private

Definition at line 59 of file RPCSynchronizer.h.

◆ lbGate

double RPCSynchronizer::lbGate
private

Definition at line 54 of file RPCSynchronizer.h.

◆ lbGateNew

double RPCSynchronizer::lbGateNew
private

Definition at line 55 of file RPCSynchronizer.h.

◆ LHCGate

double RPCSynchronizer::LHCGate
private

Definition at line 57 of file RPCSynchronizer.h.

◆ N_BX

int RPCSynchronizer::N_BX
private

Definition at line 64 of file RPCSynchronizer.h.

◆ resEle

double RPCSynchronizer::resEle
private

Definition at line 51 of file RPCSynchronizer.h.

◆ resRPC

double RPCSynchronizer::resRPC
private

Definition at line 48 of file RPCSynchronizer.h.

◆ sspeed

double RPCSynchronizer::sspeed
private

Definition at line 52 of file RPCSynchronizer.h.

◆ the_exact_time

double RPCSynchronizer::the_exact_time
private

Definition at line 61 of file RPCSynchronizer.h.

Referenced by getExactTime().

◆ the_smeared_time

double RPCSynchronizer::the_smeared_time
private

Definition at line 62 of file RPCSynchronizer.h.

Referenced by getSmearedTime().

◆ theSimSetUp

RPCSimSetUp* RPCSynchronizer::theSimSetUp
private

Definition at line 63 of file RPCSynchronizer.h.

Referenced by getRPCSimSetUp(), and setRPCSimSetUp().

◆ timOff

double RPCSynchronizer::timOff
private

Definition at line 49 of file RPCSynchronizer.h.

Generated for CMSSW Reference Manual by doxygen 1.8.14