/data/refman/pasoursint/CMSSW_4_2_9_HLT1_bphpatch4/src/SimMuon/CSCDigitizer/src/CSCConfigurableStripConditions.cc

Go to the documentation of this file.

#include "SimMuon/CSCDigitizer/src/CSCConfigurableStripConditions.h"
#include "DataFormats/MuonDetId/interface/CSCDetId.h"
#include "Geometry/CSCGeometry/interface/CSCChamberSpecs.h"
#include "FWCore/Utilities/interface/Exception.h"

CSCConfigurableStripConditions::CSCConfigurableStripConditions(const edm::ParameterSet & p)
: theGain( p.getParameter<double>("gain") ),
  theME11Gain( p.getParameter<double>("me11gain") ),
  theGainSigma( p.getParameter<double>("ampGainSigma") ),
  thePedestal( p.getParameter<double>("pedestal") ),
  thePedestalSigma( p.getParameter<double>("pedestalSigma") ),
  theCapacitiveCrosstalk(0.0167),
  theResistiveCrosstalk(0.02)
{
  theNoisifiers.resize(9);
  makeNoisifier(1, p.getParameter<std::vector<double> >("me11a") );
  makeNoisifier(2, p.getParameter<std::vector<double> >("me11") );
  makeNoisifier(3, p.getParameter<std::vector<double> >("me12") );
  makeNoisifier(4, p.getParameter<std::vector<double> >("me13") ); // not sure about this one
  makeNoisifier(5, p.getParameter<std::vector<double> >("me21") );
  makeNoisifier(6, p.getParameter<std::vector<double> >("me22") );
  makeNoisifier(7, p.getParameter<std::vector<double> >("me31") );
  makeNoisifier(8, p.getParameter<std::vector<double> >("me32") );
  makeNoisifier(9, p.getParameter<std::vector<double> >("me31") ); // for lack of a better idea
}


CSCConfigurableStripConditions::~CSCConfigurableStripConditions()
{
  for(int i = 0; i < 9; ++i)
  {
    delete theNoisifiers[i];
  }
}


float CSCConfigurableStripConditions::gain(const CSCDetId & detId, int channel) const
{
  if(detId.station() == 1 && (detId.ring() == 1 || detId.ring() == 4) )
  {
    return theME11Gain;
  }  
  else
  {
    return theGain;
  }
}


void CSCConfigurableStripConditions::fetchNoisifier(const CSCDetId & detId, int istrip)
{
  //TODO get this moved toCSCDetId
  int chamberType = CSCChamberSpecs::whatChamberType(detId.station(), detId.ring());
  theNoisifier = theNoisifiers[chamberType-1];
}


void CSCConfigurableStripConditions::makeNoisifier(int chamberType, const std::vector<double> & correlations)
{

  // format is 33, 34, 44, 35, 45, 55
  //           46, 56, 66, 57, 67, 77
  if(correlations.size() != 12)
  {
    throw cms::Exception("CSCConfigurableStripConditions")
      << "Expect 12 noise correlation coefficients, but got "
      << correlations.size();
  }

  CSCCorrelatedNoiseMatrix matrix;
  matrix(3,3) = correlations[0];
  matrix(3,4) = correlations[1];
  matrix(4,4) = correlations[2];
  matrix(3,5) = correlations[3];
  matrix(4,5) = correlations[4];
  matrix(5,5) = correlations[5];
  matrix(4,6) = correlations[6];
  matrix(5,6) = correlations[7];
  matrix(6,6) = correlations[8];
  matrix(5,7) = correlations[9];
  matrix(6,7) = correlations[10];
  matrix(7,7) = correlations[11];

  // since I don't know how to correlate the pedestal samples,
  // take as constant
  double scaVariance = 2. * thePedestalSigma * thePedestalSigma;
  matrix(0,0) = scaVariance;
  matrix(1,1) = scaVariance;
  matrix(2,2) = scaVariance;
  theNoisifiers[chamberType-1] = new CSCCorrelatedNoisifier(matrix);

}

void CSCConfigurableStripConditions::crosstalk(const CSCDetId&detId, int channel,
                 double stripLength, bool leftRight,
                 float & capacitive, float & resistive) const
{
  capacitive = theCapacitiveCrosstalk * stripLength;
  resistive = theResistiveCrosstalk;
}