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

#include <test/EmbeddingLHEProducer/plugins/EmbeddingLHEProducer.cc>

Inheritance diagram for EmbeddingLHEProducer:
edm::one::EDProducer< edm::BeginRunProducer, edm::EndRunProducer > edm::one::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

 EmbeddingLHEProducer (const edm::ParameterSet &)
 
 ~EmbeddingLHEProducer ()
 
- Public Member Functions inherited from edm::one::EDProducer< edm::BeginRunProducer, edm::EndRunProducer >
 EDProducer ()=default
 
- Public Member Functions inherited from edm::one::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
virtual ~EDProducerBase ()
 
- Public Member Functions inherited from edm::ProducerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
std::vector< edm::ProductResolverIndex > const & indiciesForPutProducts (BranchType iBranchType) const
 
 ProducerBase ()
 
void registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
 
std::function< void(BranchDescription const &)> registrationCallback () const
 used by the fwk to register list of products More...
 
void resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
 
virtual ~ProducerBase () noexcept(false)
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

void assign_4vector (TLorentzVector &Lepton, const pat::Muon *muon, std::string FSRmode)
 
virtual void beginJob () override
 
virtual void beginRunProduce (edm::Run &run, edm::EventSetup const &es) override
 
virtual void endJob () override
 
virtual void endRunProduce (edm::Run &, edm::EventSetup const &) override
 
void fill_lhe_from_mumu (TLorentzVector &positiveLepton, TLorentzVector &negativeLepton, lhef::HEPEUP &outlhe, CLHEP::HepRandomEngine *engine)
 
void fill_lhe_with_particle (lhef::HEPEUP &outlhe, TLorentzVector &particle, int pdgid, double spin, double ctau)
 
const reco::Candidate * find_original_muon (const reco::Candidate *muon)
 
LHERunInfoProduct::Header give_slha ()
 
void mirror (TLorentzVector &positiveLepton, TLorentzVector &negativeLepton)
 
virtual void produce (edm::Event &, const edm::EventSetup &) override
 
void rotate180 (TLorentzVector &positiveLepton, TLorentzVector &negativeLepton)
 
void transform_mumu_to_tautau (TLorentzVector &positiveLepton, TLorentzVector &negativeLepton)
 

Private Attributes

std::ofstream file
 
bool mirror_
 
edm::EDGetTokenT< edm::View< pat::Muon > > muonsCollection_
 
bool rotate180_
 
std::string studyFSRmode_
 
bool switchToMuonEmbedding_
 
const double tauMass_ = 1.77682
 
edm::EDGetTokenT< reco::VertexCollection > vertexCollection_
 
bool write_lheout
 

Additional Inherited Members

- Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase ModuleType
 
- Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList TypeLabelList
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 72 of file EmbeddingLHEProducer.cc.

Constructor & Destructor Documentation

EmbeddingLHEProducer::EmbeddingLHEProducer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 117 of file EmbeddingLHEProducer.cc.

References Exception, FrontierConditions_GlobalTag_cff::file, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), edm::Service< T >::isAvailable(), MillePedeFileConverter_cfg::out, AlCaHLTBitMon_QueryRunRegistry::string, and estimatePileup_makeJSON::trunc.

118 {
119  //register your products
120  produces<LHEEventProduct>();
121  produces<LHERunInfoProduct, edm::InRun>();
122  produces<math::XYZTLorentzVectorD>("vertexPosition");
123 
124  muonsCollection_ = consumes<edm::View<pat::Muon>>(iConfig.getParameter<edm::InputTag>("src"));
125  vertexCollection_ = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("vertices"));
126  switchToMuonEmbedding_ = iConfig.getParameter<bool>("switchToMuonEmbedding");
127  mirror_ = iConfig.getParameter<bool>("mirror");
128  rotate180_ = iConfig.getParameter<bool>("rotate180");
129  studyFSRmode_ = iConfig.getUntrackedParameter<std::string>("studyFSRmode","");
130 
131  write_lheout=false;
132  std::string lhe_ouputfile = iConfig.getUntrackedParameter<std::string>("lhe_outputfilename","");
133  if (lhe_ouputfile !=""){
134  write_lheout=true;
135  file.open(lhe_ouputfile, std::fstream::out | std::fstream::trunc);
136  }
137 
138 
139  edm::Service<edm::RandomNumberGenerator> rng;
140  if ( ! rng.isAvailable()) {
141  throw cms::Exception("Configuration")
142  << "The EmbeddingLHEProducer requires the RandomNumberGeneratorService\n"
143  "which is not present in the configuration file. \n"
144  "You must add the service\n"
145  "in the configuration file or remove the modules that require it.";
146  }
147 
148 
149 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:46
estimatePileup_makeJSON.trunc
trunc
Definition: estimatePileup_makeJSON.py:214
MillePedeFileConverter_cfg.out
out
Definition: MillePedeFileConverter_cfg.py:31
EmbeddingLHEProducer::muonsCollection_
edm::EDGetTokenT< edm::View< pat::Muon > > muonsCollection_
Definition: EmbeddingLHEProducer.cc:99
EmbeddingLHEProducer::vertexCollection_
edm::EDGetTokenT< reco::VertexCollection > vertexCollection_
Definition: EmbeddingLHEProducer.cc:100
EmbeddingLHEProducer::~EmbeddingLHEProducer ( )

