/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/MagneticField/Engine/interface/MagneticFieldHelpers.h

Go to the documentation of this file.

#ifndef MagneticFieldHelpers_h
#define MagneticFieldHelpers_h

// #include "CondFormats/RunInfo/interface/RunInfo.h"
// #include "CondFormats/DataRecord/interface/RunSummaryRcd.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/EventSetup.h"

namespace magneticFieldHelpers  {

  int closerNominalField(float current) {
    int zeroFieldThreshold = 1000; //fixme
    float nominalCurrents[5] = {9558,14416,16819,18268,19262} ; //FIXME: replace with correct values...
    int nominalFields[5] = {20,30,35,38,40} ; //in kGauss
    if(current < zeroFieldThreshold) return 0;
    int i=0;
    for(;i<4;i++)
      {
        if(2*current < nominalCurrents[i]+nominalCurrents[i+1] )
          return nominalFields[i];
      }
    return nominalFields[i];
  }

//   /// Return the closer nominal field value (kGauss) to the average current stored in the DB.
//   int closerNominalField(const RunInfo & runInfo) {
//     return closerNominalField(runInfo.m_avg_current);
//   }

//   /// Return the closer nominal field value (kGauss) to the average current stored in the DB.
//   int closerNominalField(const edm::EventSetup& es)
//   {
//     edm::ESHandle<RunInfo> sum;
//     es.get<RunInfoRcd>().get(sum);
//     return closerNominalField( *sum.product() );
//   }


}

#endif // MagneticFieldHelpers_h