Class for blinding data, adapted from G0 blinder class. More...

#include <QwBlinder.h>

Public Types
enum	EQwBlindingStrategy { kDisabled, kAdditive, kMultiplicative, kAdditiveMultiplicative }
	Available blinding strategies. More...

enum	EQwBlinderStatus { kIndeterminate = 0, kNotBlindable, kBlindable, kBlindableFail }
	Status of the blinding process or intermediate steps of the process. More...

Public Member Functions
	QwBlinder (const EQwBlindingStrategy blinding_strategy=kAdditive)
	Default constructor with optional database. More...

virtual	~QwBlinder ()
	Default destructor. More...

void	ProcessOptions (QwOptions &options)
	Update the status with new external information. More...

void	Update (QwParityDB *db)
	Update the status with new external information. More...

void	Update (const QwSubsystemArrayParity &detectors)
	Update the status with new external information. More...

void	Update (const QwEPICSEvent &epics)
	Update the status with new external information. More...

void	ClearEventData ()

void	WriteFinalValuesToDB (QwParityDB *db)

void	PrintFinalValues ()

void	FillDB (QwParityDB *db, TString datatype)
	Write to the database. More...

void	FillErrDB (QwParityDB *db, TString datatype)

void	ModifyThisErrorCode (UInt_t &errorcode) const

void	BlindValue (Double_t &value) const
	Asymmetry blinding. More...

void	UnBlindValue (Double_t &value) const
	Asymmetry unblinding. More...

void	BlindValue (Double_t &value, const Double_t &yield) const
	Difference blinding. More...

void	UnBlindValue (Double_t &value, const Double_t &yield) const
	Difference unblinding. More...

void	Blind (QwSubsystemArrayParity &diff)
	Blind the asymmetry of an array of subsystems. More...

void	UnBlind (QwSubsystemArrayParity &diff)
	Unblind the asymmetry of an array of subsystems. More...

void	Blind (QwSubsystemArrayParity &diff, const QwSubsystemArrayParity &yield)
	Blind the difference of an array of subsystems. More...

void	UnBlind (QwSubsystemArrayParity &diff, const QwSubsystemArrayParity &yield)
	Unblind the difference of an array of subsystems. More...

const Bool_t &	IsBlinderOkay () const

Static Public Member Functions
static void	DefineOptions (QwOptions &options)

Static Public Attributes
static const TString	fStatusName [4]

static const UInt_t	kErrorFlag_BlinderFail = 0x200
	Error flag value. More...

Private Member Functions
void	SetTargetBlindability (EQwBlinderStatus status)

void	SetWienState (EQwWienMode wienmode)

void	SetIHWPPolarity (Int_t ihwppolarity)

EQwBlinderStatus	CheckBlindability ()

	QwBlinder (const QwBlinder &__attribute__((unused)) blinder)
	Private copy constructor. More...

const QwBlinder &	operator= (const QwBlinder &__attribute__((unused)) blinder)
	Private assignment operator. More...

void	InitBlinders (const UInt_t seed_id)
	Vector of test values, after unblinding. More...

void	InitTestValues (const int n)
	Initializes fBlindingFactor from fSeed. More...

Bool_t	CheckTestValues ()

Int_t	UseMD5 (const TString &barestring)
	Recomputes fBlindTestValue to check for memory errors. More...

Int_t	UseStringManip (const TString &barestring)
	Returns an integer from a string using MD5. More...

Int_t	UsePseudorandom (const TString &barestring)

Int_t	ReadSeed (QwParityDB *db, const UInt_t seed_id)
	Reads the seed with specified id from the database object. More...

Int_t	ReadSeed (QwParityDB *db)
	Reads the seed from the database object. More...

void	WriteChecksum (QwParityDB *db)

void	WriteTestValues (QwParityDB *db)
	Writes fSeedID and fBFChecksum to DB for this analysis ID. More...

std::vector< UChar_t >	GenerateDigest (const TString &input) const
	Writes fTestNumber and fBlindTestValue to DB for this analysis ID. More...

Private Attributes
EQwBlinderStatus	fTargetBlindability_firstread
	Indicates the first value recieved of the blindability of the target. More...

EQwBlinderStatus	fTargetBlindability

Bool_t	fTargetPositionForced

EQwWienMode	fWienMode_firstread

EQwWienMode	fWienMode

Int_t	fIHWPPolarity_firstread

Int_t	fIHWPPolarity

Double_t	fBeamCurrentThreshold

Bool_t	fBeamIsPresent

Bool_t	fBlinderIsOkay

const EQwBlindingStrategy	fBlindingStrategy

Double_t	fBlindingOffset
	Blinding strategy. More...

Double_t	fBlindingOffset_Base
	The term to be added to detector asymmetries. More...

Double_t	fBlindingFactor
	The term to be added to detector asymmetries, before polarity correction. More...

UInt_t	fSeedID
	Maximum blinding factor (in % from identity) More...

TString	fSeed
	ID of seed used (seeds.seed_id) More...

std::vector< UChar_t >	fDigest
	Default seed. More...

std::string	fChecksum
	Checksum in raw hex. More...

std::vector< double >	fTestValues
	Checksum in ASCII hex. More...

std::vector< double >	fBlindTestValues
	Vector of test values, original. More...

std::vector< double >	fUnBlindTestValues
	Vector of test values, after blinding. More...

std::vector< Int_t >	fPatternCounters
	Counts the number of events in each failure mode. More...

Static Private Attributes
static const Double_t	kMaximumBlindingAsymmetry = 0.06
	The factor to be mutliplied to detector asymmetries. More...

static const Double_t	kMaximumBlindingFactor = 0.1
	Maximum blinding asymmetry (in ppm) More...

static const TString	kDefaultSeed = "Default seed, should not be used!"
	Seed string (seeds.seed) More...

Detailed Description

Class for blinding data, adapted from G0 blinder class.

Asymmetry blinding scheme:
$Asym_{blinded} = (Asym_{actual} + fBlindingOffset) * fBlindingFactor$
where $fBlindingOffset = F \times sign(\lambda/2)$ , F is an encrypted factor with |F| < 0.06 ppm. This offset fBlindingOffset will be applied on the block and blocksum of the asymmetry.
Difference blinding scheme: For blinding the helicity correlated differences of the detectors, we'd have to do:
$Diff_{blinded} = (Diff_{raw} + Yield_{raw} * fBlindingOffset) * fBlindingFactor)$
where $Asym_{raw} = Diff_{raw} / Yield_{raw}$ is the unblinded asymmetry, and $Asym_{blinded} = Diff_{blinded} / Yield_{blinded}$ the blinded asymmetry. Blinding the differences allows that the difference can be written to the output ROOT files without compromising the blinding.

Definition at line 64 of file QwBlinder.h.

Member Enumeration Documentation

enum QwBlinder::EQwBlinderStatus

Status of the blinding process or intermediate steps of the process.

Enumerator
kIndeterminate
kNotBlindable
kBlindable
kBlindableFail

Definition at line 76 of file QwBlinder.h.

                         {
     kIndeterminate = 0,
     kNotBlindable,
     kBlindable,
     kBlindableFail
   };

enum QwBlinder::EQwBlindingStrategy

Available blinding strategies.

Enumerator
kDisabled
kAdditive
kMultiplicative
kAdditiveMultiplicative

Definition at line 69 of file QwBlinder.h.

                            {
     kDisabled,
     kAdditive,
     kMultiplicative,
     kAdditiveMultiplicative
   };

Constructor & Destructor Documentation

QwBlinder::QwBlinder ( const EQwBlindingStrategy blinding_strategy = kAdditive )

Default constructor with optional database.

Default constructor using optional database connection and blinding strategy

Parameters

blinding_strategy Blinding strategy

Definition at line 68 of file QwBlinder.cc.

