CSCRecoConditions Class Reference

#include <CSCRecoConditions.h>

Public Member Functions
float	anodeBXoffset (const CSCDetId &detId) const
float	averageGain () const
	return average gain over entire CSC system More...
bool	badStrip (const CSCDetId &id, int geomStrip, int nstrips) const
	Is the strip bad? More...
const std::bitset< 112 > &	badWireWord (const CSCDetId &id) const
	Get bad wiregroup word. More...
float	chamberTimingCorrection (const CSCDetId &id) const
float	chipCorrection (const CSCDetId &detId, int channel) const
	All other functions are accessed by geometrical strip label (i.e. strip number according to local coordinates) More...
void	crossTalk (const CSCDetId &id, int centralStrip, std::vector< float > &xtalks) const
	CSCRecoConditions (const edm::ParameterSet &pset)
void	fillBadChannelWords (const CSCDetId &id)
	fill bad strip & bad wiregroup bitsets from conditions data More...
float	gain (const CSCDetId &id, int geomStrip) const
	channels and geomstrips count from 1 More...
float	gasGainCorrection (const CSCDetId &id, int strip, int wireGroup) const
	returns gas-gain correction More...
void	initializeEvent (const edm::EventSetup &es)
	fetch the cond data from the database More...
bool	nearBadStrip (const CSCDetId &id, int geomStrip, int nstrips) const
	Is a neighbour bad? More...
void	noiseMatrix (const CSCDetId &id, int centralStrip, std::vector< float > &nme) const
float	pedestal (const CSCDetId &id, int channel) const
	static pedestal in ADC counts for strip channel (e.g. 1-16 for ganged ME1a, 1-48 for unganged ME1a) More...
float	pedestalSigma (const CSCDetId &id, int channel) const
	sigma of static pedestal in ADC counts for strip channel (e.g. 1-16 for ganged ME1a, 1-48 for unganged ME1a) More...
void	stripWeights (const CSCDetId &id, short int nstrips, float *weights) const
	~CSCRecoConditions ()

Private Member Functions
float	stripWeight (const CSCDetId &id, int geomStrip) const
	return gain weight for given strip channel More...

Private Attributes
CSCConditions	theConditions

Detailed Description

Wrap CSCConditions class for use in CSC local reconstruction, in analogy with wrapper classes Rick uses in CSCDigitizer.

CSCConditions encapsulates the conditions data (e.g. calibration data) from the database and presents it as CSCRecHitD requires (and that is somewhat historical!)

All functions in public interface accept CSCDetId for ME1A (i.e. ring 4) with channel number 1-16 (and not the raw ME11 channel 65-80).

Author

Tim Cox - UC Davis

Definition at line 24 of file CSCRecoConditions.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 6 of file CSCRecoConditions.cc.

                                                                 : theConditions( ps ) {
}

CSCRecoConditions::~CSCRecoConditions

(

)

Definition at line 9 of file CSCRecoConditions.cc.

                                      {
}

Member Function Documentation

float CSCRecoConditions::anodeBXoffset

(

const CSCDetId &

detId

)

const

Definition at line 171 of file CSCRecoConditions.cc.

References CSCConditions::anodeBXoffset(), and theConditions.

Referenced by averageGain(), and CSCMake2DRecHit::findWireBx().

                                                                {
  return theConditions.anodeBXoffset(id);
}

float CSCRecoConditions::averageGain ( ) const

inline

return average gain over entire CSC system

Definition at line 41 of file CSCRecoConditions.h.

References anodeBXoffset(), CSCConditions::averageGain(), badStrip(), badWireWord(), chamberTimingCorrection(), chipCorrection(), crossTalk(), fillBadChannelWords(), gasGainCorrection(), nearBadStrip(), noiseMatrix(), pedestal(), pedestalSigma(), digi_MixPreMix_cfi::strip, stripWeight(), stripWeights(), theConditions, and create_public_lumi_plots::weights.

Referenced by stripWeight().

                            {
     return theConditions.averageGain(); }

bool CSCRecoConditions::badStrip	(	const CSCDetId &	id,
		int	geomStrip,
		int	nstrips
	)		const

