CMS 3D CMS Logo
Classes | Typedefs | Functions | Variables
externalGenerator.cc File Reference
#include "boost/program_options.hpp"
#include <atomic>
#include <csignal>
#include <iostream>
#include <string>
#include <thread>
#include <memory>
#include <filesystem>
#include "FWCore/TestProcessor/interface/TestProcessor.h"
#include "SimDataFormats/GeneratorProducts/interface/HepMCProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenRunInfoProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenLumiInfoHeader.h"
#include "SimDataFormats/GeneratorProducts/interface/GenLumiInfoProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"
#include "SimDataFormats/GeneratorProducts/interface/ExternalGeneratorEventInfo.h"
#include "SimDataFormats/GeneratorProducts/interface/ExternalGeneratorLumiInfo.h"
#include "FWCore/Services/interface/ExternalRandomNumberGeneratorService.h"
#include "FWCore/SharedMemory/interface/WriteBuffer.h"
#include "FWCore/SharedMemory/interface/ReadBuffer.h"
#include "FWCore/SharedMemory/interface/WorkerChannel.h"
#include "FWCore/SharedMemory/interface/ROOTSerializer.h"
#include "FWCore/SharedMemory/interface/ROOTDeserializer.h"
#include "FWCore/SharedMemory/interface/WorkerMonitorThread.h"
#include "FWCore/Utilities/interface/thread_safety_macros.h"

Go to the source code of this file.

Classes

class  Harness
 

Typedefs

template<typename T >
using Serializer = ROOTSerializer< T, WriteBuffer >
 

Functions

int main (int argc, char *argv[])
 

Variables

static char const *const kHelpCommandOpt = "help,h"
 
static char const *const kHelpOpt = "help"
 
static char const *const kMemoryNameCommandOpt = "memory-name,m"
 
static char const *const kMemoryNameOpt = "memory-name"
 
static char const *const kUniqueIDCommandOpt = "unique-id,i"
 
static char const *const kUniqueIDOpt = "unique-id"
 
static char const *const kVerboseCommandOpt = "verbose,v"
 
static char const *const kVerboseOpt = "verbose"
 
static std::string s_uniqueID
 

Typedef Documentation

◆ Serializer

template<typename T >
using Serializer = ROOTSerializer<T, WriteBuffer>

Definition at line 83 of file externalGenerator.cc.

Function Documentation

◆ main()

int main ( int  argc,
char *  argv[] 
)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk




















































































eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Definition at line 97 of file externalGenerator.cc.

References edm::shared_memory::WorkerChannel::accessLock(), dir2webdir::argc, cmsBatch::argv, cms::cuda::assert(), edm::BeginLuminosityBlock, edm::BeginRun, c, EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, CMS_SA_ALLOW, ctppsRawToDigi_cff::configuration, counter, gather_cfg::cout, edm::ServiceRegistry::createContaining(), submitPVResolutionJobs::desc, edm::shared_memory::ROOTDeserializer< T, READBUFFER >::deserialize(), edm::EndLuminosityBlock, edm::EndRun, relativeConstraints::error, edm::Event, cppFunctionSkipper::exception, edm::shared_memory::WorkerChannel::fromWorkerBufferInfo(), Harness::getBeginLumiValue(), Harness::getEndLumiValue(), Harness::getEndRunValue(), Harness::getEventValue(), edm::ExternalRandomNumberGeneratorService::getState(), edm::shared_memory::WorkerChannel::handleTransitions(), mps_fire::i, kHelpCommandOpt, kHelpOpt, kMemoryNameCommandOpt, kMemoryNameOpt, kUniqueIDCommandOpt, kUniqueIDOpt, kVerboseCommandOpt, kVerboseOpt, edm::ExternalRandomNumberGeneratorService::mySeed(), AlcaSiPixelAliHarvester0T_cff::options, AlCaHLTBitMon_ParallelJobs::p, castor_dqm_sourceclient_file_cfg::path, writedatasetfile::run, alignCSCRings::s, s_uniqueID, edm::shared_memory::ROOTSerializer< T, WRITEBUFFER >::serialize(), edm::shared_memory::WorkerMonitorThread::setAction(), edm::ExternalRandomNumberGeneratorService::setState(), edm::shared_memory::WorkerMonitorThread::setupSignalHandling(), edm::shared_memory::WorkerMonitorThread::startThread(), AlCaHLTBitMon_QueryRunRegistry::string, edm::shared_memory::WorkerChannel::toWorkerBufferInfo(), verbose, and edm::shared_memory::WorkerChannel::workerSetupDone().