References fPatternCounters, fSeed, fSeedID, InitBlinders(), InitTestValues(), kBlinderCount_NumCounters, and kDefaultSeed.

                                                                :
   fTargetBlindability_firstread(kIndeterminate),
   fTargetBlindability(kIndeterminate),
   fTargetPositionForced(kFALSE),
   //
   fWienMode_firstread(kWienIndeterminate),
   fWienMode(kWienIndeterminate),
   fIHWPPolarity_firstread(0),
   fIHWPPolarity(0),
   //
   fBeamCurrentThreshold(1.0),
   fBeamIsPresent(kFALSE),
   fBlindingStrategy(blinding_strategy),
   fBlindingOffset(0.0),
   fBlindingOffset_Base(0.0),
   fBlindingFactor(1.0)
 {
   // Set up the blinder with seed_id 0
   fSeed = kDefaultSeed;
   fSeedID = 0;
 
   InitBlinders(0);
 
   // Calculate set of test values
   InitTestValues(10);
   fPatternCounters.resize(kBlinderCount_NumCounters);
 }

Here is the call graph for this function:

QwBlinder::~QwBlinder ( )

virtual

Default destructor.

Destructor checks the validity of the blinding and unblinding

Definition at line 100 of file QwBlinder.cc.

References PrintFinalValues().

 {
   // Check the blinded values
   PrintFinalValues();
 }

Here is the call graph for this function:

QwBlinder::QwBlinder ( const QwBlinder &__attribute__((unused)) blinder )

inlineprivate

Private copy constructor.

Definition at line 221 of file QwBlinder.h.

221 : fBlindingStrategy(kDisabled) { };

QwBlinder::kDisabled

Definition: QwBlinder.h:70

QwBlinder::fBlindingStrategy

const EQwBlindingStrategy fBlindingStrategy

Definition: QwBlinder.h:223

Member Function Documentation

void QwBlinder::Blind ( QwSubsystemArrayParity & diff )

inline

Blind the asymmetry of an array of subsystems.

Definition at line 174 of file QwBlinder.h.

References QwSubsystemArrayParity::Blind(), CheckBlindability(), and kNotBlindable.

Referenced by QwHelicityPattern::CalculateAsymmetry().

                                               {
       if (CheckBlindability()!=kNotBlindable)
   diff.Blind(this);
     };

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::Blind	(	QwSubsystemArrayParity &	diff,
		const QwSubsystemArrayParity &	yield
	)

inline

Blind the difference of an array of subsystems.

Definition at line 184 of file QwBlinder.h.

References QwSubsystemArrayParity::Blind(), CheckBlindability(), and kNotBlindable.

                                                                                    {
       if (CheckBlindability()!=kNotBlindable)
   diff.Blind(this, yield);
     };

Here is the call graph for this function:

void QwBlinder::BlindValue ( Double_t & value ) const

inline

Asymmetry blinding.

Definition at line 124 of file QwBlinder.h.

References fBlindingFactor, fBlindingOffset, fBlindingStrategy, kAdditive, kAdditiveMultiplicative, and kMultiplicative.

Referenced by QwVQWK_Channel::Blind(), and InitTestValues().

                                             {
       switch (fBlindingStrategy) {
         case kAdditive:
           value += fBlindingOffset; break;
         case kMultiplicative:
           value *= fBlindingFactor; break;
         case kAdditiveMultiplicative:
           value = (value + fBlindingOffset) * fBlindingFactor; break;
         default: break;
       }
     };

Here is the caller graph for this function:

void QwBlinder::BlindValue	(	Double_t &	value,
		const Double_t &	yield
	)		const

inline

Difference blinding.

Definition at line 149 of file QwBlinder.h.

References fBlindingFactor, fBlindingOffset, fBlindingStrategy, kAdditive, kAdditiveMultiplicative, and kMultiplicative.

                                                                    {
       switch (fBlindingStrategy) {
       case kAdditive:
   value += yield * fBlindingOffset; break;
       case kMultiplicative:
   value *= fBlindingFactor; break;
       case kAdditiveMultiplicative:
   value = (value + fBlindingOffset * yield) * fBlindingFactor; break;
       default: break;
       }
     };

QwBlinder::EQwBlinderStatus QwBlinder::CheckBlindability ( )

private

Definition at line 1079 of file QwBlinder.cc.

References QwLog::endl(), fBeamIsPresent, fBlinderIsOkay, fBlindingStrategy, fIHWPPolarity, fIHWPPolarity_firstread, fPatternCounters, fTargetBlindability, fTargetBlindability_firstread, fTargetPositionForced, fWienMode, fWienMode_firstread, kBlindable, kBlindableFail, kBlinderCount_Blindable, kBlinderCount_ChangedIHWP, kBlinderCount_ChangedTarget, kBlinderCount_ChangedWien, kBlinderCount_Disabled, kBlinderCount_NoBeam, kBlinderCount_NonBlindable, kBlinderCount_OtherFailure, kBlinderCount_Transverse, kBlinderCount_UndefinedIHWP, kBlinderCount_UndefinedWien, kBlinderCount_UnknownTarget, kDisabled, kIndeterminate, kNotBlindable, kWienHorizTrans, kWienIndeterminate, kWienVertTrans, QwDebug, and QwError.

Referenced by Blind().

 {
   EQwBlinderStatus status = QwBlinder::kBlindableFail;
   if (fBlindingStrategy == kDisabled) {
     status = QwBlinder::kNotBlindable;
     fPatternCounters.at(kBlinderCount_Disabled)++;
   } else if (fTargetBlindability == QwBlinder::kIndeterminate) {
     QwDebug  << "QwBlinder::CheckBlindability:  The target blindability is not determined.  "
        << "Fail this pattern." << QwLog::endl;
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_UnknownTarget)++;    
   } else if (fTargetBlindability!=fTargetBlindability_firstread
        && !fTargetPositionForced) {
     QwDebug << "QwBlinder::CheckBlindability:  The target blindability has changed.  "
       << "Fail this pattern." << QwLog::endl;
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_ChangedTarget)++;
   } else if (fTargetBlindability==kNotBlindable) {
     //  This isn't a blindable target, so don't do anything.
     status = QwBlinder::kNotBlindable;
     fPatternCounters.at(kBlinderCount_NonBlindable)++;
   } else if (fTargetBlindability==kBlindable &&
        fWienMode != fWienMode_firstread) {
     //  Wien status changed.  Fail
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_ChangedWien)++;
   } else if (fTargetBlindability==kBlindable &&
        fIHWPPolarity != fIHWPPolarity_firstread) {
     //  IHWP status changed.  Fail
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_ChangedIHWP)++;
   } else if (fTargetBlindability==kBlindable &&
        fWienMode == kWienIndeterminate) {
     //  Wien status isn't determined.  Fail.
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_UndefinedWien)++;
   } else if (fTargetBlindability==kBlindable &&
        fIHWPPolarity==0) {
     //  IHWP status isn't determined.  Fail.
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_UndefinedIHWP)++;
   } else if (fTargetBlindability==kBlindable &&
       (fWienMode == kWienVertTrans || fWienMode == kWienHorizTrans)) {
     //  We don't have longitudinal beam, so don't blind.
     status = QwBlinder::kNotBlindable;
     fPatternCounters.at(kBlinderCount_Transverse)++;
   } else if (fTargetBlindability==kBlindable 
        && fBeamIsPresent) {
     //  This is a blindable target and the beam is sufficent.
     status = QwBlinder::kBlindable;
     fPatternCounters.at(kBlinderCount_Blindable)++;
   } else if (fTargetBlindability==kBlindable 
        && (! fBeamIsPresent) ) {
     //  This is a blindable target but there is insufficent beam present
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_NoBeam)++;
   } else {
     QwError << "QwBlinder::CheckBlindability:  The pattern blindability is unclear.  "
        << "Fail this pattern." << QwLog::endl;
     status = QwBlinder::kBlindableFail;
     fPatternCounters.at(kBlinderCount_OtherFailure)++;
   }
   //
   fBlinderIsOkay = (status != QwBlinder::kBlindableFail);
 
   return status;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Bool_t QwBlinder::CheckTestValues ( )

