---

LASProfileJudge Class Reference

check if a LASModuleProfile is usable for being stored and fitted More...

#include <Alignment/LaserAlignment/src/LASProfileJudge.h>

List of all members.

Public Member Functions
void	EnableZeroFilter (bool)
	toggle the zero filter (passed from cfg file)
bool	IsSignalIn (const LASModuleProfile &, int)
	Check if a LASModuleProfile indicates that the module has been hit, i.e.
bool	JudgeProfile (const LASModuleProfile &, int)
	Check if a LASModuleProfile is usable for being stored, i.e.
	LASProfileJudge ()
Private Member Functions
double	GetNegativity (int)
	In case of too high laser intensities, the APV baselines tend to drop down.
bool	IsNegativePeaksInProfile (int)
	sometimes when the laser intensity is too high the APVs get confused and a negative peak (dip) shows up.
bool	IsPeaksInProfile (int)
	If the laser intensity is too small, there's no peak at all.
Private Attributes
bool	isZeroFilter
LASModuleProfile	profile
std::pair< unsigned int, double >	thePeak

---

Detailed Description

check if a LASModuleProfile is usable for being stored and fitted

Definition at line 15 of file LASProfileJudge.h.

---

Constructor & Destructor Documentation

LASProfileJudge::LASProfileJudge

(

)

Definition at line 16 of file LASProfileJudge.cc.

References isZeroFilter.

                                 {

  // switch on the zero filter by default
  isZeroFilter = true;

}

---

Member Function Documentation

void LASProfileJudge::EnableZeroFilter

(

bool

zeroFilter

)

toggle the zero filter (passed from cfg file)

Definition at line 87 of file LASProfileJudge.cc.

References isZeroFilter.

                                                        {

  isZeroFilter = zeroFilter;

}

double LASProfileJudge::GetNegativity

(

int

offset

)

[private]

In case of too high laser intensities, the APV baselines tend to drop down.

here, the strip amplitudes in the area around the signal region are summed to return a variable which can indicate this.

Definition at line 102 of file LASProfileJudge.cc.

References LASModuleProfile::GetValue(), i, and profile.

Referenced by IsSignalIn(), and JudgeProfile().

                                                  {

  // here we could later run the sum only on the affected (pair of) APV

  // expected beam position (in strips)
  const unsigned int meanPosition = 256 + offset;
  // backplane "alignment hole" (="signal region") approx. half size 
  const unsigned int halfWindowSize = 33;
  // half size of range over which is summed (must be > halfWindowSize)
  const unsigned int sumHalfRange = 128;

  double neg = 0;
  
  for( unsigned int i = meanPosition - sumHalfRange; i < meanPosition - halfWindowSize; ++i ) {
    neg += profile.GetValue( i );
  }

  for( unsigned int i = meanPosition + halfWindowSize; i < meanPosition + sumHalfRange; ++i ) {
    neg += profile.GetValue( i );
  }

  return( neg );

}

bool LASProfileJudge::IsNegativePeaksInProfile

(

int

offset

)

[private]

sometimes when the laser intensity is too high the APVs get confused and a negative peak (dip) shows up.

this is filtered here.

Definition at line 174 of file LASProfileJudge.cc.

References average(), LASModuleProfile::GetValue(), profile, strip(), and thePeak.

Referenced by IsSignalIn(), and JudgeProfile().

                                                           {

  // expected beam position in middle of module (in strips)
  const unsigned int meanPosition = 256 + offset;
  // backplane "alignment hole" approx. half size (in strips)
  const unsigned int halfWindowSize = 33;

  bool returnValue = false;
  
  // calculate average out-of-signal
  double average = 0., counterD = 0.;
  for( unsigned int strip = 0; strip < 512; ++strip ) {
    if( strip < meanPosition - halfWindowSize || strip > meanPosition + halfWindowSize ) {
      average += profile.GetValue( strip );
      counterD += 1.;
    }
  }
  average /= counterD;
  
  // find strips with negative amplitude way above noise level
  const double noiseLevel = 2.;
  for( unsigned int strip = 0; strip < 512; ++strip ) {
    if( profile.GetValue( strip ) < ( average - 10. * noiseLevel ) ) { 
      returnValue = true;
      thePeak.first = strip; thePeak.second = profile.GetValue( strip );
      break;
    }
  }

  return( returnValue );

}

