QwHelicityPattern Class Reference

#include <QwHelicityPattern.h>

Collaboration diagram for QwHelicityPattern:

Public Member Functions
	QwHelicityPattern (QwSubsystemArrayParity &event)
	Constructor with subsystem array. More...
virtual	~QwHelicityPattern ()
	Virtual destructor. More...
void	ProcessOptions (QwOptions &options)
	Process the configuration options. More...
void	LoadEventData (QwSubsystemArrayParity &event)
Bool_t	HasDataLoaded () const
Bool_t	IsCompletePattern () const
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. More...
void	DisableAlternateAsymmetry ()
	Disable alternate asymmetry calculation. More...
Bool_t	IsAlternateAsymEnabled ()
	Status of alternate asymmetry calculation flag. More...
void	EnableBurstSum (const Bool_t flag=kTRUE)
	Enable/disable burst sum calculation. More...
void	DisableBurstSum ()
	Disable burst sum calculation. More...
Bool_t	IsBurstSumEnabled ()
	Status of burst sum calculation flag. More...
void	EnableRunningSum (const Bool_t flag=kTRUE)
	Enable/disable running sum calculation. More...
void	DisableRunningSum ()
	Disable running sum calculation. More...
Bool_t	IsRunningSumEnabled ()
	Status of running sum calculation flag. More...
void	EnableDifference (const Bool_t flag=kTRUE)
	Enable/disable storing pattern differences. More...
void	DisableDifference ()
	Disable storing pattern differences. More...
Bool_t	IsDifferenceEnabled ()
	Status of storing pattern differences flag. More...
void	UpdateBlinder (QwParityDB *db)
	Update the blinder status with new external information. More...
void	UpdateBlinder (const QwSubsystemArrayParity &detectors)
	Update the blinder status with new external information. More...
void	UpdateBlinder (const QwEPICSEvent &epics)
	Update the blinder status with new external information. More...
QwSubsystemArrayParity &	GetBurstYield ()
QwSubsystemArrayParity &	GetBurstDifference ()
QwSubsystemArrayParity &	GetBurstAsymmetry ()
void	AccumulateBurstSum ()
void	AccumulateRunningBurstSum ()
void	AccumulateRunningSum ()
void	CalculateBurstAverage ()
void	CalculateRunningBurstAverage ()
void	CalculateRunningAverage ()
void	PrintRunningBurstAverage () const
void	PrintRunningAverage () const
void	PrintBurstAverage () const
void	ConstructHistograms ()
void	ConstructHistograms (TDirectory *folder)
void	FillHistograms ()
void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
void	ConstructBranch (TTree *tree, TString &prefix)
void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &trim_tree)
void	FillTreeVector (std::vector< Double_t > &values) const
void	FillDB (QwParityDB *db)
void	FillErrDB (QwParityDB *db)
void	WritePromptSummary (QwPromptSummary *ps)
Bool_t	IsGoodAsymmetry ()
UInt_t	GetEventcutErrorFlag () const
void	ClearEventData ()
void	ClearBurstSum ()
void	ClearRunningSum ()
void	Print () const

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

Protected Member Functions
void	SetDataLoaded (Bool_t flag)

Protected Attributes
Bool_t	fDEBUG
std::vector < QwSubsystemArrayParity >	fEvents
std::vector< Bool_t >	fEventLoaded
std::vector< Int_t >	fHelicity
std::vector< Int_t >	fEventNumber
Int_t	fCurrentPatternNumber
Int_t	fPatternSize
Int_t	fQuartetNumber
QwBlinder	fBlinder
Bool_t	fHelicityIsMissing
Bool_t	fIgnoreHelicity
QwSubsystemArrayParity	fYield
QwSubsystemArrayParity	fAsymmetry
Bool_t	fEnableAlternateAsym
QwSubsystemArrayParity	fAsymmetry1
QwSubsystemArrayParity	fAsymmetry2
Int_t	fBurstLength
Bool_t	fEnableBurstSum
Bool_t	fPrintBurstSum
QwSubsystemArrayParity	fBurstYield
QwSubsystemArrayParity	fBurstDifference
QwSubsystemArrayParity	fBurstAsymmetry
QwSubsystemArrayParity	fRunningBurstYield
QwSubsystemArrayParity	fRunningBurstDifference
QwSubsystemArrayParity	fRunningBurstAsymmetry
Bool_t	fEnableRunningSum
Bool_t	fPrintRunningSum
QwSubsystemArrayParity	fRunningYield
QwSubsystemArrayParity	fRunningDifference
QwSubsystemArrayParity	fRunningAsymmetry
QwSubsystemArrayParity	fRunningAsymmetry1
QwSubsystemArrayParity	fRunningAsymmetry2
Bool_t	fEnableDifference
QwSubsystemArrayParity	fDifference
QwSubsystemArrayParity	fAlternateDiff
QwSubsystemArrayParity	fPositiveHelicitySum
QwSubsystemArrayParity	fNegativeHelicitySum
Long_t	fLastWindowNumber
Long_t	fLastPatternNumber
Int_t	fLastPhaseNumber
Bool_t	fPatternIsGood
Bool_t	fIsDataLoaded

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

