gen::PhotosppInterface Class Reference

#include <PhotosppInterface.h>

Inheritance diagram for gen::PhotosppInterface:

Public Member Functions
HepMC::GenEvent *	apply (HepMC::GenEvent *) override
void	avoidTauLeptonicDecays () override
void	configureOnlyFor (int) override
void	init () override
bool	isTauLeptonicDecay (HepMC::GenVertex *)
	PhotosppInterface (const edm::ParameterSet &pset)
void	setRandomEngine (CLHEP::HepRandomEngine *decayRandomEngine) override
const std::vector< std::string > &	specialSettings () override
void	statistics () override
	~PhotosppInterface () override
Public Member Functions inherited from gen::PhotosInterfaceBase
	PhotosInterfaceBase ()
	PhotosInterfaceBase (const edm::ParameterSet &)
virtual void	SetDecayRandomEngine (CLHEP::HepRandomEngine *decayRandomEngine)
virtual	~PhotosInterfaceBase ()

Static Public Member Functions
static double	flat ()

Private Attributes
bool	fAvoidTauLeptonicDecays
bool	fIsInitialized
int	fOnlyPDG
edm::ParameterSet *	fPSet

Static Private Attributes
static CLHEP::HepRandomEngine *	fRandomEngine = nullptr

Additional Inherited Members
Protected Attributes inherited from gen::PhotosInterfaceBase
std::vector< std::string >	fSpecialSettings

Detailed Description

Definition at line 19 of file PhotosppInterface.h.

Constructor & Destructor Documentation

PhotosppInterface()

PhotosppInterface::PhotosppInterface

(

const edm::ParameterSet &

pset

)

Definition at line 22 of file PhotosppInterface.cc.

References fPSet, gen::PhotosInterfaceBase::fSpecialSettings, edm::ParameterSet::getParameter(), muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

    : fOnlyPDG(-1), fAvoidTauLeptonicDecays(false), fIsInitialized(false), fPSet(nullptr) {
  // add ability to keep brem from hadronizer and only modify specific channels 10/27/2014
  bool UseHadronizerQEDBrem = false;
  fPSet = new ParameterSet(pset);
  std::vector<std::string> par = fPSet->getParameter<std::vector<std::string> >("parameterSets");
  for (unsigned int ip = 0; ip < par.size(); ++ip) {
    std::string curSet = par[ip];
    // Physics settings
    if (curSet == "UseHadronizerQEDBrem")
      UseHadronizerQEDBrem = true;
  }
  if (!UseHadronizerQEDBrem)
    fSpecialSettings.push_back("QED-brem-off:all");
}

~PhotosppInterface()

gen::PhotosppInterface::~PhotosppInterface ( )

inlineoverride

Definition at line 23 of file PhotosppInterface.h.

{}

Member Function Documentation

apply()

HepMC::GenEvent * PhotosppInterface::apply ( HepMC::GenEvent *  evt )

overridevirtual

Reimplemented from gen::PhotosInterfaceBase.

Definition at line 257 of file PhotosppInterface.cc.

References class-composition::children, fIsInitialized, flat(), GenParticle::GenParticle, LaserDQM_cfg::p1, and extraflags_cff::vtx.

