#include "AtlfastAlgs/ReconstructorFactory.h"
 #include "AtlfastAlgs/CompoundReconstructor.h"
 #include "AtlfastAlgs/ParticleSubsystemPartitionTool.h"
 #include "AtlfastAlgs/DefaultReconstructor.h"
 #include "AtlfastAlgs/ElectronReconstructor.h"
 #include "AtlfastAlgs/MuonReconstructor.h"
 #include "AtlfastAlgs/MuonMisalignmentReconstructor.h"
 #include "AtlfastAlgs/PhotonReconstructor.h"
 
 namespace Atlfast{
   ReconstructorFactory::ReconstructorFactory(std::string fileName,MsgStream &log):m_log(log),m_fl(fileName){
     /////////////////////////////////////////////////////////////////////////////////////////////////////////
     //To Add New Reconstructors Follow Instructions Below...
     /////////////////////////////////////////////////////////////////////////////////////////////////////////
     //step 1)
     //add lines like example below for adding mapping of PROBLEM to marginal and replacement reconstructors
     //all new reconstructors must be able to take a single construction arguement that is a ProblematicRegion
     //which will contain the smearing parameters and the acceptance criteria.
     //example:
     //    m_reconstructorTypeMap[DetectorParameters::Scaling] = &smearerForge<testReconstructor>;
     /////////////////////////////////////////////////////////////////////////////////////////////////////////
     m_reconstructorTypeMap[DetectorParameters::Misalignment] = &reconstructorForge<MuonMisalignmentReconstructor>;
                                                                
     /////////////////////////////////////////////////////////////////////////////////////////////////////////
     //step 2)
     // Below we set up the mapping of particle pdgID to SUBSYSTEMs that affect it.
     /////////////////////////////////////////////////////////////////////////////////////////////////////////
     m_particleSubsystemMap.insert(std::pair<int,DetectorParameters::SUBSYSTEM>(11,DetectorParameters::Calorimeter));
     m_particleSubsystemMap.insert(std::pair<int,DetectorParameters::SUBSYSTEM>(13,DetectorParameters::Calorimeter));
     m_particleSubsystemMap.insert(std::pair<int,DetectorParameters::SUBSYSTEM>(13,DetectorParameters::MuonSystem));
    
   }
 
   IReconstructor* ReconstructorFactory::makeReconstructor(   int particleType,
                                                              const int randSeed, 
                                                              const int lumi, 
                                                              const int muSmearKey, 
                                                              std::string& muonResFile, 
                                                              std::vector<double> smearParamArray, 
                                                              int smearParamSchema )const{
 
     // A) firstly set up the compound reconstructor with the default reconstructor
     IDefaultReconstructor *defaultReconstructor;
 
     if (particleType == 11) 
       {
         m_log << MSG::DEBUG << "Creating new ElectronReconstructor randSeed: "<<randSeed<<" lumi: "<<lumi<<endreq;
         defaultReconstructor = new ElectronReconstructor(randSeed, lumi, m_log);
       } 
     else if (particleType == 13) 
       {
         m_log << MSG::DEBUG << "Creating new MuonReconstructor randSeed: "<<randSeed<<" lumi: "<<lumi<<endreq;
         m_log << MSG::DEBUG << "                         smearKey: "<<muSmearKey<<" resFile: "<<muonResFile<<endreq;
         defaultReconstructor = new MuonReconstructor(randSeed, lumi, muSmearKey, muonResFile, m_log);
       } 
     else if (particleType == 22) 
       {
         m_log << MSG::DEBUG << "Creating new PhotonReconstructor randSeed: "<<randSeed<<" lumi: "<<lumi<<endreq;
         defaultReconstructor = new PhotonReconstructor(randSeed, lumi, m_log);
       } 
     else 
       {
         m_log << MSG::DEBUG << "Creating new DefaultReconstructor randSeed: "<<randSeed<<endreq;
         defaultReconstructor = new DefaultReconstructor(randSeed);
       }
     defaultReconstructor->setSmearParameters(smearParamArray);
     defaultReconstructor->setSmearParamSchema(smearParamSchema);
     CompoundReconstructor *compoundReconstructor = new CompoundReconstructor(defaultReconstructor);
 
     // B) now add marginal and replacement reconstructors that apply for the given particleType
    m_log << MSG::DEBUG << "Partitioning problematic regions."<<endreq;
    std::vector<const ProblematicRegion*> pbr(m_fl.GetRegions());//need to copy to a new one as partition re-orders i.e. not const
    std::vector<const ProblematicRegion*>::iterator itr, start, end, bound;
    start = pbr.begin();
    end = pbr.end();
    ParticleSubsystemPartitionTool ptp(particleType,m_particleSubsystemMap);
    bound = std::partition(start, end, ptp); //now sorted into those that affect particle 
                                             //and those that dont
    
    m_log << MSG::DEBUG << "Adding defined replacement/marginal reconstructor to compound reconstructor."<<endreq;
    std::map<DetectorParameters::PROBLEM,reconstructorForgePointer>::const_iterator mItr;
    for(itr=pbr.begin();itr!=bound;++itr){
      mItr = m_reconstructorTypeMap.find((*itr)->problem());
      if(mItr!=m_reconstructorTypeMap.end()){
        IAcceptingReconstructor *ias = mItr->second(**itr);
        if(ias->is_replacement()){ compoundReconstructor->addReplacement(ias); }
        else{ compoundReconstructor->addMarginal(ias); }
      }
      else{
        m_log<<MSG::ERROR<<"ReconstructorFactory: You shouldn't be seeing this message. If you are then no mapping has been"<<endreq;
        m_log<<MSG::ERROR<<"                      established between your chosen PROBLEM type and a concrete reconstructor implementation!"<<endreq;
        assert(false);
      }
    }
 
    m_log << MSG::DEBUG << "Returning CompoundReconstructor."<<endreq;
    return compoundReconstructor;
   }
   
 
 }//end of namespace
 

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.