QwHelicityPattern Class Reference

Manages yields/differences/asymmetries for helicity patterns. More...

#include <QwHelicityPattern.h>

Inheritance diagram for QwHelicityPattern:

Collaboration diagram for QwHelicityPattern:

Public Member Functions
	QwHelicityPattern (QwSubsystemArrayParity &event, const TString &run="0")
	Constructor with subsystem array.
	QwHelicityPattern (const QwHelicityPattern &source)
	Copy constructor by reference.
virtual	~QwHelicityPattern ()
	Virtual destructor.
void	ProcessOptions (QwOptions &options)
	Process the configuration options.
void	LoadEventData (QwSubsystemArrayParity &event)
Bool_t	HasDataLoaded () const
Bool_t	PairAsymmetryIsGood ()
Bool_t	NextPairIsComplete ()
void	CalculatePairAsymmetry ()
void	ClearPairData ()
Bool_t	IsCompletePattern () const
void	PrintIndexMapFile (Int_t runNum)
Bool_t	IsEndOfBurst ()
void	CalculateAsymmetry ()
void	GetTargetChargeStat (Double_t &asym, Double_t &error, Double_t &width)
void	EnableAlternateAsymmetry (const Bool_t flag=kTRUE)
	Enable/disable alternate asymmetry calculation.
void	DisableAlternateAsymmetry ()
	Disable alternate asymmetry calculation.
Bool_t	IsAlternateAsymEnabled ()
	Status of alternate asymmetry calculation flag.
void	EnableBurstSum (const Bool_t flag=kTRUE)
	Enable/disable burst sum calculation.
void	DisableBurstSum ()
	Disable burst sum calculation.
Bool_t	IsBurstSumEnabled ()
	Status of burst sum calculation flag.
void	EnableRunningSum (const Bool_t flag=kTRUE)
	Enable/disable running sum calculation.
void	DisableRunningSum ()
	Disable running sum calculation.
Bool_t	IsRunningSumEnabled ()
	Status of running sum calculation flag.
void	EnableDifference (const Bool_t flag=kTRUE)
	Enable/disable storing pattern differences.
void	DisableDifference ()
	Disable storing pattern differences.
Bool_t	IsDifferenceEnabled ()
	Status of storing pattern differences flag.
void	EnablePairs (const Bool_t flag=kTRUE)
	Enable/disable storing pair differences.
void	DisablePairs ()
	Disable storing pair differences.
Bool_t	IsPairsEnabled ()
	Status of storing pair differences flag.
void	UpdateBlinder ()
	Update the blinder status using a random number generator.
void	UpdateBlinder (const QwSubsystemArrayParity &detectors)
	Update the blinder status with new external information.
void	UpdateBlinder (const QwEPICSEvent &epics)
	Update the blinder status with new external information.
QwSubsystemArrayParity &	GetYield ()
QwSubsystemArrayParity &	GetDifference ()
QwSubsystemArrayParity &	GetAsymmetry ()
QwSubsystemArrayParity &	GetPairYield ()
QwSubsystemArrayParity &	GetPairDifference ()
QwSubsystemArrayParity &	GetPairAsymmetry ()
void	AccumulateRunningSum (QwHelicityPattern &entry, Int_t count=0, Int_t ErrorMask=0xFFFFFFF)
void	AccumulatePairRunningSum (QwHelicityPattern &entry)
void	CalculateRunningAverage ()
void	PrintValue () const
void	ConstructObjects ()
void	ConstructObjects (TDirectory *folder)
void	ConstructHistograms ()
void	ConstructHistograms (TDirectory *folder)
void	FillHistograms ()
void	ConstructBranchAndVector (TTree *tree, TString &prefix, QwRootTreeBranchVector &values)
void	ConstructBranch (TTree *tree, TString &prefix)
void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &trim_tree)
void	FillTreeVector (QwRootTreeBranchVector &values) const
void	WritePromptSummary (QwPromptSummary *ps)
Bool_t	IsGoodAsymmetry ()
UInt_t	GetEventcutErrorFlag () const
const UInt_t *	GetEventcutErrorFlagPointer () const
Bool_t	HasBurstData ()
void	IncrementBurstCounter ()
Short_t	GetBurstCounter () const
void	ClearEventData ()
void	Print () const

Static Public Member Functions
static void	DefineOptions (QwOptions &options)
	Define the configuration options.

Protected Member Functions
void	SetDataLoaded (Bool_t flag)

Protected Attributes
std::vector< QwSubsystemArrayParity >	fEvents
std::vector< Bool_t >	fEventLoaded
std::vector< Int_t >	fHelicity
std::vector< Int_t >	fEventNumber
Int_t	fCurrentPatternNumber
size_t	fPatternSize
Int_t	fQuartetNumber
QwBlinder	fBlinder
Bool_t	fHelicityIsMissing
Bool_t	fIgnoreHelicity
QwSubsystemArrayParity	fYield
QwSubsystemArrayParity	fDifference
QwSubsystemArrayParity	fAsymmetry
Bool_t	fEnableAlternateAsym
QwSubsystemArrayParity	fAsymmetry1
QwSubsystemArrayParity	fAsymmetry2
Bool_t	fEnablePairs
QwSubsystemArrayParity	fPairYield
QwSubsystemArrayParity	fPairDifference
QwSubsystemArrayParity	fPairAsymmetry
Int_t	fBurstLength
Int_t	fMaxBurstIndex
Bool_t	fPrintIndexFile
Int_t	fBurstMinGoodPatterns
Int_t	fGoodPatterns
Short_t	fBurstCounter
Bool_t	fEnableBurstSum
Bool_t	fPrintBurstSum
Bool_t	fEnableRunningSum
Bool_t	fPrintRunningSum
Bool_t	fEnableDifference
QwSubsystemArrayParity	fAlternateDiff
QwSubsystemArrayParity	fPositiveHelicitySum
QwSubsystemArrayParity	fNegativeHelicitySum
Long_t	fLastWindowNumber
Long_t	fLastPatternNumber
Int_t	fLastPhaseNumber
size_t	fNextPair
Bool_t	fPairIsGood
Bool_t	fPatternIsGood
TString	run_label
Bool_t	fIsDataLoaded

Private Member Functions
	QwHelicityPattern ()
	Private default constructor (not implemented, will throw linker error on use)

Friends
class	QwDataHandlerArray

Detailed Description

Manages yields/differences/asymmetries for helicity patterns.

Aggregates per-window subsystem arrays into pairs and patterns, supports burst and running sums, and integrates with the blinder strategy.

Definition at line 42 of file QwHelicityPattern.h.

Constructor & Destructor Documentation

QwHelicityPattern() [1/3]

QwHelicityPattern::QwHelicityPattern ( )

private

Private default constructor (not implemented, will throw linker error on use)

References QwHelicityPattern().

Referenced by AccumulatePairRunningSum(), AccumulateRunningSum(), QwHelicityCorrelatedFeedback::QwHelicityCorrelatedFeedback(), QwHelicityPattern(), and QwHelicityPattern().

Here is the call graph for this function:

Here is the caller graph for this function:

QwHelicityPattern() [2/3]

QwHelicityPattern::QwHelicityPattern	(	QwSubsystemArrayParity &	event,
		const TString &	run = "0" )

Constructor with subsystem array.

Definition at line 106 of file QwHelicityPattern.cc.

  : fBlinder(),
    fHelicityIsMissing(kFALSE),   
    fIgnoreHelicity(kFALSE),
    fYield(event), 
    fDifference(event), 
    fAsymmetry(event),   
    fEnableAlternateAsym(kFALSE), 
    fAsymmetry1(event), 
    fAsymmetry2(event),
    fEnablePairs(kTRUE),
    fPairYield(event), 
    fPairDifference(event), 
    fPairAsymmetry(event),
    fBurstLength(0),
    fMaxBurstIndex(0x7fffffff),
    fPrintIndexFile(kFALSE),
    fBurstMinGoodPatterns(0),
    fGoodPatterns(0),
    fBurstCounter(0),
    fEnableBurstSum(kFALSE),
    fPrintBurstSum(kFALSE),
    fEnableRunningSum(kTRUE),     
    fPrintRunningSum(kFALSE),
    fEnableDifference(kFALSE),
    fAlternateDiff(event),
    fPositiveHelicitySum(event), 
    fNegativeHelicitySum(event),
    fLastWindowNumber(0),
    fLastPatternNumber(0),
    fLastPhaseNumber(0),
    fNextPair(0),
    fPairIsGood(false),
    fPatternIsGood(false),
    fIsDataLoaded(false)
{
  // Retrieve the helicity subsystem to query for
  std::vector<VQwSubsystem*> subsys_helicity = event.GetSubsystemByType("QwHelicity");
  if (subsys_helicity.size() > 0) {
 
    // Take the first helicity subsystem
    //    fHelicitySubsystem = dynamic_cast<QwHelicity*>(subsys_helicity.at(0));
    QwHelicity* helicity = dynamic_cast<QwHelicity*>(subsys_helicity.at(0));
 
    // Get the maximum pattern phase (i.e. pattern size)
    fPatternSize = helicity->GetMaxPatternPhase();
 
    // Warn if more than one helicity subsystem defined
    if (subsys_helicity.size() > 1)
      QwWarning << "Multiple helicity subsystems defined! "
                << "Using " << helicity->GetName() << "."
                << QwLog::endl;
 
  } else {
    // We are not using any helicity subsystem
    QwError << "No helicity subsystem defined!  " 
        << "Calculate asymmetries based on (+--+) quartets!"
        << QwLog::endl;
    fHelicityIsMissing = kTRUE;
    fPatternSize = 4; // default to quartets
  }
  QwMessage << "QwHelicity::MaxPatternPhase = " << fPatternSize << QwLog::endl;
 
  try
    {
      if(fPatternSize%2 == 0)
        {
          fEvents.resize(fPatternSize,event);
          fHelicity.resize(fPatternSize,-9999);
          fEventNumber.resize(fPatternSize,-1);
          fEventLoaded.resize(fPatternSize,kFALSE);
 
          // Initialize the pattern number
          fQuartetNumber = 0;
          fCurrentPatternNumber = -1;
        }
      else
        {
          TString loc=
            "Standard exception from QwHelicityPattern : the pattern size has to be even;  right now pattern_size=";
          loc+=Form("%zu",fPatternSize);
          throw std::invalid_argument(loc.Data());
        }
    }
  catch (std::exception& e)
    {
      std::cerr << e.what() << std::endl;
    }
 
}

