CalDrift.cxx

Go to the documentation of this file.

// CalDrift.cxx
//
// Package: Dbi (Database Interface).
//
// The DbiTableRow class responsible for the drift calibration
//
//
//
//
// j.evans2@physics.ox.ac.uk

#include "DatabaseInterface/DbiOutRowStream.h"
#include "MessageService/MsgService.h"

#include "Calibrator/CalDrift.h"

ClassImp(CalDrift)

CVSID("$Id: CalDrift.cxx,v 1.13 2009/05/22 10:57:19 tinti Exp $");

#include "DatabaseInterface/DbiResultPtr.tpl"
#include "DatabaseInterface/DbiResultSet.h"
#include "DatabaseInterface/DbiWriter.tpl"
#include "DatabaseInterface/DbiValidityRec.h"
#include "Validity/VldContext.h"

template class DbiResultPtr<CalDrift>;
template class DbiWriter<CalDrift>;

//____________________________________________________________________72
CalDrift::CalDrift() : DbiTableRow()
{
  LEA_CTOR; // Leak detection macro.
  fbackMedian = 0.0;
  fbackMedianErr = 0.0;
  ffrontMedian = 0.0;
  ffrontMedianErr = 0.0;
  fnonFidMedian = 0.0;
  fnonFidMedianErr = 0.0;
  fnumMuons = 0;
  fnumPlanes = 0;
  fmean = 0.0;
  fmeanErr = 0.0;
  fmedian = 0.0;
  fmedianErr = 0.0;
  fpeak = 0.0;
  fpeakErr = 0.0;
  fquantile50 = 0.0;
  frecoVersion = "Dogwood"; //default to Dogwood
  //frecoVersion = "R1_18_2";
  fdriftQuantity1 = 1.0;
  fdriftQuantity2 = 1.0;
}

//____________________________________________________________________72
CalDrift::~CalDrift()
{
  LEA_DTOR; // Leak detection macro.
}


//____________________________________________________________________72
const Float_t
CalDrift::CalibratedPH(const Float_t ph,
                       const VldContext& cxNow,
                       const VldContext& cxRef) const
{
  // This code shouldn't be used except for test puposes.
  // The MuonDriftCalScheme does a much better job.
  // Nathaniel, Feb/06

  //set the task to 1=Dogwood
  //all previous constants where in the default task=0  
  Dbi::Task task = 1;
  
  DbiResultPtr<CalDrift> nowTable("CALDRIFT",
                                  cxNow,
                                  task,
                                  Dbi::kTableMissing,
                                  false);
  DbiResultPtr<CalDrift> t0Table("CALDRIFT",
                                 cxRef, 
                                 task,
                                 Dbi::kTableMissing,
                                 false);

  Double_t multipFactor = 1.0;
  if (nowTable.GetNumRows() && t0Table.GetNumRows()){
    const CalDrift* nowRow = nowTable.GetRow(0);
    const CalDrift* t0Row = t0Table.GetRow(0);
    multipFactor =
      t0Row->GetMedianAsR1_18_2(cxRef.GetDetector())
      /nowRow->GetMedianAsR1_18_2(cxRef.GetDetector());
//     multipFactor =
//       t0Row->GetMeanAsR1_18_2(cxRef.GetDetector())
//       /nowRow->GetMeanAsR1_18_2(cxRef.GetDetector());
//     multipFactor =
//       t0Row->GetMean()
//       /nowRow->GetMean();
  }

  return (Float_t) ((Double_t) ph)*multipFactor;
}

//____________________________________________________________________72
DbiTableRow* CalDrift::CreateTableRow() const
{
  return new CalDrift;
}

//____________________________________________________________________72
void CalDrift::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.

//Only using the dumb method for now.