Definition at line 152 of file EmbeddingLHEProducer.cc.

153 {
154 }

Member Function Documentation

void EmbeddingLHEProducer::assign_4vector ( TLorentzVector &  Lepton,
const pat::Muon *  muon,
std::string  FSRmode 
)
private

Definition at line 391 of file EmbeddingLHEProducer.cc.

References pat::PATObject< ObjectType >::genParticle(), reco::Candidate::p4(), and reco::LeafCandidate::p4().

392 {
393  if("afterFSR" == FSRmode && muon->genParticle() != 0)
394  {
395  const reco::GenParticle* afterFSRMuon = muon->genParticle();
396  Lepton.SetPxPyPzE(afterFSRMuon->p4().px(),afterFSRMuon->p4().py(),afterFSRMuon->p4().pz(), afterFSRMuon->p4().e());
397  }
398  else if ("beforeFSR" == FSRmode && muon->genParticle() != 0)
399  {
400  const reco::Candidate* beforeFSRMuon = find_original_muon(muon->genParticle());
401  Lepton.SetPxPyPzE(beforeFSRMuon->p4().px(),beforeFSRMuon->p4().py(),beforeFSRMuon->p4().pz(), beforeFSRMuon->p4().e());
402  }
403  else
404  {
405  Lepton.SetPxPyPzE(muon->p4().px(),muon->p4().py(),muon->p4().pz(), muon->p4().e());
406  }
407  return;
408 }
pat::PATObject::genParticle
const reco::GenParticle * genParticle(size_t idx=0) const
Definition: PATObject.h:244
reco::Candidate::p4
virtual const LorentzVector & p4() const =0
four-momentum Lorentz vector
EmbeddingLHEProducer::find_original_muon
const reco::Candidate * find_original_muon(const reco::Candidate *muon)
Definition: EmbeddingLHEProducer.cc:410
Lepton
Definition: Lepton.py:1
reco::LeafCandidate::p4
virtual const LorentzVector & p4() const final
four-momentum Lorentz vector
Definition: LeafCandidate.h:99
void EmbeddingLHEProducer::beginJob ( void  )
overrideprivatevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 221 of file EmbeddingLHEProducer.cc.

222 {
223 }
void EmbeddingLHEProducer::beginRunProduce ( edm::Run &  run,
edm::EventSetup const &  es 
)
overrideprivatevirtual

Definition at line 233 of file EmbeddingLHEProducer.cc.

References popcon2dropbox::copy(), FrontierConditions_GlobalTag_cff::file, lhef::HEPRUP::IDWTUP, lhef::HEPRUP::LPRUP, eostools::move(), edm::Run::put(), lhef::HEPRUP::resize(), lhef::HEPRUP::XERRUP, lhef::HEPRUP::XMAXUP, and lhef::HEPRUP::XSECUP.

234 {
235  // fill HEPRUP common block and store in edm::Run
236  lhef::HEPRUP heprup;
237 
238  // set number of processes: 1 for Z to tau tau
239  heprup.resize(1);
240 
241  //Process independent information
242 
243  //beam particles ID (two protons)
244  //heprup.IDBMUP.first = 2212;
245  //heprup.IDBMUP.second = 2212;
246 
247  //beam particles energies (both 6.5 GeV)
248  //heprup.EBMUP.first = 6500.;
249  //heprup.EBMUP.second = 6500.;
250 
251  //take default pdf group for both beamparticles
252  //heprup.PDFGUP.first = -1;
253  //heprup.PDFGUP.second = -1;
254 
255  //take certan pdf set ID (same as in officially produced DYJets LHE files)
256  //heprup.PDFSUP.first = -1;
257  //heprup.PDFSUP.second = -1;
258 
259  //master switch for event weight iterpretation (same as in officially produced DYJets LHE files)
260  heprup.IDWTUP = 3;
261 
262  //Information for first process (Z to tau tau), for now only placeholder:
263  heprup.XSECUP[0] = 1.;
264  heprup.XERRUP[0] = 0;
265  heprup.XMAXUP[0] = 1;
266  heprup.LPRUP[0]= 1;
267 
268 
269  std::unique_ptr<LHERunInfoProduct> runInfo(new LHERunInfoProduct(heprup));
270  runInfo->addHeader(give_slha());
271 
272 
273  if (write_lheout)std::copy(runInfo->begin(), runInfo->end(),std::ostream_iterator<std::string>(file));
274  run.put(std::move(runInfo));
275 
276 
277 }
lhef::HEPRUP::resize
void resize(int nrup)
Definition: LesHouches.h:52
EmbeddingLHEProducer::give_slha
LHERunInfoProduct::Header give_slha()
Definition: EmbeddingLHEProducer.cc:472
lhef::HEPRUP::XERRUP
std::vector< double > XERRUP
Definition: LesHouches.h:114
lhef::HEPRUP::XMAXUP
std::vector< double > XMAXUP
Definition: LesHouches.h:119
edm::Run::put
void put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Run.h:111
lhef::HEPRUP::XSECUP
std::vector< double > XSECUP
Definition: LesHouches.h:108
eostools.move
def move(src, dest)
Definition: eostools.py:510
lhef::HEPRUP::LPRUP
std::vector< int > LPRUP
Definition: LesHouches.h:124
void EmbeddingLHEProducer::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 227 of file EmbeddingLHEProducer.cc.

