// ================================================
 // 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/MCCuts.h"
 #include "AtlfastUtils/HepMC_helper/AscendingEta.h"
 #include "AtlfastUtils/HepMC_helper/SelectJetTag.h"
 #include "AtlfastUtils/HepMC_helper/SelectTauTag.h"
 #include "AtlfastUtils/HepMC_helper/IsCharged.h"
 #include "AtlfastUtils/HepMC_helper/BFieldCutter.h"
 #include "AtlfastUtils/HepMC_helper/RejectType.h"
 #include "AtlfastUtils/HepMC_helper/Unseen.h"
 #include "TruthHelper/GenIMCselector.h"
 #include "TruthHelper/IsGenType.h"
 #include "TruthHelper/IsGenStable.h"
 #include "TruthHelper/NCutter.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,
                                   std::string detEffectsFileName,
                                   std::string muonAlignFileName,
                                   std::string muonLinearFitAlignFileName,
                                   double muonMisalignment){
     /** The events Luninosity */
     m_lumi = lumi;
     /**On/Off switch B-field */
     m_fieldOn = fieldOn;
     /**Geometry info */
     m_barrelForwardEta = barrelForwardEta;
     /** randm seed for Smearer */
     m_randSeed = randSeed;
     /** ids of invis partcles */
     m_invisibles = invisibles;
     /** Storegate key for the mc data */
     m_mcLocation = mcLocation;
     /** Adjust missing ET for isolation? */
     m_adjustMissEtForIsolation = true;
     /** PDG ID for magnetic monopoles */
     m_monopoleIDs = monopoleIDs;
     /** Just consider the hard scatter GenEvent? */
     m_justHardScatter = justHardScatter;
     /** detector effects data file name */
     m_detEffectsFileName = detEffectsFileName;
     /** muon misalignment data file name */
     m_muonAlignFileName = muonAlignFileName;
     /** muon misalignment data file name (linear fit) */
     m_muonLinearFitAlignFileName = muonLinearFitAlignFileName;
     /** muon misalignment value */
     m_muonMisalignment = muonMisalignment;
     /** ids of invis partcles */
     std::vector<int> temp = m_invisibles;
     TruthHelper::GenIMCselector* selector;
     std::vector<TruthHelper::GenIMCselector*> selectors;
     temp.push_back(ParticleCodes::MUON);
     /** Make a selector with particles belonging to the final state*/
     selector = new TruthHelper::IsGenStable();
     selectors.push_back( selector ); 
     /** Rejects GenParticles accordimg to particle pdg id */
     selector = new HepMC_helper::RejectType(temp);
     selectors.push_back( selector ); 
     
     try{
     /** Select GenParticles above pt threshold to reach barrel calorimeter */
     /** Threshold is determined by the magnetic field and GenParticle charge */
       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 TruthHelper::NCutter(selectors);
     
     // make a selector of invisibles to Atlas
     // neutrinos, muons out of acceptance and user selected pdg ids.
     std::vector<TruthHelper::GenIMCselector*>::iterator iter;
 
     for(iter=selectors.begin(); iter!=selectors.end(); delete *iter, ++iter);
     selectors.erase(selectors.begin(), selectors.end());
 
     selector = new TruthHelper::IsGenStable();
     selectors.push_back( selector ); 
 
     selector = new HepMC_helper::Unseen(m_invisibles); 
     selectors.push_back( selector );
     /** Chain selectors together to make particles visible to Atlas */
     m_visibleToAtlas = new TruthHelper::NCutter(selectors);
 
     for(iter=selectors.begin(); iter!=selectors.end(); delete *iter, ++iter);
   }
 
   
 }  //end of namespace bracket
 
 
 
 
 
 
 
 
 
 
 
 

Due to the LXR bug, the updates fail sometimes to remove references to deleted files. The Saturday's full rebuilds fix these problems
This page was automatically generated by the LXR engine.