GlobalEventData.cxx

Go to the documentation of this file.

// ================================================
// GlobalEventData Singleton class Implementation
// ================================================
//
// This file contains all the common data variables 
// needed by Atlfast to run
// They are stored as statics so that the class needs not be instantiated
// and they can be accesed in algorithms execute() methods.
//
// Namespace Atlfast
//

#include "AtlfastAlgs/GlobalEventData.h"

#include "AtlfastUtils/HepMC_helper/IMCselector.h"
#include "AtlfastUtils/HepMC_helper/IsFinalState.h"
#include "AtlfastUtils/HepMC_helper/IsFromHardScatter.h"
#include "AtlfastUtils/HepMC_helper/RejectType.h"
#include "AtlfastUtils/HepMC_helper/BFieldCutter.h"
#include "AtlfastUtils/HepMC_helper/NCutter.h"
#include "AtlfastUtils/HepMC_helper/Unseen.h"

#include "AtlfastEvent/ParticleCodes.h"

namespace Atlfast {
  //declare static member variables
  GlobalEventData* GlobalEventData::m_instance = 0;
  //Singleton Instance method to makesure only one object instanciated
  GlobalEventData* GlobalEventData::Instance() {
    if(m_instance == 0){
      m_instance = new GlobalEventData();
    }
    return m_instance;
  }
  //set values of static member variables
  void GlobalEventData::setValues(int lumi,
                                  bool fieldOn,
                                  double barrelForwardEta,
                                  int randSeed,
                                  std::vector<int> invisibles,
                                  std::string mcLocation,
                                  std::vector<int> monopoleIDs,
                                  bool justHardScatter){
    m_lumi = lumi;
    m_fieldOn = fieldOn;
    m_barrelForwardEta = barrelForwardEta;
    m_randSeed = randSeed;
    m_invisibles = invisibles;
    m_mcLocation = mcLocation;
    m_adjustMissEtForIsolation = true;
    m_monopoleIDs = monopoleIDs;
    m_justHardScatter = justHardScatter;

    std::vector<int> temp = m_invisibles;
    HepMC_helper::IMCselector* selector;
    std::vector<HepMC_helper::IMCselector*> selectors;
    temp.push_back(ParticleCodes::MUON);

    selector = new HepMC_helper::IsFromHardScatter();
    selectors.push_back( selector ); 
    
    selector = new HepMC_helper::IsFinalState();
    selectors.push_back( selector ); 
    
    selector = new HepMC_helper::RejectType(temp);
    selectors.push_back( selector ); 
    
    
    try{
      selector = new HepMC_helper::BFieldCutter(m_fieldOn);
    }catch(std::string errMsg){
      std::cerr<<"GlobalEventData: error makign a BFieldCutter "
               <<errMsg
               <<std::endl;
    }catch(...){
      std::cerr<<"GlobalEventData: error makign a BFieldCutter "
               <<std::endl;
    }      
    selectors.push_back(selector);
    

    m_visibleToCal = new HepMC_helper::NCutter(selectors); 
    
    // make a selector of invisibles to Atlas
    // neutrinos, muons out of acceptance and user selected pdg ids.
    std::vector<HepMC_helper::IMCselector*>::iterator iter;

    for(iter=selectors.begin(); iter!=selectors.end(); delete *iter, ++iter);
    selectors.erase(selectors.begin(), selectors.end());

    selector = new HepMC_helper::IsFromHardScatter();
    selectors.push_back( selector ); 
    
    selector = new HepMC_helper::IsFinalState();
    selectors.push_back( selector ); 
    
    selector = new HepMC_helper::Unseen(m_invisibles);
    selectors.push_back( selector ); 
    
    m_visibleToAtlas = new HepMC_helper::NCutter(selectors); 

    for(iter=selectors.begin(); iter!=selectors.end(); delete *iter, ++iter);
  }

  
}  //end of namespace bracket

---