CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
BaseEvtVtxGenerator.cc
Go to the documentation of this file.
1 
2 /*
3 */
4 
6 
9 
12 
15 
17 
20 
21 //#include "HepMC/GenEvent.h"
22 // #include "CLHEP/Vector/ThreeVector.h"
23 // #include "HepMC/SimpleVector.h"
24 
25 using namespace edm;
26 using namespace CLHEP;
27 //using namespace HepMC;
28 
29 
31 {
33  if ( ! rng.isAvailable()) {
34  throw cms::Exception("Configuration")
35  << "The BaseEvtVtxGenerator requires the RandomNumberGeneratorService\n"
36  "which is not present in the configuration file. \n"
37  "You must add the service\n"
38  "in the configuration file or remove the modules that require it.";
39  }
40 
41  sourceToken=consumes<edm::HepMCProduct>(pset.getParameter<edm::InputTag>("src"));
42  produces<edm::HepMCProduct>();
43 }
44 
46 {
47 }
48 
50 {
52  CLHEP::HepRandomEngine* engine = &rng->getEngine(evt.streamID());
53 
54  Handle<HepMCProduct> HepUnsmearedMCEvt ;
55 
56  evt.getByToken( sourceToken, HepUnsmearedMCEvt ) ;
57 
58  // Copy the HepMC::GenEvent
59  HepMC::GenEvent* genevt = new HepMC::GenEvent(*HepUnsmearedMCEvt->GetEvent());
60  std::unique_ptr<edm::HepMCProduct> HepMCEvt(new edm::HepMCProduct(genevt));
61  // generate new vertex & apply the shift
62  //
63  HepMCEvt->applyVtxGen( newVertex(engine) ) ;
64 
65  //HepMCEvt->LorentzBoost( 0., 142.e-6 );
66  HepMCEvt->boostToLab( GetInvLorentzBoost(), "vertex" );
67  HepMCEvt->boostToLab( GetInvLorentzBoost(), "momentum" );
68 
69  evt.put(std::move(HepMCEvt)) ;
70 
71  return ;
72 }
T getParameter(std::string const &) const
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:127
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:508
void produce(edm::Event &, const edm::EventSetup &) override
virtual CLHEP::HepRandomEngine & getEngine(StreamID const &)=0
Use this engine in event methods.
BaseEvtVtxGenerator(const edm::ParameterSet &)
bool isAvailable() const
Definition: Service.h:46
const HepMC::GenEvent * GetEvent() const
Definition: HepMCProduct.h:38
return(e1-e2)*(e1-e2)+dp *dp
HLT enums.
StreamID streamID() const
Definition: Event.h:86
def move(src, dest)
Definition: eostools.py:510