227  {
228 }
void EmbeddingLHEProducer::endRunProduce ( edm::Run &  run,
edm::EventSetup const &  es 
)
overrideprivatevirtual

Definition at line 281 of file EmbeddingLHEProducer.cc.

References LHERunInfoProduct::endOfFile(), and FrontierConditions_GlobalTag_cff::file.

282 {
283  if (write_lheout) {
284  file << LHERunInfoProduct::endOfFile();
285  file.close();
286  }
287 }
LHERunInfoProduct::endOfFile
static const std::string & endOfFile()
Definition: LHERunInfoProduct.cc:163
void EmbeddingLHEProducer::fill_lhe_from_mumu ( TLorentzVector &  positiveLepton,
TLorentzVector &  negativeLepton,
lhef::HEPEUP &  outlhe,
CLHEP::HepRandomEngine *  engine 
)
private

Definition at line 290 of file EmbeddingLHEProducer.cc.

References Gflash::Z.

291 {
292 
293  TLorentzVector Z = positiveLepton + negativeLepton;
294  int leptonPDGID = switchToMuonEmbedding_ ? 13 : 15;
295 
296  // double tau_ctau = 0.00871100; //cm
297  double tau_ctau0 = 8.71100e-02; // mm (for Pythia)
298  double tau_ctau_p = tau_ctau0 * CLHEP::RandExponential::shoot(engine); // return -std::log(HepRandom::getTheEngine()->flat());
299  // replaces tau = process[iNow].tau0() * rndmPtr->exp(); from pythia8212/src/ProcessContainer.cc which is not initialized for ProcessLevel:all = off mode (no beam particle mode)
300  double tau_ctau_n = tau_ctau0 * CLHEP::RandExponential::shoot(engine);
301  //std::cout<<"tau_ctau P: "<<tau_ctau_p<<" tau_ctau N: "<<tau_ctau_n<<std::endl;
302 
303  fill_lhe_with_particle(outlhe, Z,23,9.0, 0);
304  fill_lhe_with_particle(outlhe, positiveLepton,-leptonPDGID,1.0, tau_ctau_p);
305  fill_lhe_with_particle(outlhe, negativeLepton,leptonPDGID,-1.0, tau_ctau_n);
306 
307  return;
308 }
Gflash::Z
const double Z[kNumberCalorimeter]
Definition: GflashNameSpace.h:58
EmbeddingLHEProducer::fill_lhe_with_particle
void fill_lhe_with_particle(lhef::HEPEUP &outlhe, TLorentzVector &particle, int pdgid, double spin, double ctau)
Definition: EmbeddingLHEProducer.cc:310
void EmbeddingLHEProducer::fill_lhe_with_particle ( lhef::HEPEUP &  outlhe,
TLorentzVector &  particle,
int  pdgid,
double  spin,
double  ctau 
)
private

Definition at line 310 of file EmbeddingLHEProducer.cc.

References funct::abs(), lhef::HEPEUP::ICOLUP, lhef::HEPEUP::IDUP, lhef::HEPEUP::ISTUP, lhef::HEPEUP::MOTHUP, lhef::HEPEUP::NUP, BPhysicsValidation_cfi::pdgid, lhef::HEPEUP::PUP, lhef::HEPEUP::resize(), spin(), lhef::HEPEUP::SPINUP, and lhef::HEPEUP::VTIMUP.