Is the strip bad?

Is strip itself a bad strip?

Definition at line 138 of file CSCRecoConditions.cc.

References CSCConditions::badStripWord(), CSCConditions::channelFromStrip(), CSCConditions::idOfBadChannelWords(), CSCConditions::rawStripChannel(), and theConditions.

Referenced by CSCRecoBadChannelsAnalyzer::analyze(), averageGain(), CSCHitFromStripOnly::isDeadStrip(), and nearBadStrip().

                                                                                       {

  // input nstrips is no. of strips in the layer (could get this from CSCChamberSpecs or CSCLayerGeometry
  // but then need a CSCLayer*)

  bool bad = true; // if geomStrip out of range, call strip bad

  if ( id != theConditions.idOfBadChannelWords() ) {
    bad = false; // if bad channel words for this id not filled, call strip good
    return bad;
  }

  if(geomStrip>0 && geomStrip<=nstrips){
    bad = false;  // default to good
    int geomChan = theConditions.channelFromStrip( id, geomStrip );
    int rawChan = theConditions.rawStripChannel( id, geomChan );
    if( rawChan>0 && rawChan<113 ){
      const std::bitset<112>& badStrips = theConditions.badStripWord();
      bad = badStrips.test(rawChan-1); // 112 bits max, labelled 0-111.
    }
  }
  return bad;
}

const std::bitset< 112 > & CSCRecoConditions::badWireWord

(

const CSCDetId &

id

)

const

Get bad wiregroup word.

Definition at line 163 of file CSCRecoConditions.cc.

References CSCConditions::badWireWord(), and theConditions.

Referenced by averageGain(), and CSCHitFromWireOnly::isDeadWG().

                                                                               {
    return theConditions.badWireWord();
}

float CSCRecoConditions::chamberTimingCorrection

(

const CSCDetId &

id

)

const

Definition at line 167 of file CSCRecoConditions.cc.

References CSCConditions::chamberTimingCorrection(), and theConditions.

Referenced by averageGain(), and CSCMake2DRecHit::hitFromStripAndWire().

                                                                          {
  return theConditions.chamberTimingCorrection(id);
}

float CSCRecoConditions::chipCorrection	(	const CSCDetId &	detId,
		int	channel
	)		const

All other functions are accessed by geometrical strip label (i.e. strip number according to local coordinates)

Definition at line 40 of file CSCRecoConditions.cc.

References CSCConditions::channelFromStrip(), CSCConditions::chipCorrection(), and theConditions.

Referenced by averageGain(), and CSCMake2DRecHit::hitFromStripAndWire().

                                                                                {
  // geometric strip to geometric channel (e.g. ME1a, 1-48->1-16 ganged or 1-48 unganged)
  int geomChannel = theConditions.channelFromStrip( id, geomStrip );
  return theConditions.chipCorrection(id, geomChannel);
}

void CSCRecoConditions::crossTalk	(	const CSCDetId &	id,
		int	centralStrip,
		std::vector< float > &	xtalks
	)		const

fill crosstalk information for 3 neighbouring strips as linear vector (must be allocated by caller) Note that centralStrip is a 'geomStrip' and e.g. always ranges 1-48 in ME1a.

Definition at line 114 of file CSCRecoConditions.cc.

References CSCConditions::channelFromStrip(), CSCConditions::crossTalk(), mps_fire::i, and theConditions.

Referenced by averageGain(), and CSCXonStrip_MatchGatti::findXOnStrip().

                                                                                                    {

  // xtalks will be filled with crosstalk for geomStrip and its immediate neighbours

  xtalks.clear();

  for ( short int i = geomStrip-1; i < geomStrip+2; ++i) {
    int geomChannel = theConditions.channelFromStrip( id, i );
    std::vector<float> ct(4);
    theConditions.crossTalk(id, geomChannel, ct);
    xtalks.push_back(ct[0]);
    xtalks.push_back(ct[1]);
    xtalks.push_back(ct[2]);
    xtalks.push_back(ct[3]);
  }
}