References QwLog::endl(), fAlternateDiff, fAsymmetry, fAsymmetry1, fAsymmetry2, fBlinder, fBurstCounter, fBurstLength, fBurstMinGoodPatterns, fDifference, fEnableAlternateAsym, fEnableBurstSum, fEnableDifference, fEnablePairs, fEnableRunningSum, fGoodPatterns, fHelicityIsMissing, fIgnoreHelicity, fIsDataLoaded, fLastPatternNumber, fLastPhaseNumber, fLastWindowNumber, fMaxBurstIndex, fNegativeHelicitySum, fNextPair, fPairAsymmetry, fPairDifference, fPairIsGood, fPairYield, fPatternIsGood, fPatternSize, fPositiveHelicitySum, fPrintBurstSum, fPrintIndexFile, fPrintRunningSum, fYield, and QwMessage.

Here is the call graph for this function:

QwHelicityPattern() [3/3]

QwHelicityPattern::QwHelicityPattern

(

const QwHelicityPattern &

source

)

Copy constructor by reference.

Definition at line 198 of file QwHelicityPattern.cc.

: fPatternSize(source.fPatternSize),
  fYield(source.fYield),
  fDifference(source.fDifference),
  fAsymmetry(source.fAsymmetry),
  fEnableAlternateAsym(source.fEnableAlternateAsym),
  fAsymmetry1(source.fAsymmetry1),
  fAsymmetry2(source.fAsymmetry2),
  fPairYield(source.fYield),
  fPairDifference(source.fYield),
  fPairAsymmetry(source.fYield),
  fBurstLength(source.fBurstLength),
  fMaxBurstIndex(source.fMaxBurstIndex),
  fPrintIndexFile(source.fPrintIndexFile),
  fBurstMinGoodPatterns(source.fBurstMinGoodPatterns),
  fGoodPatterns(source.fGoodPatterns),
  fBurstCounter(source.fBurstCounter),
  fEnableBurstSum(source.fEnableBurstSum),
  fPrintBurstSum(source.fPrintBurstSum),
  fEnableRunningSum(source.fEnableRunningSum),
  fPrintRunningSum(source.fPrintRunningSum),
  fEnableDifference(source.fEnableDifference),
  fAlternateDiff(source.fYield),
  fPositiveHelicitySum(source.fYield),
  fNegativeHelicitySum(source.fYield),
  fLastWindowNumber(0),
  fLastPatternNumber(0),
  fLastPhaseNumber(0),
  fNextPair(0),
  fPairIsGood(false),
  fPatternIsGood(false),
  fIsDataLoaded(false)
{
};

References fAlternateDiff, fAsymmetry, fAsymmetry1, fAsymmetry2, fBurstCounter, fBurstLength, fBurstMinGoodPatterns, fDifference, fEnableAlternateAsym, fEnableBurstSum, fEnableDifference, fEnableRunningSum, fGoodPatterns, fIsDataLoaded, fLastPatternNumber, fLastPhaseNumber, fLastWindowNumber, fMaxBurstIndex, fNegativeHelicitySum, fNextPair, fPairAsymmetry, fPairDifference, fPairIsGood, fPairYield, fPatternIsGood, fPatternSize, fPositiveHelicitySum, fPrintBurstSum, fPrintIndexFile, fPrintRunningSum, fYield, and QwHelicityPattern().

Here is the call graph for this function:

~QwHelicityPattern()

virtual QwHelicityPattern::~QwHelicityPattern ( )

inlinevirtual

Virtual destructor.

Definition at line 57 of file QwHelicityPattern.h.

{ };

Member Function Documentation

AccumulatePairRunningSum()

void QwHelicityPattern::AccumulatePairRunningSum

(

QwHelicityPattern &

entry

)

Accumulate the running sum for pairs by adding this helicity pattern to the running sums of yield, difference and asymmetry.

Definition at line 726 of file QwHelicityPattern.cc.

{
  if (entry.fPairIsGood){
    fPairYield.AccumulateRunningSum(entry.fPairYield);
    fPairAsymmetry.AccumulateRunningSum(entry.fPairAsymmetry);
    if (fEnableDifference){
      fPairDifference.AccumulateRunningSum(entry.fPairDifference);
      // The difference is blinded, so the running difference is also blinded.
    }
    fPairIsGood = entry.fPairIsGood;
  }
}

References fEnableDifference, fPairAsymmetry, fPairDifference, fPairIsGood, fPairYield, and QwHelicityPattern().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

AccumulateRunningSum()

void QwHelicityPattern::AccumulateRunningSum	(	QwHelicityPattern &	entry,
		Int_t	count = 0,
		Int_t	ErrorMask = 0xFFFFFFF )

Accumulate the running sum by adding this helicity pattern to the running sums of yield, difference and asymmetry.

Definition at line 698 of file QwHelicityPattern.cc.