Referenced by heavyIonTools.ConfigureHeavyIons::__call__(), coreTools.RunOnData::__call__(), trackTools.MakeAODTrackCandidates::__call__(), runJetUncertainties.RunJetUncertainties::__call__(), metTools.AddMETCollection::__call__(), heavyIonTools.ProductionDefaults::__call__(), cmsswVersionTools.PickRelValInputFiles::__call__(), coreTools.RemoveMCMatching::__call__(), trackTools.MakePATTrackCandidates::__call__(), trigTools.SwitchOnTrigger::__call__(), runMETCorrectionsAndUncertainties.RunMETCorrectionsAndUncertainties::__call__(), heavyIonTools.SelectionDefaults::__call__(), heavyIonTools.DisbaleMonteCarloDeps::__call__(), trigTools.SwitchOnTriggerStandAlone::__call__(), tauTools.AddTauCollection::__call__(), trackTools.MakeTrackCandidates::__call__(), trigTools.SwitchOnTriggerMatching::__call__(), trigTools.SwitchOnTriggerMatchingStandAlone::__call__(), trigTools.SwitchOnTriggerMatchEmbedding::__call__(), jetTools.AddJetCollection::__call__(), jetTools.SwitchJetCollection::__call__(), jetTools.UpdateJetCollection::__call__(), jetTools.AddJetID::__call__(), and jetTools.SetTagInfos::__call__().

                                                          {
  Photospp::Photos::setRandomGenerator(PhotosppInterface::flat);
  if (!fIsInitialized)
    return evt;
  int NPartBefore = evt->particles_size();
  Photospp::PhotosHepMCEvent PhotosEvt(evt);
  PhotosEvt.process();
  //Fix the vertices and barcodes based on Julia Yarba's solution from TauolaInterface
  for (HepMC::GenEvent::vertex_const_iterator vtx = evt->vertices_begin(); vtx != evt->vertices_end(); vtx++) {
    std::vector<int> BCodes;
    BCodes.clear();
    if (*vtx) {
      for (HepMC::GenVertex::particle_iterator pitr = (*vtx)->particles_begin(HepMC::children);
           pitr != (*vtx)->particles_end(HepMC::children);
           ++pitr) {
        if ((*pitr)->barcode() > 10000) {
          BCodes.push_back((*pitr)->barcode());
        }
      }
    }
    if (!BCodes.empty()) {
      for (size_t ibc = 0; ibc < BCodes.size(); ibc++) {
        HepMC::GenParticle* p1 = evt->barcode_to_particle(BCodes[ibc]);
        int nbc = p1->barcode() - 10000 + NPartBefore;
        p1->suggest_barcode(nbc);
      }
    }
  }
  return evt;
}

avoidTauLeptonicDecays()

void gen::PhotosppInterface::avoidTauLeptonicDecays ( )

inlineoverridevirtual

Implements gen::PhotosInterfaceBase.

Definition at line 29 of file PhotosppInterface.h.

References fAvoidTauLeptonicDecays.

                                           {
      fAvoidTauLeptonicDecays = true;
      return;
    }

configureOnlyFor()

void PhotosppInterface::configureOnlyFor ( int  ipdg )

overridevirtual

Implements gen::PhotosInterfaceBase.

Definition at line 42 of file PhotosppInterface.cc.

References fOnlyPDG, and gen::PhotosInterfaceBase::fSpecialSettings.

                                                 {
  fOnlyPDG = ipdg;
  fSpecialSettings.clear();
  return;
}

flat()

double PhotosppInterface::flat ( )

static

Definition at line 288 of file PhotosppInterface.cc.

References Exception, and fRandomEngine.

Referenced by apply().

                               {
  if (!fRandomEngine) {
    throw cms::Exception("LogicError")
        << "PhotosppInterface::flat: Attempt to generate random number when engine pointer is null\n"
        << "This might mean that the code was modified to generate a random number outside the\n"
        << "event and beginLuminosityBlock methods, which is not allowed.\n";
  }
  return fRandomEngine->flat();
}

init()

void PhotosppInterface::init ( void  )

overridevirtual

Implements gen::PhotosInterfaceBase.

Definition at line 48 of file PhotosppInterface.cc.