private

Initializes the test values: fTestNumber, fTestValue, fBlindTestValue, if fBlindingFactor is set.

This routines checks to see if the stored fBlindTestValues[i] match a recomputed blinded test value. The values are cast into floats, and their difference must be less than a change

of the least-significant-bit of fBlindTestValues[i].

First test: compare a blinded value with a second computation

Second test: compare the unblinded value with the original value

Definition at line 807 of file QwBlinder.cc.

References QwLog::endl(), fBlindingOffset, fBlindingOffset_Base, fBlindTestValues, fTestValues, fUnBlindTestValues, QwError, and UnBlindValue().

Referenced by FillDB(), PrintFinalValues(), and WriteFinalValuesToDB().

 {
   Bool_t status = kTRUE;
 
   Double_t tmp_offset = fBlindingOffset;
   fBlindingOffset = fBlindingOffset_Base;
   double epsilon = std::numeric_limits<double>::epsilon();
   for (size_t i = 0; i < fTestValues.size(); i++) {
 
     /// First test: compare a blinded value with a second computation
     double checkval = fBlindTestValues[i];
     UnBlindValue(checkval);
 
     double test1 = fTestValues[i];
     double test2 = checkval;
     if ((test1 - test2) <= -epsilon || (test1 - test2) >= epsilon) {
       QwError << "QwBlinder::CheckTestValues():  Unblinded test value "
               << i
               << " does not agree with original test value, "
               << "with a difference of "
               << (test1 - test2) 
         << " (epsilon==" << epsilon << ")"
         << "." << QwLog::endl;
       status = kFALSE;
     }
 
     /// Second test: compare the unblinded value with the original value
     test1 = fUnBlindTestValues[i];
     test2 = fTestValues[i];
     if ((test1 - test2) <= -epsilon || (test1 - test2) >= epsilon) {
       QwError << "QwBlinder::CheckTestValues():  Unblinded test value "
               << i
               << " does not agree with original test value, "
               << "with a difference of "
               << (test1 - test2) << "." << QwLog::endl;
       status = kFALSE;
     }
   }
   fBlindingOffset = tmp_offset;
   return status;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::ClearEventData ( )

inline

Definition at line 105 of file QwBlinder.h.

References fBeamIsPresent.

Referenced by QwHelicityPattern::ClearEventData().

                          {
       fBeamIsPresent = kTRUE;
     };

Here is the caller graph for this function:

void QwBlinder::DefineOptions ( QwOptions & options )

static

Definition at line 54 of file QwBlinder.cc.

References QwOptions::AddOptions(), and default_bool_value.

Referenced by QwHelicityPattern::DefineOptions().

                                                {
   options.AddOptions("Blinder")("blinder.force-target-lh2", po::value<bool>()->default_bool_value(false),
            "Forces the blinder to interpret the target position as LH2");
   options.AddOptions("Blinder")("blinder.force-target-out", po::value<bool>()->default_bool_value(false),
            "Forces the blinder to interpret the target position as target-out");
   options.AddOptions("Blinder")("blinder.beam-current-threshold", po::value<double>()->default_value(1.0),
            "Beam current in microamps below which data will not be blinded");
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::FillDB	(	QwParityDB *	db,
		TString	datatype
	)

Write to the database.

Write the blinding parameters to the database

Parameters

db	Database connection
datatype	Datatype

For each analyzed run the database contains a digest of the blinding parameters and a number of blinded test entries.

Definition at line 937 of file QwBlinder.cc.

References CheckTestValues(), QwDatabase::Connect(), QwDatabase::Disconnect(), QwLog::endl(), fBlindTestValues, fChecksum, fSeedID, fTestValues, QwParityDB::GetAnalysisID(), QwDatabase::Query(), QwDebug, QwError, and QwMessage.

Referenced by QwHelicityPattern::FillDB().

 {
   QwDebug << " --------------------------------------------------------------- " << QwLog::endl;
   QwDebug << "                         QwBlinder::FillDB                       " << QwLog::endl;
   QwDebug << " --------------------------------------------------------------- " << QwLog::endl;
 
   // Get the analysis ID
   UInt_t analysis_id = db->GetAnalysisID();
 
   // Fill the rows of the QwParitySSQLS::bf_test table
   if (! CheckTestValues()) {
     QwError << "QwBlinder::FillDB():  "
             << "Blinded test values have changed; "
       << "may be a problem in the analysis!!!"
             << QwLog::endl;
   }
 
   QwParitySSQLS::bf_test bf_test_row(0);
   std::vector<QwParitySSQLS::bf_test> bf_test_list;
   for (size_t i = 0; i < fTestValues.size(); i++) {
     bf_test_row.bf_test_id = 0;
     bf_test_row.analysis_id = analysis_id;
     bf_test_row.test_number = i;
     bf_test_row.test_value  = fBlindTestValues[i];
     bf_test_list.push_back(bf_test_row);
   }
 
 
   // Connect to the database
   db->Connect();
 
   // Modify the seed_id and bf_checksum in the analysis table
   try {
     // Get the rows of the QwParitySSQLS::analysis table
     mysqlpp::Query query = db->Query();
     query << "select * from analysis where analysis_id = " 
     << analysis_id;
     std::vector<QwParitySSQLS::analysis> analysis_res;
     QwDebug << "Query: " << query.str() << QwLog::endl;
     query.storein(analysis_res);
     if (analysis_res.size() == 1) {
       QwParitySSQLS::analysis analysis_row_new  = analysis_res[0];
       // Modify the seed_id and bf_checksum
       analysis_row_new.seed_id = fSeedID;
       analysis_row_new.bf_checksum = fChecksum;
       // Update the analysis table
       query.update(analysis_res[0], analysis_row_new);
       QwDebug << "Query: " << query.str() << QwLog::endl;
       query.execute();
     } else {
       QwError << "Unable to update analysis table with blinder information. "
         << QwLog::endl;
     }
   } catch (const mysqlpp::Exception& err) {
     QwError << err.what() << QwLog::endl;
   }
 
   // Add the bf_test rows
   try {
     if (bf_test_list.size()) {
       mysqlpp::Query query = db->Query();
       query.insert(bf_test_list.begin(), bf_test_list.end());
       QwDebug << "Query: " << query.str() << QwLog::endl;
       query.execute();
     } else {
       QwMessage << "QwBlinder::FillDB(): No bf_test entries to write." 
     << QwLog::endl;
     }
   } catch (const mysqlpp::Exception& err) {
     QwError << err.what() << QwLog::endl;
   }
 
   // Disconnect from database
   db->Disconnect();
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::FillErrDB	(	QwParityDB *	db,
		TString	datatype
	)

Definition at line 1014 of file QwBlinder.cc.

References QwErrDBInterface::AddThisEntryToList(), QwDatabase::Connect(), QwDatabase::Disconnect(), QwLog::endl(), fPatternCounters, QwParityDB::GetAnalysisID(), kBlinderCount_NumCounters, QwDatabase::Query(), QwDebug, QwErrDBInterface::SetAnalysisID(), QwErrDBInterface::SetDeviceID(), QwErrDBInterface::SetErrorCodeId(), and QwErrDBInterface::SetN().

Referenced by QwHelicityPattern::FillErrDB().

 {
   QwDebug << " --------------------------------------------------------------- " << QwLog::endl;
   QwDebug << "                     QwBlinder::FillErrDB                        " << QwLog::endl;
   QwDebug << " --------------------------------------------------------------- " << QwLog::endl;
   QwErrDBInterface row;
   std::vector<QwParitySSQLS::general_errors> entrylist;
 
   UInt_t analysis_id = db->GetAnalysisID();
 
   row.SetAnalysisID(analysis_id);
   row.SetDeviceID(kMaxUInt);
   for (size_t index=0; index<kBlinderCount_NumCounters; index++){
     row.SetErrorCodeId(index+20); 
     row.SetN(fPatternCounters.at(index));
     row.AddThisEntryToList( entrylist );
   }
 
   db->Connect();
   // Check the entrylist size, if it isn't zero, start to query..
   if( entrylist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(entrylist.begin(), entrylist.end());
     query.execute();
   }
   db->Disconnect();
 
   return;
 };

Here is the call graph for this function:

Here is the caller graph for this function:

std::vector< UChar_t > QwBlinder::GenerateDigest ( const TString & input ) const

private

Writes fTestNumber and fBlindTestValue to DB for this analysis ID.

Generate an MD5 digest of the blinding parameters

Parameters

input Input string

Returns: MD5 digest of the input string

Definition at line 855 of file QwBlinder.cc.

Referenced by InitBlinders(), and UseMD5().

 {
   // Initialize MD5 checksum array
   const UInt_t length = 16;
   UChar_t value[length];
   for (UInt_t i = 0; i < length; i++)
     value[i] = 0;
 
   // Calculate MD5 checksum from input and store in md5_value
   TMD5 md5;
   md5.Update((UChar_t*) input.Data(), input.Length());
   md5.Final(value);
 
   // Copy the MD5 checksum in an STL vector
   std::vector<UChar_t> output;
   for (UInt_t i = 0; i < length; i++)
     output.push_back(value[i]);
 
   return output;
 }

Here is the caller graph for this function:

void QwBlinder::InitBlinders ( const UInt_t seed_id )

private

Vector of test values, after unblinding.

Initialize the blinder parameters

The blinding constants are determined in two steps.

First, the blinding asymmetry (offset) is determined. It is generated from a signed number between +/- 0.244948974 that is squared to get a number between +/- 0.06 ppm.

Secondly, the multiplicative blinding factor is determined. This number is generated from the blinding asymmetry between, say, 0.9 and 1.1 by an oscillating but uniformly distributed sawtooth function.

TODO: This section of InitBlinders doesn't calculate a reasonable fBlindingFactor but we don't use it for anything.

Definition at line 441 of file QwBlinder.cc.

References QwLog::endl(), fBlindingFactor, fBlindingOffset, fBlindingOffset_Base, fBlindingStrategy, fChecksum, fDigest, fSeed, fSeedID, GenerateDigest(), kDefaultSeed, kDisabled, kMaximumBlindingAsymmetry, QwMessage, QwWarning, and UseMD5().

Referenced by QwBlinder(), and Update().

 {
   // If the blinding strategy is disabled
   if (fBlindingStrategy == kDisabled) {
 
       fSeed = kDefaultSeed;
       fSeedID = 0;
       fBlindingFactor = 1.0;
       fBlindingOffset = 0.0;
       fBlindingOffset_Base = 0.0;
       QwWarning << "Blinding parameters have been disabled!"<< QwLog::endl;
 
   // Else blinding is enabled
   } else {
 
     Int_t finalseed = UseMD5(fSeed);
 
     Double_t newtempout;
     if ((finalseed & 0x80000000) == 0x80000000) {
       newtempout = -1.0 * (finalseed & 0x7FFFFFFF);
     } else {
       newtempout =  1.0 * (finalseed & 0x7FFFFFFF);
     }
 
 
     /// The blinding constants are determined in two steps.
     ///
     /// First, the blinding asymmetry (offset) is determined.  It is
     /// generated from a signed number between +/- 0.244948974 that
     /// is squared to get a number between +/- 0.06 ppm.
     static Double_t maximum_asymmetry_sqrt = sqrt(kMaximumBlindingAsymmetry);
     Double_t tmp1 = maximum_asymmetry_sqrt * (newtempout / Int_t(0x7FFFFFFF));
     fBlindingOffset = tmp1 * fabs(tmp1) * 0.000001;
 
     //  Do another little calulation to round off the blinding asymmetry
     Double_t tmp2;
     tmp1 = fBlindingOffset * 4;    // Exactly shifts by two binary places
     tmp2 = tmp1 + fBlindingOffset; // Rounds 5*fBlindingOffset
     fBlindingOffset = tmp2 - tmp1; // fBlindingOffset has been rounded.
     
     //  Set the base blinding offset.
     fBlindingOffset_Base = fBlindingOffset;
 
     /// Secondly, the multiplicative blinding factor is determined.  This
     /// number is generated from the blinding asymmetry between, say, 0.9 and 1.1
     /// by an oscillating but uniformly distributed sawtooth function.
     fBlindingFactor = 1.0;
     if (kMaximumBlindingAsymmetry > 0.0) {
       /// TODO:  This section of InitBlinders doesn't calculate a reasonable fBlindingFactor but we don't use it for anything.
       fBlindingFactor  = 1.0 + fmod(30.0 * fBlindingOffset, kMaximumBlindingAsymmetry);
       fBlindingFactor /= (kMaximumBlindingAsymmetry > 0.0 ? kMaximumBlindingAsymmetry : 1.0);
     }
 
     QwMessage << "Blinding parameters have been calculated."<< QwLog::endl;
 
   }
 
   // Generate checksum
   TString hex_string;
   hex_string.Form("%.16llx%.16llx", *(ULong64_t*)(&fBlindingFactor), *(ULong64_t*)(&fBlindingOffset));
   fDigest = GenerateDigest(hex_string);
   fChecksum = "";
   for (size_t i = 0; i < fDigest.size(); i++)
     fChecksum += string(Form("%.2x",fDigest[i]));
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::InitTestValues ( const int n )

private

Initializes fBlindingFactor from fSeed.

Generate a set of test values of similar size as measured asymmetries

Parameters

n	Number of test values

Definition at line 525 of file QwBlinder.cc.

References BlindValue(), QwLog::endl(), fBlindingOffset, fBlindingOffset_Base, fBlindTestValues, fSeed, fTestValues, fUnBlindTestValues, QwMessage, UnBlindValue(), and UsePseudorandom().

Referenced by QwBlinder(), and Update().

 {
   // Use the stored seed to get a pseudorandom number
   Int_t finalseed = UsePseudorandom(fSeed);
 
   fTestValues.clear();
   fBlindTestValues.clear();
   fUnBlindTestValues.clear();
 
   Double_t tmp_offset = fBlindingOffset;
   fBlindingOffset = fBlindingOffset_Base;
   // For each test case
   for (int i = 0; i < n; i++) {
 
     // Generate a pseudorandom number
     for (Int_t j = 0; j < 16; j++) {
       finalseed &= 0x7FFFFFFF;
       if ((finalseed & 0x800000) == 0x800000) {
         finalseed = ((finalseed ^ 0x00000d) << 1) | 0x1;
       } else {
         finalseed <<= 1;
       }
     }
 
     // Mask out the low digits of the finalseed, multiply by two,
     // divide by the mask value, subtract from 1, and divide result by
     // 1.0e6 to get a range of about -1000 to +1000 ppb.
     Int_t mask = 0xFFFFFF;
     Double_t tempval = (1.0 - 2.0*(finalseed&mask)/mask) / (1.0e6);
 
     // Store the test values
     fTestValues.push_back(tempval);
     BlindValue(tempval);
     fBlindTestValues.push_back(tempval);
     UnBlindValue(tempval);
     fUnBlindTestValues.push_back(tempval);
   }
   fBlindingOffset = tmp_offset;
   QwMessage << "QwBlinder::InitTestValues(): A total of " << fTestValues.size()
             << " test values have been calculated successfully." << QwLog::endl;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

const Bool_t& QwBlinder::IsBlinderOkay ( ) const

inline

Definition at line 193 of file QwBlinder.h.

References fBlinderIsOkay.

Referenced by QwVQWK_Channel::Blind().

193 {return fBlinderIsOkay;};

QwBlinder::fBlinderIsOkay

Bool_t fBlinderIsOkay

Definition: QwBlinder.h:213

Here is the caller graph for this function:

void QwBlinder::ModifyThisErrorCode ( UInt_t & errorcode ) const

inline

Modifies the device error code variable passed to it, if the blinder is not okay.

Definition at line 119 of file QwBlinder.h.

References kErrorFlag_BlinderFail.

Referenced by QwVQWK_Channel::Blind().

                                                       {
       errorcode |= kErrorFlag_BlinderFail;
     };

Here is the caller graph for this function:

const QwBlinder& QwBlinder::operator= ( const QwBlinder &__attribute__((unused)) blinder )

inlineprivate

Private assignment operator.

Definition at line 223 of file QwBlinder.h.

223 { return *this; };

void QwBlinder::PrintFinalValues ( )

Print a summary of the blinding/unblinding test

Definition at line 880 of file QwBlinder.cc.

References CheckTestValues(), QwLog::endl(), fBlindTestValues, fChecksum, fPatternCounters, fSeedID, fTestValues, fUnBlindTestValues, kBlinderCount_Blindable, kBlinderCount_ChangedIHWP, kBlinderCount_ChangedTarget, kBlinderCount_ChangedWien, kBlinderCount_Disabled, kBlinderCount_NoBeam, kBlinderCount_NonBlindable, kBlinderCount_OtherFailure, kBlinderCount_Transverse, kBlinderCount_UndefinedIHWP, kBlinderCount_UndefinedWien, kBlinderCount_UnknownTarget, and QwMessage.

Referenced by ~QwBlinder().

 {
   QwMessage << "QwBlinder::PrintFinalValues():  Begin summary"    << QwLog::endl;
   QwMessage << "================================================" << QwLog::endl;
   QwMessage << "Blinder Passed Patterns"  << QwLog::endl;
   QwMessage << "\tPatterns with blinding disabled: " << fPatternCounters.at(kBlinderCount_Disabled) << QwLog::endl;
   QwMessage << "\tPatterns on a non-blindable target: " << fPatternCounters.at(kBlinderCount_NonBlindable) << QwLog::endl;
   QwMessage << "\tPatterns with transverse beam: " << fPatternCounters.at(kBlinderCount_Transverse) << QwLog::endl;
   QwMessage << "\tPatterns on blindable target with beam present: " << fPatternCounters.at(kBlinderCount_Blindable) << QwLog::endl;
   QwMessage << "Blinder Failed Patterns"  << QwLog::endl;
   QwMessage << "\tPatterns with unknown target position: " << fPatternCounters.at(kBlinderCount_UnknownTarget) << QwLog::endl;
   QwMessage << "\tPatterns with changed target position: " << fPatternCounters.at(kBlinderCount_ChangedTarget) << QwLog::endl;
   QwMessage << "\tPatterns with an undefined Wien setting: " << fPatternCounters.at(kBlinderCount_UndefinedWien) << QwLog::endl;
   QwMessage << "\tPatterns with a changed Wien setting: " << fPatternCounters.at(kBlinderCount_ChangedWien) << QwLog::endl;
   QwMessage << "\tPatterns with an undefined IHWP setting: " << fPatternCounters.at(kBlinderCount_UndefinedIHWP) << QwLog::endl;
   QwMessage << "\tPatterns with a changed IHWP setting: " << fPatternCounters.at(kBlinderCount_ChangedIHWP) << QwLog::endl;
   QwMessage << "\tPatterns on blindable target with no beam: " << fPatternCounters.at(kBlinderCount_NoBeam) << QwLog::endl;
   QwMessage << "\tPatterns with other blinding failure: " << fPatternCounters.at(kBlinderCount_OtherFailure) << QwLog::endl;
   QwMessage << "================================================" << QwLog::endl;
   QwMessage << "The blinding parameters checksum for seed ID "
             << fSeedID << " is:" << QwLog::endl;
   QwMessage << fChecksum << QwLog::endl;
   QwMessage << "================================================" << QwLog::endl;
   CheckTestValues();
   double epsilon = std::numeric_limits<double>::epsilon();
   TString diff;
   QwMessage << "The test results are:" << QwLog::endl;
   QwMessage << std::setw(8)  << "Index"
             << std::setw(16) << "Original value"
             << std::setw(16) << "Blinded value"
             << std::setw(22) << "Orig.-Unblind value"
             << QwLog::endl;
   for (size_t i = 0; i < fTestValues.size(); i++) {
     if ((fTestValues[i]-fUnBlindTestValues[i]) < -epsilon 
   || (fTestValues[i]-fUnBlindTestValues[i]) > epsilon )
       diff = Form("% 9g ppb", fTestValues[i]-fUnBlindTestValues[i]*1e9);
     else
       diff = "epsilon";
     QwMessage << std::setw(8)  << i
               << std::setw(16) << Form(" [CENSORED]")
         << std::setw(16) << Form("% 9.3f ppb",fBlindTestValues[i]*1e9)
               << std::setw(22) << diff
               << QwLog::endl;
   }
   QwMessage << "================================================" << QwLog::endl;
   QwMessage << "QwBlinder::PrintFinalValues():  End of summary"   << QwLog::endl;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::ProcessOptions ( QwOptions & options )

Update the status with new external information.

Update the blinder status with new external information

Parameters

options Qweak option handler

Definition at line 112 of file QwBlinder.cc.

References QwLog::endl(), fBeamCurrentThreshold, fTargetPositionForced, QwOptions::GetValue(), QwOptions::HasValue(), kBlindable, kNotBlindable, QwError, and SetTargetBlindability().

Referenced by QwHelicityPattern::ProcessOptions().

 {
   if (options.GetValue<bool>("blinder.force-target-out")
       && options.GetValue<bool>("blinder.force-target-lh2")){
     QwError << "QwBlinder::Update:  Both blinder.force-target-lh2 and blinder.force-target-out are set.  "
       << "Only one can be in force at one time.  Exit and choose one option."
       << QwLog::endl;
     exit(10);
   } else if (options.GetValue<bool>("blinder.force-target-lh2")){
     fTargetPositionForced = kTRUE;
     SetTargetBlindability(QwBlinder::kBlindable);
   } else if (options.GetValue<bool>("blinder.force-target-out")){
     fTargetPositionForced = kTRUE;
     SetTargetBlindability(QwBlinder::kNotBlindable);
   }
   if (options.HasValue("blinder.beam-current-threshold")){
     fBeamCurrentThreshold = options.GetValue<double>("blinder.beam-current-threshold");
   }
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwBlinder::ReadSeed	(	QwParityDB *	db,
		const UInt_t	seed_id
	)

private

Reads the seed with specified id from the database object.

Returns an integer from a string using a version of the helicity bit pseudorandom algorithm.

--------------------------------------------------------—

Function to read the seed in from the database.

Parameters: seed_id = ID number of seed to use (0 = most recent seed)

Return: Int_t

Definition at line 360 of file QwBlinder.cc.

References QwDatabase::Connect(), QwDatabase::Disconnect(), QwLog::endl(), fSeed, fSeedID, kDefaultSeed, QwDatabase::Query(), QwError, and QwWarning.

Referenced by Update().

 {
   // Return unchanged if no database specified
   if (! db) {
     QwWarning << "QwBlinder::ReadSeed(): No database specified" << QwLog::endl;
     fSeedID = 0;
     fSeed   = "Default seed, No database specified";
     return 0;
   }
 
   // Try to connect to the database
   if (db->Connect()) {
 
     // Build up query
     string s_sql = "SELECT * FROM seeds ";
     if (fSeedID > 0) {
       // Use specified seed
       Char_t s_seed_id[10];
       sprintf(s_seed_id, "%d", seed_id);
       s_sql += "WHERE seed_id = ";
       s_sql += string(s_seed_id);
     } else {
       // Use most recent seed
       s_sql += "ORDER BY seed_id DESC LIMIT 1";
     }
 
     // Send query
     mysqlpp::Query query = db->Query();
     query << s_sql;
     std::vector<QwParitySSQLS::seeds> res;
     query.storein(res);
 
     // Store seed_id and seed value in fSeedID and fSeed (want to store actual seed_id in those
     // cases where only the most recent was requested (fSeedID = 0))
     // There should be one and only one seed_id for each seed.
     if (res.size() == 1) {
       // Store seed_id in fSeedID (want to store actual seed_id in those
       // cases where only the most recent was requested (fSeedID = 0))
       fSeedID = res[0].seed_id;
 
       // Store seed value
       if (!res[0].seed.is_null) {
         fSeed = res[0].seed.data;
       } else {
         QwError << "QwBlinder::ReadSeed(): Seed value came back NULL from the database." << QwLog::endl;
         fSeedID = 0;
         fSeed = kDefaultSeed;
       }
 
       std::cout << "QwBlinder::ReadSeed():  Successfully read "
       << Form("the fSeed with ID %d from the database.", fSeedID)
       << std::endl;
 
     } else {
       // Error Condition.
       // There should be one and only one seed_id for each seed.
       fSeedID = 0;
       fSeed   = Form("ERROR:  There should be one and only one seed_id for each seed, but this had %zu.",
                      res.size());
       std::cerr << "QwBlinder::ReadSeed(): "<<fSeed<<std::endl;
     }
 
     // Disconnect from database
     db->Disconnect();
 
   } else {
 
     //  We were unable to open the connection.
     fSeedID = 0;
     fSeed   = "ERROR:  Unable to open the connection to the database.";
     QwError << "QwBlinder::ReadSeed(): Unable to open connection to database." << QwLog::endl;
   }
 
   return fSeedID;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwBlinder::ReadSeed ( QwParityDB * db )

private

Reads the seed from the database object.

--------------------------------------------------------—

Function to read the seed in from the database.

Parameters:

Return: Int_t

Definition at line 273 of file QwBlinder.cc.

References QwDatabase::Connect(), QwDatabase::Disconnect(), QwLog::endl(), fSeed, fSeedID, QwParityDB::GetRunNumber(), kDefaultSeed, QwDatabase::Query(), QwError, and QwWarning.

 {
   // Return unchanged if no database specified
   if (! db) {
     QwWarning << "QwBlinder::ReadSeed(): No database specified" << QwLog::endl;
     fSeedID = 0;
     fSeed   = "Default seed, No database specified";
     return 0;
   }
 
   // Try to connect to the database
   if (db->Connect()) {
 
     QwError << "QwBlinder::ReadSeed  db->GetRunNumber() returns "
       <<  db->GetRunNumber() << QwLog::endl;
     // Build up query
     string s_sql = "SELECT seed_id, seed FROM seeds, run as rf, run as rl WHERE seeds.first_run_id = rf.run_id AND seeds.last_run_id = rl.run_id ";
     s_sql += "AND rf.run_number <= ";
     s_sql += Form("%d",db->GetRunNumber());
     s_sql += " AND rl.run_number >= ";
     s_sql += Form("%d",db->GetRunNumber());
     s_sql += " AND seed_id>2";
 
     QwError << "QwBlinder::ReadSeed s_sql contains \'"
       << s_sql  << "\'" << QwLog::endl;
 
     // Send query
     mysqlpp::Query query = db->Query();
     query << s_sql;
     mysqlpp::StoreQueryResult res = query.store();
 
     // Store seed_id and seed value in fSeedID and fSeed (want to store actual seed_id in those
     // cases where only the most recent was requested (fSeedID = 0))
     // There should be one and only one seed_id for each seed.
     if (res.size() == 1) {
       // Store seed_id in fSeedID (want to store actual seed_id in those
       // cases where only the most recent was requested (fSeedID = 0))
       fSeedID = res[0]["seed_id"];
 
       // Store seed value
       if (!res[0]["seed"].is_null()) {
         fSeed = res[0]["seed"].data();
       } else {
         QwError << "QwBlinder::ReadSeed(): Seed value came back NULL from the database." << QwLog::endl;
         fSeedID = 0;
         fSeed = kDefaultSeed;
       }
 
       std::cout << "QwBlinder::ReadSeed():  Successfully read "
       << Form("the fSeed with ID %d from the database.", fSeedID)
       << std::endl;
 
     } else {
       // Error Condition.
       // There should be one and only one seed_id for each seed.
       fSeedID = 0;
       fSeed   = Form("ERROR:  There should be one and only one seed_id for each seed, but this had %zu.",
                      res.size());
       std::cerr << "QwBlinder::ReadSeed(): "<<fSeed<<std::endl;
     }
 
     // Disconnect from database
     db->Disconnect();
 
   } else {
 
     //  We were unable to open the connection.
     fSeedID = 0;
     fSeed   = "ERROR:  Unable to open the connection to the database.";
     QwError << "QwBlinder::ReadSeed(): Unable to open connection to database." << QwLog::endl;
   }
 
   return fSeedID;
 }

Here is the call graph for this function:

void QwBlinder::SetIHWPPolarity ( Int_t ihwppolarity )

private

Definition at line 1068 of file QwBlinder.cc.

References QwLog::endl(), fIHWPPolarity, fIHWPPolarity_firstread, and QwMessage.

Referenced by Update().

 {
   fIHWPPolarity = ihwppolarity;
   if (fIHWPPolarity_firstread == 0 && fIHWPPolarity != 0){
     fIHWPPolarity_firstread = fIHWPPolarity;
     QwMessage << "QwBlinder:  First set IHWP state to " 
         << fIHWPPolarity << QwLog::endl;
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::SetTargetBlindability ( QwBlinder::EQwBlinderStatus status )

private

Definition at line 1046 of file QwBlinder.cc.

References QwLog::endl(), fStatusName, fTargetBlindability, fTargetBlindability_firstread, kIndeterminate, and QwMessage.

Referenced by ProcessOptions(), and Update().

 {
   fTargetBlindability = status;
   if (fTargetBlindability_firstread == QwBlinder::kIndeterminate
       && fTargetBlindability != QwBlinder::kIndeterminate){
     fTargetBlindability_firstread = fTargetBlindability;
     QwMessage << "QwBlinder:  First set target blindability to " 
         << fStatusName[fTargetBlindability] << QwLog::endl;
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::SetWienState ( EQwWienMode wienmode )

private

Definition at line 1057 of file QwBlinder.cc.

References QwLog::endl(), fWienMode, fWienMode_firstread, kWienIndeterminate, QwMessage, and WienModeName().

Referenced by Update().

 {
   fWienMode = wienmode;
   if (fWienMode_firstread == kWienIndeterminate
       && fWienMode != kWienIndeterminate){
     fWienMode_firstread = fWienMode;
     QwMessage << "QwBlinder:  First set Wien state to " 
         << WienModeName(fWienMode) << QwLog::endl;
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::UnBlind ( QwSubsystemArrayParity & diff )

inline

Unblind the asymmetry of an array of subsystems.

Definition at line 179 of file QwBlinder.h.

References QwSubsystemArrayParity::UnBlind().

                                                 {
       diff.UnBlind(this);
     };

Here is the call graph for this function:

void QwBlinder::UnBlind	(	QwSubsystemArrayParity &	diff,
		const QwSubsystemArrayParity &	yield
	)

inline

Unblind the difference of an array of subsystems.

Definition at line 189 of file QwBlinder.h.

References QwSubsystemArrayParity::UnBlind().

                                                                                      {
       diff.UnBlind(this, yield);
     };

Here is the call graph for this function:

void QwBlinder::UnBlindValue ( Double_t & value ) const

inline

Asymmetry unblinding.

Definition at line 136 of file QwBlinder.h.

References fBlindingFactor, fBlindingOffset, fBlindingStrategy, kAdditive, kAdditiveMultiplicative, and kMultiplicative.

Referenced by CheckTestValues(), and InitTestValues().

                                               {
       switch (fBlindingStrategy) {
         case kAdditive:
           value -= fBlindingOffset; break;
         case kMultiplicative:
           value /= fBlindingFactor; break;
         case kAdditiveMultiplicative:
           value = value / fBlindingFactor - fBlindingOffset; break;
         default: break;
       }
     };

Here is the caller graph for this function:

void QwBlinder::UnBlindValue	(	Double_t &	value,
		const Double_t &	yield
	)		const

inline

Difference unblinding.

Definition at line 161 of file QwBlinder.h.

References fBlindingFactor, fBlindingOffset, fBlindingStrategy, kAdditive, kAdditiveMultiplicative, and kMultiplicative.

                                                                      {
       switch (fBlindingStrategy) {
         case kAdditive:
           value -= yield * fBlindingOffset; break;
         case kMultiplicative:
           value /= fBlindingFactor; break;
         case kAdditiveMultiplicative:
           value = value / fBlindingFactor - yield * fBlindingOffset; break;
         default: break;
       }
     };

void QwBlinder::Update ( QwParityDB * db )

Update the status with new external information.

Update the blinder status with new external information

Parameters

db	Database connection

Definition at line 138 of file QwBlinder.cc.

References QwLog::endl(), fSeed, fSeedID, InitBlinders(), InitTestValues(), kDefaultSeed, QwWarning, and ReadSeed().

Referenced by QwHelicityPattern::UpdateBlinder().

 {
   //  Update the seed ID then tell us if it has changed.
   UInt_t old_seed_id = fSeedID;
   ReadSeed(db);
   // If the blinding seed has changed, re-initialize the blinder
   if (fSeedID != old_seed_id ||
       (fSeedID==0 && fSeed!=kDefaultSeed) ) {
     QwWarning << "Changing blinder seed to " << fSeedID
               << " from " << old_seed_id << "." << QwLog::endl;
     InitBlinders(fSeedID);
     InitTestValues(10);
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::Update ( const QwSubsystemArrayParity & detectors )

Update the status with new external information.

Update the blinder status with new external information

Parameters

detectors Current subsystem array

Definition at line 158 of file QwBlinder.cc.

References fBeamCurrentThreshold, fBeamIsPresent, fBlindingStrategy, fTargetBlindability, QwVQWK_Channel::GetValue(), kBlindable, kDisabled, and QwSubsystemArray::ReturnInternalValue().

 {
   static QwVQWK_Channel q_targ("q_targ");
   if (fBlindingStrategy != kDisabled && fTargetBlindability==kBlindable) {
     // Check for the target blindability flag
     
 
     // Check that the current on target is above acceptable limit
     Bool_t tmp_beam = kFALSE;
     //    if (detectors.ReturnInternalValue(q_targ.GetElementName(), &q_targ)) {
     if (detectors.ReturnInternalValue("q_targ", &q_targ)) {
       if (q_targ.GetValue() > fBeamCurrentThreshold){
   //  std::cerr << "q_targ.GetValue()==" 
   //      << q_targ.GetValue() << std::endl;
   tmp_beam = kTRUE;
       }
     }
     fBeamIsPresent &= tmp_beam;
   }
 }

Here is the call graph for this function:

void QwBlinder::Update ( const QwEPICSEvent & epics )

Update the status with new external information.

Update the blinder status with new external information

Parameters

epics Current EPICS event

Definition at line 184 of file QwBlinder.cc.

References QwEPICSEvent::DetermineIHWPPolarity(), QwEPICSEvent::DetermineWienMode(), QwLog::endl(), fBlindingOffset, fBlindingOffset_Base, fBlindingStrategy, fIHWPPolarity, fTargetBlindability, fTargetPositionForced, fWienMode, QwEPICSEvent::GetDataString(), QwEPICSEvent::GetDataValue(), kBlindable, kDisabled, kIndeterminate, kNotBlindable, kWienBackward, kWienForward, QwDebug, QwWarning, SetIHWPPolarity(), SetTargetBlindability(), and SetWienState().

 {
   // Check for the target information
   // Position:
   //     QWtgt_name == "HYDROGEN-CELL" || QWTGTPOS > 350
   // Temperature:
   //     QWT_miA < 22.0 && QWT_moA < 22.0
   // Pressure:
   //     QW_PT3 in [20,35] && QW_PT4 in [20,35]
   if (fBlindingStrategy != kDisabled && !(fTargetPositionForced) ) {
     TString  position  = epics.GetDataString("QWtgt_name");
     Double_t tgt_pos   = epics.GetDataValue("QWTGTPOS");
     Double_t tgt_temperture  = epics.GetDataValue("QWT_miB");
     Double_t tgt_temperture2 = epics.GetDataValue("QWT_moB");
     Double_t tgt_pressure    = epics.GetDataValue("QW_PT3");
     Double_t tgt_pressure2   = epics.GetDataValue("QW_PT4");
     QwDebug << "Target parameters used by the blinder: "
       << "QWtgt_name=" << position << " "
       << "QWTGTPOS=" << tgt_pos << " "
       << "QWT_miB=" << tgt_temperture << " "
       << "QWT_moB=" << tgt_temperture2 << " "
       << "QW_PT3=" << tgt_pressure << " "
       << "QW_PT4=" << tgt_pressure2 << " "
       << QwLog::endl;
     if (position == "HYDROGEN-CELL"
   && tgt_pos > 350. 
   && (tgt_temperture>18.0 && tgt_temperture<22.0)
   && (tgt_temperture2>18.0 && tgt_temperture2<22.0)
   && (tgt_pressure>20.0 && tgt_pressure < 35.0)
   && (tgt_pressure2>20.0 && tgt_pressure2 < 35.0)){
       SetTargetBlindability(QwBlinder::kBlindable);
     } else if ((position == "HYDROGEN-CELL"
     && tgt_pos > 350.)
          && (tgt_temperture > 22.0)
          && (tgt_temperture2 > 22.0)
          && (tgt_pressure > 35.0)
          && (tgt_pressure2 > 35.0)){
       //  Hydrogen cell is in, but temperature and pressures
       //  are all higher than they should be for LH2.
       SetTargetBlindability(QwBlinder::kNotBlindable);
     } else if ((position != "HYDROGEN-CELL"
     && tgt_pos < 350.)){
       //  Name and position agree that this isn't the hydrogen
       //  cell.
       SetTargetBlindability(QwBlinder::kNotBlindable);
     } else {
       SetTargetBlindability(QwBlinder::kIndeterminate);
       QwWarning << "Target parameters used by the blinder are indeterminate: "
     << "QWtgt_name=" << position << " "
     << "QWTGTPOS=" << tgt_pos << " "
     << "QWT_miB=" << tgt_temperture << " "
     << "QWT_moB=" << tgt_temperture2 << " "
     << "QW_PT3=" << tgt_pressure << " "
     << "QW_PT4=" << tgt_pressure2 << " "
     << QwLog::endl;
 
 
     }
   }
   // Check for the beam polarity information
   //     IGL1I00DI24_24M         Beam Half-wave plate Read(off=out)
   //
   if (fBlindingStrategy != kDisabled &&
       (fTargetBlindability == QwBlinder::kBlindable) ) {
     //  Use the EPICS class functions to determine the
     //  Wien mode and IHWP polarity.
     SetWienState(epics.DetermineWienMode());
     SetIHWPPolarity(epics.DetermineIHWPPolarity());
     
     if (fWienMode == kWienForward){
       fBlindingOffset = fBlindingOffset_Base * fIHWPPolarity;
     } else if (fWienMode == kWienBackward){
       fBlindingOffset = -1 * fBlindingOffset_Base * fIHWPPolarity;
     } else {
       fBlindingOffset = 0.0;
     }
   }
 }

Here is the call graph for this function:

Int_t QwBlinder::UseMD5 ( const TString & barestring )

private

Recomputes fBlindTestValue to check for memory errors.

Use an MD5 checksum to get a number from the seed string

Parameters

barestring Seed string

Returns: Integer number

Definition at line 682 of file QwBlinder.cc.

References GenerateDigest().

Referenced by InitBlinders().

 {
   Int_t temp = 0;
   Int_t tempout = 0;
 
   std::vector<UChar_t> digest = GenerateDigest(barestring);
   for (size_t i = 0; i < digest.size(); i++)
     {
       Int_t j = i%4;
       if (j == 0)
         {
           temp = 0;
         }
       temp |= (digest[i])<<(24-(j*8));
       if ( (j==3) || (i==(digest.size()-1)))
         {
           tempout ^= temp;
         }
     }
 
   if ((tempout & 0x80000000) == 0x80000000) {
       tempout = -1 * (tempout & 0x7FFFFFFF);
   } else {
       tempout = (tempout & 0x7FFFFFFF);
   }
 
   return tempout;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwBlinder::UsePseudorandom ( const TString & barestring )

private

Returns an integer from a string using a character manipulation algorithm

Use pseudo-random number generator to get a number from the seed string

Parameters

barestring Seed string

Returns: Integer number

Definition at line 617 of file QwBlinder.cc.

Referenced by InitTestValues().

 {
   ULong64_t finalseed;
   Int_t bitcount;
   Int_t tempout = 0;
 
   //  This is an attempt to build a portable 64-bit constant
   ULong64_t longmask = (0x7FFFFFFF);
   longmask<<=32;
   longmask|=0xFFFFFFFF;
 
   finalseed = 0;
   bitcount  = 0;
   for (Int_t i=0; i<barestring.Length(); i++)
     {
       if ( ((barestring[i])&0xC0)!=0 && ((barestring[i])&0xC0)!=0xC0)
         {
           finalseed = ((finalseed&longmask)<<1) | (((barestring[i])&0x40)>>6);
           bitcount++;
         }
       if ( ((barestring[i])&0x30)!=0 && ((barestring[i])&0x30)!=0x30)
         {
           finalseed = ((finalseed&longmask)<<1) | (((barestring[i])&0x10)>>4);
           bitcount++;
         }
       if ( ((barestring[i])&0xC)!=0 && ((barestring[i])&0xC)!=0xC)
         {
           finalseed = ((finalseed&longmask)<<1) | (((barestring[i])&0x4)>>2);
           bitcount++;
         }
       if ( ((barestring[i])&0x3)!=0 && ((barestring[i])&0x3)!=0x3)
         {
           finalseed = ((finalseed&longmask)<<1) | ((barestring[i])&0x1);
           bitcount++;
         }
     }
   for (Int_t i=0; i<(192-bitcount); i++)
     {
       if ((finalseed & 0x800000) == 0x800000)
         {
           finalseed = ((finalseed^0x00000d)<<1) | 0x1;
         }
       else
         {
           finalseed<<=1;
         }
     }
   tempout = (finalseed&0xFFFFFFFF) ^ ((finalseed>>32)&0xFFFFFFFF);
   if ((tempout&0x80000000) == 0x80000000)
     {
       tempout = -1 * (tempout&0x7FFFFFFF);
     }
   else
     {
       tempout =  1 * (tempout&0x7FFFFFFF);
     }
   return tempout;
 }

Here is the caller graph for this function:

Int_t QwBlinder::UseStringManip ( const TString & barestring )

private

Returns an integer from a string using MD5.

Use string manipulation to get a number from the seed string

Parameters

barestring Seed string

Returns: Integer number

Definition at line 572 of file QwBlinder.cc.

 {
   std::vector<UInt_t> choppedwords;
   UInt_t tmpword;
   Int_t finalseed = 0;
 
   for (Int_t i = 0; i < barestring.Length(); i++)
     {
       if (i % 4 == 0) tmpword = 0;
       tmpword |= (char(barestring[i]))<<(24-8*(i%4));
       if (i%4 == 3 || i == barestring.Length()-1)
         {
           choppedwords.push_back(tmpword);
           finalseed ^= (tmpword);
         }
     }
   for (Int_t i=0; i<64; i++)
     {
       finalseed &= 0x7FFFFFFF;
       if ((finalseed & 0x800000) == 0x800000)
         {
           finalseed = ((finalseed^0x00000d)<<1) | 0x1;
         }
       else
         {
           finalseed<<=1;
         }
     }
   if ((finalseed&0x80000000) == 0x80000000)
     {
       finalseed = -1 * (finalseed&0x7FFFFFFF);
     }
   else
     {
       finalseed = (finalseed&0x7FFFFFFF);
     }
   return finalseed;
 }

void QwBlinder::WriteChecksum ( QwParityDB * db )

private

--------------------------------------------------------—

Function to write the checksum into the analysis table

Parameters: void

Return: void

Note: This function assumes that the analysis table has already

been filled for the run.

Definition at line 725 of file QwBlinder.cc.

References QwDatabase::Connect(), QwDatabase::Disconnect(), fChecksum, fSeedID, QwParityDB::GetAnalysisID(), and QwDatabase::Query().

Referenced by WriteFinalValuesToDB().

 {
   //----------------------------------------------------------
   // Construct SQL
   //----------------------------------------------------------
   Char_t s_number[20];
   string s_sql = "UPDATE analysis SET seed_id = ";
   sprintf(s_number, "%d", fSeedID);
   s_sql += string(s_number);
   s_sql += ", bf_checksum = ";
   s_sql += "\'" + fChecksum + "\'";
   s_sql += " WHERE analysis_id = ";
   sprintf(s_number, "%d", db->GetAnalysisID());
   s_sql += string(s_number);
 
   //----------------------------------------------------------
   // Execute SQL
   //----------------------------------------------------------
   db->Connect();
   mysqlpp::Query query = db->Query();
   query <<s_sql;
   query.execute();
   db->Disconnect();
 } //End QwBlinder::WriteChecksum

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBlinder::WriteFinalValuesToDB ( QwParityDB * db )

Definition at line 508 of file QwBlinder.cc.

References CheckTestValues(), QwLog::endl(), QwError, WriteChecksum(), and WriteTestValues().

 {
   WriteChecksum(db);
   if (! CheckTestValues()) {
     QwError << "QwBlinder::WriteFinalValuesToDB():  "
             << "Blinded test values have changed; may be a problem in the analysis!!!"
             << QwLog::endl;
   }
   WriteTestValues(db);
 }

Here is the call graph for this function:

void QwBlinder::WriteTestValues ( QwParityDB * db )

private

Writes fSeedID and fBFChecksum to DB for this analysis ID.

---------------------------------------------------------—

Function to write the test values into the database

Parameters: void

Return: void

Definition at line 759 of file QwBlinder.cc.

References QwDatabase::Connect(), QwDatabase::Disconnect(), fBlindTestValues, fTestValues, QwParityDB::GetAnalysisID(), and QwDatabase::Query().

Referenced by WriteFinalValuesToDB().

 {
   //----------------------------------------------------------
   // Construct Initial SQL
   //----------------------------------------------------------
   Char_t s_number[20];
 
   string s_sql_pre = "INSERT INTO bf_test (analysis_id, test_number, test_value) VALUES (";
   // analysis_id
   sprintf(s_number, "%d", db->GetAnalysisID());
   s_sql_pre += string(s_number);
   s_sql_pre += ", ";
 
   //----------------------------------------------------------
   // Construct Individual SQL and Execute
   //----------------------------------------------------------
   // Loop over all test values
   for (size_t i = 0; i < fTestValues.size(); i++)
     {
       string s_sql = s_sql_pre;
 
       // test_number
       sprintf(s_number, "%d", (int) i);
       s_sql += string(s_number);
       s_sql += ", ";
 
       // test_value
       sprintf(s_number, "%9g", fBlindTestValues[i]);
       s_sql += s_number;
       s_sql += ")";
 
       // Execute SQL
       db->Connect();
       mysqlpp::Query query = db->Query();
       query <<s_sql;
       query.execute();
       db->Disconnect();
     }
 }