AiresUtilities.cc

Go to the documentation of this file.

#include <io/AiresUtilities.h>

#include <utl/Particle.h>
using utl::Particle;
#include <utl/GeometryUtilities.h>
#include <utl/AugerUnits.h>
#include <utl/ErrorLogger.h>

#include <map>
using std::map;
#include <utility>
#include <sstream>
using std::ostringstream;

#include <iostream>
using namespace std;

namespace {

  map<int, int> airesToPDGMap;

  void InsertAiresToPDG(int theAiresCode, int thePDGCode) {

    if (!airesToPDGMap.insert(std::make_pair(theAiresCode, 
                                             thePDGCode)).second) {
      ostringstream msg;
      msg << "Cannot insert pair ("
          << theAiresCode << ", " << thePDGCode
          << ") into AiresToPDG map.";
      ERROR(msg);
    
    }
  
  }

  void InitAiresToPDGMap() {
    
    InsertAiresToPDG(1, Particle::ePhoton);
    InsertAiresToPDG(-2, Particle::eElectron);
    InsertAiresToPDG(2, Particle::ePositron);
    InsertAiresToPDG(3, Particle::eAntiMuon);
    InsertAiresToPDG(-3, Particle::eMuon);
    InsertAiresToPDG(4, Particle::eAntiTau);
    InsertAiresToPDG(-4, Particle::eTau);
    InsertAiresToPDG(6, Particle::eNuElectron);
    InsertAiresToPDG(-6, Particle::eAntiNuElectron);
    InsertAiresToPDG(7, Particle::eNuMuon);
    InsertAiresToPDG(-7, Particle::eAntiNuMuon);
    InsertAiresToPDG(8, Particle::eNuTau);
    InsertAiresToPDG(-8, Particle::eAntiNuTau);
    InsertAiresToPDG(10, Particle::ePiZero);
    InsertAiresToPDG(-11, Particle::ePiPlus);
    InsertAiresToPDG(11, Particle::ePiMinus);
    InsertAiresToPDG(12, Particle::eKaon0S);
    InsertAiresToPDG(13, Particle::eKaon0L);
    InsertAiresToPDG(14, Particle::eKaonPlus);
    InsertAiresToPDG(-14, Particle::eKaonMinus);
    InsertAiresToPDG(15, Particle::eEta);
    InsertAiresToPDG(20, Particle::eLambda);
    InsertAiresToPDG(-20, Particle::eAntiLambda);
    InsertAiresToPDG(-30, Particle::eAntiNeutron);
    InsertAiresToPDG(30, Particle::eNeutron);
    InsertAiresToPDG(31, Particle::eProton);
    InsertAiresToPDG(-31, Particle::eAntiProton);

  }

} // namespace

int io::Aires::AiresToPDG(int theAiresCode) {

  if (theAiresCode < 100) {

    if (!airesToPDGMap.size())
      InitAiresToPDGMap();
    
    map<int, int>::const_iterator index =
      airesToPDGMap.find(theAiresCode);
    
    if (index != airesToPDGMap.end())
      return index->second;
  
  } else if (theAiresCode < 100 * 32*26 + 31) {

    // From Aires manual :
    // code = 100 + 32*Z + (N - Z +8)
    // 0 <= N - Z + 8 <= 31
    // 1 <= z <= 26 
    //
    unsigned int Z = (theAiresCode-100) / 32;
    unsigned int A = Z + theAiresCode-108-31*Z;
    return Particle::NucleusCode(Z,A);
  } 
  return Particle::eUndefined;
} 

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: