32 #include "CLHEP/Random/RandFlat.h"
33 #include "CLHEP/Random/RandPoissonQ.h"
54 std::cout <<
"Average Cluster Size = " << aveCls <<
" strips" << std::endl;
55 std::cout <<
"RPC Time Resolution = " << resRPC <<
" ns" << std::endl;
56 std::cout <<
"RPC Signal formation time = " << timOff <<
" ns" << std::endl;
57 std::cout <<
"RPC adjacent strip delay = " << dtimCs <<
" ns" << std::endl;
58 std::cout <<
"Electronic Jitter = " << resEle <<
" ns" << std::endl;
59 std::cout <<
"Signal propagation time = " << sspeed <<
" x c" << std::endl;
60 std::cout <<
"Link Board Gate Width = " << lbGate <<
" ns" << std::endl;
76 double rr_cl = CLHEP::RandFlat::shoot(engine);
77 if (0.0 <= posX && posX < 0.2) {
78 func = (clsMap[1])[(clsMap[1]).size() - 1] * (rr_cl);
79 sum_clsize = clsMap[1];
81 if (0.2 <= posX && posX < 0.4) {
82 func = (clsMap[2])[(clsMap[2]).size() - 1] * (rr_cl);
83 sum_clsize = clsMap[2];
85 if (0.4 <= posX && posX < 0.6) {
86 func = (clsMap[3])[(clsMap[3]).size() - 1] * (rr_cl);
87 sum_clsize = clsMap[3];
89 if (0.6 <= posX && posX < 0.8) {
90 func = (clsMap[4])[(clsMap[4]).size() - 1] * (rr_cl);
91 sum_clsize = clsMap[4];
93 if (0.8 <= posX && posX < 1.0) {
94 func = (clsMap[5])[(clsMap[5]).size() - 1] * (rr_cl);
95 sum_clsize = clsMap[5];
98 for (vector<double>::iterator iter = sum_clsize.begin(); iter != sum_clsize.end(); ++iter) {
100 if (func > (*iter)) {
102 }
else if (func < (*iter)) {
111 CLHEP::HepRandomEngine* engine) {
119 for (edm::PSimHitContainer::const_iterator _hit = rpcHits.begin(); _hit != rpcHits.end(); ++_hit) {
123 float posX = roll->
strip(_hit->localPosition()) - static_cast<int>(roll->
strip(_hit->localPosition()));
127 if (CLHEP::RandFlat::shoot(engine) <
aveEff) {
128 int centralStrip = topology.channel(entr) + 1;
129 int fstrip = centralStrip;
130 int lstrip = centralStrip;
134 int clsize = this->
getClSize(posX, engine);
136 std::vector<int> cls;
137 cls.push_back(centralStrip);
139 for (
int cl = 0;
cl < (clsize - 1) / 2;
cl++) {
140 if (centralStrip -
cl - 1 >= 1) {
141 fstrip = centralStrip -
cl - 1;
142 cls.push_back(fstrip);
144 if (centralStrip +
cl + 1 <= roll->
nstrips()) {
145 lstrip = centralStrip +
cl + 1;
146 cls.push_back(lstrip);
149 if (clsize % 2 == 0) {
154 if (lstrip < roll->nstrips()) {
156 cls.push_back(lstrip);
161 cls.push_back(fstrip);
167 for (std::vector<int>::iterator
i = cls.begin();
i != cls.end();
i++) {
169 std::pair<int, int> digi(*
i, time_hit);
181 unsigned int nstrips = roll->
nstrips();
185 if (rpcId.
region() == 0) {
190 area = striplength * (xmax -
xmin);
196 area = striplength * (xmax -
xmin);
198 for (
unsigned int j = 0;
j < vnoise.size(); ++
j) {
203 CLHEP::RandPoissonQ randPoissonQ(*engine, ave);
204 N_hits = randPoissonQ.fire();
207 int time_hit = (
static_cast<int>(CLHEP::RandFlat::shoot(engine, (
nbxing *
gate) / gate))) -
nbxing / 2;
208 std::pair<int, int> digi(
j + 1, time_hit);
float strip(const LocalPoint &lp) const
LocalPoint centreOfStrip(int strip) const
LocalPoint localPosition(float strip) const override
DetectorHitMap theDetectorHitMap
void setRPCSimSetUp(RPCSimSetUp *simsetup)
constexpr uint32_t rawId() const
get the raw id
RPCSimSetUp * getRPCSimSetUp()
uint32_t T const *__restrict__ uint32_t const *__restrict__ int32_t int Histo::index_type cudaStream_t Func __host__ __device__ V int Func func
std::map< int, std::vector< double > > clsMap
const RPCRollSpecs * specs() const
int getSimHitBx(const PSimHit *, CLHEP::HepRandomEngine *)
int getClSize(float posX, CLHEP::HepRandomEngine *)
float stripLength() const override
edm::DetSet< RPCDigiSimLink > RPCDigiSimLinks
Container::value_type value_type
std::set< std::pair< int, int > > strips
void simulateNoise(const RPCRoll *, CLHEP::HepRandomEngine *) override
const std::map< int, std::vector< double > > & getClsMap()
void simulate(const RPCRoll *roll, const edm::PSimHitContainer &rpcHits, CLHEP::HepRandomEngine *) override
RPCSimAverageNoise(const edm::ParameterSet &config)
const Topology & topology() const override
~RPCSimAverageNoise() override
T getParameter(std::string const &) const
RPCSynchronizer * _rpcSync
tuple config
parse the configuration file
const Topology & topology() const override
std::vector< PSimHit > PSimHitContainer
float stripLength() const override
det heigth (strip length in the middle)
RPCDigiSimLinks theRpcDigiSimLinks
LocalPoint localPosition(float strip) const override
std::vector< double > sum_clsize
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.