References looper::cfg, fAvoidTauLeptonicDecays, fIsInitialized, fOnlyPDG, PhotosppPythia8TauSettings_cfi::forceBremForDecay, fPSet, edm::ParameterSet::getParameter(), mps_fire::i, initialize(), PhotosppPythia8TauSettings_cfi::setMomentumConservationThreshold, AlCaHLTBitMon_QueryRunRegistry::string, PhotosppPythia8TauSettings_cfi::suppressAll, and gen::v.

                             {
  if (fIsInitialized)
    return;  // do init only once
  Photospp::Photos::initialize();
  Photospp::Photos::createHistoryEntries(true, 746);  // P-H-O
  std::vector<std::string> par = fPSet->getParameter<std::vector<std::string> >("parameterSets");
  for (unsigned int ip = 0; ip < par.size(); ++ip) {
    std::string curSet = par[ip];

    // Physics settings
    if (curSet == "maxWtInterference")
      Photospp::Photos::maxWtInterference(fPSet->getParameter<double>(curSet));
    if (curSet == "setInfraredCutOff")
      Photospp::Photos::setInfraredCutOff(fPSet->getParameter<double>(curSet));
    if (curSet == "setAlphaQED")
      Photospp::Photos::setAlphaQED(fPSet->getParameter<double>(curSet));
    if (curSet == "setInterference")
      Photospp::Photos::setInterference(fPSet->getParameter<bool>(curSet));
    if (curSet == "setDoubleBrem")
      Photospp::Photos::setDoubleBrem(fPSet->getParameter<bool>(curSet));
    if (curSet == "setQuatroBrem")
      Photospp::Photos::setQuatroBrem(fPSet->getParameter<bool>(curSet));
    if (curSet == "setExponentiation")
      Photospp::Photos::setExponentiation(fPSet->getParameter<bool>(curSet));
    if (curSet == "setCorrectionWtForW")
      Photospp::Photos::setCorrectionWtForW(fPSet->getParameter<bool>(curSet));
    if (curSet == "setMeCorrectionWtForScalar")
      Photospp::Photos::setMeCorrectionWtForScalar(fPSet->getParameter<bool>(curSet));
    if (curSet == "setMeCorrectionWtForW")
      Photospp::Photos::setMeCorrectionWtForW(fPSet->getParameter<bool>(curSet));
    if (curSet == "setMeCorrectionWtForZ")
      Photospp::Photos::setMeCorrectionWtForZ(fPSet->getParameter<bool>(curSet));
    if (curSet == "initializeKinematicCorrections")
      Photospp::Photos::initializeKinematicCorrections(fPSet->getParameter<int>(curSet));
    if (curSet == "forceMassFrom4Vector")
      Photospp::Photos::forceMassFrom4Vector(fPSet->getParameter<bool>(curSet));
    if (curSet == "forceMassFromEventRecord")
      Photospp::Photos::forceMassFromEventRecord(fPSet->getParameter<int>(curSet));
    if (curSet == "ignoreParticlesOfStatus")
      Photospp::Photos::ignoreParticlesOfStatus(fPSet->getParameter<int>(curSet));
    if (curSet == "deIgnoreParticlesOfStatus")
      Photospp::Photos::deIgnoreParticlesOfStatus(fPSet->getParameter<int>(curSet));
    if (curSet == "setMomentumConservationThreshold")
      Photospp::Photos::setMomentumConservationThreshold(fPSet->getParameter<double>(curSet));
    if (curSet == "suppressAll")
      if (fPSet->getParameter<bool>(curSet) == true)
        Photospp::Photos::suppressAll();
    if (curSet == "setPairEmission")
      Photospp::Photos::setPairEmission(fPSet->getParameter<bool>(curSet));
    if (curSet == "setPhotonEmission")
      Photospp::Photos::setPhotonEmission(fPSet->getParameter<bool>(curSet));
    if (curSet == "setStopAtCriticalError")
      Photospp::Photos::setStopAtCriticalError(fPSet->getParameter<bool>(curSet));
    if (curSet == "createHistoryEntries")
      Photospp::Photos::createHistoryEntries(fPSet->getParameter<bool>(curSet), 746);

    // Now setup more complicated radiation/mass supression and forcing.
    if (curSet == "suppressBremForBranch") {
      edm::ParameterSet cfg = fPSet->getParameter<edm::ParameterSet>(curSet);
      std::vector<std::string> v = cfg.getParameter<std::vector<std::string> >("parameterSets");
      for (unsigned int i = 0; i < v.size(); i++) {
        std::string vs = v[i];
        std::vector<int> vpar = cfg.getParameter<std::vector<int> >(vs);
        if (vpar.size() == 1)
          Photospp::Photos::suppressBremForBranch(0, vpar[0]);
        if (vpar.size() == 2)
          Photospp::Photos::suppressBremForBranch(0 /*vpar[0]*/, vpar[1]);
        if (vpar.size() == 3)
          Photospp::Photos::suppressBremForBranch(vpar[0], vpar[1], vpar[2]);
        if (vpar.size() == 4)
          Photospp::Photos::suppressBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3]);
        if (vpar.size() == 5)
          Photospp::Photos::suppressBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4]);
        if (vpar.size() == 6)
          Photospp::Photos::suppressBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5]);
        if (vpar.size() == 7)
          Photospp::Photos::suppressBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6]);
        if (vpar.size() == 8)
          Photospp::Photos::suppressBremForBranch(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7]);
        if (vpar.size() == 9)
          Photospp::Photos::suppressBremForBranch(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8]);
        if (vpar.size() == 10)
          Photospp::Photos::suppressBremForBranch(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8], vpar[9]);
      }
    }
    if (curSet == "suppressBremForDecay") {
      edm::ParameterSet cfg = fPSet->getParameter<edm::ParameterSet>(curSet);
      std::vector<std::string> v = cfg.getParameter<std::vector<std::string> >("parameterSets");
      for (unsigned int i = 0; i < v.size(); i++) {
        std::string vs = v[i];
        std::vector<int> vpar = cfg.getParameter<std::vector<int> >(vs);
        if (vpar.size() == 1)
          Photospp::Photos::suppressBremForDecay(0, vpar[0]);
        if (vpar.size() == 2)
          Photospp::Photos::suppressBremForDecay(0 /*vpar[0]*/, vpar[1]);
        if (vpar.size() == 3)
          Photospp::Photos::suppressBremForDecay(vpar[0], vpar[1], vpar[2]);
        if (vpar.size() == 4)
          Photospp::Photos::suppressBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3]);
        if (vpar.size() == 5)
          Photospp::Photos::suppressBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4]);
        if (vpar.size() == 6)
          Photospp::Photos::suppressBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5]);
        if (vpar.size() == 7)
          Photospp::Photos::suppressBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6]);
        if (vpar.size() == 8)
          Photospp::Photos::suppressBremForDecay(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7]);
        if (vpar.size() == 9)
          Photospp::Photos::suppressBremForDecay(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8]);
        if (vpar.size() == 10)
          Photospp::Photos::suppressBremForDecay(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8], vpar[9]);
      }
    }

    if (curSet == "forceBremForBranch") {
      edm::ParameterSet cfg = fPSet->getParameter<edm::ParameterSet>(curSet);
      std::vector<std::string> v = cfg.getParameter<std::vector<std::string> >("parameterSets");
      for (unsigned int i = 0; i < v.size(); i++) {
        std::string vs = v[i];
        std::vector<int> vpar = cfg.getParameter<std::vector<int> >(vs);
        if (vpar.size() == 1)
          Photospp::Photos::forceBremForBranch(0, vpar[0]);
        if (vpar.size() == 2)
          Photospp::Photos::forceBremForBranch(0 /*vpar[0]*/, vpar[1]);
        if (vpar.size() == 3)
          Photospp::Photos::forceBremForBranch(vpar[0], vpar[1], vpar[2]);
        if (vpar.size() == 4)
          Photospp::Photos::forceBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3]);
        if (vpar.size() == 5)
          Photospp::Photos::forceBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4]);
        if (vpar.size() == 6)
          Photospp::Photos::forceBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5]);
        if (vpar.size() == 7)
          Photospp::Photos::forceBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6]);
        if (vpar.size() == 8)
          Photospp::Photos::forceBremForBranch(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7]);
        if (vpar.size() == 9)
          Photospp::Photos::forceBremForBranch(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8]);
        if (vpar.size() == 10)
          Photospp::Photos::forceBremForBranch(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8], vpar[9]);
      }
    }
    if (curSet == "forceBremForDecay") {
      edm::ParameterSet cfg = fPSet->getParameter<edm::ParameterSet>(curSet);
      std::vector<std::string> v = cfg.getParameter<std::vector<std::string> >("parameterSets");
      for (unsigned int i = 0; i < v.size(); i++) {
        std::string vs = v[i];
        std::vector<int> vpar = cfg.getParameter<std::vector<int> >(vs);
        if (vpar.size() == 1)
          Photospp::Photos::forceBremForDecay(0, vpar[0]);
        if (vpar.size() == 2)
          Photospp::Photos::forceBremForDecay(0 /*vpar[0]*/, vpar[1]);
        if (vpar.size() == 3)
          Photospp::Photos::forceBremForDecay(vpar[0], vpar[1], vpar[2]);
        if (vpar.size() == 4)
          Photospp::Photos::forceBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3]);
        if (vpar.size() == 5)
          Photospp::Photos::forceBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4]);
        if (vpar.size() == 6)
          Photospp::Photos::forceBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5]);
        if (vpar.size() == 7)
          Photospp::Photos::forceBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6]);
        if (vpar.size() == 8)
          Photospp::Photos::forceBremForDecay(vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7]);
        if (vpar.size() == 9)
          Photospp::Photos::forceBremForDecay(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8]);
        if (vpar.size() == 10)
          Photospp::Photos::forceBremForDecay(
              vpar[0], vpar[1], vpar[2], vpar[3], vpar[4], vpar[5], vpar[6], vpar[7], vpar[8], vpar[9]);
      }
    }

    if (curSet == "forceMass") {
      edm::ParameterSet cfg = fPSet->getParameter<edm::ParameterSet>(curSet);
      std::vector<std::string> v = cfg.getParameter<std::vector<std::string> >("parameterSets");
      for (unsigned int i = 0; i < v.size(); i++) {
        std::string vs = v[i];
        std::vector<double> vpar = cfg.getParameter<std::vector<double> >(vs);
        if (vpar.size() == 2)
          Photospp::Photos::forceMass((int)vpar[0], vpar[1]);
      }
    }
  }
  // implement options set via c++
  if (fOnlyPDG >= 0) {
    Photospp::Photos::suppressAll();
    Photospp::Photos::forceBremForBranch(0, fOnlyPDG);
    Photospp::Photos::forceBremForBranch(0, -1 * fOnlyPDG);
  }
  if (fAvoidTauLeptonicDecays) {
    Photospp::Photos::suppressBremForDecay(3, 15, 16, 11, -12);
    Photospp::Photos::suppressBremForDecay(3, -15, -16, -11, 12);
    Photospp::Photos::suppressBremForDecay(3, 15, 16, 13, -14);
    Photospp::Photos::suppressBremForDecay(3, -15, -16, -13, -14);
  }
  Photospp::Photos::iniInfo();
  fIsInitialized = true;
  return;
}

