---

PionBinData.cxx

Go to the documentation of this file.

#include "AtlfastAlgs/PionBinData.h"

namespace Atlfast {

  //-----------------------------------------------
  // PUBLIC - Constructor
  //-----------------------------------------------
  
  PionBinData::PionBinData(
                           BinID& id,
                           vector< ParameterResolutions* > core,
                           vector< ParameterResolutions* > tails,
                           vector< ParameterResolutions* > fractions,
                           vector< ParameterResolutions* > correlations,
                           int randSeed
                           ):
    m_id(id),
    m_cores(core),    
    m_tails(tails),
    m_fractions(fractions),
    m_correlations(correlations)

{
    m_randomEngine = new HepJamesRandom(randSeed);
}    

  
  //-----------------------------------------------
  // PUBLIC - HepMatrix getMatrix(random)
  // returns appropriate sigma matrix using random
  // number 
  //-----------------------------------------------
  
  HepMatrix PionBinData::getMatrix(const TrackTrajectory& traj) const{

    HepMatrix sigma(5,5,0);

    double resolution, random[5];

    random[0] = m_randomEngine->flat() ;
    random[1] = m_randomEngine->flat() ;

    random[2] = random[0];
    random[3] = random[1];
    random[4] = random[0];

    // Diagonals
    for (int param = 0; param<5; param++)
      {
        resolution = m_fractions[param]->resolution( traj );

        //frac << random[param] << endl;
        //frac << resolution << endl;

        if (random[param] < resolution )
          {
            sigma[param][param] = 
              pow( m_cores[param]->resolution( traj ), 2 );
          }
        else
          {
            sigma[param][param] = 
              pow( m_tails[param]->resolution( traj ), 2 );
          }
      }
    // off diagonals
    // (1,3)

    sigma(3,1) = sigma(1,3) =
      m_correlations[0]->resolution( traj ) *
      pow( (sigma(1,1)*sigma(3,3)) , 0.5 );
    // (1,5)
    sigma(5,1) = sigma(1,5) =
      m_correlations[1]->resolution( traj ) *
      pow( (sigma(1,1)*sigma(5,5)) , 0.5 );
    // (2,4)
    sigma(2,4) = sigma(4,2) =
      m_correlations[2]->resolution( traj ) *
      pow( (sigma(2,2)*sigma(4,4)) , 0.5 );
    // (3,5)
    sigma(3,5) = sigma(5,3) =
      m_correlations[3]->resolution( traj ) * 
      pow( (sigma(3,3)*sigma(5,5)) , 0.5 );

    // DONE!
    return sigma;
    
  }
  
}//namespace bracket

---