TT6CommonModeCalculator Class Reference

#include <TT6CommonModeCalculator.h>

Inheritance diagram for TT6CommonModeCalculator:

Public Member Functions
TkCommonMode *	commonMode () override
ApvAnalysis::PedestalType	doIt (const ApvAnalysis::PedestalType &) override
float	getCMSlope () override
void	newEvent () override
void	setCM (TkCommonMode *in) override
void	setCM (const std::vector< float > &in) override
	TT6CommonModeCalculator (TkNoiseCalculator *noise_calc, TkApvMask *mask_calc, float sig_cut)
	~TT6CommonModeCalculator () override
Public Member Functions inherited from TkCommonModeCalculator
virtual	~TkCommonModeCalculator ()

Protected Member Functions
void	calculateCMSlope (ApvAnalysis::PedestalType &)
void	calculateCommonMode (ApvAnalysis::PedestalType &)

Protected Attributes
bool	alreadyUsedEvent
float	cutToAvoidSignal
float	slope
TkApvMask *	theApvMask
std::vector< float >	theCommonModeValues
TkNoiseCalculator *	theNoiseCalculator
TkCommonMode *	theTkCommonMode

Detailed Description

Concrete implementation of TkCommonModeCalculator for TT6.

Definition at line 11 of file TT6CommonModeCalculator.h.

Constructor & Destructor Documentation

TT6CommonModeCalculator()

TT6CommonModeCalculator::TT6CommonModeCalculator	(	TkNoiseCalculator *	noise_calc,
		TkApvMask *	mask_calc,
		float	sig_cut
	)

Definition at line 5 of file TT6CommonModeCalculator.cc.

References gather_cfg::cout, and cutToAvoidSignal.

    : theNoiseCalculator(noise_calc), theApvMask(mask_calc), alreadyUsedEvent(false) {
  if (false)
    cout << "Constructing TT6CommonMode Calculator ..." << endl;
  cutToAvoidSignal = sig_cut;
}

~TT6CommonModeCalculator()

TT6CommonModeCalculator::~TT6CommonModeCalculator ( )

override

Definition at line 14 of file TT6CommonModeCalculator.cc.

References gather_cfg::cout.

                                                  {
  if (false)
    cout << "Destructing TT6CommonModeCalculator " << endl;
}

Member Function Documentation

calculateCMSlope()

void TT6CommonModeCalculator::calculateCMSlope ( ApvAnalysis::PedestalType &  indat )

protected

Definition at line 79 of file TT6CommonModeCalculator.cc.

References alignCSCRings::s, slope, and jetUpdater_cfi::sort.

Referenced by calculateCommonMode().

                                                                             {
  if (indat.size() != 128) {
    slope = -100.0;
    return;
  }
  ApvAnalysis::PedestalType diffVec;
  diffVec.clear();
  for (int s = 0; s < 64; s++)
    diffVec.push_back(indat[s + 64] - indat[s]);
  std::sort(diffVec.begin(), diffVec.end());
  slope = (diffVec[31] + diffVec[32]) / 2. / 64.;
}

calculateCommonMode()

void TT6CommonModeCalculator::calculateCommonMode ( ApvAnalysis::PedestalType &  indat )

protected

Definition at line 39 of file TT6CommonModeCalculator.cc.

References alreadyUsedEvent, calculateCMSlope(), cutToAvoidSignal, TkCommonModeTopology::finalStrips(), mps_fire::i, TkCommonModeTopology::initialStrips(), dqmiolumiharvest::j, TkApvMask::mask(), TkNoiseCalculator::noise(), TkCommonModeTopology::numberOfSets(), TkApvMask::ok, setCM(), theApvMask, theCommonModeValues, theNoiseCalculator, theTkCommonMode, and TkCommonMode::topology().

Referenced by doIt().

                                                                                {
  if (alreadyUsedEvent == false) {
    alreadyUsedEvent = true;
    //  cout<< "I am inside the calculateCommonMode"<<endl;
    TkApvMask::MaskType strip_mask = theApvMask->mask();
    ApvAnalysis::PedestalType strip_noise = theNoiseCalculator->noise();
    theCommonModeValues.clear();

    if (!strip_noise.empty()) {
      int nSet = theTkCommonMode->topology().numberOfSets();
      for (int i = 0; i < nSet; i++) {
        int initial = theTkCommonMode->topology().initialStrips()[i];
        int final = theTkCommonMode->topology().finalStrips()[i];
        double sumVal = 0.0;
        double sumWt = 0.0;
        for (int j = initial; j <= final; j++) {
          if (strip_mask[j] == TkApvMask::ok) {
            if (fabs(indat[j]) < cutToAvoidSignal * strip_noise[j]) {
              double nWeight = 1 / (strip_noise[j] * strip_noise[j]);
              sumVal += (indat[j] * nWeight);
              sumWt += nWeight;
            }
          }
        }
        double avVal = (sumWt) ? sumVal / sumWt : 0.0;
        theCommonModeValues.push_back(static_cast<float>(avVal));
        //cout <<"Setting CM values"<<endl;
      }
    }
  }
  TT6CommonModeCalculator::setCM(theCommonModeValues);
  calculateCMSlope(indat);
}