311 {
312  // Pay attention to different index conventions:
313  // 'particleindex' follows usual C++ index conventions starting at 0 for a list.
314  // 'motherindex' follows the LHE index conventions: 0 is for 'not defined', so the listing starts at 1.
315  // That means: LHE index 1 == C++ index 0.
316  int particleindex = outlhe.NUP;
317  outlhe.resize(outlhe.NUP+1);
318 
319  outlhe.PUP[particleindex][0] = particle.Px();
320  outlhe.PUP[particleindex][1] = particle.Py();
321  outlhe.PUP[particleindex][2] = particle.Pz();
322  outlhe.PUP[particleindex][3] = particle.E();
323  outlhe.PUP[particleindex][4] = particle.M();
324  outlhe.IDUP[particleindex] = pdgid;
325  outlhe.SPINUP[particleindex] = spin;
326  outlhe.VTIMUP[particleindex] = ctau;
327 
328  outlhe.ICOLUP[particleindex].first = 0;
329  outlhe.ICOLUP[particleindex].second = 0;
330 
331 
332  if (std::abs(pdgid) == 23){
333  outlhe.MOTHUP[particleindex].first = 0; // No Mother
334  outlhe.MOTHUP[particleindex].second = 0;
335  outlhe.ISTUP[particleindex] = 2; // status
336 
337  }
338 
339  if (std::abs(pdgid) == 15){
340  outlhe.MOTHUP[particleindex].first = 1; // Mother is the Z (first partile)
341  outlhe.MOTHUP[particleindex].second = 1; // Mother is the Z (first partile)
342 
343  outlhe.ISTUP[particleindex] = 1;//status
344 
345  }
346 
347 
348 
349  return;
350 }
lhef::HEPEUP::VTIMUP
std::vector< double > VTIMUP
Definition: LesHouches.h:254
lhef::HEPEUP::MOTHUP
std::vector< std::pair< int, int > > MOTHUP
Definition: LesHouches.h:236
lhef::HEPEUP::resize
void resize(int nup)
Definition: LesHouches.h:161
lhef::HEPEUP::PUP
std::vector< FiveVector > PUP
Definition: LesHouches.h:248
lhef::HEPEUP::SPINUP
std::vector< double > SPINUP
Definition: LesHouches.h:261
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
BPhysicsValidation_cfi.pdgid
pdgid
Definition: BPhysicsValidation_cfi.py:7
lhef::HEPEUP::ISTUP
std::vector< int > ISTUP
Definition: LesHouches.h:230
lhef::HEPEUP::IDUP
std::vector< int > IDUP
Definition: LesHouches.h:225
lhef::HEPEUP::ICOLUP
std::vector< std::pair< int, int > > ICOLUP
Definition: LesHouches.h:242
spin
float spin(float ph)
Definition: ValidHitPairFilter.cc:44
void EmbeddingLHEProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 569 of file EmbeddingLHEProducer.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

569  {
570  //The following says we do not know what parameters are allowed so do no validation
571  // Please change this to state exactly what you do use, even if it is no parameters
572  edm::ParameterSetDescription desc;
573  desc.setUnknown();
574  descriptions.addDefault(desc);
575 }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99
const reco::Candidate * EmbeddingLHEProducer::find_original_muon ( const reco::Candidate *  muon)
private

Definition at line 410 of file EmbeddingLHEProducer.cc.

References reco::Candidate::mother(), metsig::muon, and reco::Candidate::pdgId().

411 {
412  if(muon->mother(0) == 0) return muon;
413  if(muon->pdgId() == muon->mother(0)->pdgId()) return find_original_muon(muon->mother(0));
414  else return muon;
415 }
reco::Candidate::mother
virtual const Candidate * mother(size_type i=0) const =0
return pointer to mother
reco::Candidate::pdgId
virtual int pdgId() const =0
PDG identifier.
EmbeddingLHEProducer::find_original_muon
const reco::Candidate * find_original_muon(const reco::Candidate *muon)
Definition: EmbeddingLHEProducer.cc:410
LHERunInfoProduct::Header EmbeddingLHEProducer::give_slha ( )
private

Definition at line 472 of file EmbeddingLHEProducer.cc.

References LHERunInfoProduct::Header::addLine().

