/Users/jmonk/Physics/ForIA/src/AnalysisTools/MatchToTruth.cxx

Go to the documentation of this file.

#include "ForIA/AnalysisTools/MatchToTruth.hh"
#include "ForIA/Utils.hh"
#include "ForIA/ConfigurationFile.hh"
namespace ForIA {
  // constructor
  MatchToTruth::MatchToTruth() {
  }
  
  void MatchToTruth::init(){
    
    // grab the truth partciles
    ConfigurationFile configFile;
    configFile.init("analysis_config.dat");
    
    double TruthE,TruthET,TruthEtaMin, TruthEtaMax;
    configFile.readInto( TruthE, "TruthE", 0.);
    configFile.readInto( TruthET, "TruthET", 0.);
    configFile.readInto( TruthEtaMin, "TruthEtaMin", 0.0);
    configFile.readInto( TruthEtaMax, "TruthEtaMax", 4.8);
    
    
    m_truthSelection.setEMin( TruthE );
    m_truthSelection.setPTMin( TruthET );
    m_truthSelection.setEtaMin( TruthEtaMin );
    m_truthSelection.setEtaMax( TruthEtaMax );
    m_truthSelection.setIsPrimary(true);
    
    m_truthSelection.setPdgId(22);    
    
    return;
  }
  
  bool MatchToTruth::setEvent(const Event& evt) {
    m_truthSelection.setEvent(evt);
    return true;
  }

  float MatchToTruth::MinEta() {

    TruthParticleVector truthVectors = m_truthSelection.particles();
    float minEta = 99.0F;
    for(TruthParticleVector::const_iterator truth = truthVectors.begin();
        truth != truthVectors.end(); ++truth){
      float feta = fabs((*truth)->eta());
      if (feta < minEta) minEta = feta;
    }
    return minEta;
  }


  bool MatchToTruth::Match(const MomentumVector& fourVectors) {
    
    TruthParticleVector truthVectors = m_truthSelection.particles();
    std::cout<<"Matching "<<fourVectors.size() << " " << truthVectors.size() << std::endl;
    for(MomentumVector::const_iterator fourVector = fourVectors.begin();
        fourVector != fourVectors.end(); ++fourVector){
//      float minDR = 99.0F;
      for(TruthParticleVector::const_iterator truth = truthVectors.begin();
          truth != truthVectors.end(); ++truth){
        // std::cout<<(*fourVector)->eta() << " " << (*truth)->eta()<<std::endl;
        // float DR = Kinematics::deltaR((*fourVector),(*truth));
        //if (DR < minDR) minDR = DR;
      }
    }
    return true;
  }
  
 float MatchToTruth::MatchDR(IFourMomentumConstPtr p1) {

     TruthParticleVector truthVectors = m_truthSelection.particles();
     float minDR = 99.0F;
     for(TruthParticleVector::const_iterator truth = truthVectors.begin();
           truth != truthVectors.end(); ++truth){
       //       std::cout<<p1->eta() << " " << (*truth)->eta()<<std::endl;
         float DR = deltaR(p1,(*truth));
         if (DR < minDR) {
           minDR = DR;
         }
     }
     return minDR;
  }



IFourMomentumConstPtr  MatchToTruth::MatchedParticle(IFourMomentumConstPtr p1) {

     TruthParticleVector truthVectors = m_truthSelection.particles();
     IFourMomentumConstPtr t1;
     float minDR = 99.0F;
     for(TruthParticleVector::const_iterator truth = truthVectors.begin();
           truth != truthVectors.end(); ++truth){
       //       std::cout<<p1->eta() << " " << (*truth)->eta()<<std::endl;
         float DR = deltaR(p1,(*truth));
         if (DR < minDR) {
           minDR = DR;
           t1 = (*truth);
         }
     }
     return t1;
  }

}