commonMode()

TkCommonMode* TT6CommonModeCalculator::commonMode ( )

inlineoverridevirtual

Get CM value

Implements TkCommonModeCalculator.

Definition at line 21 of file TT6CommonModeCalculator.h.

References theTkCommonMode.

{ return theTkCommonMode; }

doIt()

ApvAnalysis::PedestalType TT6CommonModeCalculator::doIt ( const ApvAnalysis::PedestalType &  )

overridevirtual

Return CM-subtracted data in APV

Implements TkCommonModeCalculator.

Definition at line 21 of file TT6CommonModeCalculator.cc.

References calculateCommonMode(), mps_fire::i, MillePedeFileConverter_cfg::out, TkCommonModeTopology::setOfStrip(), theCommonModeValues, theTkCommonMode, and TkCommonMode::topology().

Referenced by MassReplace.MassSearchReplaceAnyInputTagVisitor::doIt(), ConfigBuilder.ConfigBuilder.MassSearchReplaceProcessNameVisitor::doIt(), MassReplace.MassSearchReplaceAnyInputTagVisitor::enter(), MassReplace.MassSearchReplaceParamVisitor::enter(), and ConfigBuilder.ConfigBuilder.MassSearchReplaceProcessNameVisitor::enter().

                                                                                           {
  ApvAnalysis::PedestalType indat = _indat;
  ApvAnalysis::PedestalType out;
  calculateCommonMode(indat);
  int setNumber;
  if (!theCommonModeValues.empty()) {
    for (unsigned int i = 0; i < indat.size(); i++) {
      setNumber = theTkCommonMode->topology().setOfStrip(i);
      out.push_back(indat[i] - theCommonModeValues[setNumber]);
    }
  } else {
    out = indat;
  }
  return out;
}

getCMSlope()

float TT6CommonModeCalculator::getCMSlope ( )

inlineoverridevirtual

Get Slope

Implements TkCommonModeCalculator.

Definition at line 24 of file TT6CommonModeCalculator.h.

References slope.

{ return slope; }

newEvent()

void TT6CommonModeCalculator::newEvent ( )

overridevirtual

Tell CM calculator that a new event is available

Reimplemented from TkCommonModeCalculator.

Definition at line 75 of file TT6CommonModeCalculator.cc.

References alreadyUsedEvent.

{ alreadyUsedEvent = false; }

setCM() [1/2]

void TT6CommonModeCalculator::setCM ( TkCommonMode *  in )

inlineoverridevirtual

Implements TkCommonModeCalculator.

Definition at line 19 of file TT6CommonModeCalculator.h.

References recoMuon::in, and theTkCommonMode.

Referenced by calculateCommonMode().

{ theTkCommonMode = in; }

setCM() [2/2]

void TT6CommonModeCalculator::setCM ( const std::vector< float > &  in )

inlineoverridevirtual

Implements TkCommonModeCalculator.

Definition at line 20 of file TT6CommonModeCalculator.h.

References recoMuon::in, TkCommonMode::setCommonMode(), and theTkCommonMode.

{ theTkCommonMode->setCommonMode(in); }

Member Data Documentation

alreadyUsedEvent

bool TT6CommonModeCalculator::alreadyUsedEvent

protected

Definition at line 34 of file TT6CommonModeCalculator.h.

Referenced by calculateCommonMode(), and newEvent().

cutToAvoidSignal

float TT6CommonModeCalculator::cutToAvoidSignal

protected

Definition at line 37 of file TT6CommonModeCalculator.h.

Referenced by calculateCommonMode(), and TT6CommonModeCalculator().

slope

float TT6CommonModeCalculator::slope

protected

Definition at line 35 of file TT6CommonModeCalculator.h.

Referenced by calculateCMSlope(), and getCMSlope().

theApvMask

TkApvMask* TT6CommonModeCalculator::theApvMask

protected

Definition at line 33 of file TT6CommonModeCalculator.h.

Referenced by calculateCommonMode().

theCommonModeValues

std::vector<float> TT6CommonModeCalculator::theCommonModeValues

protected

Definition at line 31 of file TT6CommonModeCalculator.h.

Referenced by calculateCommonMode(), and doIt().

theNoiseCalculator

TkNoiseCalculator* TT6CommonModeCalculator::theNoiseCalculator

protected

Definition at line 32 of file TT6CommonModeCalculator.h.

Referenced by calculateCommonMode().

theTkCommonMode

TkCommonMode* TT6CommonModeCalculator::theTkCommonMode

protected

Definition at line 30 of file TT6CommonModeCalculator.h.

Referenced by calculateCommonMode(), commonMode(), doIt(), and setCM().