CalStripAtten.cxx

Go to the documentation of this file.

// $Id: CalStripAtten.cxx,v 1.7 2007/01/15 19:52:01 rhatcher Exp $
//
// CalStripAtten
//
// DB Row class corresponding to the double-exponential attenuation curve along a strip.
//
// Nathaniel Tagg n.tagg1@physics.ox.ac.uk

#include "Calibrator/CalStripAtten.h"
#include "MessageService/MsgService.h"
#include "DatabaseInterface/DbiOutRowStream.h"
#include "DatabaseInterface/DbiResultSet.h"
#include "DatabaseInterface/DbiValidityRec.h"
#include <cmath>

ClassImp(CalStripAtten)

//   Definition of static data members
//   *********************************

CVSID("$Id: CalStripAtten.cxx,v 1.7 2007/01/15 19:52:01 rhatcher Exp $\n  \
      CVSID_DBIRESULTPTR ");

//  Instantiate associated Result Pointer class.
//  *******************************************

#include "DatabaseInterface/DbiResultPtr.tpl"
template class  DbiResultPtr<CalStripAtten>;

#include "DatabaseInterface/DbiWriter.tpl"
template class  DbiWriter<CalStripAtten>;

// Definition of member functions (alphabetical order)
// ***************************************************


//.....................................................................

void CalStripAtten::Fill(DbiResultSet& rs, 
                              const DbiValidityRec* /* vrec */) {

//
//
//  Purpose:  Fill object from Result Set
//
//  Arguments: 
//    rs           in    Result Set used to fill object
//    vrec         in    Associated validity record (or 0 if filling
//                                                    DbiValidityRec)
//
//  o Fill object from current row of Result Set.

//  Program Notes:-
//  =============

//  This method demonstrates both the "dumb" fill method (just
//  load the data as it comes) and the smart method (check column
//  name and load according to column order). 

   if ( rs.TableName() == "CALSTRIPATTEN" ) {
      // Dumb method.
     rs  >> fSEIdEncoded >> fLambda1 >> fLambda2 >> fFrac1 >> fLambda1Err >> fLambda2Err >> fFrac1Err; 
     fSEIdKey = PlexStripEndId(fSEIdEncoded).BuildPlnStripEndKey();
   }

   else {
      
      // Smart method
      Int_t numCol = rs.NumCols();
      //  The first column (SeqNo) has already been processed.
      for (Int_t curCol = rs.HasRowCounter() ? 3 : 2; curCol <= numCol; ++curCol) {
         string colName = rs.CurColName();
         if ( colName == "SEIDENCODED" ) {
           rs >> fSEIdEncoded; 
           fSEIdKey =  PlexStripEndId(fSEIdEncoded).BuildPlnStripEndKey();       
         }
         else if( colName == "LAMBDA1" )     rs >> fLambda1;
         else if( colName == "LAMBDA2" )     rs >> fLambda2;
         else if( colName == "FRAC1"   )     rs >> fFrac1;
         else if( colName == "LAMBDA1ERR" )     rs >> fLambda1Err;
         else if( colName == "LAMBDA2ERR" )     rs >> fLambda2Err;
         else if( colName == "FRAC1ERR"   )     rs >> fFrac1Err;

         else {
           MSG("Dbi",Msg::kDebug) << "Ignoring column " << curCol 
                                  << "(" << colName << ")"
                                  << "; not part of CalStripAtten" 
                                  << endl;
           rs.IncrementCurCol();
         }
      }
   }
   
   
}


//.....................................................................
void CalStripAtten::Store(DbiOutRowStream& ors,
                               const DbiValidityRec* /* vrec */) const {
  //
  //
  //  Purpose:  Stream object to output row stream
  //
  //  Arguments: 
  //    ors          in     Output row stream.
  //    vrec         in    Associated validity record (or 0 if filling
  
  ors << fSEIdEncoded  << fLambda1 << fLambda2 << fFrac1
      << fLambda1Err << fLambda2Err << fFrac1Err; 
}
//.....................................................................




Float_t  CalStripAtten::GetAttenuation(const Float_t wlsLen) const
{

  if(wlsLen<0) return 1; // Disallow nonphysical queries.

  Float_t atten =  fFrac1  * exp(-wlsLen/fLambda1) // Short atten length
            + (1.0-fFrac1) * exp(-wlsLen/fLambda2); // Long atten length

  return  atten;
}


Float_t  CalStripAtten::GetAttenuation(const Float_t wlsLen, 
                                       const Float_t* sigma_err) const
{

  if(wlsLen<0) return 1;

  Float_t frac = fFrac1 + fFrac1Err*sigma_err[0];
  Float_t l1  = fLambda1 + fLambda1Err*sigma_err[1];
  Float_t l2  = fLambda2 + fLambda2Err*sigma_err[2];

  Float_t atten =  frac  * exp(-wlsLen/l1) // Short atten length
            + (1.0-frac) * exp(-wlsLen/l2); // Long atten length

  return atten;


}

---