#include <TT6ApvMask.h>

Inheritance diagram for TT6ApvMask:

Public Member Functions
void	calculateMask (const ApvAnalysis::PedestalType &)

MaskType	mask ()

void	setMask (const MaskType &in)

	TT6ApvMask (int ctype, float ncut, float dcut, float tcut)

virtual	~TT6ApvMask ()

Public Member Functions inherited from TkApvMask
virtual	~TkApvMask ()

Protected Member Functions
bool	defineNoisy (float avrg, float rms, float noise)

Private Attributes
int	theCalculationFlag_

float	theDeadCut_

MaskType	theMask_

float	theNoiseCut_

float	theTruncationCut_

Additional Inherited Members
Public Types inherited from TkApvMask
typedef std::vector< StripMaskType >	MaskType

enum	StripMaskType { ok =0, dead =1, noisy =2 }

Detailed Description

Concrete implementation of TkApvMask for TT6.

Definition at line 10 of file TT6ApvMask.h.

Constructor & Destructor Documentation

TT6ApvMask::TT6ApvMask	(	int	ctype,
		float	ncut,
		float	dcut,
		float	tcut
	)

Definition at line 9 of file TT6ApvMask.cc.

                                                                     {
   theCalculationFlag_ = ctype;
   theNoiseCut_        = ncut;
   theDeadCut_         = dcut;
   theTruncationCut_   = tcut;
 }

TT6ApvMask::~TT6ApvMask ( )

virtual

Definition at line 19 of file TT6ApvMask.cc.

References gather_cfg::cout.

                        {
   if (0) cout << "Destructing TT6ApvMask " << endl;
 }

Member Function Documentation

void TT6ApvMask::calculateMask ( const ApvAnalysis::PedestalType & in )

virtual

Implements TkApvMask.

Definition at line 25 of file TT6ApvMask.cc.

References gather_cfg::cout, mps_fire::i, createfilelist::int, convertSQLiteXML::ok, and mathSSE::sqrt().

Referenced by mask().

                                                                {
 
   theMask_.clear();
   ApvAnalysis::PedestalType temp_in(in);
   double sumVal,sqSumVal,avVal,sqAvVal,rmsVal; 
   sort(temp_in.begin(), temp_in.end());
   int nSize    = in.size();
   int cutLow   = int(nSize * theTruncationCut_);
   int cutHigh  = int(nSize * theTruncationCut_);
   int effSize  = nSize - cutLow - cutHigh;
   sumVal = 0.0;
   sqSumVal = 0.0;
   sumVal   = accumulate((temp_in.begin()+cutLow), (temp_in.end()-cutHigh), 0.0);
   sqSumVal = inner_product((temp_in.begin()+cutLow), (temp_in.end()-cutHigh), 
                            (temp_in.begin()+cutLow), 0.0);
   
   avVal    = (effSize) ? sumVal/float(effSize):0.0;
   sqAvVal  = (effSize) ? sqSumVal/float(effSize):0.0;
   rmsVal   = (sqAvVal - avVal*avVal > 0.0) ? sqrt(sqAvVal - avVal*avVal):0.0; 
   if (0) cout << " TT6ApvMask::calculateMask  Mean " << avVal <<
        " RMS " << rmsVal << " " <<  effSize << endl;       
   for (unsigned int i=0; i<in.size(); i++){
     if (defineNoisy( static_cast<float>(avVal),
                      static_cast<float>(rmsVal),in[i])) {
       theMask_.push_back(noisy);
     } else if (in[i] < theDeadCut_*avVal) {
            theMask_.push_back(dead);
     } else {
       theMask_.push_back(ok);
     }     
   }
 }

bool TT6ApvMask::defineNoisy	(	float	avrg,
		float	rms,
		float	noise
	)

protected

Definition at line 62 of file TT6ApvMask.cc.

References gather_cfg::cout, and groupFilesInBlocks::temp.

Referenced by mask().

                                                             {
   bool temp;
   temp=false;
   if (theCalculationFlag_ == 1){
    if ((noise-avrg) > theNoiseCut_*rms) {
      temp=true;
      if (0) cout << " Mean " << avrg << " rms " << rms << " Noise " << noise << endl;
    }
   } else if (theCalculationFlag_ == 2){
     if ((noise-avrg) > avrg*theNoiseCut_) temp=true;
   } else if (theCalculationFlag_ == 3){
     if (noise > theNoiseCut_) temp=true;
   } 
   return temp;
 }