void CSCRecoConditions::fillBadChannelWords

(

const CSCDetId &

id

)

fill bad strip & bad wiregroup bitsets from conditions data

Fill badStrip word and badWire word for given CSC layer, id.

Definition at line 17 of file CSCRecoConditions.cc.

References CSCConditions::fillBadChannelWords(), and theConditions.

Referenced by CSCRecoBadChannelsAnalyzer::analyze(), averageGain(), and CSCRecHitDBuilder::build().

                                                                {
  theConditions.fillBadChannelWords ( id );
}

float CSCRecoConditions::gain	(	const CSCDetId &	id,
		int	geomStrip
	)		const

channels and geomstrips count from 1

return gain for given strip

Definition at line 33 of file CSCRecoConditions.cc.

References CSCConditions::gain(), LogTrace, and theConditions.

Referenced by stripWeight().

                                                                       {
  LogTrace("CSCRecoConditions") << id << " geomChannel " <<  geomChannel << " gain " << theConditions.gain(id, geomChannel);
  return theConditions.gain(id, geomChannel);
}

float CSCRecoConditions::gasGainCorrection	(	const CSCDetId &	id,
		int	strip,
		int	wireGroup
	)		const

returns gas-gain correction

Definition at line 175 of file CSCRecoConditions.cc.

References CSCConditions::channelFromStrip(), CSCConditions::gasGainCorrection(), and theConditions.

Referenced by averageGain(), and CSCMake2DRecHit::hitFromStripAndWire().

                                                                                                   {
  int geomChannel = theConditions.channelFromStrip( id, geomStrip);
  return theConditions.gasGainCorrection(id, geomChannel, wiregroup);
}

void CSCRecoConditions::initializeEvent

(

const edm::EventSetup &

es

)

fetch the cond data from the database

Definition at line 12 of file CSCRecoConditions.cc.

References CSCConditions::initializeEvent(), and theConditions.

Referenced by CSCRecoBadChannelsAnalyzer::analyze(), and CSCRecHitDProducer::produce().

                                                                 {
  theConditions.initializeEvent( es );
}

bool CSCRecoConditions::nearBadStrip	(	const CSCDetId &	id,
		int	geomStrip,
		int	nstrips
	)		const

Is a neighbour bad?

Is an immediate neighbour a bad strip?

Definition at line 132 of file CSCRecoConditions.cc.

References badStrip().

Referenced by CSCRecoBadChannelsAnalyzer::analyze(), averageGain(), and CSCHitFromStripOnly::isNearDeadStrip().

                                                                                           {
  bool nearBad = ( badStrip(id,geomStrip-1, nstrips) || badStrip(id,geomStrip+1, nstrips) );
  return nearBad;
}

void CSCRecoConditions::noiseMatrix	(	const CSCDetId &	id,
		int	centralStrip,
		std::vector< float > &	nme
	)		const

fill expanded noise matrix for 3 neighbouring strips as linear vector (must be allocated by caller) Note that centralStrip is a 'geomStrip' and ranges 1-48 in ME1a.

Definition at line 71 of file CSCRecoConditions.cc.

References CSCConditions::channelFromStrip(), HTMLExport::elem(), mps_fire::i, gen::k, funct::m, CSCConditions::noiseMatrixElements(), stripWeight(), theConditions, and w.