Friends
class	QwRegression

Detailed Description

Definition at line 29 of file QwHelicityPattern.h.

Constructor & Destructor Documentation

QwHelicityPattern::QwHelicityPattern ( )

private

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

QwHelicityPattern::QwHelicityPattern

(

QwSubsystemArrayParity &

event

)

Constructor with subsystem array.

Definition at line 80 of file QwHelicityPattern.cc.

References ClearBurstSum(), ClearRunningSum(), Qw::e, QwLog::endl(), fCurrentPatternNumber, fEventLoaded, fEventNumber, fEvents, fHelicity, fHelicityIsMissing, fPatternSize, fQuartetNumber, QwHelicity::GetMaxPatternPhase(), VQwSubsystem::GetSubsystemName(), QwError, QwMessage, and QwWarning.

  : fBlinder(),
    fHelicityIsMissing(kFALSE),   
    fIgnoreHelicity(kFALSE),
    fYield(event), 
    fAsymmetry(event),   
    fEnableAlternateAsym(kFALSE), 
    fAsymmetry1(event), 
    fAsymmetry2(event),
    fEnableBurstSum(kFALSE),      
    fPrintBurstSum(kFALSE),
    fBurstYield(event), 
    fBurstDifference(event), 
    fBurstAsymmetry(event),
    fRunningBurstYield(event), 
    fRunningBurstDifference(event), 
    fRunningBurstAsymmetry(event),
    fEnableRunningSum(kTRUE),     
    fPrintRunningSum(kFALSE),
    fRunningYield(event), 
    fRunningDifference(event), 
    fRunningAsymmetry(event),
    fRunningAsymmetry1(event), 
    fRunningAsymmetry2(event),
    fEnableDifference(kFALSE),
    fDifference(event), 
    fAlternateDiff(event),
    fPositiveHelicitySum(event), 
    fNegativeHelicitySum(event),
    fLastWindowNumber(0),
    fLastPatternNumber(0),
    fLastPhaseNumber(0)
{
  // 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->GetSubsystemName() << "."
                << 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;

          // Clear the burst sum
          ClearBurstSum();

          // Clear the running sum
          ClearRunningSum();

        }
      else
        {
          TString loc=
            "Standard exception from QwHelicityPattern : the pattern size has to be even;  right now pattern_size=";
          loc+=Form("%d",fPatternSize);
          throw std::invalid_argument(loc.Data());
        }
    }
  catch (std::exception& e)
    {
      std::cerr << e.what() << std::endl;
    }

}

Here is the call graph for this function:

virtual QwHelicityPattern::~QwHelicityPattern ( )

inlinevirtual

Virtual destructor.

Definition at line 42 of file QwHelicityPattern.h.

{ };

Member Function Documentation

void QwHelicityPattern::AccumulateBurstSum

(

)

Accumulate the burst sum by adding this helicity pattern to the burst sums of yield and difference. There is no burst sum of asymmetry, because that can only be calculated with meaningful moments at the end of a burst.

Definition at line 562 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::AccumulateRunningSum(), fBurstDifference, fBurstYield, fDifference, fPatternIsGood, and fYield.

Referenced by CalculateAsymmetry().