{
  if (entry.fPatternIsGood){
    //    if ( (*(entry.fAsymmetry.GetEventcutErrorFlagPointer()) & ErrorMask) == 0 )  {
    if (entry.fAsymmetry.GetEventcutErrorFlag()==0){
      fGoodPatterns++;
    }
    fBurstCounter = entry.fBurstCounter;
    fYield.AccumulateRunningSum(entry.fYield, count, ErrorMask);
    fAsymmetry.AccumulateRunningSum(entry.fAsymmetry, count, ErrorMask);
    if (fEnableDifference){
      fDifference.AccumulateRunningSum(entry.fDifference, count, ErrorMask);
      // The difference is blinded, so the running difference is also blinded.
    }
    if (fEnableAlternateAsym) {
      fAsymmetry1.AccumulateRunningSum(entry.fAsymmetry1, count, ErrorMask);
      fAsymmetry2.AccumulateRunningSum(entry.fAsymmetry2, count, ErrorMask);
    }
    fPatternIsGood = entry.fPatternIsGood;
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fBurstCounter, fDifference, fEnableAlternateAsym, fEnableDifference, fGoodPatterns, fPatternIsGood, fYield, QwSubsystemArrayParity::GetEventcutErrorFlag(), and QwHelicityPattern().

Referenced by QwHelicityCorrelatedFeedback::AccumulateRunningSum(), and main().

Here is the call graph for this function:

Here is the caller graph for this function:

CalculateAsymmetry()

void QwHelicityPattern::CalculateAsymmetry

(

)

Calculate asymmetries for the current pattern.

Definition at line 481 of file QwHelicityPattern.cc.

{
 
  Bool_t localdebug=kFALSE;
 
  if(localdebug)  std::cout<<"Entering QwHelicityPattern::CalculateAsymmetry \n";
 
  Int_t plushel  = 1;
  Int_t minushel = 0;
  Int_t checkhel = 0;
  Bool_t firstplushel=kTRUE;
  Bool_t firstminushel=kTRUE;
 
  fPositiveHelicitySum.ClearEventData();
  fNegativeHelicitySum.ClearEventData();
 
  if (fIgnoreHelicity){
    //  Don't check to see if we have equal numbers of even and odd helicity states in this pattern.
    //  Build an asymmetry with even-parity phases as "+" and odd-parity phases as "-"
    for (size_t i = 0; i < fPatternSize; i++) {
      Int_t localhel = 1;
      for (size_t j = 0; j < fPatternSize/2; j++) {
    localhel ^= ((i >> j)&0x1);
      }
      if (localhel == plushel) {
    if (firstplushel) {
      fPositiveHelicitySum = fEvents.at(i);
      firstplushel = kFALSE;
    } else {
      fPositiveHelicitySum += fEvents.at(i);
    }
      } else if (localhel == minushel) {
    if (firstminushel) {
      fNegativeHelicitySum = fEvents.at(i);
      firstminushel = kFALSE;
    } else {
      fNegativeHelicitySum += fEvents.at(i);
    }
      }
    }
  } else {
    //  
    for (size_t i = 0; i < fPatternSize; i++) {
      if (fHelicity[i] == plushel) {
    if (localdebug) std::cout<<"QwHelicityPattern::CalculateAsymmetry:  here filling fPositiveHelicitySum \n";
    if (firstplushel) {
      if (localdebug) std::cout<<"QwHelicityPattern::CalculateAsymmetry:  with = \n";
      fPositiveHelicitySum = fEvents.at(i);
      firstplushel = kFALSE;
    } else {
      if (localdebug) std::cout<<"QwHelicityPattern::CalculateAsymmetry:  with += \n";
      fPositiveHelicitySum += fEvents.at(i);
    }
    checkhel += 1;
      } else if (fHelicity[i] == minushel) {
    if (localdebug) std::cout<<"QwHelicityPattern::CalculateAsymmetry:  here filling fNegativeHelicitySum \n";
    if (firstminushel) {
      if (localdebug) std::cout<<"QwHelicityPattern::CalculateAsymmetry:  with = \n";
      fNegativeHelicitySum = fEvents.at(i);
      firstminushel = kFALSE;
    } else {
      if (localdebug) std::cout<<"QwHelicityPattern::CalculateAsymmetry:  with += \n";
      fNegativeHelicitySum += fEvents.at(i);
    }
    checkhel -= 1;
      } else {
    QwDebug << "QwHelicityPattern::CalculateAsymmetry:  "
        << "Helicity should be "<<plushel<<" or "<<minushel
        <<" but is "<< fHelicity[i]
        << "; Asymmetry computation aborted!"<<QwLog::endl;
    ClearEventData();
    checkhel = -9999;
  break;
    // This is an unknown helicity event.
      }
    }
  }
 
  if (checkhel == -9999) {
    //do nothing the asymmetry computation has been aborted earlier in this function
    fPatternIsGood = kFALSE;
  } else if (checkhel!=0) {
    fPatternIsGood = kFALSE;
    // there is a different number of plus and minus helicity window.
    QwError<<" QwHelicityPattern::CalculateAsymmetry == \n"
       <<" you do not have the same number of positive and negative \n"
       <<" impossible to compute asymmetry \n"
       <<" dropping every thing -- pattern number ="<<fCurrentPatternNumber<<QwLog::endl;
  } else {
    //  This is a good pattern.
    //  Calculate the asymmetry.
    fPatternIsGood = kTRUE;
    fQuartetNumber++;//Then increment the quartet number
    //std::cout<<" quartet count ="<<fQuartetNumber<<"\n";
 
    fYield.Sum(fPositiveHelicitySum,fNegativeHelicitySum);
    fYield.Scale(1.0/fPatternSize);
    fDifference.Difference(fPositiveHelicitySum,fNegativeHelicitySum);
    fDifference.Scale(1.0/fPatternSize);
 
    
    
    if (! fIgnoreHelicity){
      // Update the blinder if conditions have changed
      UpdateBlinder(fYield);
      //  Only blind the difference if we're using the real helicity.
      fBlinder.Blind(fDifference,fYield);
      //  Update the global error code in fDifference, and use it
      //  to update the errors in fYield, in case blinder errors
      //  can propagate to the global error.
      fDifference.UpdateErrorFlag();
      fYield.UpdateErrorFlag(fDifference);
      if (! fBlinder.IsBlinderOkay()){
    fYield.UpdateErrorFlag(QwBlinder::kErrorFlag_BlinderFail);
    fDifference.UpdateErrorFlag(QwBlinder::kErrorFlag_BlinderFail);
      }
    }
    fAsymmetry.Ratio(fDifference,fYield);
    fAsymmetry.IncrementErrorCounters();
 
    /*
      With additional two asymmetry calculations
      Don't blind them!
 
      quartet pattern + - - +
                      1 2 3 4
                      fAsymmetry  = (1+4)-(2+3)/(1+2+3+4)
                      fAsymmetry1 = (1+2)-(3+4)/(1+2+3+4)
                      fAsymmetry2 = (1+3)-(2+4)/(1+2+3+4)
    */
    if (fEnableAlternateAsym) {
      //  fAsymmetry1:  (first 1/2 pattern - second 1/2 pattern)/fYield
      fPositiveHelicitySum.ClearEventData();
      fNegativeHelicitySum.ClearEventData();
      fPositiveHelicitySum = fEvents.at(0);
      fNegativeHelicitySum = fEvents.at(fPatternSize/2);
      if (fPatternSize/2 > 1){
    for (size_t i = 1; i < fPatternSize/2 ; i++){
      fPositiveHelicitySum += fEvents.at(i);
      fNegativeHelicitySum += fEvents.at(fPatternSize/2 +i);
    }
      }
      fAlternateDiff.ClearEventData();
      fAlternateDiff.Difference(fPositiveHelicitySum,fNegativeHelicitySum);
      fAlternateDiff.Scale(1.0/fPatternSize);
      //  Do not blind this helicity-uncorrelated difference.
      fAsymmetry1.Ratio(fAlternateDiff,fYield);
      //  fAsymmetry2:  (even events - odd events)/fYield
      //  Only build fAsymmetry2 if fPatternSize>2.
      if (fPatternSize > 2) {
    fPositiveHelicitySum.ClearEventData();
    fNegativeHelicitySum.ClearEventData();
    fPositiveHelicitySum = fEvents.at(0);
    fNegativeHelicitySum = fEvents.at(1);
    if (fPatternSize/2 > 1){
      for (size_t i = 1; i < fPatternSize/2 ; i++){
        fPositiveHelicitySum += fEvents.at(2*i);
        fNegativeHelicitySum += fEvents.at(2*i + 1);
      }
    }
    fAlternateDiff.ClearEventData();
    fAlternateDiff.Difference(fPositiveHelicitySum,fNegativeHelicitySum);
    fAlternateDiff.Scale(1.0/fPatternSize);
    //  Do not blind this helicity-uncorrelated difference.
    fAsymmetry2.Ratio(fAlternateDiff,fYield);
      }
    }
 
    if (localdebug) QwDebug << " pattern number =" << fQuartetNumber << QwLog::endl;
  }
}

References ClearEventData(), QwLog::endl(), fAlternateDiff, fAsymmetry, fAsymmetry1, fAsymmetry2, fBlinder, fCurrentPatternNumber, fDifference, fEnableAlternateAsym, fEvents, fHelicity, fIgnoreHelicity, fNegativeHelicitySum, fPatternIsGood, fPatternSize, fPositiveHelicitySum, fQuartetNumber, fYield, QwBlinder::kErrorFlag_BlinderFail, QwDebug, QwError, and UpdateBlinder().

Referenced by IsGoodAsymmetry().

Here is the call graph for this function:

Here is the caller graph for this function:

CalculatePairAsymmetry()

void QwHelicityPattern::CalculatePairAsymmetry

(

)

Definition at line 363 of file QwHelicityPattern.cc.

{
 
  Bool_t localdebug=kFALSE;
 
  if(localdebug)  std::cout<<"Entering QwHelicityPattern::CalculateAsymmetry \n";
 
  Int_t plushel  = 1;
  Int_t minushel = 0;
  [[maybe_unused]]
  Int_t checkhel = 0;
 
  if (fNextPair<fPatternSize/2){
    size_t firstevt  = fNextPair*2;
    size_t secondevt = firstevt + 1;
    fPairIsGood = kTRUE;
    fNextPair++;
 
    fPairYield.Sum(fEvents.at(firstevt), fEvents.at(secondevt));
    fPairYield.Scale(0.5);
  
    if (fIgnoreHelicity){
      fPairDifference.Difference(fEvents.at(firstevt), fEvents.at(secondevt));
      fPairDifference.Scale(0.5);
    } else {
      if (fHelicity[firstevt] == plushel && fHelicity[firstevt]!=fHelicity[secondevt]) {
    fPairDifference.Difference(fEvents.at(firstevt), fEvents.at(secondevt));
    fPairDifference.Scale(0.5);
      } else if (fHelicity[firstevt] == minushel && fHelicity[firstevt]!=fHelicity[secondevt]) {
    fPairDifference.Difference(fEvents.at(secondevt),fEvents.at(firstevt));
    fPairDifference.Scale(0.5);
      } else if (fHelicity[firstevt] == -9999 || fHelicity[secondevt]==-9999) {
    checkhel= -9999;
    // Helicity polarity is undefined.
    QwDebug<<" QwHelicityPattern::CalculatePairAsymmetry == \n"
           <<" undefined local helicity (-9999) \n"
           <<" impossible to compute asymmetry \n"
           <<" dropping every thing -- pattern number ="<<fCurrentPatternNumber<<QwLog::endl;
    // This is an unknown helicity event.
    fPairIsGood = kFALSE;
      }else {
    QwDebug << "QwHelicityPattern::CalculatePairAsymmetry:  "
        << "Helicity should be "<<plushel<<" or "<<minushel
        <<" but is "<< fHelicity[firstevt] << " and " << fHelicity[secondevt]
        << "; Asymmetry computation aborted!"<<QwLog::endl;
    fPairIsGood = kFALSE;
    // there is a different number of plus and minus helicity window.
    QwError<<" QwHelicityPattern::CalculatePairAsymmetry == \n"
           <<" you do not have the same number of positive and negative \n"
           <<" impossible to compute asymmetry \n"
           <<" dropping every thing -- pattern number ="<<fCurrentPatternNumber<<QwLog::endl;
    // This is an unknown helicity event.
      }
    }
  }
 
  if (fPairIsGood){
    if (! fIgnoreHelicity){
      //      // Update the blinder if conditions have changed
      UpdateBlinder(fPairYield);
      //  Only blind the difference if we're using the real helicity.
      fBlinder.BlindPair(fPairDifference,fPairYield);
      //  Update the global error code in fDifference, and use it
      //  to update the errors in fYield, in case blinder errors
      //  can propagate to the global error.
      fPairDifference.UpdateErrorFlag();
      fPairYield.UpdateErrorFlag(fPairDifference);
      if (! fBlinder.IsBlinderOkay()){
    fPairYield.UpdateErrorFlag(QwBlinder::kErrorFlag_BlinderFail);
    fPairDifference.UpdateErrorFlag(QwBlinder::kErrorFlag_BlinderFail);
      }
    }
    fPairAsymmetry.Ratio(fPairDifference,fPairYield);
    fPairAsymmetry.IncrementErrorCounters();
  }
}

References QwLog::endl(), fBlinder, fCurrentPatternNumber, fEvents, fHelicity, fIgnoreHelicity, fNextPair, fPairAsymmetry, fPairDifference, fPairIsGood, fPairYield, fPatternSize, QwBlinder::kErrorFlag_BlinderFail, QwDebug, QwError, and UpdateBlinder().

Referenced by PairAsymmetryIsGood().

Here is the call graph for this function:

Here is the caller graph for this function:

CalculateRunningAverage()

void QwHelicityPattern::CalculateRunningAverage

(

)

Definition at line 741 of file QwHelicityPattern.cc.

{
  fAsymmetry.CalculateRunningAverage();
  if (fEnableDifference){
    fDifference.CalculateRunningAverage();
  }
  fYield.CalculateRunningAverage();
  if (fEnableAlternateAsym) {
    fAsymmetry1.CalculateRunningAverage();
    fAsymmetry2.CalculateRunningAverage();
  }
  if (fEnablePairs) {
    //  Pair averages
    fPairYield.CalculateRunningAverage();
    fPairAsymmetry.CalculateRunningAverage();
    if (fEnableDifference){
      fPairDifference.CalculateRunningAverage();
    }
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, fEnablePairs, fPairAsymmetry, fPairDifference, fPairYield, and fYield.

Referenced by QwHelicityCorrelatedFeedback::CalculateRunningAverage(), and main().

Here is the caller graph for this function:

ClearEventData()

void QwHelicityPattern::ClearEventData

(

)

Clear event data and the vectors used for the calculation of. yields and asymmetries.

Definition at line 658 of file QwHelicityPattern.cc.

{
  fIgnoreHelicity = kFALSE;
  for(size_t i=0; i<fEvents.size(); i++)
    {
      fEvents[i].ClearEventData();
      fEventLoaded[i]=kFALSE;
      fHelicity[i]=-999;
    }
  fBlinder.ClearEventData();
 
  // Primary yield and asymmetry
  fYield.ClearEventData();
  fAsymmetry.ClearEventData();
  // Alternate asymmetries
  if (fEnableAlternateAsym){
    fAsymmetry1.ClearEventData();
    fAsymmetry2.ClearEventData();
  }
 
  fPositiveHelicitySum.ClearEventData();
  fNegativeHelicitySum.ClearEventData();
  fDifference.ClearEventData();
  fAlternateDiff.ClearEventData();
 
  fPairIsGood = kFALSE;
  fNextPair   = 0;
 
 
  fGoodPatterns = 0;
  fPatternIsGood = kFALSE;
  SetDataLoaded(kFALSE);
}

References fAlternateDiff, fAsymmetry, fAsymmetry1, fAsymmetry2, fBlinder, fDifference, fEnableAlternateAsym, fEventLoaded, fEvents, fGoodPatterns, fHelicity, fIgnoreHelicity, fNegativeHelicitySum, fNextPair, fPairIsGood, fPatternIsGood, fPositiveHelicitySum, fYield, and SetDataLoaded().

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), LoadEventData(), and main().

Here is the call graph for this function:

Here is the caller graph for this function:

ClearPairData()

void QwHelicityPattern::ClearPairData ( )

inline

Definition at line 70 of file QwHelicityPattern.h.

                        {
    fPairYield.ClearEventData();
    fPairDifference.ClearEventData();
    fPairAsymmetry.ClearEventData();
  }

References fPairAsymmetry, fPairDifference, and fPairYield.

Referenced by main().

Here is the caller graph for this function:

ConstructBranch() [1/2]

void QwHelicityPattern::ConstructBranch	(	TTree *	tree,
		TString &	prefix )

Definition at line 881 of file QwHelicityPattern.cc.

{
  TString basename = prefix(0, (prefix.First("|") >= 0)? prefix.First("|"): prefix.Length())+"BurstCounter";
  tree->Branch(basename,&fBurstCounter,basename+"/S");
 
  TString newprefix = "yield_" + prefix;
  fYield.ConstructBranch(tree, newprefix);
  newprefix = "asym_" + prefix;
  fAsymmetry.ConstructBranch(tree, newprefix);
 
  if (fEnableDifference) {
    newprefix = "diff_" + prefix;
    fDifference.ConstructBranch(tree, newprefix);
  }
  if (fEnableAlternateAsym) {
    newprefix = "asym1_" + prefix;
    fAsymmetry1.ConstructBranch(tree, newprefix);
    newprefix = "asym2_" + prefix;
    fAsymmetry2.ConstructBranch(tree, newprefix);
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fBurstCounter, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

ConstructBranch() [2/2]

void QwHelicityPattern::ConstructBranch	(	TTree *	tree,
		TString &	prefix,
		QwParameterFile &	trim_tree )

Definition at line 903 of file QwHelicityPattern.cc.

{
  TString basename = prefix(0, (prefix.First("|") >= 0)? prefix.First("|"): prefix.Length())+"BurstCounter";
  tree->Branch(basename,&fBurstCounter,basename+"/S");
  TString newprefix = "yield_" + prefix;
  fYield.ConstructBranch(tree, newprefix, trim_tree);
  newprefix = "asym_" + prefix;
  fAsymmetry.ConstructBranch(tree, newprefix, trim_tree);
 
  if (fEnableDifference) {
    newprefix = "diff_" + prefix;
    fDifference.ConstructBranch(tree, newprefix);
  }
  if (fEnableAlternateAsym) {
    newprefix = "asym1_" + prefix;
    fAsymmetry1.ConstructBranch(tree, newprefix, trim_tree);
    newprefix = "asym2_" + prefix;
    fAsymmetry2.ConstructBranch(tree, newprefix, trim_tree);
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fBurstCounter, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

ConstructBranchAndVector()

void QwHelicityPattern::ConstructBranchAndVector	(	TTree *	tree,
		TString &	prefix,
		QwRootTreeBranchVector &	values )

Definition at line 860 of file QwHelicityPattern.cc.

{
TString basename = prefix(0, (prefix.First("|") >= 0)? prefix.First("|"): prefix.Length())+"BurstCounter";
  tree->Branch(basename,&fBurstCounter,basename+"/S");
  TString newprefix = "yield_" + prefix;
  fYield.ConstructBranchAndVector(tree, newprefix, values);
  newprefix = "asym_" + prefix;
  fAsymmetry.ConstructBranchAndVector(tree, newprefix, values);
 
  if (fEnableDifference) {
    newprefix = "diff_" + prefix;
    fDifference.ConstructBranchAndVector(tree, newprefix, values);
  }
  if (fEnableAlternateAsym) {
    newprefix = "asym1_" + prefix;
    fAsymmetry1.ConstructBranchAndVector(tree, newprefix, values);
    newprefix = "asym2_" + prefix;
    fAsymmetry2.ConstructBranchAndVector(tree, newprefix, values);
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fBurstCounter, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

Referenced by QwHelicityCorrelatedFeedback::ConstructBranchAndVector().

Here is the caller graph for this function:

ConstructHistograms() [1/2]

void QwHelicityPattern::ConstructHistograms ( )

inline

Definition at line 180 of file QwHelicityPattern.h.

{ConstructHistograms((TDirectory*)NULL);};

References ConstructHistograms().

Referenced by ConstructHistograms().

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructHistograms() [2/2]

void QwHelicityPattern::ConstructHistograms

(

TDirectory *

folder

)

Definition at line 825 of file QwHelicityPattern.cc.

{
  TString prefix = "yield_";
  fYield.ConstructHistograms(folder,prefix);
  prefix = "asym_";
  fAsymmetry.ConstructHistograms(folder,prefix);
 
  if (fEnableDifference) {
    prefix = "diff_";
    fDifference.ConstructHistograms(folder,prefix);
  }
  if (fEnableAlternateAsym) {
    prefix = "asym1_";
    fAsymmetry1.ConstructHistograms(folder,prefix);
    prefix = "asym2_";
    fAsymmetry2.ConstructHistograms(folder,prefix);
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

ConstructObjects() [1/2]

void QwHelicityPattern::ConstructObjects ( )

inline

Definition at line 177 of file QwHelicityPattern.h.

{ConstructObjects((TDirectory*)NULL);};

References ConstructObjects().

Referenced by ConstructObjects().

Here is the call graph for this function:

Here is the caller graph for this function:

ConstructObjects() [2/2]

void QwHelicityPattern::ConstructObjects

(

TDirectory *

folder

)

Definition at line 802 of file QwHelicityPattern.cc.

{
  TString prefix = "blinder_";
  fBlinder.ConstructObjects(folder,prefix);
 
  prefix = "yield_";
  fYield.ConstructObjects(folder,prefix);
  prefix = "asym_";
  fAsymmetry.ConstructObjects(folder,prefix);
 
  if (fEnableDifference) {
    prefix = "diff_";
    fDifference.ConstructObjects(folder,prefix);
  }
  if (fEnableAlternateAsym) {
    prefix = "asym1_";
    fAsymmetry1.ConstructObjects(folder,prefix);
    prefix = "asym2_";
    fAsymmetry2.ConstructObjects(folder,prefix);
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fBlinder, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

DefineOptions()

void QwHelicityPattern::DefineOptions ( QwOptions & options )

static

Define the configuration options.

Defines configuration options using QwOptions functionality.

Parameters

options

Options object

Definition at line 33 of file QwHelicityPattern.cc.

{
  options.AddOptions("Helicity pattern")
    ("enable-burstsum", po::value<bool>()->default_bool_value(false),
     "enable burst sum calculation");
  options.AddOptions("Helicity pattern")
    ("enable-runningsum", po::value<bool>()->default_bool_value(true),
     "enable running sum calculation");
  options.AddOptions("Helicity pattern")
    ("enable-differences", po::value<bool>()->default_bool_value(false),
     "store pattern differences in tree");
  options.AddOptions("Helicity pattern")
    ("enable-alternateasym", po::value<bool>()->default_bool_value(false),
     "enable alternate asymmetries");
 
  options.AddOptions("Helicity pattern")
    ("print-burstsum", po::value<bool>()->default_bool_value(false),
     "print burst sum of subsystems");
  options.AddOptions("Helicity pattern")
    ("print-patternsum", po::value<bool>()->default_bool_value(false),
     "print pattern sum of subsystems");
  options.AddOptions("Helicity pattern")
    ("print-runningsum", po::value<bool>()->default_bool_value(false),
     "print running sum of subsystems");
 
  options.AddOptions("Helicity pattern")
    ("burstlength", po::value<int>()->default_value(9000),
     "number of patterns per burst");
 
  options.AddOptions("Helicity pattern")
    ("max-burst-index", po::value<int>()->default_value(0x7fffffff), // Default to max signed int
     "max number of bursts for a run");
 
  options.AddOptions("Helicity pattern")
    ("print-burst-index-map", po::value<bool>()->default_value(kFALSE),
     "whether to print the shorter max burst index if the final is shorter than min-burstlength");
 
  options.AddOptions("Helicity pattern")
    ("min-burstlength", po::value<int>()->default_value(2333), // Default 1/4 of 9000
     "minimum acceptable burst length");
 
  QwBlinder::DefineOptions(options);
}

References QwOptions::AddOptions(), and QwBlinder::DefineOptions().

Referenced by DefineOptionsParity().

Here is the call graph for this function:

Here is the caller graph for this function:

DisableAlternateAsymmetry()

void QwHelicityPattern::DisableAlternateAsymmetry ( )

inline

Disable alternate asymmetry calculation.

Definition at line 109 of file QwHelicityPattern.h.

{ fEnableAlternateAsym = kFALSE; };

References fEnableAlternateAsym.

DisableBurstSum()

void QwHelicityPattern::DisableBurstSum ( )

inline

Disable burst sum calculation.

Definition at line 116 of file QwHelicityPattern.h.

{ fEnableBurstSum = kFALSE; };

References fEnableBurstSum.

Referenced by ProcessOptions().

Here is the caller graph for this function:

DisableDifference()

void QwHelicityPattern::DisableDifference ( )

inline

Disable storing pattern differences.

Definition at line 130 of file QwHelicityPattern.h.

{ fEnableDifference = kFALSE; };

References fEnableDifference.

DisablePairs()

void QwHelicityPattern::DisablePairs ( )

inline

Disable storing pair differences.

Definition at line 137 of file QwHelicityPattern.h.

{ fEnablePairs = kFALSE; };

References fEnablePairs.

Referenced by main().

Here is the caller graph for this function:

DisableRunningSum()

void QwHelicityPattern::DisableRunningSum ( )

inline

Disable running sum calculation.

Definition at line 123 of file QwHelicityPattern.h.

{ fEnableRunningSum = kFALSE; };

References fEnableRunningSum.

EnableAlternateAsymmetry()

void QwHelicityPattern::EnableAlternateAsymmetry ( const Bool_t flag = kTRUE )

inline

Enable/disable alternate asymmetry calculation.

Definition at line 107 of file QwHelicityPattern.h.

{ fEnableAlternateAsym = flag; };

References fEnableAlternateAsym.

EnableBurstSum()

void QwHelicityPattern::EnableBurstSum ( const Bool_t flag = kTRUE )

inline

Enable/disable burst sum calculation.

Definition at line 114 of file QwHelicityPattern.h.

{ fEnableBurstSum = flag; };

References fEnableBurstSum.

EnableDifference()

void QwHelicityPattern::EnableDifference ( const Bool_t flag = kTRUE )

inline

Enable/disable storing pattern differences.

Definition at line 128 of file QwHelicityPattern.h.

{ fEnableDifference = flag; };

References fEnableDifference.

EnablePairs()

void QwHelicityPattern::EnablePairs ( const Bool_t flag = kTRUE )

inline

Enable/disable storing pair differences.

Definition at line 135 of file QwHelicityPattern.h.

{ fEnablePairs = flag; };

References fEnablePairs.

EnableRunningSum()

void QwHelicityPattern::EnableRunningSum ( const Bool_t flag = kTRUE )

inline

Enable/disable running sum calculation.

Definition at line 121 of file QwHelicityPattern.h.

{ fEnableRunningSum = flag; };

References fEnableRunningSum.

FillHistograms()

void QwHelicityPattern::FillHistograms

(

)

Definition at line 844 of file QwHelicityPattern.cc.

{
  if (fPatternIsGood) {
    fYield.FillHistograms();
    fAsymmetry.FillHistograms();
    if (fEnableDifference) {
      fDifference.FillHistograms();
    }
    if (fEnableAlternateAsym){
      fAsymmetry1.FillHistograms();
      fAsymmetry2.FillHistograms();
    }
  } else {
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, fPatternIsGood, and fYield.

FillTreeVector()

void QwHelicityPattern::FillTreeVector

(

QwRootTreeBranchVector &

values

)

const

Definition at line 924 of file QwHelicityPattern.cc.

{
  if (fPatternIsGood) {
    fYield.FillTreeVector(values);
    fAsymmetry.FillTreeVector(values);
    if (fEnableDifference) {
      fDifference.FillTreeVector(values);
    }
    if (fEnableAlternateAsym) {
      fAsymmetry1.FillTreeVector(values);
      fAsymmetry2.FillTreeVector(values);
    }
  }
}

References fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, fPatternIsGood, and fYield.

Referenced by QwHelicityCorrelatedFeedback::FillTreeVector().

Here is the caller graph for this function:

GetAsymmetry()

QwSubsystemArrayParity & QwHelicityPattern::GetAsymmetry ( )

inline

Definition at line 163 of file QwHelicityPattern.h.

{ return fAsymmetry; };

References fAsymmetry.

Referenced by VQwDataHandler::ConnectChannels().

Here is the caller graph for this function:

GetBurstCounter()

Short_t QwHelicityPattern::GetBurstCounter ( ) const

inline

Definition at line 217 of file QwHelicityPattern.h.

{return fBurstCounter;}

References fBurstCounter.

Referenced by main().

Here is the caller graph for this function:

GetDifference()

QwSubsystemArrayParity & QwHelicityPattern::GetDifference ( )

inline

Definition at line 162 of file QwHelicityPattern.h.

{ return fDifference; };

References fDifference.

Referenced by VQwDataHandler::ConnectChannels().

Here is the caller graph for this function:

GetEventcutErrorFlag()

UInt_t QwHelicityPattern::GetEventcutErrorFlag ( ) const

inline

Definition at line 203 of file QwHelicityPattern.h.

                                      {
    return fAsymmetry.GetEventcutErrorFlag();
  };

References fAsymmetry.

GetEventcutErrorFlagPointer()

const UInt_t * QwHelicityPattern::GetEventcutErrorFlagPointer ( ) const

inline

Definition at line 206 of file QwHelicityPattern.h.

                                                    {
    return fAsymmetry.GetEventcutErrorFlagPointer();
  };

References fAsymmetry.

Referenced by VQwDataHandler::SetPointer().

Here is the caller graph for this function:

GetPairAsymmetry()

QwSubsystemArrayParity & QwHelicityPattern::GetPairAsymmetry ( )

inline

Definition at line 168 of file QwHelicityPattern.h.

{ return fPairAsymmetry; };

References fPairAsymmetry.

Referenced by main().

Here is the caller graph for this function:

GetPairDifference()

QwSubsystemArrayParity & QwHelicityPattern::GetPairDifference ( )

inline

Definition at line 167 of file QwHelicityPattern.h.

{ return fPairDifference; };

References fPairDifference.

Referenced by main().

Here is the caller graph for this function:

GetPairYield()

QwSubsystemArrayParity & QwHelicityPattern::GetPairYield ( )

inline

Definition at line 166 of file QwHelicityPattern.h.

{ return fPairYield; };

References fPairYield.

Referenced by main().

Here is the caller graph for this function:

GetTargetChargeStat()

void QwHelicityPattern::GetTargetChargeStat	(	Double_t &	asym,
		Double_t &	error,
		Double_t &	width )

GetYield()

QwSubsystemArrayParity & QwHelicityPattern::GetYield ( )

inline

Definition at line 161 of file QwHelicityPattern.h.

{ return fYield; };

References fYield.

Referenced by VQwDataHandler::ConnectChannels().

Here is the caller graph for this function:

HasBurstData()

Bool_t QwHelicityPattern::HasBurstData ( )

inline

Definition at line 210 of file QwHelicityPattern.h.

{return fGoodPatterns>0;};

References fGoodPatterns.

Referenced by main().

Here is the caller graph for this function:

HasDataLoaded()

Bool_t QwHelicityPattern::HasDataLoaded ( ) const

inline

Definition at line 65 of file QwHelicityPattern.h.

{ return fIsDataLoaded; };

References fIsDataLoaded.

IncrementBurstCounter()

void QwHelicityPattern::IncrementBurstCounter ( )

inline

Definition at line 211 of file QwHelicityPattern.h.

                               {
    if (fBurstCounter < fMaxBurstIndex) {
      fBurstCounter++;
    }
    // Else we just park here and don't try to increment any more. This is a parameter from command line or map file
  }

References fBurstCounter, and fMaxBurstIndex.

Referenced by main().

Here is the caller graph for this function:

IsAlternateAsymEnabled()

Bool_t QwHelicityPattern::IsAlternateAsymEnabled ( )

inline

Status of alternate asymmetry calculation flag.

Definition at line 111 of file QwHelicityPattern.h.

{ return fEnableAlternateAsym; };

References fEnableAlternateAsym.

IsBurstSumEnabled()

Bool_t QwHelicityPattern::IsBurstSumEnabled ( )

inline

Status of burst sum calculation flag.

Definition at line 118 of file QwHelicityPattern.h.

{ return fEnableBurstSum; };

References fEnableBurstSum.

IsCompletePattern()

Bool_t QwHelicityPattern::IsCompletePattern

(

)

const

Check to see if the pattern is complete.

Definition at line 455 of file QwHelicityPattern.cc.

{
  Bool_t localdebug=kFALSE;
  Bool_t filled=kTRUE;
  Int_t i=fPatternSize-1;
  while(filled && i>-1)
    {
      if (localdebug){
        std::cout<<" i="<<i<<" is loaded ?"
                 <<fEventLoaded[fEvents.size()-i-1]<<"\n";
      }
      if(!fEventLoaded[i])
        filled=kFALSE;
      i--;
    }
 
 
 
  return filled;
}

References fEventLoaded, fEvents, and fPatternSize.

Referenced by IsGoodAsymmetry(), main(), and Print().

Here is the caller graph for this function:

IsDifferenceEnabled()

Bool_t QwHelicityPattern::IsDifferenceEnabled ( )

inline

Status of storing pattern differences flag.

Definition at line 132 of file QwHelicityPattern.h.

{ return fEnableDifference; };

References fEnableDifference.

IsEndOfBurst()

Bool_t QwHelicityPattern::IsEndOfBurst ( )

inline

Definition at line 98 of file QwHelicityPattern.h.

                       {
    //  Is this the end of a burst? And is this not the final burst?
    return (( fBurstLength > 0 && fGoodPatterns >= fBurstLength ) && ( fBurstCounter<fMaxBurstIndex ));
  }

References fBurstCounter, fBurstLength, fGoodPatterns, and fMaxBurstIndex.

Referenced by main().

Here is the caller graph for this function:

IsGoodAsymmetry()

Bool_t QwHelicityPattern::IsGoodAsymmetry

(

)

Definition at line 441 of file QwHelicityPattern.cc.

{
  Bool_t complete_and_good = kFALSE;
  if (IsCompletePattern()){
    CalculateAsymmetry();
    complete_and_good = fPatternIsGood;
  }
  return complete_and_good;
};

References CalculateAsymmetry(), fPatternIsGood, and IsCompletePattern().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

IsPairsEnabled()

Bool_t QwHelicityPattern::IsPairsEnabled ( )

inline

Status of storing pair differences flag.

Definition at line 139 of file QwHelicityPattern.h.

{ return fEnablePairs; };

References fEnablePairs.

IsRunningSumEnabled()

Bool_t QwHelicityPattern::IsRunningSumEnabled ( )

inline

Status of running sum calculation flag.

Definition at line 125 of file QwHelicityPattern.h.

{ return fEnableRunningSum; };

References fEnableRunningSum.

LoadEventData()

void QwHelicityPattern::LoadEventData

(

QwSubsystemArrayParity &

event

)

Load event data corresponding to the current pattern from the subsystems.

Definition at line 240 of file QwHelicityPattern.cc.

{
 
  Bool_t localdebug = kFALSE;
  fPatternIsGood = kFALSE;
 
  Long_t localEventNumber = -1;
  Long_t localPatternNumber = -1;
  Int_t  localPhaseNumber = -1;
  Int_t  localHelicityActual = -1;
  Bool_t localIgnoreHelicity = kFALSE;
  
  
  // Get the list of helicity subsystems
  if (! fHelicityIsMissing){
    std::vector<VQwSubsystem*> subsys_helicity = event.GetSubsystemByType("QwHelicity");
    QwHelicity* helicity = 0;
    
    if (subsys_helicity.size() > 0) {
      // Take the first helicity subsystem
      helicity = dynamic_cast<QwHelicity*>(subsys_helicity.at(0));
      if (helicity->HasDataLoaded()){
    localIgnoreHelicity = helicity->IsHelicityIgnored();
    // Get the event, pattern, phase number and helicity
    localEventNumber    = helicity->GetEventNumber();
    localPatternNumber  = helicity->GetPatternNumber();
    localPhaseNumber    = helicity->GetPhaseNumber() - helicity->GetMinPatternPhase(); // Use "reduced pattern phase", counts from 0.
    localHelicityActual = helicity->GetHelicityActual();
      } else {
    QwError << "QwHelicityPattern::LoadEventData:  The helicity subsystem does not have valid data!"
        << QwLog::endl;
      }
    } else {
      // We are not using any helicity subsystem
      static Bool_t user_has_been_warned = kFALSE;
      if (! user_has_been_warned) {
    QwError << "No helicity subsystem found!  Dropping to \"Missing Helicity\" mode!" << QwLog::endl;
    user_has_been_warned = kTRUE;
    fHelicityIsMissing = kTRUE;
      }
    }
  }
  if (fHelicityIsMissing){
    localIgnoreHelicity = kTRUE;
    localPatternNumber = fLastPatternNumber;
    localEventNumber   = fLastWindowNumber + 1;
    localPhaseNumber   = fLastPhaseNumber + 1;
    if(localPhaseNumber>=fPatternSize){
      localPatternNumber++;
      localPhaseNumber = 0; // Use "reduced pattern phase", counts from 0.
    }
    fLastWindowNumber  = localEventNumber;
    fLastPhaseNumber   = localPhaseNumber;
    fLastPatternNumber = localPatternNumber; 
  }
  if(localdebug) {
    std::cout<<"\n ###################################\n";
    std::cout<<"QwHelicityPattern::LoadEventData :: ";
    std::cout<<" event, pattern, phase # "<<localEventNumber<<" "<<localPatternNumber<<" "<<localPhaseNumber<<"\n";
    std::cout<<" helicity ="<< localHelicityActual<<"\n";
    for(size_t i=0; i<fEvents.size(); i++)
      std::cout<<i<<":"<<fEventLoaded[i]<<"  ";
    std::cout<<"\n";
  }
  if(fCurrentPatternNumber!=localPatternNumber){
    // new pattern
    ClearEventData();
    fCurrentPatternNumber=localPatternNumber;
  }
  if(localPhaseNumber<0){
    QwWarning << "QwHelicityPattern::LoadEventData:  "
        << "Reduced event phase number is less than zero; ignore this event."
        << QwLog::endl;
    ClearEventData();
  } else if(localPhaseNumber>=fPatternSize){
    QwWarning<<" In QwHelicityPattern::LoadEventData trying upload an event with a phase larger than expected \n"
       <<" phase ="<<localPhaseNumber+1<<" maximum expected phase="<<fPatternSize<<"\n"
       <<" operation impossible, pattern reset to 0: no asymmetries will be computed "<<QwLog::endl;
    ClearEventData();
  } else {
    if(localdebug){
      std::cout<<"QwHelicityPattern::LoadEventData local i="
           <<localPhaseNumber<<"\n";
    }
    fEvents[localPhaseNumber]      = event;
    fEventLoaded[localPhaseNumber] = kTRUE;
    fHelicity[localPhaseNumber]    = localHelicityActual;
    fEventNumber[localPhaseNumber] = localEventNumber;
    // Check to see if we should ignore the helicity; this is
    // reset to false in ClearEventData.
    fIgnoreHelicity |= localIgnoreHelicity;
    SetDataLoaded(kTRUE);
  }
  if(localdebug){
    Print();
  }
  return;
}

References ClearEventData(), QwLog::endl(), fCurrentPatternNumber, fEventLoaded, fEventNumber, fEvents, fHelicity, fHelicityIsMissing, fIgnoreHelicity, fLastPatternNumber, fLastPhaseNumber, fLastWindowNumber, fPatternIsGood, fPatternSize, QwHelicity::GetEventNumber(), QwHelicity::GetHelicityActual(), QwHelicity::GetMinPatternPhase(), QwHelicity::GetPatternNumber(), QwHelicity::GetPhaseNumber(), VQwSubsystem::HasDataLoaded(), QwHelicity::IsHelicityIgnored(), Print(), QwError, QwWarning, and SetDataLoaded().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

NextPairIsComplete()

Bool_t QwHelicityPattern::NextPairIsComplete

(

)

Definition at line 349 of file QwHelicityPattern.cc.

{
  Bool_t filled=kFALSE;
  // Potentially this function could check all possible pairs to see which are complete, if the initial
  // fNextPair is not complete.  The function must end with fNextPair equal to the index of the lowest pair
  // that is complete.
  if (fNextPair<fPatternSize/2){
    size_t firstevt  = fNextPair*2;
    size_t secondevt = firstevt + 1;
    filled = fEventLoaded.at(firstevt) && fEventLoaded.at(secondevt);
  }
  return (filled);
}

References fEventLoaded, fNextPair, and fPatternSize.

Referenced by PairAsymmetryIsGood().

Here is the caller graph for this function:

PairAsymmetryIsGood()

Bool_t QwHelicityPattern::PairAsymmetryIsGood

(

)

Definition at line 339 of file QwHelicityPattern.cc.

{
  Bool_t complete_and_good = kFALSE;
  if (NextPairIsComplete()){
    CalculatePairAsymmetry();  /*  Uses the same calculational variables as the pattern */
    complete_and_good = fPairIsGood;
  }
  return complete_and_good;
};

References CalculatePairAsymmetry(), fPairIsGood, and NextPairIsComplete().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

Print()

void QwHelicityPattern::Print

(

)

const

Definition at line 1013 of file QwHelicityPattern.cc.

{
  QwOut << "Pattern number = " << fCurrentPatternNumber << QwLog::endl;
  for (size_t i = 0; i < fPatternSize; i++)
    QwOut << "Event " << fEventNumber[i] << ": "
          << fEventLoaded[i] << ", " << fHelicity[i] << QwLog::endl;
  QwOut << "Is a complete pattern? (n/y:0/1) " << IsCompletePattern() << QwLog::endl;
}

References QwLog::endl(), fCurrentPatternNumber, fEventLoaded, fEventNumber, fHelicity, fPatternSize, IsCompletePattern(), and QwOut.

Referenced by LoadEventData().

Here is the call graph for this function:

Here is the caller graph for this function:

PrintIndexMapFile()

void QwHelicityPattern::PrintIndexMapFile ( Int_t runNum )

inline

Definition at line 78 of file QwHelicityPattern.h.

                                      {
    if ( fPrintIndexFile && (fGoodPatterns < fBurstMinGoodPatterns) ) {
      // Print a text map file to max_burst_index.####.conf
      Int_t maxBurst = 0;
      if (fBurstCounter==0) {
        // It's a single burst run anyway
        maxBurst = 0;
      }
      else {
        // It's a multi burst run and we want to merge the final two bursts
        maxBurst = fBurstCounter-1;
      }
      QwWarning << "Printing max_burst_index." << runNum << ".conf file with " << maxBurst << " max burst number" << QwLog::endl;
      std::ofstream output;
      output.open(Form("max_burst_index.%d.conf",runNum));
      output<< "max-burst-index=" << maxBurst << std::endl; // Print the current index before incrementing further, this will be the max index and the next pass will overflow this one instead of having another
      output.close();
    }
  }

References QwLog::endl(), fBurstCounter, fBurstMinGoodPatterns, fGoodPatterns, fPrintIndexFile, and QwWarning.

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

PrintValue()

void QwHelicityPattern::PrintValue

(

)

const

Definition at line 764 of file QwHelicityPattern.cc.

{
  // Pattern
  QwMessage << " Pattern yields                " << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fYield.PrintValue();
  QwMessage << " Pattern asymmetry " << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fAsymmetry.PrintValue();
  if (fEnableAlternateAsym) {
    QwMessage << " First/second half asymmetry "  << QwLog::endl;
    QwMessage << " ==============================" << QwLog::endl;
    fAsymmetry1.PrintValue();
    QwMessage << " Even/odd asymmetry " << QwLog::endl;
    QwMessage << " ==============================" << QwLog::endl;
    fAsymmetry2.PrintValue();
  }
  if (fEnableDifference){
    QwMessage << " Pattern difference " << QwLog::endl;
    QwMessage << " ==============================" << QwLog::endl;
    fDifference.PrintValue();
  }
  if (fEnablePairs){
    QwMessage << " Pair yield " << QwLog::endl;
    QwMessage << " ==============================" << QwLog::endl;
    fPairYield.PrintValue();
    QwMessage << " Pair asymmetry " << QwLog::endl;
    QwMessage << " ==============================" << QwLog::endl;
    fPairAsymmetry.PrintValue();
    if (fEnableDifference) {
      QwMessage << " Pair difference " << QwLog::endl;
      QwMessage << " ==============================" << QwLog::endl;
      fPairDifference.PrintValue();
    }
  }
}

References QwLog::endl(), fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, fEnablePairs, fPairAsymmetry, fPairDifference, fPairYield, fYield, and QwMessage.

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

ProcessOptions()

void QwHelicityPattern::ProcessOptions

(

QwOptions &

options

)

Process the configuration options.

Definition at line 78 of file QwHelicityPattern.cc.

{
  fEnableBurstSum   = options.GetValue<bool>("enable-burstsum");
  fEnableRunningSum = options.GetValue<bool>("enable-runningsum");
  fPrintBurstSum    = options.GetValue<bool>("print-burstsum");
  fPrintRunningSum  = options.GetValue<bool>("print-runningsum");
 
  fEnableDifference    = options.GetValue<bool>("enable-differences");
  fEnableAlternateAsym = options.GetValue<bool>("enable-alternateasym");
 
  fBurstLength = options.GetValue<int>("burstlength");
  if (fBurstLength == 0) DisableBurstSum();
 
  fMaxBurstIndex        = options.GetValue<int>("max-burst-index");
  fPrintIndexFile       = options.GetValue<bool>("print-burst-index-map");
  fBurstMinGoodPatterns = options.GetValue<int>("min-burstlength");
 
  if (fEnableAlternateAsym && fPatternSize <= 2){
    QwWarning << "QwHelicityPattern::ProcessOptions: "
          << "The 'enable-alternateasym' flag is disabled for pair analysis."
          << QwLog::endl;
    fEnableAlternateAsym = kFALSE;
  }
 
  fBlinder.ProcessOptions(options);
}

References DisableBurstSum(), QwLog::endl(), fBlinder, fBurstLength, fBurstMinGoodPatterns, fEnableAlternateAsym, fEnableBurstSum, fEnableDifference, fEnableRunningSum, fMaxBurstIndex, fPatternSize, fPrintBurstSum, fPrintIndexFile, fPrintRunningSum, QwOptions::GetValue(), and QwWarning.

Referenced by main(), and QwHelicityCorrelatedFeedback::ProcessOptions().

Here is the call graph for this function:

Here is the caller graph for this function:

SetDataLoaded()

void QwHelicityPattern::SetDataLoaded ( Bool_t flag )

inlineprotected

Definition at line 289 of file QwHelicityPattern.h.

{ fIsDataLoaded = flag; };

References fIsDataLoaded.

Referenced by ClearEventData(), and LoadEventData().

Here is the caller graph for this function:

UpdateBlinder() [1/3]

void QwHelicityPattern::UpdateBlinder ( )

inline

Update the blinder status using a random number generator.

Definition at line 148 of file QwHelicityPattern.h.

                      {
    fBlinder.Update();
  };

References fBlinder.

Referenced by CalculateAsymmetry(), CalculatePairAsymmetry(), and main().

Here is the caller graph for this function:

UpdateBlinder() [2/3]

void QwHelicityPattern::UpdateBlinder ( const QwEPICSEvent & epics )

inline

Update the blinder status with new external information.

Definition at line 156 of file QwHelicityPattern.h.

                                                {
    fBlinder.Update(epics);
  };

References fBlinder.

UpdateBlinder() [3/3]

void QwHelicityPattern::UpdateBlinder ( const QwSubsystemArrayParity & detectors )

inline

Update the blinder status with new external information.

Definition at line 152 of file QwHelicityPattern.h.

                                                              {
    fBlinder.Update(detectors);
  };

References fBlinder.

WritePromptSummary()

void QwHelicityPattern::WritePromptSummary

(

QwPromptSummary *

ps

)

Definition at line 1002 of file QwHelicityPattern.cc.

{
  ps->SetPatternSize(fPatternSize);
  fYield.WritePromptSummary(ps, "yield");
  fAsymmetry.WritePromptSummary(ps, "asymmetry");
  return;
};

References fAsymmetry, fPatternSize, fYield, and QwPromptSummary::SetPatternSize().

Referenced by main().

Here is the call graph for this function:

Here is the caller graph for this function:

Friends And Related Symbol Documentation

QwDataHandlerArray

friend class QwDataHandlerArray

friend

Definition at line 291 of file QwHelicityPattern.h.

References QwDataHandlerArray.

Referenced by QwDataHandlerArray.

Field Documentation

fAlternateDiff

QwSubsystemArrayParity QwHelicityPattern::fAlternateDiff

protected

Definition at line 272 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), ClearEventData(), QwHelicityPattern(), and QwHelicityPattern().

fAsymmetry

QwSubsystemArrayParity QwHelicityPattern::fAsymmetry

protected

Definition at line 245 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::AccumulateRunningSum(), AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), FillHistograms(), FillTreeVector(), GetAsymmetry(), GetEventcutErrorFlag(), GetEventcutErrorFlagPointer(), PrintValue(), QwHelicityPattern(), QwHelicityPattern(), and WritePromptSummary().

fAsymmetry1

QwSubsystemArrayParity QwHelicityPattern::fAsymmetry1

protected

Definition at line 248 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), FillHistograms(), FillTreeVector(), PrintValue(), QwHelicityPattern(), and QwHelicityPattern().

fAsymmetry2

QwSubsystemArrayParity QwHelicityPattern::fAsymmetry2

protected

Definition at line 249 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), FillHistograms(), FillTreeVector(), PrintValue(), QwHelicityPattern(), and QwHelicityPattern().

fBlinder

QwBlinder QwHelicityPattern::fBlinder

protected

Definition at line 233 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculatePairAsymmetry(), ClearEventData(), ConstructObjects(), ProcessOptions(), QwHelicityPattern(), UpdateBlinder(), UpdateBlinder(), and UpdateBlinder().

fBurstCounter

Short_t QwHelicityPattern::fBurstCounter

protected

Definition at line 263 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), GetBurstCounter(), IncrementBurstCounter(), IsEndOfBurst(), PrintIndexMapFile(), QwHelicityPattern(), and QwHelicityPattern().

fBurstLength

Int_t QwHelicityPattern::fBurstLength

protected

Definition at line 258 of file QwHelicityPattern.h.

Referenced by IsEndOfBurst(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fBurstMinGoodPatterns

Int_t QwHelicityPattern::fBurstMinGoodPatterns

protected

Definition at line 261 of file QwHelicityPattern.h.

Referenced by PrintIndexMapFile(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fCurrentPatternNumber

Int_t QwHelicityPattern::fCurrentPatternNumber

protected

Definition at line 228 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculatePairAsymmetry(), LoadEventData(), and Print().

fDifference

QwSubsystemArrayParity QwHelicityPattern::fDifference

protected

Definition at line 244 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), FillHistograms(), FillTreeVector(), GetDifference(), PrintValue(), QwHelicityPattern(), and QwHelicityPattern().

fEnableAlternateAsym

Bool_t QwHelicityPattern::fEnableAlternateAsym

protected

Definition at line 247 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), DisableAlternateAsymmetry(), EnableAlternateAsymmetry(), FillHistograms(), FillTreeVector(), IsAlternateAsymEnabled(), PrintValue(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fEnableBurstSum

Bool_t QwHelicityPattern::fEnableBurstSum

protected

Definition at line 264 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), DisableBurstSum(), EnableBurstSum(), IsBurstSumEnabled(), ProcessOptions(), QwHelicityCorrelatedFeedback::QwHelicityCorrelatedFeedback(), QwHelicityPattern(), and QwHelicityPattern().

fEnableDifference

Bool_t QwHelicityPattern::fEnableDifference

protected

Definition at line 271 of file QwHelicityPattern.h.

Referenced by AccumulatePairRunningSum(), AccumulateRunningSum(), CalculateRunningAverage(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), DisableDifference(), EnableDifference(), FillHistograms(), FillTreeVector(), IsDifferenceEnabled(), PrintValue(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fEnablePairs

Bool_t QwHelicityPattern::fEnablePairs

protected

Definition at line 252 of file QwHelicityPattern.h.

Referenced by CalculateRunningAverage(), DisablePairs(), EnablePairs(), IsPairsEnabled(), PrintValue(), and QwHelicityPattern().

fEnableRunningSum

Bool_t QwHelicityPattern::fEnableRunningSum

protected

Definition at line 268 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), DisableRunningSum(), EnableRunningSum(), IsRunningSumEnabled(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fEventLoaded

std::vector<Bool_t> QwHelicityPattern::fEventLoaded

protected

Definition at line 225 of file QwHelicityPattern.h.

Referenced by ClearEventData(), IsCompletePattern(), LoadEventData(), NextPairIsComplete(), and Print().

fEventNumber

std::vector<Int_t> QwHelicityPattern::fEventNumber

protected

Definition at line 227 of file QwHelicityPattern.h.

Referenced by LoadEventData(), and Print().

fEvents

std::vector<QwSubsystemArrayParity> QwHelicityPattern::fEvents

protected

Definition at line 224 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculatePairAsymmetry(), ClearEventData(), IsCompletePattern(), and LoadEventData().

fGoodPatterns

Int_t QwHelicityPattern::fGoodPatterns

protected

Definition at line 262 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), ClearEventData(), HasBurstData(), IsEndOfBurst(), PrintIndexMapFile(), QwHelicityPattern(), and QwHelicityPattern().

fHelicity

std::vector<Int_t> QwHelicityPattern::fHelicity

protected

Definition at line 226 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculatePairAsymmetry(), ClearEventData(), LoadEventData(), and Print().

fHelicityIsMissing

Bool_t QwHelicityPattern::fHelicityIsMissing

protected

This indicates if the subsystem arrays are missing the helicity object. It is updated once during initialization and once when processing the first event

Definition at line 237 of file QwHelicityPattern.h.

Referenced by LoadEventData(), and QwHelicityPattern().

fIgnoreHelicity

Bool_t QwHelicityPattern::fIgnoreHelicity

protected

This is true if any of the helicity objects of this pattern have indicated that we should ignore the helicity. It is updated every event and reset by ClearEventData.

Definition at line 240 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculatePairAsymmetry(), ClearEventData(), LoadEventData(), and QwHelicityPattern().

fIsDataLoaded

Bool_t QwHelicityPattern::fIsDataLoaded

protected

Definition at line 288 of file QwHelicityPattern.h.

Referenced by HasDataLoaded(), QwHelicityPattern(), QwHelicityPattern(), and SetDataLoaded().

fLastPatternNumber

Long_t QwHelicityPattern::fLastPatternNumber

protected

Definition at line 277 of file QwHelicityPattern.h.

Referenced by LoadEventData(), QwHelicityPattern(), and QwHelicityPattern().

fLastPhaseNumber

Int_t QwHelicityPattern::fLastPhaseNumber

protected

Definition at line 278 of file QwHelicityPattern.h.

Referenced by LoadEventData(), QwHelicityPattern(), and QwHelicityPattern().

fLastWindowNumber

Long_t QwHelicityPattern::fLastWindowNumber

protected

Definition at line 276 of file QwHelicityPattern.h.

Referenced by LoadEventData(), QwHelicityPattern(), and QwHelicityPattern().

fMaxBurstIndex

Int_t QwHelicityPattern::fMaxBurstIndex

protected

Definition at line 259 of file QwHelicityPattern.h.

Referenced by IncrementBurstCounter(), IsEndOfBurst(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fNegativeHelicitySum

QwSubsystemArrayParity QwHelicityPattern::fNegativeHelicitySum

protected

Definition at line 274 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), ClearEventData(), QwHelicityPattern(), and QwHelicityPattern().

fNextPair

size_t QwHelicityPattern::fNextPair

protected

Definition at line 280 of file QwHelicityPattern.h.

Referenced by CalculatePairAsymmetry(), ClearEventData(), NextPairIsComplete(), QwHelicityPattern(), and QwHelicityPattern().

fPairAsymmetry

QwSubsystemArrayParity QwHelicityPattern::fPairAsymmetry

protected

Definition at line 255 of file QwHelicityPattern.h.

Referenced by AccumulatePairRunningSum(), CalculatePairAsymmetry(), CalculateRunningAverage(), ClearPairData(), GetPairAsymmetry(), PrintValue(), QwHelicityPattern(), and QwHelicityPattern().

fPairDifference

QwSubsystemArrayParity QwHelicityPattern::fPairDifference

protected

Definition at line 254 of file QwHelicityPattern.h.

Referenced by AccumulatePairRunningSum(), CalculatePairAsymmetry(), CalculateRunningAverage(), ClearPairData(), GetPairDifference(), PrintValue(), QwHelicityPattern(), and QwHelicityPattern().

fPairIsGood

Bool_t QwHelicityPattern::fPairIsGood

protected

Definition at line 281 of file QwHelicityPattern.h.

Referenced by AccumulatePairRunningSum(), CalculatePairAsymmetry(), ClearEventData(), PairAsymmetryIsGood(), QwHelicityPattern(), and QwHelicityPattern().

fPairYield

QwSubsystemArrayParity QwHelicityPattern::fPairYield

protected

Definition at line 253 of file QwHelicityPattern.h.

Referenced by AccumulatePairRunningSum(), CalculatePairAsymmetry(), CalculateRunningAverage(), ClearPairData(), GetPairYield(), PrintValue(), QwHelicityPattern(), and QwHelicityPattern().

fPatternIsGood

Bool_t QwHelicityPattern::fPatternIsGood

protected

Definition at line 283 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), ClearEventData(), FillHistograms(), FillTreeVector(), IsGoodAsymmetry(), LoadEventData(), QwHelicityPattern(), and QwHelicityPattern().

fPatternSize

size_t QwHelicityPattern::fPatternSize

protected

Definition at line 229 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculatePairAsymmetry(), IsCompletePattern(), LoadEventData(), NextPairIsComplete(), Print(), ProcessOptions(), QwHelicityPattern(), QwHelicityPattern(), and WritePromptSummary().

fPositiveHelicitySum

QwSubsystemArrayParity QwHelicityPattern::fPositiveHelicitySum

protected

Definition at line 273 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), ClearEventData(), QwHelicityPattern(), and QwHelicityPattern().

fPrintBurstSum

Bool_t QwHelicityPattern::fPrintBurstSum

protected

Definition at line 265 of file QwHelicityPattern.h.

Referenced by ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fPrintIndexFile

Bool_t QwHelicityPattern::fPrintIndexFile

protected

Definition at line 260 of file QwHelicityPattern.h.

Referenced by PrintIndexMapFile(), ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fPrintRunningSum

Bool_t QwHelicityPattern::fPrintRunningSum

protected

Definition at line 269 of file QwHelicityPattern.h.

Referenced by ProcessOptions(), QwHelicityPattern(), and QwHelicityPattern().

fQuartetNumber

Int_t QwHelicityPattern::fQuartetNumber

protected

Definition at line 230 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), QwHelicityCorrelatedFeedback::LogHAParameters(), QwHelicityCorrelatedFeedback::LogParameters(), and QwHelicityCorrelatedFeedback::LogParameters().

fYield

QwSubsystemArrayParity QwHelicityPattern::fYield

protected

Definition at line 243 of file QwHelicityPattern.h.

Referenced by QwHelicityCorrelatedFeedback::AccumulateRunningSum(), AccumulateRunningSum(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ConstructBranch(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), ConstructObjects(), FillHistograms(), FillTreeVector(), GetYield(), PrintValue(), QwHelicityPattern(), QwHelicityPattern(), and WritePromptSummary().

run_label

TString QwHelicityPattern::run_label

protected

Definition at line 285 of file QwHelicityPattern.h.

---

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

• QwHelicityPattern.h
• QwHelicityPattern.cc