#include <CSCDbStripConditions.h>

Inheritance diagram for CSCDbStripConditions:

Public Member Functions
virtual void	crosstalk (const CSCDetId &detId, int channel, double stripLength, bool leftRight, float &capacitive, float &resistive) const

	CSCDbStripConditions (const edm::ParameterSet &pset)

virtual float	gain (const CSCDetId &detId, int channel) const
	channels count from 1 More...

virtual float	gainSigma (const CSCDetId &detId, int channel) const
	total calibration precision More...

virtual void	initializeEvent (const edm::EventSetup &es)
	fetch the maps from the database More...

virtual bool	isInBadChamber (const CSCDetId &id) const
	check list of bad chambers from db More...

virtual float	pedestal (const CSCDetId &detId, int channel) const
	in ADC counts More...

virtual float	pedestalSigma (const CSCDetId &detId, int channel) const

virtual	~CSCDbStripConditions ()

Public Member Functions inherited from CSCStripConditions
float	analogNoise (const CSCDetId &detId, int channel) const
	calculated from pedestalSigma & gain More...

	CSCStripConditions ()

void	noisify (const CSCDetId &detId, CSCAnalogSignal &signal)
	superimposes noise, in fC, on the signal More...

void	setRandomEngine (CLHEP::HepRandomEngine &engine)

virtual float	smearedGain (const CSCDetId &detId, int channel) const

virtual	~CSCStripConditions ()

Private Member Functions
virtual void	fetchNoisifier (const CSCDetId &detId, int istrip)

Private Attributes
bool	doCorrelatedNoise_

float	theCapacitiveCrosstalk

CSCConditions	theConditions

float	theGainsConstant

float	theResistiveCrosstalkScaling

Additional Inherited Members
Public Types inherited from CSCStripConditions
typedef math::ErrorD< 8 >::type	CSCCorrelatedNoiseMatrix

typedef CorrelatedNoisifier < CSCCorrelatedNoiseMatrix >	CSCCorrelatedNoisifier

Protected Attributes inherited from CSCStripConditions
CSCCorrelatedNoisifier *	theNoisifier

CLHEP::RandGaussQ *	theRandGaussQ

Detailed Description

Definition at line 8 of file CSCDbStripConditions.h.

Constructor & Destructor Documentation

CSCDbStripConditions::CSCDbStripConditions ( const edm::ParameterSet & pset )

explicit

Definition at line 12 of file CSCDbStripConditions.cc.

 : CSCStripConditions(),
   theConditions( pset ),
   theCapacitiveCrosstalk(pset.getParameter<double>("capacativeCrosstalk")),
   theResistiveCrosstalkScaling(pset.getParameter<double>("resistiveCrosstalkScaling")),
   theGainsConstant(pset.getParameter<double>("gainsConstant")),
   doCorrelatedNoise_(pset.getParameter<bool>("doCorrelatedNoise"))
 {
 //  theCapacitiveCrosstalk = = 1/maxslope/maxsignal) = 1/ (0.00231/0.143);
 // Howoever, need a bit more.  Maybe the slope gets smeared?
 }

CSCDbStripConditions::~CSCDbStripConditions ( )

virtual

Definition at line 25 of file CSCDbStripConditions.cc.

References CSCStripConditions::theNoisifier.

 {
   if(theNoisifier != 0) delete theNoisifier;
 }

Member Function Documentation

void CSCDbStripConditions::crosstalk	(	const CSCDetId &	detId,
		int	channel,
		double	stripLength,
		bool	leftRight,
		float &	capacitive,
		float &	resistive
	)		const

virtual

Implements CSCStripConditions.

Definition at line 65 of file CSCDbStripConditions.cc.

References CSCConditions::crosstalkIntercept(), CSCConditions::crosstalkSlope(), CSCChannelTranslator::rawCSCDetId(), CSCChannelTranslator::rawStripChannel(), slope, theCapacitiveCrosstalk, theConditions, and theResistiveCrosstalkScaling.

 {
   CSCChannelTranslator translate;
   CSCDetId idraw = translate.rawCSCDetId( detId );
   int iraw = translate.rawStripChannel( detId, channel );
     
   resistive = theConditions.crosstalkIntercept(idraw, iraw, leftRight)
              * theResistiveCrosstalkScaling;
   float slope = theConditions.crosstalkSlope(idraw, iraw, leftRight);
   // ns before the peak where slope is max
   float maxSlopeTime = 60.; 
   // some confusion about +/-
   float capacitiveFraction = fabs(slope)*maxSlopeTime;
   // theCapacitiveCrosstalk is the number needed for 100% xtalk, so
   capacitive = theCapacitiveCrosstalk * capacitiveFraction;
 }

void CSCDbStripConditions::fetchNoisifier	(	const CSCDetId &	detId,
		int	istrip
	)

privatevirtual

Implements CSCStripConditions.

Definition at line 86 of file CSCDbStripConditions.cc.

