TestTrkParametersIdentificationHelpers.cxx File Reference

#include "TrkParametersIdentificationHelpers/TrackParametersIdentificationHelper.h"
#include <iostream>
#include <iomanip>
#include <bitset>

Include dependency graph for TestTrkParametersIdentificationHelpers.cxx:

Go to the source code of this file.

Functions
void	testBitField ()
void	testCalo ()
int	main ()

Function Documentation

main()

int main

(

)

getCoefficientMap(label, EigenIdxList)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light
2. EigenIdxList is user defined vector containing all eigenvector index that user interested in. output: Map of format map<string, map<string, float>> containing decomposition coefficient of the list of eigenvectors defined by EigenIdxList.

getCoefficients(label, evIdx)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light
2. evIdx: The index of eigenvector user interested in. output value: vector of coefficient values. The order is the same as output given by getListOfOriginalNuisanceParameters()

getListOfOriginalNuisanceParameters(label)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light output value: List of original nuisance parameter names.

getNumEigenVectors(label)

input value:

1. label: falvor label in std::string format, could be one of B, C, T, Light return value: number of eigen vectors used for chosen label. Return 0 if error occured.

Definition at line 80 of file TestTrkParametersIdentificationHelpers.cxx.

           {
  
  testBitField();
  testCalo();
}

testBitField()

void testBitField

(

)

Definition at line 11 of file TestTrkParametersIdentificationHelpers.cxx.

                   {
  Trk::BitField<unsigned int> field(31,1);
  TrackParametersIdentifier id=0;
  unsigned int max = -1;
  std::cout << " starting loop " << max << std::endl;
  for( unsigned int i=0;i<max;++i ){
    if( !field.encode(i,id) ){
      std::cout << " encoding failed for " << i << std::endl; 
      break;
    }
    unsigned int value = field.decode(id);
    std::cout << "  i " << std::setw(3) << i << " value " << std::setw(3) << value << "  " << std::bitset<32>(i) << " id " << std::bitset<32>(id) << std::endl;
  }
} 

testCalo()

void testCalo

(

)

Definition at line 26 of file TestTrkParametersIdentificationHelpers.cxx.

                {
  
  Trk::TrackParametersIdHelper helper;
  unsigned int good = 0;
  unsigned int bad  = 0;
  unsigned int goodCalo = 0;
  unsigned int badCalo  = 0;
  
  // loop over technologies
  for( unsigned int tech = AtlasDetDescr::fFirstAtlasDetTechnology; tech < AtlasDetDescr::fNumAtlasDetTechnologies; ++tech ){
    
    bool isCalo = (tech >= AtlasDetDescr::fFirstAtlasCaloTechnology) && (tech <= AtlasDetDescr::fLastAtlasCaloTechnology);
    
    // loop over samplings
    for( unsigned int sample = CaloSampling::PreSamplerB; sample < CaloSampling::Unknown; ++sample ){
      
      for( unsigned int isEntry = 0; isEntry < 2 ; ++isEntry ){
        bool entry = isEntry == 1;
        TrackParametersIdentifier id = helper.encode( static_cast<AtlasDetDescr::AtlasDetTechnology>(tech), 
                                                      static_cast<CaloSampling::CaloSample>(sample), entry );
        AtlasDetDescr::AtlasDetTechnology itech = helper.technology(id);
        CaloSampling::CaloSample isample = helper.caloSample(id);
        bool ientry = helper.isEntryToVolume(id);
        bool ivalid = helper.isValid(id);
        bool printDetails = false;
        if( !ivalid || sample != isample || entry != ientry ) {
          ++bad;
          if( isCalo ) {
            std::cout << "WARNING bad decoding/encoding " << std::endl;
            ++badCalo;
            printDetails = true;
          }
        }else{ 
          ++good;
          if( isCalo ) { ++goodCalo;
          } else {
            std::cout << " good but no calo!!! " << std::endl;
            printDetails = true;
          }
        }
        if( printDetails ){
          std::cout << " tech  " << std::setw(4) << tech  << " sample  " << std::setw(4) << sample  << " isEntry " << entry << std::endl
                    << " itech " << std::setw(4) << itech << " isample " << std::setw(4) << isample << " isEntry " << ientry
                    << " valid " << ivalid << " id " << std::bitset<32>(id) << std::endl;
        }
      }
    }
  }
  std::cout << "all:  good " << good << " bad " << bad << std::endl;
  std::cout << "calo: good " << goodCalo << " bad " << badCalo << std::endl;
    
}