gen::PhotosppInterface Class Reference

#include <PhotosppInterface.h>

Inheritance diagram for gen::PhotosppInterface:

Public Member Functions
HepMC::GenEvent *	apply (HepMC::GenEvent *)
void	avoidTauLeptonicDecays ()
void	configureOnlyFor (int)
void	init ()
bool	isTauLeptonicDecay (HepMC::GenVertex *)
	PhotosppInterface (const edm::ParameterSet &pset)
void	setRandomEngine (CLHEP::HepRandomEngine *decayRandomEngine)
const std::vector< std::string > &	specialSettings ()
void	statistics ()
	~PhotosppInterface ()
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 = 0

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

Detailed Description

Definition at line 21 of file PhotosppInterface.h.

Constructor & Destructor Documentation

PhotosppInterface::PhotosppInterface

(

const edm::ParameterSet &

pset

)

Definition at line 22 of file PhotosppInterface.cc.

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

  : fOnlyPDG(-1),
    fAvoidTauLeptonicDecays(false),
    fIsInitialized(false),
    fPSet(0)
{
  // 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");
}

gen::PhotosppInterface::~PhotosppInterface ( )

inline

Definition at line 26 of file PhotosppInterface.h.

References init.

{}

Member Function Documentation

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

virtual

Reimplemented from gen::PhotosInterfaceBase.

Definition at line 182 of file PhotosppInterface.cc.

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

Referenced by heavyIonTools.ConfigureHeavyIons::__call__(), editorTools.UserCodeTool::__call__(), HiCoreTools.RestrictInputToAOD::__call__(), coreTools.RunOnData::__call__(), trackTools.MakeAODTrackCandidates::__call__(), runJetUncertainties.RunJetUncertainties::__call__(), metTools.AddMETCollection::__call__(), heavyIonTools.ProductionDefaults::__call__(), editorTools.ChangeSource::__call__(), HiCoreTools.RemoveMCMatching::__call__(), cmsswVersionTools.PickRelValInputFiles::__call__(), coreTools.RemoveMCMatching::__call__(), trackTools.MakePATTrackCandidates::__call__(), trigTools.SwitchOnTrigger::__call__(), runMETCorrectionsAndUncertainties.RunMETCorrectionsAndUncertainties::__call__(), heavyIonTools.SelectionDefaults::__call__(), HiCoreTools.RemoveAllPATObjectsBut::__call__(), heavyIonTools.DisbaleMonteCarloDeps::__call__(), HiCoreTools.RemoveSpecificPATObjects::__call__(), trigTools.SwitchOnTriggerStandAlone::__call__(), trackTools.MakeTrackCandidates::__call__(), tauTools.AddTauCollection::__call__(), trigTools.SwitchOnTriggerMatching::__call__(), HiCoreTools.RemoveCleaning::__call__(), HiCoreTools.AddCleaning::__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.size() > 0){
      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;
}

void gen::PhotosppInterface::avoidTauLeptonicDecays ( )

inlinevirtual

Implements gen::PhotosInterfaceBase.

Definition at line 32 of file PhotosppInterface.h.

References decayRandomEngine, and myMessageLogger_cff::statistics.

{ fAvoidTauLeptonicDecays=true; return; }

void PhotosppInterface::configureOnlyFor ( int  ipdg )

virtual

Implements gen::PhotosInterfaceBase.

Definition at line 42 of file PhotosppInterface.cc.

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

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

double PhotosppInterface::flat ( )

static

Definition at line 210 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();
}

void PhotosppInterface::init ( void  )

virtual

Implements gen::PhotosInterfaceBase.

Definition at line 48 of file PhotosppInterface.cc.

References looper::cfg, fAvoidTauLeptonicDecays, fIsInitialized, fOnlyPDG, fPSet, edm::ParameterSet::getParameter(), mps_fire::i, initialize(), AlCaHLTBitMon_QueryRunRegistry::string, 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;
}

bool gen::PhotosppInterface::isTauLeptonicDecay

(

HepMC::GenVertex *

)

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

virtual

Implements gen::PhotosInterfaceBase.

Definition at line 40 of file PhotosppInterface.cc.

References decayRandomEngine, and fRandomEngine.

{fRandomEngine=decayRandomEngine;}

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

inlinevirtual

Reimplemented from gen::PhotosInterfaceBase.

Definition at line 29 of file PhotosppInterface.h.

References apply().

{ return fSpecialSettings; }

void PhotosppInterface::statistics ( )

virtual

Reimplemented from gen::PhotosInterfaceBase.

Definition at line 220 of file PhotosppInterface.cc.

References DEFINE_EDM_PLUGIN.

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

Member Data Documentation

bool gen::PhotosppInterface::fAvoidTauLeptonicDecays

private

Definition at line 40 of file PhotosppInterface.h.

Referenced by init().

bool gen::PhotosppInterface::fIsInitialized

private

Definition at line 41 of file PhotosppInterface.h.

Referenced by apply(), and init().

int gen::PhotosppInterface::fOnlyPDG

private

Definition at line 39 of file PhotosppInterface.h.

Referenced by configureOnlyFor(), and init().

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

private

Definition at line 42 of file PhotosppInterface.h.

Referenced by init(), and PhotosppInterface().

CLHEP::HepRandomEngine * PhotosppInterface::fRandomEngine = 0

staticprivate

Definition at line 44 of file PhotosppInterface.h.

Referenced by flat(), and setRandomEngine().