472  {
473  LHERunInfoProduct::Header slhah("slha");
474 
475  slhah.addLine("######################################################################\n");
476  slhah.addLine("## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####\n");
477  slhah.addLine("######################################################################\n");
478  slhah.addLine("## ##\n");
479  slhah.addLine("## Width set on Auto will be computed following the information ##\n");
480  slhah.addLine("## present in the decay.py files of the model. ##\n");
481  slhah.addLine("## See arXiv:1402.1178 for more details. ##\n");
482  slhah.addLine("## ##\n");
483  slhah.addLine("######################################################################\n");
484  slhah.addLine("\n");
485  slhah.addLine("###################################\n");
486  slhah.addLine("## INFORMATION FOR MASS\n");
487  slhah.addLine("###################################\n");
488  slhah.addLine("Block mass \n");
489  slhah.addLine(" 6 1.730000e+02 # MT \n");
490  slhah.addLine(" 15 1.777000e+00 # MTA \n");
491  slhah.addLine(" 23 9.118800e+01 # MZ \n");
492  slhah.addLine(" 25 1.250000e+02 # MH \n");
493  slhah.addLine("## Dependent parameters, given by model restrictions.\n");
494  slhah.addLine("## Those values should be edited following the \n");
495  slhah.addLine("## analytical expression. MG5 ignores those values \n");
496  slhah.addLine("## but they are important for interfacing the output of MG5\n");
497  slhah.addLine("## to external program such as Pythia.\n");
498  slhah.addLine(" 1 0.000000 # d : 0.0 \n");
499  slhah.addLine(" 2 0.000000 # u : 0.0 \n");
500  slhah.addLine(" 3 0.000000 # s : 0.0 \n");
501  slhah.addLine(" 4 0.000000 # c : 0.0 \n");
502  slhah.addLine(" 5 0.000000 # b : 0.0 \n");
503  slhah.addLine(" 11 0.000000 # e- : 0.0 \n");
504  slhah.addLine(" 12 0.000000 # ve : 0.0 \n");
505  slhah.addLine(" 13 0.000000 # mu- : 0.0 \n");
506  slhah.addLine(" 14 0.000000 # vm : 0.0 \n");
507  slhah.addLine(" 16 0.000000 # vt : 0.0 \n");
508  slhah.addLine(" 21 0.000000 # g : 0.0 \n");
509  slhah.addLine(" 22 0.000000 # a : 0.0 \n");
510  slhah.addLine(" 24 80.419002 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2))) \n");
511  slhah.addLine("\n");
512  slhah.addLine("###################################\n");
513  slhah.addLine("## INFORMATION FOR SMINPUTS\n");
514  slhah.addLine("###################################\n");
515  slhah.addLine("Block sminputs \n");
516  slhah.addLine(" 1 1.325070e+02 # aEWM1 \n");
517  slhah.addLine(" 2 1.166390e-05 # Gf \n");
518  slhah.addLine(" 3 1.180000e-01 # aS \n");
519  slhah.addLine("\n");
520  slhah.addLine("###################################\n");
521  slhah.addLine("## INFORMATION FOR WOLFENSTEIN\n");
522  slhah.addLine("###################################\n");
523  slhah.addLine("Block wolfenstein \n");
524  slhah.addLine(" 1 2.253000e-01 # lamWS \n");
525  slhah.addLine(" 2 8.080000e-01 # AWS \n");
526  slhah.addLine(" 3 1.320000e-01 # rhoWS \n");
527  slhah.addLine(" 4 3.410000e-01 # etaWS \n");
528  slhah.addLine("\n");
529  slhah.addLine("###################################\n");
530  slhah.addLine("## INFORMATION FOR YUKAWA\n");
531  slhah.addLine("###################################\n");
532  slhah.addLine("Block yukawa \n");
533  slhah.addLine(" 6 1.730000e+02 # ymt \n");
534  slhah.addLine(" 15 1.777000e+00 # ymtau \n");
535  slhah.addLine("\n");
536  slhah.addLine("###################################\n");
537  slhah.addLine("## INFORMATION FOR DECAY\n");
538  slhah.addLine("###################################\n");
539  slhah.addLine("DECAY 6 1.491500e+00 # WT \n");
540  slhah.addLine("DECAY 15 2.270000e-12 # WTau \n");
541  slhah.addLine("DECAY 23 2.441404e+00 # WZ \n");
542  slhah.addLine("DECAY 24 2.047600e+00 # WW \n");
543  slhah.addLine("DECAY 25 6.382339e-03 # WH \n");
544  slhah.addLine("## Dependent parameters, given by model restrictions.\n");
545  slhah.addLine("## Those values should be edited following the \n");
546  slhah.addLine("## analytical expression. MG5 ignores those values \n");
547  slhah.addLine("## but they are important for interfacing the output of MG5\n");
548  slhah.addLine("## to external program such as Pythia.\n");
549  slhah.addLine("DECAY 1 0.000000 # d : 0.0 \n");
550  slhah.addLine("DECAY 2 0.000000 # u : 0.0 \n");
551  slhah.addLine("DECAY 3 0.000000 # s : 0.0 \n");
552  slhah.addLine("DECAY 4 0.000000 # c : 0.0 \n");
553  slhah.addLine("DECAY 5 0.000000 # b : 0.0 \n");
554  slhah.addLine("DECAY 11 0.000000 # e- : 0.0 \n");
555  slhah.addLine("DECAY 12 0.000000 # ve : 0.0 \n");
556  slhah.addLine("DECAY 13 0.000000 # mu- : 0.0 \n");
557  slhah.addLine("DECAY 14 0.000000 # vm : 0.0 \n");
558  slhah.addLine("DECAY 16 0.000000 # vt : 0.0 \n");
559  slhah.addLine("DECAY 21 0.000000 # g : 0.0 \n");
560  slhah.addLine("DECAY 22 0.000000 # a : 0.0\n");
561 
562  return slhah;
563 }
void EmbeddingLHEProducer::mirror ( TLorentzVector &  positiveLepton,
TLorentzVector &  negativeLepton 
)
private