Referenced by averageGain(), and CSCXonStrip_MatchGatti::findXOnStrip().

                                                                                                        {

  // nMatrix will be filled with expanded noise matrix elements for strip 'geomStrip' and its immediate neighbours

  nMatrix.clear();

  // These are ME1/2 constants as fall-back
  const float fakeme12[15] = {8.64, 3.47, 2.45, 8.60, 3.28, 1.88, 8.61, 3.18, 1.99, 7.67, 2.64, 0., 7.71, 0., 0.};

  float elem[15];

  for ( short int i = geomStrip-1; i < geomStrip+2; ++i) {

    std::vector<float> me(12);

    float w = stripWeight(id, i);
    w = w*w;
    int geomChannel = theConditions.channelFromStrip( id, i );
    theConditions.noiseMatrixElements(id, geomChannel, me);
    for ( short int j=0; j<11; ++j ) {
      elem[j] = me[j] * w;
    }
    elem[11]= 0.;
    elem[12]= me[11] * w;
    elem[13]= 0.;
    elem[14]= 0.;

    // Test that elements make sense:
    bool isFlawed = false;
    for ( short int k = 0; k < 15; ++k) {
      if (elem[k] < 0.001) elem[k] = 0.001; // fix if too small...
      if (elem[k] > 50.) isFlawed = true; // fail if too big...
    }

    if ( isFlawed ) {
      // These are fake ME1/2:
      for ( short int m = 0; m < 15; ++m ) { elem[m] = fakeme12[m]; }
    }

    for (int k = 0; k < 15; ++k) { nMatrix.push_back( elem[k] ); }
  }
}

float CSCRecoConditions::pedestal	(	const CSCDetId &	id,
		int	geomChannel
	)		const

static pedestal in ADC counts for strip channel (e.g. 1-16 for ganged ME1a, 1-48 for unganged ME1a)

gains & pedestals are requested by geometric channel (as in CSCStripDigi- e.g. 1-16 for ganged ME1a, and with any readout flips already removed)

Definition at line 24 of file CSCRecoConditions.cc.

References LogTrace, CSCConditions::pedestal(), and theConditions.

Referenced by averageGain(), and CSCStaticPedestal::pedestal().

                                                                           {
  LogTrace("CSCRecoConditions") << id << " geomChannel " << geomChannel << " pedestal " << theConditions.pedestal(id, geomChannel);
  return theConditions.pedestal(id, geomChannel);
}

float CSCRecoConditions::pedestalSigma	(	const CSCDetId &	id,
		int	channel
	)		const

sigma of static pedestal in ADC counts for strip channel (e.g. 1-16 for ganged ME1a, 1-48 for unganged ME1a)

Definition at line 29 of file CSCRecoConditions.cc.

References CSCConditions::pedestalSigma(), and theConditions.

Referenced by averageGain().

                                                                                {
  return theConditions.pedestalSigma(id, geomChannel);
}

float CSCRecoConditions::stripWeight ( const CSCDetId &  id, int  geomStrip  ) const

private

return gain weight for given strip channel

Definition at line 61 of file CSCRecoConditions.cc.

References averageGain(), CSCConditions::channelFromStrip(), gain(), LogTrace, theConditions, and w.

Referenced by averageGain(), noiseMatrix(), and stripWeights().

                                                                              {
   int geomChannel = theConditions.channelFromStrip( id, geomStrip );
   float w = averageGain() / gain(id, geomChannel); // averageGain() from CSCConditions
   // Weights are forced to lie within 0.5 and 1.5
   if (w > 1.5) w = 1.5;
   if (w < 0.5) w = 0.5;
   LogTrace("CSCRecoConditions") << id << " geomStrip " << geomStrip << " stripWeight " << w;
   return w;
}

void CSCRecoConditions::stripWeights	(	const CSCDetId &	id,
		short int	nstrips,
		float *	weights
	)		const

calculate gain weights for all strips in a CSC layer (total in layer = nstrips) this is averageGain()/gain for each strip filled into a C-array which caller must have allocated. values are constrained to be in [0.5, 1.5]

Definition at line 51 of file CSCRecoConditions.cc.

References mps_fire::i, and stripWeight().

Referenced by averageGain(), and CSCHitFromStripOnly::runStrip().

                                                                                                  {

  for ( short int i = 1; i < nstrips+1; ++i) {
      weights[i-1] = stripWeight(id, i) ;
  }
}

Member Data Documentation

CSCConditions CSCRecoConditions::theConditions

private

Definition at line 93 of file CSCRecoConditions.h.

Referenced by anodeBXoffset(), averageGain(), badStrip(), badWireWord(), chamberTimingCorrection(), chipCorrection(), crossTalk(), fillBadChannelWords(), gain(), gasGainCorrection(), initializeEvent(), noiseMatrix(), pedestal(), pedestalSigma(), and stripWeight().