#include "GaudiKernel/ToolFactory.h"
 #include "GaudiKernel/SmartDataPtr.h"
 #include "GaudiKernel/IDataProviderSvc.h"
 
 //----------------------------------------------------------------//
 #include "TrkTrack/Track.h"
 
 //----------------------------------------------------------------//
 #include "muonEvent/Muon.h"
 #include "Particle/TrackParticle.h"
 #include "Particle/TrackParticleContainer.h"
 #include "IsolationTool/IIsolationTool.h"
 
 /////////////////////////////////////////////////////////
 #include "STACOTools/StacoToCombinedMuonTool.h"
 #include "StacoMuonQualityWordHelper.h"
 
 //----------------------------------------------------------------//
 #include "STACOEvent/StacoObject.h"
 
 StacoToCombinedMuonTool::StacoToCombinedMuonTool(const std::string& t, 
                                             const std::string& n,
                                             const IInterface*  p ):AthAlgTool(t,n,p),
  p_IIsolationTool("MuonIsolationTool/MuonIsolationTool")                                            
 {
 
   declareInterface<IStacoToCombinedMuonTool>(this);
 
   declareProperty("IIsolationTool" ,p_IIsolationTool );
 
 }
 
 StacoToCombinedMuonTool::~StacoToCombinedMuonTool(){}
 
 // Initialize
 StatusCode StacoToCombinedMuonTool::initialize() {
 
   StatusCode sc = StatusCode::SUCCESS;
 
   sc = AthAlgTool::initialize(); 
   if ( sc.isFailure() ) {
     msg(MSG::FATAL) << " AthAlgTool::initialize() failed" << endreq;
     return( StatusCode::FAILURE );
   }
 
 
 //Retrieve p_IIsolationTool
   if ( p_IIsolationTool.retrieve().isFailure() ) {
     msg(MSG::WARNING) << "Failed to retrieve tool " << p_IIsolationTool << endreq;
     p_IIsolationTool = 0 ;
   }else{
     msg(MSG::INFO) << "Retrieved tool " << p_IIsolationTool << endreq;
   }
 
   msg(MSG::INFO) << "Initialisation ended     " << endreq;
 
   return StatusCode::SUCCESS;
 
 }
 
 // Finalize
 StatusCode StacoToCombinedMuonTool::finalize(){return StatusCode::SUCCESS;}
 
 // Convert
 Analysis::Muon*  StacoToCombinedMuonTool::convert(
                          const StacoObject*                 pStacoObject,
                          const Rec::TrackParticleContainer* pTPIDColl  ,
                          const Rec::TrackParticleContainer* pTPMSColl ,
                          const Rec::TrackParticleContainer* pTPMSOnlyColl ,
                          const Rec::TrackParticleContainer* pTPStacoColl    ,
                          const Rec::TrackParticle*          pTrackPartStaco  
 ){
 
 
   const Rec::TrackParticle* pTrackPartID     = pStacoObject->MuIdet();
 
   const Rec::TrackParticle* pTrackPartMS     = pStacoObject->MuSpectro();
     
   const Rec::TrackParticle* pTrackPartMSOnly = pStacoObject->MuSpectroMSOnly();
     
   //Now create a MS Standalone for eta > 2.5
   if(pTrackPartID==0) {
     Analysis::Muon* MuonSA = new Analysis::Muon(
        MuonParameters::STACO,
        pTrackPartMS,
        true);
        
     MuonSA->set_muonExtrapolatedTrackParticleContainer(pTPMSColl); 
 
 //Hopefully temporary 
     if (pTrackPartMSOnly) {
       MuonSA->set_muonSpectrometerTrackParticle(pTPMSOnlyColl,pTrackPartMSOnly);
     }else{
       MuonSA->set_muonSpectrometerTrackParticle(pTPMSColl,pTrackPartMS);
     }
 
     //Calculate parametrized energy loss     
     double paramEnergy = pStacoObject->Eloss();
     double paramEnergyErr = 0.20 * paramEnergy ;
     Trk::EnergyLoss aEnergyLoss (paramEnergy, paramEnergyErr);
     
     MuonSA->set_bestMatch(false);
   //calculate track isolation energy  
     if (p_IIsolationTool){
       StatusCode sc = p_IIsolationTool->dressing( MuonSA, &aEnergyLoss  );
       if ( sc.isFailure() )  msg(MSG::WARNING) << "Fail to dress muon " << endreq;
     }
 
     if ( MuonSA ) {
       StacoMuonQualityWordHelper aStacoMuonQualityWordHelper ;
       MuonQualityWord aMuonQualityWord ;
       aStacoMuonQualityWordHelper.Set_isTight(aMuonQualityWord) ;
       MuonSA->set_isMuonBits(aMuonQualityWord);
     }
 
     return MuonSA ;
 
 
   }else{
   //Now create a CombinedMuon 0.< eta <2.5 
  
       Analysis::Muon* CBMuon = new Analysis::Muon(
         MuonParameters::STACO,
         pTrackPartID, 
         pTrackPartMS, 
         pTrackPartStaco);
 
     CBMuon->set_muonExtrapolatedTrackParticleContainer(pTPMSColl);
     CBMuon->set_inDetTrackParticleContainer(pTPIDColl);
     CBMuon->set_combinedMuonTrackParticleContainer(pTPStacoColl);
     
 //Hopefully temporary 
     if (pTrackPartMSOnly) {
       CBMuon->set_muonSpectrometerTrackParticle(pTPMSOnlyColl,pTrackPartMSOnly);
     }else{
       CBMuon->set_muonSpectrometerTrackParticle(pTPMSColl,pTrackPartMS);
     }
 
     //Calculate parametrized energy loss
     double paramEnergy = pStacoObject->Eloss();
     double paramEnergyErr = 0.20 * paramEnergy ;
     Trk::EnergyLoss aEnergyLoss (paramEnergy, paramEnergyErr);
     
     CBMuon->set_bestMatch(true);
     
     double chi2CB = pTrackPartStaco->fitQuality()->chiSquared()-
                     pTrackPartID->fitQuality()->chiSquared()-
                     pTrackPartMS->fitQuality()->chiSquared();
 
     CBMuon->set_matchChi2(chi2CB);
     int Nodf = 5;
     CBMuon->set_matchNumberDoF(Nodf);
     
     if (p_IIsolationTool){
       StatusCode sc = p_IIsolationTool->dressing( CBMuon, &aEnergyLoss  );
       if ( sc.isFailure() )  msg(MSG::WARNING) << "Fail to dress muon " << endreq;
     }
 
     if ( CBMuon ) {
       StacoMuonQualityWordHelper aStacoMuonQualityWordHelper ;
       MuonQualityWord aMuonQualityWord ;
       aStacoMuonQualityWordHelper.Set_isTight(aMuonQualityWord) ;
       CBMuon->set_isMuonBits(aMuonQualityWord);
     }
 
     return CBMuon ;
 
   }     
           
   return 0;
 
 }

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.