Definition at line 445 of file EmbeddingLHEProducer.cc.

References Gflash::Z.

446 {
447  if(!mirror_) return;
448  edm::LogInfo("TauEmbedding")<< "Applying mirroring" ;
449  TLorentzVector Z = positiveLepton + negativeLepton;
450 
451  edm::LogInfo("TauEmbedding") << "MuMinus before. Pt: " << negativeLepton.Pt() << " Eta: " << negativeLepton.Eta() << " Phi: " << negativeLepton.Phi() << " Mass: " << negativeLepton.M() ;
452 
453  TVector3 Z3 = Z.Vect();
454  TVector3 positiveLepton3 = positiveLepton.Vect();
455  TVector3 negativeLepton3 = negativeLepton.Vect();
456 
457  TVector3 beam(0.,0.,1.);
458  TVector3 perpToZandBeam = Z3.Cross(beam).Unit();
459 
460  positiveLepton3 -= 2*positiveLepton3.Dot(perpToZandBeam)*perpToZandBeam;
461  negativeLepton3 -= 2*negativeLepton3.Dot(perpToZandBeam)*perpToZandBeam;
462 
463  positiveLepton.SetVect(positiveLepton3);
464  negativeLepton.SetVect(negativeLepton3);
465 
466  edm::LogInfo("TauEmbedding") << "MuMinus after. Pt: " << negativeLepton.Pt() << " Eta: " << negativeLepton.Eta() << " Phi: " << negativeLepton.Phi() << " Mass: " << negativeLepton.M() ;
467 
468  return;
469 }
Gflash::Z
const double Z[kNumberCalorimeter]
Definition: GflashNameSpace.h:58
void EmbeddingLHEProducer::produce ( edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivatevirtual

Definition at line 163 of file EmbeddingLHEProducer.cc.

References lhef::HEPEUP::AQCDUP, lhef::HEPEUP::AQEDUP, edm::View< T >::begin(), popcon2dropbox::copy(), edm::View< T >::end(), FrontierConditions_GlobalTag_cff::file, edm::Event::getByToken(), edm::RandomNumberGenerator::getEngine(), lhef::HEPEUP::IDPRUP, EmbeddingLHEProducer_cfi::mirror, eostools::move(), metsig::muon, edm::Event::put(), EmbeddingLHEProducer_cfi::rotate180, lhef::HEPEUP::SCALUP, edm::Event::streamID(), and lhef::HEPEUP::XWGTUP.

164 {
165  using namespace edm;
166 
167  edm::Service<edm::RandomNumberGenerator> rng;
168  CLHEP::HepRandomEngine* engine = &rng->getEngine(iEvent.streamID());
169 
170 
171  edm::Handle< edm::View<pat::Muon> > muonHandle;
172  iEvent.getByToken(muonsCollection_, muonHandle);
173  edm::View<pat::Muon> coll_muons = *muonHandle;
174 
175  Handle<std::vector<reco::Vertex>> coll_vertices;
176  iEvent.getByToken(vertexCollection_ , coll_vertices);
177 
178  TLorentzVector positiveLepton, negativeLepton;
179  bool mu_plus_found = false;
180  bool mu_minus_found = false;
181  lhef::HEPEUP hepeup;
182  hepeup.IDPRUP = 0;
183  hepeup.XWGTUP = 1;
184  hepeup.SCALUP = -1;
185  hepeup.AQEDUP = -1;
186  hepeup.AQCDUP = -1;
187  // Assuming Pt-Order
188  for (edm::View<pat::Muon>::const_iterator muon= coll_muons.begin(); muon!= coll_muons.end(); ++muon)
189  {
190  if (muon->charge() == 1 && !mu_plus_found)
191  {
192  assign_4vector(positiveLepton, &(*muon), studyFSRmode_);
193  mu_plus_found = true;
194  }
195  else if (muon->charge() == -1 && !mu_minus_found)
196  {
197  assign_4vector(negativeLepton, &(*muon), studyFSRmode_);
198  mu_minus_found = true;
199  }
200  else if (mu_minus_found && mu_plus_found) break;
201  }
202  mirror(positiveLepton,negativeLepton); // if no mirror, function does nothing.
203  rotate180(positiveLepton,negativeLepton); // if no rotate180, function does nothing
204  transform_mumu_to_tautau(positiveLepton,negativeLepton); // if MuonEmbedding, function does nothing.
205  fill_lhe_from_mumu(positiveLepton,negativeLepton,hepeup,engine);
206 
207  double originalXWGTUP_ = 1.;
208  std::unique_ptr<LHEEventProduct> product( new LHEEventProduct(hepeup,originalXWGTUP_) );
209 
210  if (write_lheout) std::copy(product->begin(), product->end(), std::ostream_iterator<std::string>(file));
211 
212  iEvent.put(std::move(product));
213  // Saving vertex position
214  std::unique_ptr<math::XYZTLorentzVectorD> vertex_position (new math::XYZTLorentzVectorD(coll_vertices->at(0).x(),coll_vertices->at(0).y(),coll_vertices->at(0).z(),0.0));
215  iEvent.put(std::move(vertex_position), "vertexPosition");
216 
217 }
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:122
EmbeddingLHEProducer::mirror
void mirror(TLorentzVector &positiveLepton, TLorentzVector &negativeLepton)
Definition: EmbeddingLHEProducer.cc:445
math::XYZTLorentzVectorD
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:14
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
EmbeddingLHEProducer::transform_mumu_to_tautau
void transform_mumu_to_tautau(TLorentzVector &positiveLepton, TLorentzVector &negativeLepton)
Definition: EmbeddingLHEProducer.cc:355
edm::RandomNumberGenerator::getEngine
virtual CLHEP::HepRandomEngine & getEngine(StreamID const &)=0
Use this engine in event methods.
edm::View::begin
const_iterator begin() const
EmbeddingLHEProducer::assign_4vector
void assign_4vector(TLorentzVector &Lepton, const pat::Muon *muon, std::string FSRmode)
Definition: EmbeddingLHEProducer.cc:391
lhef::HEPEUP::AQCDUP
double AQCDUP
Definition: LesHouches.h:220
EmbeddingLHEProducer::rotate180
void rotate180(TLorentzVector &positiveLepton, TLorentzVector &negativeLepton)
Definition: EmbeddingLHEProducer.cc:417
EmbeddingLHEProducer::muonsCollection_
edm::EDGetTokenT< edm::View< pat::Muon > > muonsCollection_
Definition: EmbeddingLHEProducer.cc:99
edm
HLT enums.
Definition: AlignableModifier.h:17
edm::View::const_iterator
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
Definition: View.h:81
edm::Event::streamID
StreamID streamID() const
Definition: Event.h:81
EmbeddingLHEProducer::vertexCollection_
edm::EDGetTokenT< reco::VertexCollection > vertexCollection_
Definition: EmbeddingLHEProducer.cc:100
lhef::HEPEUP::AQEDUP
double AQEDUP
Definition: LesHouches.h:215
edm::View::end
const_iterator end() const
EmbeddingLHEProducer::fill_lhe_from_mumu
void fill_lhe_from_mumu(TLorentzVector &positiveLepton, TLorentzVector &negativeLepton, lhef::HEPEUP &outlhe, CLHEP::HepRandomEngine *engine)
Definition: EmbeddingLHEProducer.cc:290
lhef::HEPEUP::XWGTUP
double XWGTUP
Definition: LesHouches.h:196
eostools.move
def move(src, dest)
Definition: eostools.py:510
lhef::HEPEUP::SCALUP
double SCALUP
Definition: LesHouches.h:210
void EmbeddingLHEProducer::rotate180 ( TLorentzVector &  positiveLepton,
TLorentzVector &  negativeLepton 
)
private

Definition at line 417 of file EmbeddingLHEProducer.cc.

References Gflash::Z.

418 {
419  if (!rotate180_) return;
420  edm::LogInfo("TauEmbedding") << "Applying 180<C2><B0> rotation" ;
421  // By construction, the 3-momenta of mu-, mu+ and Z are in one plane.
422  // That means, one vector for perpendicular projection can be used for both leptons.
423  TLorentzVector Z = positiveLepton + negativeLepton;
424 
425  edm::LogInfo("TauEmbedding") << "MuMinus before. Pt: " << negativeLepton.Pt() << " Eta: " << negativeLepton.Eta() << " Phi: " << negativeLepton.Phi() << " Mass: " << negativeLepton.M();
426 
427  TVector3 Z3 = Z.Vect();
428  TVector3 positiveLepton3 = positiveLepton.Vect();
429  TVector3 negativeLepton3 = negativeLepton.Vect();
430 
431  TVector3 p3_perp = positiveLepton3 - positiveLepton3.Dot(Z3)/Z3.Dot(Z3)*Z3;
432  p3_perp = p3_perp.Unit();
433 
434  positiveLepton3 -= 2*positiveLepton3.Dot(p3_perp)*p3_perp;
435  negativeLepton3 -= 2*negativeLepton3.Dot(p3_perp)*p3_perp;
436 
437  positiveLepton.SetVect(positiveLepton3);
438  negativeLepton.SetVect(negativeLepton3);
439 
440  edm::LogInfo("TauEmbedding") << "MuMinus after. Pt: " << negativeLepton.Pt() << " Eta: " << negativeLepton.Eta() << " Phi: " << negativeLepton.Phi() << " Mass: " << negativeLepton.M();
441 
442  return;
443 }
Gflash::Z
const double Z[kNumberCalorimeter]
Definition: GflashNameSpace.h:58
void EmbeddingLHEProducer::transform_mumu_to_tautau ( TLorentzVector &  positiveLepton,
TLorentzVector &  negativeLepton 
)
private

Definition at line 355 of file EmbeddingLHEProducer.cc.

References Scenarios_cff::scale, mathSSE::sqrt(), and Gflash::Z.

356 {
357  // No corrections applied for muon embedding
358  if (switchToMuonEmbedding_) return;
359 
360  TLorentzVector Z = positiveLepton + negativeLepton;
361 
362  TVector3 boost_from_Z_to_LAB = Z.BoostVector();
363  TVector3 boost_from_LAB_to_Z = -Z.BoostVector();
364 
365  // Boosting the two leptons to Z restframe, then both are back to back. This means, same 3-momentum squared
366  positiveLepton.Boost(boost_from_LAB_to_Z);
367  negativeLepton.Boost(boost_from_LAB_to_Z);
368 
369  // Energy of tau = 0.5*Z-mass
370  double tau_mass_squared = tauMass_*tauMass_;
371  double tau_energy_squared = 0.25*Z.M2();
372  double tau_3momentum_squared = tau_energy_squared - tau_mass_squared;
373  if (tau_3momentum_squared < 0)
374  {
375  edm::LogWarning("TauEmbedding") << "3-Momentum squared is negative";
376  return;
377  }
378 
379  //Computing scale, applying it on the 3-momenta and building new 4 momenta of the taus
380  double scale = std::sqrt(tau_3momentum_squared/positiveLepton.Vect().Mag2());
381  positiveLepton.SetPxPyPzE(scale*positiveLepton.Px(),scale*positiveLepton.Py(),scale*positiveLepton.Pz(),std::sqrt(tau_energy_squared));
382  negativeLepton.SetPxPyPzE(scale*negativeLepton.Px(),scale*negativeLepton.Py(),scale*negativeLepton.Pz(),std::sqrt(tau_energy_squared));
383 
384  //Boosting the new taus back to LAB frame
385  positiveLepton.Boost(boost_from_Z_to_LAB);
386  negativeLepton.Boost(boost_from_Z_to_LAB);
387 
388  return;
389 }
Gflash::Z
const double Z[kNumberCalorimeter]
Definition: GflashNameSpace.h:58
Scenarios_cff.scale
scale
Definition: Scenarios_cff.py:2190
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18

Member Data Documentation

std::ofstream EmbeddingLHEProducer::file
private

Definition at line 105 of file EmbeddingLHEProducer.cc.

bool EmbeddingLHEProducer::mirror_
private

Definition at line 102 of file EmbeddingLHEProducer.cc.

edm::EDGetTokenT<edm::View<pat::Muon> > EmbeddingLHEProducer::muonsCollection_
private

Definition at line 99 of file EmbeddingLHEProducer.cc.

bool EmbeddingLHEProducer::rotate180_
private

Definition at line 102 of file EmbeddingLHEProducer.cc.

std::string EmbeddingLHEProducer::studyFSRmode_
private

Definition at line 111 of file EmbeddingLHEProducer.cc.

bool EmbeddingLHEProducer::switchToMuonEmbedding_
private

Definition at line 101 of file EmbeddingLHEProducer.cc.

const double EmbeddingLHEProducer::tauMass_ = 1.77682
private

Definition at line 103 of file EmbeddingLHEProducer.cc.

edm::EDGetTokenT<reco::VertexCollection> EmbeddingLHEProducer::vertexCollection_
private

Definition at line 100 of file EmbeddingLHEProducer.cc.

bool EmbeddingLHEProducer::write_lheout
private

Definition at line 106 of file EmbeddingLHEProducer.cc.

Generated for CMSSW Reference Manual by doxygen 1.8.11