{
  if (fPatternIsGood){
    fBurstYield.AccumulateRunningSum(fYield);
    fBurstDifference.AccumulateRunningSum(fDifference);
    // The difference is blinded, so the burst difference is also blinded.
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::AccumulateRunningBurstSum

(

)

Accumulate the running burst sum by adding the current burst sum to the running sums of burst yield, difference and asymmetry.

Definition at line 597 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::AccumulateRunningSum(), fBurstAsymmetry, fBurstDifference, fBurstYield, fEnableRunningSum, fRunningBurstAsymmetry, fRunningBurstDifference, fRunningBurstYield, and QwSubsystemArrayParity::Ratio().

Referenced by main().

{
  // Accumulate the burst yield and difference
  if (fEnableRunningSum) {
    fRunningBurstYield.AccumulateRunningSum(fBurstYield);
    fRunningBurstDifference.AccumulateRunningSum(fBurstDifference);
    // The burst difference is blinded, so the running burst difference is also blinded.
  }

  // Calculate asymmetry over this entire burst
  fBurstAsymmetry.Ratio(fBurstDifference, fBurstYield);
  // Accumulate this burst asymmetry
  fRunningBurstAsymmetry.AccumulateRunningSum(fBurstAsymmetry);

  // Be sure to clear the burst sums after this function!
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::AccumulateRunningSum

(

)

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

Definition at line 576 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::AccumulateRunningSum(), fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, fPatternIsGood, fRunningAsymmetry, fRunningAsymmetry1, fRunningAsymmetry2, fRunningDifference, fRunningYield, and fYield.

Referenced by CalculateAsymmetry().

{
  if (fPatternIsGood){
    fRunningYield.AccumulateRunningSum(fYield);
    if (fEnableDifference){
      fRunningDifference.AccumulateRunningSum(fDifference);
      // The difference is blinded, so the running difference is also blinded.
    }
    fRunningAsymmetry.AccumulateRunningSum(fAsymmetry);
    if (fEnableAlternateAsym) {
      fRunningAsymmetry1.AccumulateRunningSum(fAsymmetry1);
      fRunningAsymmetry2.AccumulateRunningSum(fAsymmetry2);
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::CalculateAsymmetry

(

)

Calculate asymmetries for the current pattern.

Definition at line 310 of file QwHelicityPattern.cc.

References AccumulateBurstSum(), AccumulateRunningSum(), QwBlinder::Blind(), QwSubsystemArray::ClearEventData(), ClearEventData(), QwSubsystemArrayParity::Difference(), QwLog::endl(), fAlternateDiff, fAsymmetry, fAsymmetry1, fAsymmetry2, fBlinder, fCurrentPatternNumber, fDifference, fEnableAlternateAsym, fEnableBurstSum, fEnableRunningSum, fEvents, fHelicity, fIgnoreHelicity, fNegativeHelicitySum, fPatternIsGood, fPatternSize, fPositiveHelicitySum, fQuartetNumber, fYield, QwSubsystemArrayParity::IncrementErrorCounters(), QwDebug, QwError, QwSubsystemArrayParity::Ratio(), QwSubsystemArrayParity::Scale(), QwSubsystemArrayParity::Sum(), and QwSubsystemArrayParity::UpdateErrorFlag().

Referenced by main().

{

  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 < (size_t) fPatternSize; i++) {
      Int_t localhel = 1;
      for (size_t j = 0; j < (size_t) 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 < (size_t) 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();
  i = fPatternSize;
  checkhel = -9999;
  // 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 assymetry \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){
      //  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);
    }
    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 < (size_t) 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 < (size_t) 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);
      }
    }

    // Accumulate the burst and running sums
    if (fEnableBurstSum)   AccumulateBurstSum();
    if (fEnableRunningSum) AccumulateRunningSum();
    
    if (localdebug) QwDebug << " pattern number =" << fQuartetNumber << QwLog::endl;
  }

  return;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::CalculateBurstAverage

(

)

Calculate the average burst yield, difference and asymmetry.

Definition at line 618 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::CalculateRunningAverage(), fBurstAsymmetry, fBurstDifference, fBurstYield, fPrintBurstSum, and PrintBurstAverage().

Referenced by main().

{
  fBurstAsymmetry.CalculateRunningAverage();
  fBurstDifference.CalculateRunningAverage();
  fBurstYield.CalculateRunningAverage();

  if (fPrintBurstSum) PrintBurstAverage();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::CalculateRunningAverage

(

)

Definition at line 641 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::CalculateRunningAverage(), fEnableAlternateAsym, fEnableDifference, fPrintRunningSum, fRunningAsymmetry, fRunningAsymmetry1, fRunningAsymmetry2, fRunningDifference, fRunningYield, and PrintRunningAverage().

Referenced by main().

{
  fRunningAsymmetry.CalculateRunningAverage();
  if (fEnableDifference){
    fRunningDifference.CalculateRunningAverage();
  }
  fRunningYield.CalculateRunningAverage();
  if (fEnableAlternateAsym) {
    fRunningAsymmetry1.CalculateRunningAverage();
    fRunningAsymmetry2.CalculateRunningAverage();
  }
  if (fPrintRunningSum) PrintRunningAverage();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::CalculateRunningBurstAverage

(

)

Calculate the average running burst yield, difference and asymmetry.

Definition at line 631 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::CalculateRunningAverage(), fPrintBurstSum, fRunningBurstAsymmetry, fRunningBurstDifference, fRunningBurstYield, and PrintRunningBurstAverage().

Referenced by main().

{
  fRunningBurstAsymmetry.CalculateRunningAverage();
  fRunningBurstDifference.CalculateRunningAverage();
  fRunningBurstYield.CalculateRunningAverage();

  if (fPrintBurstSum) PrintRunningBurstAverage();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::ClearBurstSum

(

)

Clear the burst sums of yield and difference. No asymmetry burst sum is used.

Definition at line 546 of file QwHelicityPattern.cc.

References QwSubsystemArray::ClearEventData(), fBurstAsymmetry, fBurstDifference, fBurstYield, and fEnableBurstSum.

Referenced by main(), and QwHelicityPattern().

{
  if (fEnableBurstSum) {
    fBurstYield.ClearEventData();
    fBurstDifference.ClearEventData();
    fBurstAsymmetry.ClearEventData();
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::ClearEventData

(

)

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

Definition at line 487 of file QwHelicityPattern.cc.

References QwBlinder::ClearEventData(), QwSubsystemArray::ClearEventData(), fAlternateDiff, fAsymmetry, fAsymmetry1, fAsymmetry2, fBlinder, fDifference, fEnableAlternateAsym, fEventLoaded, fEvents, fHelicity, fIgnoreHelicity, fNegativeHelicitySum, fPatternIsGood, fPositiveHelicitySum, fYield, and SetDataLoaded().

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

{
  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();

  fPatternIsGood = kFALSE;
  SetDataLoaded(kFALSE);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::ClearRunningSum

(

)

Clear the running sums of yield, difference and asymmetry. Also clear the running burst sums if enabled.

Definition at line 521 of file QwHelicityPattern.cc.

References QwSubsystemArray::ClearEventData(), fEnableAlternateAsym, fEnableBurstSum, fEnableRunningSum, fRunningAsymmetry, fRunningAsymmetry1, fRunningAsymmetry2, fRunningBurstAsymmetry, fRunningBurstDifference, fRunningBurstYield, fRunningDifference, and fRunningYield.

Referenced by main(), and QwHelicityPattern().

{
  if (fEnableRunningSum) {
    fRunningYield.ClearEventData();
    fRunningDifference.ClearEventData();
    fRunningAsymmetry.ClearEventData();
    // Running alternate asymmetries
    if (fEnableAlternateAsym) {
      fRunningAsymmetry1.ClearEventData();
      fRunningAsymmetry2.ClearEventData();
    }
    // Running burst sums
    if (fEnableBurstSum) {
      fRunningBurstYield.ClearEventData();
      fRunningBurstDifference.ClearEventData();
      fRunningBurstAsymmetry.ClearEventData();
    }
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 762 of file QwHelicityPattern.cc.

References QwSubsystemArray::ConstructBranch(), fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

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

Here is the call graph for this function:

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

Definition at line 781 of file QwHelicityPattern.cc.

References QwSubsystemArray::ConstructBranch(), fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

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

Here is the call graph for this function:

void QwHelicityPattern::ConstructBranchAndVector	(	TTree *	tree,
		TString &	prefix,
		std::vector< Double_t > &	values
	)

Definition at line 743 of file QwHelicityPattern.cc.

References QwSubsystemArrayParity::ConstructBranchAndVector(), fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

Referenced by main().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::ConstructHistograms ( )

inline

Definition at line 120 of file QwHelicityPattern.h.

References ConstructHistograms().

Referenced by ConstructHistograms(), and main().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::ConstructHistograms

(

TDirectory *

folder

)

Definition at line 709 of file QwHelicityPattern.cc.

References QwSubsystemArray::ConstructHistograms(), fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, and fYield.

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

Here is the call graph for this function:

void QwHelicityPattern::DefineOptions ( QwOptions &  options )

static

Define the configuration options.

Defines configuration options using QwOptions functionality.

Parameters

options

Options object

Definition at line 26 of file QwHelicityPattern.cc.

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

Referenced by DefineOptionsParity().

{
  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-runningsum", po::value<bool>()->default_bool_value(false),
     "print running sum of subsystems");

  options.AddOptions("Helicity pattern")
    ("burstlength", po::value<int>()->default_value(240),
     "number of patterns per burst");

  QwBlinder::DefineOptions(options);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::DisableAlternateAsymmetry ( )

inline

Disable alternate asymmetry calculation.

Definition at line 65 of file QwHelicityPattern.h.

References fEnableAlternateAsym.

{ fEnableAlternateAsym = kFALSE; };

void QwHelicityPattern::DisableBurstSum ( )

inline

Disable burst sum calculation.

Definition at line 72 of file QwHelicityPattern.h.

References fEnableBurstSum.

Referenced by ProcessOptions().

{ fEnableBurstSum = kFALSE; };

Here is the caller graph for this function:

void QwHelicityPattern::DisableDifference ( )

inline

Disable storing pattern differences.

Definition at line 86 of file QwHelicityPattern.h.

References fEnableDifference.

{ fEnableDifference = kFALSE; };

void QwHelicityPattern::DisableRunningSum ( )

inline

Disable running sum calculation.

Definition at line 79 of file QwHelicityPattern.h.

References fEnableRunningSum.

{ fEnableRunningSum = kFALSE; };

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

inline

Enable/disable alternate asymmetry calculation.

Definition at line 63 of file QwHelicityPattern.h.

References fEnableAlternateAsym.

{ fEnableAlternateAsym = flag; };

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

inline

Enable/disable burst sum calculation.

Definition at line 70 of file QwHelicityPattern.h.

References fEnableBurstSum.

{ fEnableBurstSum = flag; };

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

inline

Enable/disable storing pattern differences.

Definition at line 84 of file QwHelicityPattern.h.

References fEnableDifference.

{ fEnableDifference = flag; };

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

inline

Enable/disable running sum calculation.

Definition at line 77 of file QwHelicityPattern.h.

References fEnableRunningSum.

{ fEnableRunningSum = flag; };

void QwHelicityPattern::FillDB

(

QwParityDB *

db

)

Definition at line 816 of file QwHelicityPattern.cc.

References fBlinder, fEnableAlternateAsym, fEnableDifference, QwSubsystemArrayParity::FillDB(), QwBlinder::FillDB(), fRunningAsymmetry, fRunningAsymmetry1, fRunningAsymmetry2, fRunningDifference, and fRunningYield.

Referenced by main().

{
  fBlinder.FillDB(db,"");

  fRunningYield.FillDB(db, "yield");
  fRunningAsymmetry.FillDB(db, "asymmetry");
  if (fEnableDifference) {
    fRunningDifference.FillDB(db, "difference");
  }
  if (fEnableAlternateAsym) {
    fRunningAsymmetry1.FillDB(db, "asymmetry1");
    fRunningAsymmetry2.FillDB(db, "asymmetry2");
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::FillErrDB

(

QwParityDB *

db

)

Definition at line 832 of file QwHelicityPattern.cc.

References fAsymmetry, fBlinder, QwSubsystemArrayParity::FillErrDB(), and QwBlinder::FillErrDB().

Referenced by main().

{
  fBlinder.FillErrDB(db,"");
  fAsymmetry.FillErrDB(db, "");
  return;
};

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::FillHistograms

(

)

Definition at line 728 of file QwHelicityPattern.cc.

References fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, QwSubsystemArrayParity::FillHistograms(), fPatternIsGood, and fYield.

Referenced by main().

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::FillTreeVector

(

std::vector< Double_t > &

values

)

const

Definition at line 800 of file QwHelicityPattern.cc.

References fAsymmetry, fAsymmetry1, fAsymmetry2, fDifference, fEnableAlternateAsym, fEnableDifference, QwSubsystemArrayParity::FillTreeVector(), fPatternIsGood, and fYield.

Referenced by main().

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

Here is the call graph for this function:

Here is the caller graph for this function:

QwSubsystemArrayParity& QwHelicityPattern::GetBurstAsymmetry ( )

inline

Definition at line 106 of file QwHelicityPattern.h.

References fBurstAsymmetry.

Referenced by main().

{ return fBurstAsymmetry; };

Here is the caller graph for this function:

QwSubsystemArrayParity& QwHelicityPattern::GetBurstDifference ( )

inline

Definition at line 105 of file QwHelicityPattern.h.

References fBurstDifference.

Referenced by main().

{ return fBurstDifference; };

Here is the caller graph for this function:

QwSubsystemArrayParity& QwHelicityPattern::GetBurstYield ( )

inline

Definition at line 104 of file QwHelicityPattern.h.

References fBurstYield.

Referenced by main().

{ return fBurstYield; };

Here is the caller graph for this function:

UInt_t QwHelicityPattern::GetEventcutErrorFlag ( ) const

inline

Definition at line 133 of file QwHelicityPattern.h.

References fAsymmetry, and QwSubsystemArrayParity::GetEventcutErrorFlag().

Referenced by QwRegression::LinearRegression().

                                     {
    return fAsymmetry.GetEventcutErrorFlag();
  };

Here is the call graph for this function:

Here is the caller graph for this function:

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

Bool_t QwHelicityPattern::HasDataLoaded ( ) const

inline

Definition at line 50 of file QwHelicityPattern.h.

References fIsDataLoaded.

{ return fIsDataLoaded; };

Bool_t QwHelicityPattern::IsAlternateAsymEnabled ( )

inline

Status of alternate asymmetry calculation flag.

Definition at line 67 of file QwHelicityPattern.h.

References fEnableAlternateAsym.

{ return fEnableAlternateAsym; };

Bool_t QwHelicityPattern::IsBurstSumEnabled ( )

inline

Status of burst sum calculation flag.

Definition at line 74 of file QwHelicityPattern.h.

References fEnableBurstSum.

Referenced by main().

{ return fEnableBurstSum; };

Here is the caller graph for this function:

Bool_t QwHelicityPattern::IsCompletePattern

(

)

const

Check to see if the pattern is complete.

Definition at line 284 of file QwHelicityPattern.cc.

References fEventLoaded, fEvents, and fPatternSize.

Referenced by main(), and Print().

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

Here is the caller graph for this function:

Bool_t QwHelicityPattern::IsDifferenceEnabled ( )

inline

Status of storing pattern differences flag.

Definition at line 88 of file QwHelicityPattern.h.

References fEnableDifference.

{ return fEnableDifference; };

Bool_t QwHelicityPattern::IsEndOfBurst ( )

inline

Definition at line 54 of file QwHelicityPattern.h.

References fBurstLength, and fCurrentPatternNumber.

Referenced by main().

                       {
    //  Is this the end of a burst?
    return (fBurstLength > 0 && fCurrentPatternNumber % fBurstLength == 0);
  }

Here is the caller graph for this function:

Bool_t QwHelicityPattern::IsGoodAsymmetry ( )

inline

Definition at line 132 of file QwHelicityPattern.h.

References fPatternIsGood.

Referenced by main().

{ return fPatternIsGood;};

Here is the caller graph for this function:

Bool_t QwHelicityPattern::IsRunningSumEnabled ( )

inline

Status of running sum calculation flag.

Definition at line 81 of file QwHelicityPattern.h.

References fEnableRunningSum.

Referenced by main().

{ return fEnableRunningSum; };

Here is the caller graph for this function:

void QwHelicityPattern::LoadEventData

(

QwSubsystemArrayParity &

event

)

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

Definition at line 181 of file QwHelicityPattern.cc.

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

{

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

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::Print

(

)

const

Definition at line 851 of file QwHelicityPattern.cc.

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

Referenced by LoadEventData().

{
  QwOut << "Pattern number = " << fCurrentPatternNumber << QwLog::endl;
  for (Int_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;
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::PrintBurstAverage

(

)

const

Definition at line 693 of file QwHelicityPattern.cc.

References QwLog::endl(), fBurstAsymmetry, fBurstDifference, fBurstYield, QwSubsystemArrayParity::PrintValue(), and QwMessage.

Referenced by CalculateBurstAverage().

{
  QwMessage << " Burst average of asymmetry  " << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fBurstAsymmetry.PrintValue();

  QwMessage << " Burst average of difference " << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fBurstDifference.PrintValue();

  QwMessage << " Burst average of yields     " << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fBurstYield.PrintValue();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::PrintRunningAverage

(

)

const

Definition at line 672 of file QwHelicityPattern.cc.

References QwLog::endl(), fEnableAlternateAsym, fEnableDifference, fRunningAsymmetry, fRunningAsymmetry1, fRunningAsymmetry2, fRunningDifference, fRunningYield, QwSubsystemArrayParity::PrintValue(), and QwMessage.

Referenced by CalculateRunningAverage(), and main().

{
  QwMessage << "QwHelicityPattern::PrintRunningAverage() const " << QwLog::endl;
  //  QwMessage << " Running average of yields     " << QwLog::endl;
  fRunningYield.PrintValue();
  //  QwMessage << " Running average of asymmetry  " << QwLog::endl;
  fRunningAsymmetry.PrintValue();
  if (fEnableAlternateAsym) {
    //    QwMessage << " Running average of first half/second half asymmetry  "  << QwLog::endl;
    fRunningAsymmetry1.PrintValue();
    //    QwMessage << " Running average of even/odd asymmetry  " << QwLog::endl;
    fRunningAsymmetry2.PrintValue();
  }
  if (fEnableDifference){
    //    QwMessage << " Running average of difference " << QwLog::endl;
    fRunningDifference.PrintValue();
  }

}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::PrintRunningBurstAverage

(

)

const

Definition at line 656 of file QwHelicityPattern.cc.

References QwLog::endl(), fRunningBurstAsymmetry, fRunningBurstDifference, fRunningBurstYield, QwSubsystemArrayParity::PrintValue(), and QwMessage.

Referenced by CalculateRunningBurstAverage(), and main().

{
  QwMessage << " Running burst average of asymmetry" << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fRunningBurstAsymmetry.PrintValue();

  QwMessage << " Running burst average of difference" << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fRunningBurstDifference.PrintValue();

  QwMessage << " Running burst average of yields" << QwLog::endl;
  QwMessage << " ==============================" << QwLog::endl;
  fRunningBurstYield.PrintValue();
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::ProcessOptions

(

QwOptions &

options

)

Process the configuration options.

Definition at line 56 of file QwHelicityPattern.cc.

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

Referenced by main().

{
  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();

  if (fEnableAlternateAsym && fPatternSize <= 2){
    QwWarning << "QwHelicityPattern::ProcessOptions: "
        << "The 'enable-alternateasym' flag is disabled for pair analysis."
        << QwLog::endl;
    fEnableAlternateAsym = kFALSE;
  }

  fBlinder.ProcessOptions(options);
}

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::SetDataLoaded ( Bool_t  flag )

inlineprotected

Definition at line 208 of file QwHelicityPattern.h.

References fIsDataLoaded.

Referenced by ClearEventData(), and LoadEventData().

{ fIsDataLoaded = flag; };

Here is the caller graph for this function:

void QwHelicityPattern::UpdateBlinder ( QwParityDB *  db )

inline

Update the blinder status with new external information.

Definition at line 91 of file QwHelicityPattern.h.

References fBlinder, and QwBlinder::Update().

Referenced by main().

                                    {
    fBlinder.Update(db);
  };

Here is the call graph for this function:

Here is the caller graph for this function:

void QwHelicityPattern::UpdateBlinder ( const QwSubsystemArrayParity &  detectors )

inline

Update the blinder status with new external information.

Definition at line 95 of file QwHelicityPattern.h.

References fBlinder, and QwBlinder::Update().

                                                              {
    fBlinder.Update(detectors);
  };

Here is the call graph for this function:

void QwHelicityPattern::UpdateBlinder ( const QwEPICSEvent &  epics )

inline

Update the blinder status with new external information.

Definition at line 99 of file QwHelicityPattern.h.

References fBlinder, and QwBlinder::Update().

                                                {
    fBlinder.Update(epics);
  };

Here is the call graph for this function:

void QwHelicityPattern::WritePromptSummary

(

QwPromptSummary *

ps

)

Definition at line 840 of file QwHelicityPattern.cc.

References fRunningAsymmetry, fRunningYield, and QwSubsystemArrayParity::WritePromptSummary().

Referenced by main().

{
  fRunningYield.WritePromptSummary(ps, "yield");
  fRunningAsymmetry.WritePromptSummary(ps, "asymmetry");
  
  return;
};

Here is the call graph for this function:

Here is the caller graph for this function:

Friends And Related Function Documentation

friend class QwRegression

friend

Definition at line 208 of file QwHelicityPattern.h.

Field Documentation

QwSubsystemArrayParity QwHelicityPattern::fAlternateDiff

protected

Definition at line 196 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), and ClearEventData().

QwSubsystemArrayParity QwHelicityPattern::fAsymmetry

protected

Definition at line 168 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateAsymmetry(), ClearEventData(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), FillErrDB(), FillHistograms(), FillTreeVector(), GetEventcutErrorFlag(), and QwRegression::QwRegression().

QwSubsystemArrayParity QwHelicityPattern::fAsymmetry1

protected

Definition at line 171 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateAsymmetry(), ClearEventData(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), and FillTreeVector().

QwSubsystemArrayParity QwHelicityPattern::fAsymmetry2

protected

Definition at line 172 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateAsymmetry(), ClearEventData(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), and FillTreeVector().

QwBlinder QwHelicityPattern::fBlinder

protected

Definition at line 157 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), ClearEventData(), FillDB(), FillErrDB(), ProcessOptions(), and UpdateBlinder().

QwSubsystemArrayParity QwHelicityPattern::fBurstAsymmetry

protected

Definition at line 180 of file QwHelicityPattern.h.

Referenced by AccumulateRunningBurstSum(), CalculateBurstAverage(), ClearBurstSum(), GetBurstAsymmetry(), and PrintBurstAverage().

QwSubsystemArrayParity QwHelicityPattern::fBurstDifference

protected

Definition at line 179 of file QwHelicityPattern.h.

Referenced by AccumulateBurstSum(), AccumulateRunningBurstSum(), CalculateBurstAverage(), ClearBurstSum(), GetBurstDifference(), and PrintBurstAverage().

Int_t QwHelicityPattern::fBurstLength

protected

Definition at line 175 of file QwHelicityPattern.h.

Referenced by IsEndOfBurst(), and ProcessOptions().

QwSubsystemArrayParity QwHelicityPattern::fBurstYield

protected

Definition at line 178 of file QwHelicityPattern.h.

Referenced by AccumulateBurstSum(), AccumulateRunningBurstSum(), CalculateBurstAverage(), ClearBurstSum(), GetBurstYield(), and PrintBurstAverage().

Int_t QwHelicityPattern::fCurrentPatternNumber

protected

Definition at line 152 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), IsEndOfBurst(), LoadEventData(), Print(), and QwHelicityPattern().

Bool_t QwHelicityPattern::fDEBUG

protected

Definition at line 144 of file QwHelicityPattern.h.

QwSubsystemArrayParity QwHelicityPattern::fDifference

protected

Definition at line 195 of file QwHelicityPattern.h.

Referenced by AccumulateBurstSum(), AccumulateRunningSum(), CalculateAsymmetry(), ClearEventData(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), and QwRegression::QwRegression().

Bool_t QwHelicityPattern::fEnableAlternateAsym

protected

Definition at line 170 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateAsymmetry(), CalculateRunningAverage(), ClearEventData(), ClearRunningSum(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), DisableAlternateAsymmetry(), EnableAlternateAsymmetry(), FillDB(), FillHistograms(), FillTreeVector(), IsAlternateAsymEnabled(), PrintRunningAverage(), and ProcessOptions().

Bool_t QwHelicityPattern::fEnableBurstSum

protected

Definition at line 176 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), ClearBurstSum(), ClearRunningSum(), DisableBurstSum(), EnableBurstSum(), IsBurstSumEnabled(), and ProcessOptions().

Bool_t QwHelicityPattern::fEnableDifference

protected

Definition at line 194 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), DisableDifference(), EnableDifference(), FillDB(), FillHistograms(), FillTreeVector(), IsDifferenceEnabled(), PrintRunningAverage(), and ProcessOptions().

Bool_t QwHelicityPattern::fEnableRunningSum

protected

Definition at line 186 of file QwHelicityPattern.h.

Referenced by AccumulateRunningBurstSum(), CalculateAsymmetry(), ClearRunningSum(), DisableRunningSum(), EnableRunningSum(), IsRunningSumEnabled(), and ProcessOptions().

std::vector<Bool_t> QwHelicityPattern::fEventLoaded

protected

Definition at line 149 of file QwHelicityPattern.h.

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

std::vector<Int_t> QwHelicityPattern::fEventNumber

protected

Definition at line 151 of file QwHelicityPattern.h.

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

std::vector<QwSubsystemArrayParity> QwHelicityPattern::fEvents

protected

Definition at line 148 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), ClearEventData(), IsCompletePattern(), LoadEventData(), and QwHelicityPattern().

std::vector<Int_t> QwHelicityPattern::fHelicity

protected

Definition at line 150 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), ClearEventData(), LoadEventData(), Print(), and QwHelicityPattern().

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 161 of file QwHelicityPattern.h.

Referenced by LoadEventData(), and QwHelicityPattern().

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 164 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), ClearEventData(), and LoadEventData().

Bool_t QwHelicityPattern::fIsDataLoaded

protected

Definition at line 207 of file QwHelicityPattern.h.

Referenced by HasDataLoaded(), and SetDataLoaded().

Long_t QwHelicityPattern::fLastPatternNumber

protected

Definition at line 201 of file QwHelicityPattern.h.

Referenced by LoadEventData().

Int_t QwHelicityPattern::fLastPhaseNumber

protected

Definition at line 202 of file QwHelicityPattern.h.

Referenced by LoadEventData().

Long_t QwHelicityPattern::fLastWindowNumber

protected

Definition at line 200 of file QwHelicityPattern.h.

Referenced by LoadEventData().

QwSubsystemArrayParity QwHelicityPattern::fNegativeHelicitySum

protected

Definition at line 198 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), and ClearEventData().

Bool_t QwHelicityPattern::fPatternIsGood

protected

Definition at line 204 of file QwHelicityPattern.h.

Referenced by AccumulateBurstSum(), AccumulateRunningSum(), CalculateAsymmetry(), ClearEventData(), FillHistograms(), FillTreeVector(), IsGoodAsymmetry(), and LoadEventData().

Int_t QwHelicityPattern::fPatternSize

protected

Definition at line 153 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), IsCompletePattern(), LoadEventData(), Print(), ProcessOptions(), and QwHelicityPattern().

QwSubsystemArrayParity QwHelicityPattern::fPositiveHelicitySum

protected

Definition at line 197 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), and ClearEventData().

Bool_t QwHelicityPattern::fPrintBurstSum

protected

Definition at line 177 of file QwHelicityPattern.h.

Referenced by CalculateBurstAverage(), CalculateRunningBurstAverage(), and ProcessOptions().

Bool_t QwHelicityPattern::fPrintRunningSum

protected

Definition at line 187 of file QwHelicityPattern.h.

Referenced by CalculateRunningAverage(), and ProcessOptions().

Int_t QwHelicityPattern::fQuartetNumber

protected

Definition at line 154 of file QwHelicityPattern.h.

Referenced by CalculateAsymmetry(), and QwHelicityPattern().

QwSubsystemArrayParity QwHelicityPattern::fRunningAsymmetry

protected

Definition at line 190 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), ClearRunningSum(), FillDB(), PrintRunningAverage(), and WritePromptSummary().

QwSubsystemArrayParity QwHelicityPattern::fRunningAsymmetry1

protected

Definition at line 191 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), ClearRunningSum(), FillDB(), and PrintRunningAverage().

QwSubsystemArrayParity QwHelicityPattern::fRunningAsymmetry2

protected

Definition at line 192 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), ClearRunningSum(), FillDB(), and PrintRunningAverage().

QwSubsystemArrayParity QwHelicityPattern::fRunningBurstAsymmetry

protected

Definition at line 183 of file QwHelicityPattern.h.

Referenced by AccumulateRunningBurstSum(), CalculateRunningBurstAverage(), ClearRunningSum(), and PrintRunningBurstAverage().

QwSubsystemArrayParity QwHelicityPattern::fRunningBurstDifference

protected

Definition at line 182 of file QwHelicityPattern.h.

Referenced by AccumulateRunningBurstSum(), CalculateRunningBurstAverage(), ClearRunningSum(), and PrintRunningBurstAverage().

QwSubsystemArrayParity QwHelicityPattern::fRunningBurstYield

protected

Definition at line 181 of file QwHelicityPattern.h.

Referenced by AccumulateRunningBurstSum(), CalculateRunningBurstAverage(), ClearRunningSum(), and PrintRunningBurstAverage().

QwSubsystemArrayParity QwHelicityPattern::fRunningDifference

protected

Definition at line 189 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), ClearRunningSum(), FillDB(), and PrintRunningAverage().

QwSubsystemArrayParity QwHelicityPattern::fRunningYield

protected

Definition at line 188 of file QwHelicityPattern.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), ClearRunningSum(), FillDB(), PrintRunningAverage(), and WritePromptSummary().

QwSubsystemArrayParity QwHelicityPattern::fYield

protected

Definition at line 167 of file QwHelicityPattern.h.

Referenced by AccumulateBurstSum(), AccumulateRunningSum(), CalculateAsymmetry(), ClearEventData(), ConstructBranch(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), and FillTreeVector().

---

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

• QwHelicityPattern.h
• QwHelicityPattern.cc