isTauLeptonicDecay()

bool gen::PhotosppInterface::isTauLeptonicDecay

(

HepMC::GenVertex *

)

setRandomEngine()

void PhotosppInterface::setRandomEngine ( CLHEP::HepRandomEngine *  decayRandomEngine )

overridevirtual

Implements gen::PhotosInterfaceBase.

Definition at line 38 of file PhotosppInterface.cc.

References fRandomEngine.

                                                                               {
  fRandomEngine = decayRandomEngine;
}

specialSettings()

const std::vector<std::string>& gen::PhotosppInterface::specialSettings ( )

inlineoverridevirtual

Reimplemented from gen::PhotosInterfaceBase.

Definition at line 26 of file PhotosppInterface.h.

References gen::PhotosInterfaceBase::fSpecialSettings.

{ return fSpecialSettings; }

statistics()

void PhotosppInterface::statistics ( )

overridevirtual

Reimplemented from gen::PhotosInterfaceBase.

Definition at line 298 of file PhotosppInterface.cc.

{ Photospp::Photos::iniInfo(); }

Member Data Documentation

fAvoidTauLeptonicDecays

bool gen::PhotosppInterface::fAvoidTauLeptonicDecays

private

Definition at line 40 of file PhotosppInterface.h.

Referenced by avoidTauLeptonicDecays(), and init().

fIsInitialized

bool gen::PhotosppInterface::fIsInitialized

private

Definition at line 41 of file PhotosppInterface.h.

Referenced by apply(), and init().

fOnlyPDG

int gen::PhotosppInterface::fOnlyPDG

private

Definition at line 39 of file PhotosppInterface.h.

Referenced by configureOnlyFor(), and init().

fPSet

edm::ParameterSet* gen::PhotosppInterface::fPSet

private

Definition at line 42 of file PhotosppInterface.h.

Referenced by init(), and PhotosppInterface().

fRandomEngine

CLHEP::HepRandomEngine * PhotosppInterface::fRandomEngine = nullptr

staticprivate

Definition at line 44 of file PhotosppInterface.h.

Referenced by flat(), and setRandomEngine().