AtlfastC.h

Go to the documentation of this file.

// ================================================
// AtlfastC class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
//
// This version....
//
//
// Namespace Atlfast::
//
// class: AtlfastC
// Author: S. Allwood, C.Collins-Tooth, S. Ferrag, C. Wright
//
// Description: 
//
//  Algorithm for particle and jet efficiencies.
//  Output should be new Atlfast collections with efficiencies applied.
//
// ................................................................
//

#ifndef ATLFAST_C
#define ATLFAST__C

// STL
#ifndef STD_VECTOR_H
#include <vector>
#define STD_VECTOR_H
#endif

#ifndef STD_STRING_H
#include <string>
#define STD_STRING_H
#endif

// Gaudi 
#ifndef GAUDIKERNEL_ALGORITHM_H
#include "GaudiKernel/Algorithm.h"
#define GAUDIKERNEL_ALGORITHM_H
#endif 

#include "AtlfastEvent/ReconstructedParticle.h"
#include "CLHEP/Random/RandomEngine.h"
#include "CLHEP/Random/RandFlat.h"

#include "AtlfastEvent/TypeVisitor.h"
#include "AtlfastEvent/CollectionDefs.h"
#include "AtlfastEvent/Interpolator.h"
#include "AtlfastAlgs/ContainerDispatcher.h"

#include "AtlfastEvent/ICluster.h"
#include "CLHEP/Vector/LorentzVector.h"

class MsgStream;
class ISvcLocator;


namespace Atlfast {

  class TesIO;

  using std::string;

  class AtlfastC : public Algorithm
{
 private:
  
  
  //switch for the different methods


  bool  m_ApplyElecElecEff;
  bool  m_ApplyElecPhotEff;
  bool  m_ApplyElecMuonEff;
  bool  m_ApplyPhotPhotEff;
  bool  m_ApplyPhotElecEff;
  bool  m_ApplyPhotMuonEff;
  bool  m_ApplyMuonElecEff;
  bool  m_ApplyMuonPhotEff;
  bool  m_ApplyMuonMuonEff;
  bool  m_ApplyJetJetEff;
  bool  m_ApplyJetElecEff;
  bool  m_ApplyJetPhotEff;
  bool  m_ApplyJetMuonEff;
  bool  m_ApplyElecJetEff;
  bool  m_ApplyPhotJetEff;
  bool  m_ApplyMuonJetEff;


  Jet* createJet(const ReconstructedParticle* rp);



  //other variables
  
  bool m_notInit;
  TesIO* m_tesIO;

  HepRandomEngine* m_pRandomEngine;
  RandFlat* m_pRandFlatGenerator;
  
  //  std::vector<Jet*>  m_Jets;

  std::vector<ReconstructedParticle*>  originalElectrons;
  std::vector<ReconstructedParticle*>  originalPhotons;
  std::vector<ReconstructedParticle*>  originalMuons;
  std::vector<Jet*>                    originalJets;

  ReconstructedParticleCollection* m_AtlfastCElectrons;
  ReconstructedParticleCollection* m_AtlfastCPhotons;
  ReconstructedParticleCollection* m_AtlfastCMuons;
  JetCollection* m_AtlfastCJets;

  int ptnum, etanum;
  std::vector<ReconstructedParticle*>::const_iterator m_it;
  std::vector<Jet*>::const_iterator j_it;
  
  //--------------------------------------------------
  // Paths in the Transient Event store to get/put things
  // These are currently set in member variables so that
  // they can be overwritten by the job options service
  //--------------------------------------------------
  //  ICluster  * dummy_cluster ;

  std::string m_inputeLocation ;
  std::string m_inputphotonLocation ;
  std::string m_inputmuonLocation ;
  std::string m_inputjetLocation ;

  std::string m_outputeLocation ;
  std::string m_outputphotonLocation ;
  std::string m_outputmuonLocation ;
  std::string m_outputjetLocation ;

  //holds the Stotegate MC location
  std::string m_mcLocation;

  std::string  m_elecelecEfile;
  std::string  m_elecphotEfile;
  std::string  m_elecmuonEfile;
  std::string  m_photelecEfile;
  std::string  m_photphotEfile;
  std::string  m_photmuonEfile;
  std::string  m_muonelecEfile;
  std::string  m_muonphotEfile;
  std::string  m_muonmuonEfile;
  std::string  m_jetjetEfile;
  std::string  m_jetelecEfile;
  std::string  m_jetphotEfile;
  std::string  m_jetmuonEfile;
  std::string  m_elecjetEfile;
  std::string  m_photjetEfile;
  std::string  m_muonjetEfile;
  
  Interpolator* m_electron_interpolator;
  Interpolator* m_photon_interpolator;
  Interpolator* m_muon_interpolator;
  Interpolator* m_jet_interpolator ;
  Interpolator* m_elecphot_interpolator;
  Interpolator* m_elecmuon_interpolator;
  Interpolator* m_elecjet_interpolator;
  Interpolator* m_photelec_interpolator;
  Interpolator* m_photmuon_interpolator;
  Interpolator* m_photjet_interpolator;
  Interpolator* m_muonelec_interpolator;
  Interpolator* m_muonphot_interpolator;
  Interpolator* m_muonjet_interpolator;
  Interpolator* m_jetelec_interpolator;
  Interpolator* m_jetphot_interpolator;
  Interpolator* m_jetmuon_interpolator;


 public:

  //-------------------------
  // Constructors/Destructors
  //
  // Gaudi requires that the constructor takes certain arguments
  // (and passes them directly to the constructor of the base class)
  //-------------------------

  AtlfastC( const std::string& name, ISvcLocator* pSvcLocator ) ; 

  virtual ~AtlfastC();

  StatusCode ApplyEfficiency(std::vector<ReconstructedParticle*> , ReconstructedParticleCollection* , Interpolator* ,int ) ;
  StatusCode ApplyEfficiency(std::vector<Jet*> , ReconstructedParticleCollection* , Interpolator* ,int ) ;
  StatusCode ApplyEfficiency(std::vector<ReconstructedParticle*> , JetCollection* , Interpolator* ) ;
  StatusCode ApplyEfficiency(std::vector<Jet*> , JetCollection* , Interpolator* ) ;

  double getEfficiencyMatrixElement(double , double , Interpolator* ) ;


  StatusCode initialize() ;
  StatusCode execute() ;
  StatusCode finalize() ;

};

}  //end of namespace bracket

#endif

---