References doCorrelatedNoise_, makeMuonMisalignmentScenario::matrix, CSCConditions::noiseMatrixElements(), pedestalSigma(), CSCChannelTranslator::rawCSCDetId(), CSCChannelTranslator::rawStripChannel(), theConditions, and CSCStripConditions::theNoisifier.

 {
   CSCChannelTranslator translate;
   CSCDetId idraw = translate.rawCSCDetId( detId );
   int iraw = translate.rawStripChannel( detId, istrip );
     
   std::vector<float> me(12); // buffer for matrix elements
   theConditions.noiseMatrixElements( idraw, iraw, me ); // fill it
 
   CSCCorrelatedNoiseMatrix matrix;
   //TODO get the pedestals right
   matrix(2,2) = me[0]; // item.elem33;
   matrix(3,3) = me[3]; // item.elem44;
   matrix(4,4) = me[6]; // item.elem55;
   matrix(5,5) = me[9]; // item.elem66;
   matrix(6,6) = me[11]; // item.elem77;
   
   if(doCorrelatedNoise_)
   {
     matrix(2,3) = me[1]; // item.elem34;
     matrix(2,4) = me[2]; // item.elem35;
     matrix(3,4) = me[4]; // item.elem45;
     matrix(3,5) = me[5]; // item.elem46;
     matrix(4,5) = me[7]; // item.elem56;
     matrix(4,6) = me[8]; // item.elem57;
     matrix(5,6) = me[10]; // item.elem67;
   }
 
   // the other diagonal elements can just come from the pedestal sigma
   float sigma = pedestalSigma(detId, istrip);
   //@@  float scaVariance = 2 * sigma * sigma;
   //@@ The '2 *' IS strictly correct, but currently the value in the cond db is 2x too large since
   //@@ it is the rms of the distribution of pedestals of all 8 time samples rather than the rms of
   //@@ the average of the first two time samples
   float scaVariance = sigma * sigma;
   matrix(0,0) = matrix(1,1) = matrix(7,7) = scaVariance;
 
   // unknown neighbors can be the average of the known neighbors
   //float avgNeighbor = (matrix(2,3)+matrix(3,4)+matrix(4,5)+matrix(5,6))/4.;
   //float avg2away = (matrix(2,4)+matrix(3,5)+matrix(4,6))/3.;
   //matrix(0,1) = matrix(1,2) = matrix(6,7) = avgNeighbor;
   //matrix(0,2) = matrix(1,3) = matrix(5,7) = avg2away;
 
   if(theNoisifier != 0) delete theNoisifier;
   theNoisifier = new CSCCorrelatedNoisifier(matrix);
 }

float CSCDbStripConditions::gain	(	const CSCDetId &	detId,
		int	channel
	)		const

virtual

channels count from 1

Implements CSCStripConditions.

Definition at line 37 of file CSCDbStripConditions.cc.

References CSCConditions::gain(), CSCChannelTranslator::rawCSCDetId(), CSCChannelTranslator::rawStripChannel(), theConditions, and theGainsConstant.

 {
   CSCChannelTranslator translate;
   CSCDetId idraw = translate.rawCSCDetId( detId );
   int iraw = translate.rawStripChannel( detId, channel );
   return theConditions.gain(idraw, iraw)  * theGainsConstant;
 }

virtual float CSCDbStripConditions::gainSigma	(	const CSCDetId &	detId,
		int	channel
	)		const

inlinevirtual

total calibration precision

Implements CSCStripConditions.

Definition at line 20 of file CSCDbStripConditions.h.

20 {return 0.005;}

void CSCDbStripConditions::initializeEvent ( const edm::EventSetup & es )

virtual

fetch the maps from the database

Reimplemented from CSCStripConditions.

Definition at line 31 of file CSCDbStripConditions.cc.

References CSCConditions::initializeEvent(), and theConditions.

 {
   theConditions.initializeEvent(es);
 }

bool CSCDbStripConditions::isInBadChamber ( const CSCDetId & id ) const

virtual

check list of bad chambers from db

Reimplemented from CSCStripConditions.

Definition at line 133 of file CSCDbStripConditions.cc.

References CSCConditions::isInBadChamber(), and theConditions.

 {
   return theConditions.isInBadChamber( detId );
 }

float CSCDbStripConditions::pedestal	(	const CSCDetId &	detId,
		int	channel
	)		const

virtual

in ADC counts

Implements CSCStripConditions.

Definition at line 47 of file CSCDbStripConditions.cc.

References CSCConditions::pedestal(), CSCChannelTranslator::rawCSCDetId(), CSCChannelTranslator::rawStripChannel(), and theConditions.

 {
   CSCChannelTranslator translate;
   CSCDetId idraw = translate.rawCSCDetId( detId );
   int iraw = translate.rawStripChannel( detId, channel );
   return theConditions.pedestal(idraw, iraw);
 }

float CSCDbStripConditions::pedestalSigma	(	const CSCDetId &	detId,
		int	channel
	)		const

virtual

Implements CSCStripConditions.

Definition at line 56 of file CSCDbStripConditions.cc.

References CSCConditions::pedestalSigma(), CSCChannelTranslator::rawCSCDetId(), CSCChannelTranslator::rawStripChannel(), and theConditions.

Referenced by fetchNoisifier().

 {
   CSCChannelTranslator translate;
   CSCDetId idraw = translate.rawCSCDetId( detId );
   int iraw = translate.rawStripChannel( detId, channel );
   return theConditions.pedestalSigma(idraw, iraw);
 }

Member Data Documentation

bool CSCDbStripConditions::doCorrelatedNoise_

private

Definition at line 45 of file CSCDbStripConditions.h.

Referenced by fetchNoisifier().

float CSCDbStripConditions::theCapacitiveCrosstalk

private

Definition at line 39 of file CSCDbStripConditions.h.

Referenced by crosstalk().

CSCConditions CSCDbStripConditions::theConditions

private

Definition at line 36 of file CSCDbStripConditions.h.

Referenced by crosstalk(), fetchNoisifier(), gain(), initializeEvent(), isInBadChamber(), pedestal(), and pedestalSigma().

float CSCDbStripConditions::theGainsConstant

private

Definition at line 44 of file CSCDbStripConditions.h.

Referenced by gain().

float CSCDbStripConditions::theResistiveCrosstalkScaling

private

Definition at line 42 of file CSCDbStripConditions.h.

Referenced by crosstalk().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation