ResonanceDecayFilterHook Class Reference

#include <ResonanceDecayFilterHook.h>

Inheritance diagram for ResonanceDecayFilterHook:

Public Member Functions
bool	canVetoResonanceDecays () override
bool	checkVetoResonanceDecays (const Pythia8::Event &process)
bool	doVetoResonanceDecays (Pythia8::Event &process) override
bool	initAfterBeams () override
	ResonanceDecayFilterHook ()

Private Attributes
bool	allNuAsEquivalent_
std::vector< int >	daughters_
bool	eMuAsEquivalent_
bool	eMuTauAsEquivalent_
bool	exclusive_
bool	filter_
std::set< int >	mothers_
std::map< int, int >	observedDaughters_
std::map< int, int >	requestedDaughters_
bool	udscAsEquivalent_
bool	udscbAsEquivalent_
bool	wzAsEquivalent_

Detailed Description

Definition at line 5 of file ResonanceDecayFilterHook.h.

Constructor & Destructor Documentation

ResonanceDecayFilterHook::ResonanceDecayFilterHook ( )

inline

Definition at line 10 of file ResonanceDecayFilterHook.h.

References initAfterBeams().

{}

Member Function Documentation

bool ResonanceDecayFilterHook::canVetoResonanceDecays ( )

inlineoverride

Definition at line 15 of file ResonanceDecayFilterHook.h.

{ return true; }

bool ResonanceDecayFilterHook::checkVetoResonanceDecays

(

const Pythia8::Event &

process

)

Definition at line 54 of file ResonanceDecayFilterHook.cc.

References funct::abs(), mps_fire::i, and AlCaHLTBitMon_ParallelJobs::p.

Referenced by doVetoResonanceDecays().

                                                                            {
  
  if (!filter_) return false;
    
  observedDaughters_.clear();
  
  //count decay products
  for (int i=0; i<process.size(); ++i) {
    const Particle &p = process[i];
    
    int did = std::abs(p.id());
    
    if ( did == 13 && (eMuAsEquivalent_ || eMuTauAsEquivalent_)) {
      did = 11;
    }
    if ( did == 15 && eMuTauAsEquivalent_) {
      did = 11;
    }
    if ( (did == 14 || did == 16) && allNuAsEquivalent_) {
      did = 12;
    }   
    if ( (did == 2 || did==3 || did==4 ) && udscAsEquivalent_) {
      did = 1;
    }
    if ( (did == 2 || did==3 || did==4 || did==5 ) && udscbAsEquivalent_) {
      did = 1;
    } 
    if( (did == 23 || did == 24) && wzAsEquivalent_) {
      did = 23;
    }
    
    int mid = p.mother1()>0 ? std::abs(process[p.mother1()].id()) : 0;
    
    //if no list of mothers is provided, then all particles
    //in hard process and resonance decays are counted together
    if(mothers_.empty() || mothers_.count(mid) || mothers_.count(-mid))
      ++observedDaughters_[did];
  }
    
  //check if criteria is satisfied
  //inclusive mode: at least as many decay products as requested
  //exclusive mode: exactly as many decay products as requested
  //(but additional particle types not appearing in the list of requested daughter id's are ignored)
  for (const auto &reqpair : requestedDaughters_) {
    int reqid = reqpair.first;
    int reqcount = reqpair.second;
    
    int obscount = 0;
    for (const auto &obspair : observedDaughters_) {
      int obsid = obspair.first;
      
      if (obsid == reqid) {
        obscount = obspair.second;
        break;
      }
    }
    
    //inclusive criteria not satisfied, veto event
    if (obscount < reqcount) return true;

    //exclusive criteria not satisfied, veto event
    if (exclusive_ && obscount > reqcount) return true;
    
  }
  
  //all criteria satisfied, don't veto
  return false;
}

bool ResonanceDecayFilterHook::doVetoResonanceDecays ( Pythia8::Event &  process )

inlineoverride

Definition at line 16 of file ResonanceDecayFilterHook.h.

References checkVetoResonanceDecays(), and LaserDQM_cfg::process.

{ return checkVetoResonanceDecays(process); }

bool ResonanceDecayFilterHook::initAfterBeams ( )

override

Definition at line 8 of file ResonanceDecayFilterHook.cc.

References funct::abs().

Referenced by ResonanceDecayFilterHook().

                                              {
  filter_ = settingsPtr->flag("ResonanceDecayFilter:filter");
  exclusive_ = settingsPtr->flag("ResonanceDecayFilter:exclusive");
  eMuAsEquivalent_ = settingsPtr->flag("ResonanceDecayFilter:eMuAsEquivalent");
  eMuTauAsEquivalent_ = settingsPtr->flag("ResonanceDecayFilter:eMuTauAsEquivalent");
  allNuAsEquivalent_ = settingsPtr->flag("ResonanceDecayFilter:allNuAsEquivalent");
  udscAsEquivalent_ = settingsPtr->flag("ResonanceDecayFilter:udscAsEquivalent");
  udscbAsEquivalent_ = settingsPtr->flag("ResonanceDecayFilter:udscbAsEquivalent");
  wzAsEquivalent_ = settingsPtr->flag("ResonanceDecayFilter:wzAsEquivalent");
  auto mothers = settingsPtr->mvec("ResonanceDecayFilter:mothers");
  mothers_.clear();
  mothers_.insert(mothers.begin(), mothers.end());
  daughters_ = settingsPtr->mvec("ResonanceDecayFilter:daughters");
  
  
  requestedDaughters_.clear();
  
  for (int id : daughters_) {
    int did = std::abs(id);
    if ( did == 13 && (eMuAsEquivalent_ || eMuTauAsEquivalent_)) {
      did = 11;
    }
    if ( did == 15 && eMuTauAsEquivalent_) {
      did = 11;
    }
    if ( (did == 14 || did == 16) && allNuAsEquivalent_) {
      did = 12;
    }
    if ( (did == 2 || did==3 || did==4 ) && udscAsEquivalent_) {
      did = 1;
    }
    if ( (did == 2 || did==3 || did==4 || did==5 ) && udscbAsEquivalent_) {
      did = 1;
    }
    if( (did == 23 || did == 24) && wzAsEquivalent_) {
      did = 23;
    }

    ++requestedDaughters_[std::abs(did)];
  }
  
  return true;
  
}

Member Data Documentation

bool ResonanceDecayFilterHook::allNuAsEquivalent_

private

Definition at line 26 of file ResonanceDecayFilterHook.h.

std::vector<int> ResonanceDecayFilterHook::daughters_

private

Definition at line 31 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::eMuAsEquivalent_

private

Definition at line 24 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::eMuTauAsEquivalent_

private

Definition at line 25 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::exclusive_

private

Definition at line 23 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::filter_

private

Definition at line 22 of file ResonanceDecayFilterHook.h.

std::set<int> ResonanceDecayFilterHook::mothers_

private

Definition at line 30 of file ResonanceDecayFilterHook.h.

std::map<int,int> ResonanceDecayFilterHook::observedDaughters_

private

Definition at line 34 of file ResonanceDecayFilterHook.h.

std::map<int,int> ResonanceDecayFilterHook::requestedDaughters_

private

Definition at line 33 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::udscAsEquivalent_

private

Definition at line 27 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::udscbAsEquivalent_

private

Definition at line 28 of file ResonanceDecayFilterHook.h.

bool ResonanceDecayFilterHook::wzAsEquivalent_

private

Definition at line 29 of file ResonanceDecayFilterHook.h.