97  {
98  std::string descString(argv[0]);
99  descString += " [--";
100  descString += kMemoryNameOpt;
101  descString += "] memory_name";
102  boost::program_options::options_description desc(descString);
103 
104  desc.add_options()(kHelpCommandOpt, "produce help message")(
105  kMemoryNameCommandOpt, boost::program_options::value<std::string>(), "memory name")(
106  kUniqueIDCommandOpt, boost::program_options::value<std::string>(), "unique id")(kVerboseCommandOpt,
107  "verbose output");
108 
109  boost::program_options::positional_options_description p;
110  p.add(kMemoryNameOpt, 1);
111  p.add(kUniqueIDOpt, 2);
112 
113  boost::program_options::options_description all_options("All Options");
114  all_options.add(desc);
115 
116  boost::program_options::variables_map vm;
117  try {
118  store(boost::program_options::command_line_parser(argc, argv).options(all_options).positional(p).run(), vm);
119  notify(vm);
120  } catch (boost::program_options::error const& iException) {
121  std::cout << argv[0] << ": Error while trying to process command line arguments:\n"
122  << iException.what() << "\nFor usage and an options list, please do 'cmsRun --help'.";
123  return 1;
124  }
125 
126  if (vm.count(kHelpOpt)) {
127  std::cout << desc << std::endl;
128  return 0;
129  }
130 
131  bool verbose = false;
132  if (vm.count(kVerboseOpt)) {
133  verbose = true;
134  }
135 
136  if (!vm.count(kMemoryNameOpt)) {
137  std::cout << " no argument given" << std::endl;
138  return 1;
139  }
140 
141  if (!vm.count(kUniqueIDOpt)) {
142  std::cout << " no second argument given" << std::endl;
143  return 1;
144  }
145 
146  using namespace std::string_literals;
147  using namespace std::filesystem;
148 
149  auto newDir = path("thread"s + vm[kUniqueIDOpt].as<std::string>());
150  create_directory(newDir);
151  current_path(newDir);
152 
153  WorkerMonitorThread monitorThread;
154 
155  monitorThread.startThread();
156 
157  std::string presentState = "setting up communicationChannel";
158 
159  CMS_SA_ALLOW try {
160  std::string const memoryName(vm[kMemoryNameOpt].as<std::string>());
161  std::string const uniqueID(vm[kUniqueIDOpt].as<std::string>());
162  s_uniqueID = uniqueID;
163  {
164  //This class is holding the lock
165  WorkerChannel communicationChannel(memoryName, uniqueID);
166 
167  presentState = "setting up read/write buffers";
168  WriteBuffer sm_buffer{memoryName, communicationChannel.fromWorkerBufferInfo()};
169  ReadBuffer sm_readbuffer{std::string("Rand") + memoryName, communicationChannel.toWorkerBufferInfo()};
170  int counter = 0;
171 
172  presentState = "setting up monitor thread";
173  //The lock must be released if there is a catastrophic signal
174  auto lockPtr = communicationChannel.accessLock();
175 
176  monitorThread.setAction([lockPtr]() {
177  if (lockPtr) {
178  std::cerr << s_uniqueID << " process: SIGNAL CAUGHT: unlock\n";
179  lockPtr->unlock();
180  }
181  });
182 
183  presentState = "setting up termination handler";
184  //be sure to unset the address of the shared lock before the lock goes away
185  s_sharedLock = lockPtr;
186  auto unsetLockPtr = [](void*) { s_sharedLock = nullptr; };
187  std::unique_ptr<decltype(s_sharedLock), decltype(unsetLockPtr)> sharedLockGuard{&s_sharedLock, unsetLockPtr};
188  std::atexit(atexit_handler);
189  auto releaseLock = []() {
190  if (s_sharedLock) {
191  std::cerr << s_uniqueID << " process: terminate called: unlock\n";
192  s_sharedLock->unlock();
193  s_sharedLock = nullptr;
194  //deactivate the abort signal
195 
196  struct sigaction act;
197  act.sa_sigaction = nullptr;
198  act.sa_flags = SA_SIGINFO;
199  sigemptyset(&act.sa_mask);
200  sigaction(SIGABRT, &act, nullptr);
201  std::abort();
202  }
203  };
204  std::set_terminate(releaseLock);
205 
206  presentState = "setting up serializers";
207  Serializer<ExternalGeneratorEventInfo> serializer(sm_buffer);
208  Serializer<ExternalGeneratorLumiInfo> bl_serializer(sm_buffer);
209  Serializer<GenLumiInfoProduct> el_serializer(sm_buffer);
210  Serializer<GenRunInfoProduct> er_serializer(sm_buffer);
211 
212  ROOTDeserializer<edm::RandomNumberGeneratorState, ReadBuffer> random_deserializer(sm_readbuffer);
213 
214  presentState = "reading configuration";
215  std::cerr << uniqueID << " process: initializing " << std::endl;
216  int nlines;
217  std::cin >> nlines;
218 
219  std::string configuration;
220  for (int i = 0; i < nlines; ++i) {
221  std::string c;
222  std::getline(std::cin, c);
223  if (verbose) {
224  std::cerr << c << "\n";
225  }
226  configuration += c + "\n";
227  }
228 
229  presentState = "setting up random number generator";
230  edm::ExternalRandomNumberGeneratorService* randomService = new edm::ExternalRandomNumberGeneratorService;
231  auto serviceToken =
232  edm::ServiceRegistry::createContaining(std::unique_ptr<edm::RandomNumberGenerator>(randomService));
233  Harness harness(configuration, serviceToken);
234 
235  //Some generator libraries override the signal handlers
236  monitorThread.setupSignalHandling();
237  std::set_terminate(releaseLock);
238 
239  if (verbose) {
240  std::cerr << uniqueID << " process: done initializing" << std::endl;
241  }
242  presentState = "finished initialization";
243  communicationChannel.workerSetupDone();
244 
245  presentState = "waiting for transition";
246  if (verbose)
247  std::cerr << uniqueID << " process: waiting " << counter << std::endl;
248  communicationChannel.handleTransitions([&](edm::Transition iTransition, unsigned long long iTransitionID) {
249  ++counter;
250  switch (iTransition) {
251  case edm::Transition::BeginRun: {
252  presentState = "beginRun transition";
253  if (verbose)
254  std::cerr << uniqueID << " process: start beginRun " << std::endl;
255  if (verbose)
256  std::cerr << uniqueID << " process: end beginRun " << std::endl;
257 
258  break;
259  }
260  case edm::Transition::BeginLuminosityBlock: {
261  presentState = "begin lumi";
262  if (verbose)
263  std::cerr << uniqueID << " process: start beginLumi " << std::endl;
264  auto randState = random_deserializer.deserialize();
265  presentState = "deserialized random state in begin lumi";
266  if (verbose)
267  std::cerr << uniqueID << " random " << randState.state_.size() << " " << randState.seed_ << std::endl;
268  randomService->setState(randState.state_, randState.seed_);
269  presentState = "processing begin lumi";
270  auto value = harness.getBeginLumiValue(iTransitionID);
271  value.randomState_.state_ = randomService->getState();
272  value.randomState_.seed_ = randomService->mySeed();
273 
274  presentState = "serialize lumi";
275  bl_serializer.serialize(value);
276  if (verbose)
277  std::cerr << uniqueID << " process: end beginLumi " << std::endl;
278  if (verbose)
279  std::cerr << uniqueID << " rand " << value.randomState_.state_.size() << " " << value.randomState_.seed_
280  << std::endl;
281  break;
282  }
283  case edm::Transition::Event: {
284  presentState = "begin event";
285  if (verbose)
286  std::cerr << uniqueID << " process: event " << counter << std::endl;
287  presentState = "deserialized random state in event";
288  auto randState = random_deserializer.deserialize();
289  randomService->setState(randState.state_, randState.seed_);
290  presentState = "processing event";
291  auto value = harness.getEventValue();
292  value.randomState_.state_ = randomService->getState();
293  value.randomState_.seed_ = randomService->mySeed();
294 
295  if (verbose)
296  std::cerr << uniqueID << " process: event " << counter << std::endl;
297 
298  presentState = "serialize event";
299  serializer.serialize(value);
300  if (verbose)
301  std::cerr << uniqueID << " process: "
302  << " " << counter << std::endl;
303  //usleep(10000000);
304  break;
305  }
306  case edm::Transition::EndLuminosityBlock: {
307  presentState = "begin end lumi";
308  if (verbose)
309  std::cerr << uniqueID << " process: start endLumi " << std::endl;
310  presentState = "processing end lumi";
311  auto value = harness.getEndLumiValue();
312 
313  presentState = "serialize end lumi";
314  el_serializer.serialize(value);
315  if (verbose)
316  std::cerr << uniqueID << " process: end endLumi " << std::endl;
317 
318  break;
319  }
320  case edm::Transition::EndRun: {
321  presentState = "begin end run";
322  if (verbose)
323  std::cerr << uniqueID << " process: start endRun " << std::endl;
324  presentState = "process end run";
325  auto value = harness.getEndRunValue();
326 
327  presentState = "serialize end run";
328  er_serializer.serialize(value);
329  if (verbose)
330  std::cerr << uniqueID << " process: end endRun " << std::endl;
331 
332  break;
333  }
334  default: {
335  assert(false);
336  }
337  }
338  presentState = "notifying and waiting after " + presentState;
339  if (verbose)
340  std::cerr << uniqueID << " process: notifying and waiting " << counter << std::endl;
341  });
342  }
343  } catch (std::exception const& iExcept) {
344  std::cerr << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n"
345  << s_uniqueID << " process: caught exception \n"
346  << iExcept.what() << "\n"
347  << " while " << presentState << "\n"
348  << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n";
349  return 1;
350  } catch (...) {
351  std::cerr << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n"
352  << s_uniqueID << " process: caught unknown exception\n while " << presentState << "\n"
353  << "%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%\n";
354  return 1;
355  }
356  return 0;
357 }
edm::shared_memory::WorkerChannel
Definition: WorkerChannel.h:35
CMS_SA_ALLOW
#define CMS_SA_ALLOW
mps_fire.i
i
Definition: mps_fire.py:428
castor_dqm_sourceclient_file_cfg.path
path
Definition: castor_dqm_sourceclient_file_cfg.py:37
verbose
bool verbose
Definition: beamSpotDipStandalone.cc:87
kHelpOpt
static char const *const kHelpOpt
Definition: externalGenerator.cc:36
kMemoryNameCommandOpt
static char const *const kMemoryNameCommandOpt
Definition: externalGenerator.cc:33
edm::shared_memory::WriteBuffer
Definition: WriteBuffer.h:37
edm::ExternalRandomNumberGeneratorService::setState
void setState(std::vector< unsigned long > const &, long seed)
Definition: ExternalRandomNumberGeneratorService.cc:24
edm::shared_memory::WorkerMonitorThread::setupSignalHandling
void setupSignalHandling()
Sets the unix signal handler which communicates with the thread.
Definition: WorkerMonitorThread.cc:102
kMemoryNameOpt
static char const *const kMemoryNameOpt
Definition: externalGenerator.cc:32
s_uniqueID
static std::string s_uniqueID
Definition: externalGenerator.cc:42
edm::shared_memory::WorkerMonitorThread::setAction
void setAction(std::function< void()> iFunc)
Definition: WorkerMonitorThread.h:48
cms::cuda::assert
assert(be >=bs)
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
kHelpCommandOpt
static char const *const kHelpCommandOpt
Definition: externalGenerator.cc:37
string_literals
edm::shared_memory::WorkerMonitorThread::startThread
void startThread()
Definition: WorkerMonitorThread.cc:78
EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0.cerr
cerr
Definition: EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0.py:8
edm::Transition
Transition
Definition: Transition.h:12
edm::ServiceRegistry::createContaining
static ServiceToken createContaining(std::unique_ptr< T > iService)
create a service token that holds the service defined by iService
Definition: ServiceRegistry.h:99
kVerboseOpt
static char const *const kVerboseOpt
Definition: externalGenerator.cc:38
edm::shared_memory::ReadBuffer
Definition: ReadBuffer.h:34
value
Definition: value.py:1
kVerboseCommandOpt
static char const *const kVerboseCommandOpt
Definition: externalGenerator.cc:39
ctppsRawToDigi_cff.configuration
configuration
Definition: ctppsRawToDigi_cff.py:12
kUniqueIDCommandOpt
static char const *const kUniqueIDCommandOpt
Definition: externalGenerator.cc:35
kUniqueIDOpt
static char const *const kUniqueIDOpt
Definition: externalGenerator.cc:34
edm::shared_memory::ROOTDeserializer
Definition: ROOTDeserializer.h:31
counter
static std::atomic< unsigned int > counter
Definition: SharedResourceNames.cc:17
gather_cfg.cout
cout
Definition: gather_cfg.py:144
edm::shared_memory::WorkerMonitorThread
Definition: WorkerMonitorThread.h:32
AlcaSiPixelAliHarvester0T_cff.options
options
Definition: AlcaSiPixelAliHarvester0T_cff.py:42
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
edm::shared_memory::ROOTSerializer
Definition: ROOTSerializer.h:31