//WARNING*******************************************
  //There's a nasty, nasty hack to force a RecoVersion flag in.
  //Make sure the version number in the CalDrift constructor is set to
  //the one you want before filling the database.
  //I'll fix this at the next reprocessing.
  //But I should be the only person who has to use this method for now.

  rs  >> fbackMedian >> fbackMedianErr >> ffrontMedian
      >> ffrontMedianErr >> fnonFidMedian >> fnonFidMedianErr
      >> fnumMuons >> fnumPlanes >> fmean >> fmeanErr >> fmedian
      >> fmedianErr >> fquantile50 >> frecoVersion
      >> fdriftQuantity1 >> fdriftQuantity2;
}

//____________________________________________________________________72
Float_t CalDrift::GetMeanAsR1_18_2(const Detector::Detector_t det)
  const
{
  Double_t recoFactor = 1.0;
  if (Detector::kNear == det){
    if ("R1_18"==frecoVersion){
      recoFactor = 1.00129; //R1_18_2 = r1_18Factor*R1_18;
    }
    if ("R1_18_2"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("R1_18_4"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("Cedar"==frecoVersion){
      recoFactor = 0.9719; //R1_18_2 = cedarFactor*cedar;
    }
    if ("Dogwood"==frecoVersion){
      recoFactor = 0.972141; //R1_18_2 = dogwoodFactor*Dogwood;
    }

  }
  if (Detector::kFar == det){
    if ("R1_18"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("R1_18_2"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("R1_18_4"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("Cedar"==frecoVersion){
      recoFactor = 0.92322; //R1_18_2 = cedarFactor*cedar;
    }
    if ("Dogwood"==frecoVersion){
      recoFactor = 1.01134; //R1_18_2 = dogwoodFactor*Dogwood;
    }
  }
  return fmean*recoFactor;
}

//____________________________________________________________________72
Float_t CalDrift::GetMedianRecoFactor(const Detector::Detector_t det)
  const
{
  //factors to scale to R1_18_2
  
  Double_t recoFactor = 1.0;
  if (Detector::kNear == det){
    if ("R1_18"==frecoVersion){
      recoFactor = 1.00129; //R1_18_2 = r1_18Factor*R1_18;
    }
    if ("R1_18_2"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("R1_18_4"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("Cedar"==frecoVersion){
      recoFactor = 1.0007; //R1_18_2 = cedarFactor*cedar;
    }
    if ("Dogwood"==frecoVersion){
      recoFactor = 0.950574; //R1_18_2 = dogwoodFactor*Dogwood;
    }
  }
  if (Detector::kFar == det){
    if ("R1_18"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("R1_18_2"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("R1_18_4"==frecoVersion){
      recoFactor = 1.0;
    }
    if ("Cedar"==frecoVersion){
      recoFactor = 0.96688; //R1_18_2 = cedarFactor*cedar;
    }
    if ("Dogwood"==frecoVersion){
      recoFactor = 0.964280; //R1_18_2 = dogwoodFactor*Dogwood;
    } 
  }
  return recoFactor;
}
//____________________________________________________________________72
Float_t CalDrift::GetMedianAsR1_18_2(const Detector::Detector_t det)
  const
{
  Float_t recoFactor=GetMedianRecoFactor(det);
  Float_t calc=fmedian;
  if("Dogwood"==frecoVersion && det==Detector::kNear)
    calc=ffrontMedian;
  return calc*recoFactor;
  
}
//____________________________________________________________________72
Float_t CalDrift::GetMedianErrAsR1_18_2(const Detector::Detector_t det)
  const
{
  Float_t recoFactor=GetMedianRecoFactor(det);
  Float_t calc=fmedianErr;
  if("Dogwood"==frecoVersion && det==Detector::kNear)
    calc=ffrontMedianErr;
  return calc*recoFactor;
}
//____________________________________________________________________72
void CalDrift::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 << fbackMedian << fbackMedianErr << ffrontMedian
      << ffrontMedianErr << fnonFidMedian << fnonFidMedianErr
      << fnumMuons << fnumPlanes << fmean << fmeanErr << fmedian
      << fmedianErr << fquantile50 << frecoVersion
      << fdriftQuantity1 << fdriftQuantity2;
}

---