bool LASProfileJudge::IsPeaksInProfile

(

int

offset

)

[private]

If the laser intensity is too small, there's no peak at all.

Here we look if any strip is well above noise level.

Definition at line 134 of file LASProfileJudge.cc.

References average(), LASModuleProfile::GetValue(), profile, strip(), and thePeak.

Referenced by IsSignalIn(), and JudgeProfile().

                                                   {

  // expected beam position (in strips)
  const unsigned int meanPosition = 256 + offset;
  // backplane "alignment hole" approx. half size (in strips)
  const unsigned int halfWindowSize = 33;

  bool returnValue = false;
  
  // calculate average out-of-signal
  double average = 0., counterD = 0.;
  for( unsigned int strip = 0; strip < 512; ++strip ) {
    if( strip < meanPosition - halfWindowSize || strip > meanPosition + halfWindowSize ) {
      average += profile.GetValue( strip );
      counterD += 1.;
    }
  }
  average /= counterD;

  // find peaks well above noise level
  const double noiseLevel = 2.; // to be softcoded..
  for( unsigned int strip = meanPosition - halfWindowSize; strip < meanPosition + halfWindowSize; ++strip ) {
    if( profile.GetValue( strip ) > ( average + 10. * noiseLevel ) ) { 
      returnValue = true;
      thePeak.first = strip; thePeak.second = profile.GetValue( strip );
      break;
    }
  }

  return( returnValue );

}

bool LASProfileJudge::IsSignalIn	(	const LASModuleProfile &	aProfile,
		int	offset
	)

Check if a LASModuleProfile indicates that the module has been hit, i.e.

contains a visible signal or is even distorted by too high laser amplitude. This method doesn't care if the profile is usable for analysis.

Definition at line 32 of file LASProfileJudge.cc.

References GetNegativity(), IsNegativePeaksInProfile(), IsPeaksInProfile(), profile, and HLT_VtxMuL3::result.

Referenced by LaserAlignment::produce().

                                                                               {

  profile = aProfile;
  
  double negativity = GetNegativity( offset );
  bool isPeaks = IsPeaksInProfile( offset );
  bool isNegativePeaks = IsNegativePeaksInProfile( offset );
  
  bool result = 
    ( negativity < -1000. ) ||  // if we see negativity, there was laser..
    ( isPeaks )             ||  // if we see a peak, " " "
    ( isNegativePeaks );    // same here
  
  return( result );


}

bool LASProfileJudge::JudgeProfile	(	const LASModuleProfile &	aProfile,
		int	offset = 0
	)

Check if a LASModuleProfile is usable for being stored, i.e.

contains a visible signal & no baseline distortions

Definition at line 58 of file LASProfileJudge.cc.

References GetNegativity(), IsNegativePeaksInProfile(), IsPeaksInProfile(), isZeroFilter, offset, profile, and HLT_VtxMuL3::result.

Referenced by LaserAlignment::produce().

                                                                                     {

  profile = aProfile;
  
  // run the tests
  double negativity = GetNegativity( offset );

  bool isPeaks;
  if( isZeroFilter ) isPeaks = true; // disable this test if set in cfg
  else isPeaks = IsPeaksInProfile( offset );

  bool isNegativePeaks = IsNegativePeaksInProfile( offset );

  bool result = 
    ( negativity > -1000. ) &&  // < 1000. = distorted profile
    ( isPeaks )             &&  // want to see a peak (zero filter)
    !( isNegativePeaks ); // no negative peaks

  return( result );

}

---

Member Data Documentation

bool LASProfileJudge::isZeroFilter [private]

Definition at line 29 of file LASProfileJudge.h.

Referenced by EnableZeroFilter(), JudgeProfile(), and LASProfileJudge().

LASModuleProfile LASProfileJudge::profile [private]

Definition at line 27 of file LASProfileJudge.h.

Referenced by GetNegativity(), IsNegativePeaksInProfile(), IsPeaksInProfile(), IsSignalIn(), and JudgeProfile().

std::pair<unsigned int, double> LASProfileJudge::thePeak [private]

Definition at line 28 of file LASProfileJudge.h.

Referenced by IsNegativePeaksInProfile(), and IsPeaksInProfile().

---

The documentation for this class was generated from the following files:

• Alignment/LaserAlignment/src/LASProfileJudge.h
• Alignment/LaserAlignment/src/LASProfileJudge.cc

---