Variable Documentation

◆ kHelpCommandOpt

char const* const kHelpCommandOpt = "help,h"
static

Definition at line 37 of file externalGenerator.cc.

Referenced by main().

◆ kHelpOpt

char const* const kHelpOpt = "help"
static

Definition at line 36 of file externalGenerator.cc.

Referenced by main().

◆ kMemoryNameCommandOpt

char const* const kMemoryNameCommandOpt = "memory-name,m"
static

Definition at line 33 of file externalGenerator.cc.

Referenced by main().

◆ kMemoryNameOpt

char const* const kMemoryNameOpt = "memory-name"
static

Definition at line 32 of file externalGenerator.cc.

Referenced by main().

◆ kUniqueIDCommandOpt

char const* const kUniqueIDCommandOpt = "unique-id,i"
static

Definition at line 35 of file externalGenerator.cc.

Referenced by main().

◆ kUniqueIDOpt

char const* const kUniqueIDOpt = "unique-id"
static

Definition at line 34 of file externalGenerator.cc.

Referenced by main().

◆ kVerboseCommandOpt

char const* const kVerboseCommandOpt = "verbose,v"
static

Definition at line 39 of file externalGenerator.cc.

Referenced by main().

◆ kVerboseOpt

char const* const kVerboseOpt = "verbose"
static

Definition at line 38 of file externalGenerator.cc.

Referenced by main().

◆ s_uniqueID

std::string s_uniqueID
static

Definition at line 42 of file externalGenerator.cc.

Referenced by main().

Generated for CMSSW Reference Manual by doxygen 1.8.14