#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	(	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 ( )

override

Definition at line 14 of file TT6CommonModeCalculator.cc.

References gather_cfg::cout.

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

Member Function Documentation

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

protected

Definition at line 79 of file TT6CommonModeCalculator.cc.

References alignCSCRings::s, and slope.

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.;
 }

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);
 }

TkCommonMode* TT6CommonModeCalculator::commonMode ( )

inlineoverridevirtual

Get CM value

Implements TkCommonModeCalculator.

Definition at line 21 of file TT6CommonModeCalculator.h.

References theTkCommonMode.

21 { return theTkCommonMode; }

TT6CommonModeCalculator::theTkCommonMode

TkCommonMode * theTkCommonMode

Definition: TT6CommonModeCalculator.h:30

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, submitPVResolutionJobs::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;
 }