QwMollerDetector Class Reference

Implementation of the analysis of Moller data (copied from QwComptonElectronDetector.h) More...

#include <QwMollerDetector.h>

Inheritance diagram for QwMollerDetector:

Collaboration diagram for QwMollerDetector:

Public Member Functions
	QwMollerDetector (const TString &name)
	Constructor with name. More...
	QwMollerDetector (const QwMollerDetector &source)
	Copy constructor. More...
virtual	~QwMollerDetector ()
	Virtual destructor. More...
void	ProcessOptions (QwOptions &options)
	Process the command line options. More...
Int_t	LoadChannelMap (TString mapfile)
Int_t	LoadInputParameters (TString pedestalfile)
	Mandatory parameter file definition. More...
Int_t	LoadEventCuts (TString &filename)
Bool_t	SingleEventCuts ()
Int_t	ProcessConfigurationBuffer (const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
Int_t	ProcessConfigurationBuffer (UInt_t ev_type, const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
Int_t	ProcessEvBuffer (UInt_t roc_id, UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
	TODO: The non-event-type-aware ProcessEvBuffer routine should be replaced with the event-type-aware version. More...
Int_t	ProcessEvBuffer (UInt_t ev_type, UInt_t roc_id, UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
void	ClearEventData ()
void	ProcessEvent ()
void	ConstructHistograms (TDirectory *folder, TString &prefix)
	Construct the histograms for this subsystem in a folder with a prefix. More...
void	FillHistograms ()
	Fill the histograms for this subsystem. More...
VQwSubsystem &	operator= (VQwSubsystem *value)
	Assignment Note: Must be called at the beginning of all subsystems routine call to operator=(VQwSubsystem *value) by VQwSubsystem::operator=(value) More...
VQwSubsystem &	operator+= (VQwSubsystem *value)
VQwSubsystem &	operator-= (VQwSubsystem *value)
void	Sum (VQwSubsystem *value1, VQwSubsystem *value2)
void	Difference (VQwSubsystem *value1, VQwSubsystem *value2)
void	Ratio (VQwSubsystem *value1, VQwSubsystem *value2)
void	Scale (Double_t)
void	AccumulateRunningSum (VQwSubsystem *value)
	Update the running sums for devices. More...
void	DeaccumulateRunningSum (VQwSubsystem *value)
	remove one entry from the running sums for devices More...
void	CalculateRunningAverage ()
	Calculate the average for all good events. More...
Int_t	LoadEventCuts (TString filename)
	Load the event cuts file. More...
Bool_t	ApplySingleEventCuts ()
	Apply the single event cuts. More...
void	IncrementErrorCounters ()
	Increment the error counters. More...
void	PrintErrorCounters () const
	Report the number of events failed due to HW and event cut failures. More...
UInt_t	GetEventcutErrorFlag ()
	Return the error flag to the top level routines related to stability checks and ErrorFlag updates. More...
void	UpdateErrorFlag (const VQwSubsystem *ev_error)
	update the error flag in the subsystem level from the top level routines related to stability checks. This will uniquely update the errorflag at each channel based on the error flag in the corresponding channel in the ev_error subsystem More...
void	ConstructBranchAndVector (TTree *, TString &, std::vector< double, std::allocator< double > > &)
void	ConstructBranch (TTree *tree, TString &prefix)
	Construct the branch and tree vector. More...
void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &trim_file)
	Construct the branch and tree vector based on the trim file. More...
void	FillTreeVector (std::vector< Double_t > &values) const
	Fill the tree vector. More...
Bool_t	Compare (VQwSubsystem *source)
void	print ()
void	PrintValue () const
	Print values of all channels. More...
float *	GetRawChannelArray ()
Int_t	GetChannelIndex (TString channelName, UInt_t module_number)
float	GetDataForChannelInModule (Int_t module_number, Int_t channel_index)
float	GetDataForChannelInModule (Int_t module_number, TString channel_name)
Public Member Functions inherited from VQwSubsystemParity
	VQwSubsystemParity (const TString &name)
	Constructor with name. More...
	VQwSubsystemParity (const VQwSubsystemParity &source)
	Copy constructor. More...
virtual	~VQwSubsystemParity ()
	Default destructor. More...
virtual void	FillDB_MPS (QwParityDB *db, TString type)
	Fill the database with MPS-based variables Note that most subsystems don't need to do this. More...
virtual void	FillDB (QwParityDB *db, TString type)
	Fill the database. More...
virtual void	FillErrDB (QwParityDB *db, TString type)
virtual UInt_t	UpdateErrorFlag ()
	Uses the error flags of contained data elements to update Returns the error flag to the top level routines related to stability checks and ErrorFlag updates. More...
virtual void	Blind (const QwBlinder *blinder)
	Blind the asymmetry of this subsystem. More...
virtual void	Blind (const QwBlinder *blinder, const VQwSubsystemParity *subsys)
	Blind the difference of this subsystem. More...
virtual void	WritePromptSummary (QwPromptSummary *ps, TString type)
virtual Bool_t	CheckForEndOfBurst () const
Public Member Functions inherited from VQwSubsystem
	VQwSubsystem (const TString &name)
	Constructor with name. More...
	VQwSubsystem (const VQwSubsystem &orig)
	Copy constructor by object. More...
virtual	~VQwSubsystem ()
	Default destructor. More...
TString	GetSubsystemName () const
Bool_t	HasDataLoaded () const
void	SetParent (QwSubsystemArray *parent)
	Set the parent of this subsystem to the specified array. More...
QwSubsystemArray *	GetParent (const unsigned int parent=0) const
	Get the parent of this subsystem. More...
VQwSubsystem *	GetSibling (const std::string &name) const
	Get the sibling with specified name. More...
Bool_t	PublishInternalValue (const TString &name, const TString &desc, const VQwHardwareChannel *value) const
	Publish a variable name to the parent subsystem array. More...
virtual Bool_t	PublishInternalValues () const
	Publish all variables of the subsystem. More...
virtual Bool_t	PublishByRequest (TString device_name)
	Try to publish an internal variable matching the submitted name. More...
Bool_t	RequestExternalValue (const TString &name, VQwHardwareChannel *value) const
	Request a named value which is owned by an external subsystem; the request will be handled by the parent subsystem array. More...
virtual const VQwHardwareChannel *	ReturnInternalValue (const TString &name) const
	Return a pointer to a varialbe to the parent subsystem array to be delivered to a different subsystem. More...
virtual Bool_t	ReturnInternalValue (const TString &name, VQwHardwareChannel *value) const
	Return a named value to the parent subsystem array to be delivered to a different subsystem. More...
virtual std::vector< TString >	GetParamFileNameList ()
virtual std::map< TString, TString >	GetDetectorMaps ()
virtual Int_t	LoadDetectorMaps (QwParameterFile &file)
	Parse parameter file to find the map files. More...
virtual Int_t	LoadGeometryDefinition (TString mapfile)
	Optional geometry definition. More...
virtual Int_t	LoadCrosstalkDefinition (TString mapfile)
	Optional crosstalk definition. More...
void	SetEventTypeMask (const UInt_t mask)
	Set event type mask. More...
UInt_t	GetEventTypeMask () const
	Get event type mask. More...
virtual void	ExchangeProcessedData ()
	Request processed data from other subsystems for internal use in the second event processing stage. Not all derived classes will require data from other subsystems. More...
virtual void	ProcessEvent_2 ()
	Process the event data again, including data from other subsystems. Not all derived classes will require a second stage of event data processing. More...
virtual void	AtEndOfEventLoop ()
	Perform actions at the end of the event loop. More...
virtual void	RandomizeEventData (int helicity=0, double time=0.0)
virtual void	EncodeEventData (std::vector< UInt_t > &buffer)
virtual void	PrintInfo () const
	Print some information about the subsystem. More...
virtual void	PrintDetectorMaps (Bool_t status) const
virtual void	ConstructHistograms ()
	Construct the histograms for this subsystem. More...
virtual void	ConstructHistograms (TDirectory *folder)
	Construct the histograms for this subsystem in a folder. More...
virtual void	ConstructHistograms (TString &prefix)
	Construct the histograms for this subsystem with a prefix. More...
virtual void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)=0
	Construct the branch and tree vector. More...
virtual void	ConstructBranchAndVector (TTree *tree, std::vector< Double_t > &values)
	Construct the branch and tree vector. More...
virtual void	ConstructTree ()
	Construct the tree for this subsystem. More...
virtual void	ConstructTree (TDirectory *folder)
	Construct the tree for this subsystem in a folder. More...
virtual void	ConstructTree (TString &prefix)
	Construct the tree for this subsystem with a prefix. More...
virtual void	ConstructTree (TDirectory *folder, TString &prefix)
	Construct the tree for this subsystem in a folder with a prefix. More...
virtual void	FillTree ()
	Fill the tree for this subsystem. More...
virtual void	DeleteTree ()
	Delete the tree for this subsystem. More...
Public Member Functions inherited from MQwHistograms
void	ShareHistograms (const MQwHistograms *source)
	Share histogram pointers between objects. More...
Public Member Functions inherited from MQwCloneable< VQwSubsystem, QwMollerDetector >
virtual	~MQwCloneable ()
	Virtual destructor. More...
virtual VQwSubsystem *	Clone () const
	Concrete clone method. More...
const VQwFactory< VQwSubsystem > *	Factory () const
	Factory getter. More...
Public Member Functions inherited from VQwCloneable< VQwSubsystem >
virtual	~VQwCloneable ()
	Virtual destructor. More...
std::string	GetClassName () const
	Get demangled name of this class. More...

Protected Attributes
std::vector< QwMollerChannelID >	fMollerChannelID
std::vector< std::vector < QwSTR7200_Channel > >	fSTR7200_Channel
std::vector< std::vector < QwSTR7200_Channel > >	fPrevious_STR7200_Channel
Protected Attributes inherited from VQwSubsystem
std::map< TString, VQwHardwareChannel * >	fPublishedInternalValues
	Map of published internal values. More...
std::vector< std::vector < TString > >	fPublishList
	List of parameters to be published (loaded at the channel map) More...
TString	fSystemName
	Name of this subsystem. More...
UInt_t	fEventTypeMask
	Mask of event types. More...
Bool_t	fIsDataLoaded
	Has this subsystem gotten data to be processed? More...
std::vector< TString >	fDetectorMapsNames
std::map< TString, TString >	fDetectorMaps
Int_t	fCurrentROC_ID
	ROC ID that is currently being processed. More...
Int_t	fCurrentBank_ID
	Bank ID that is currently being processed. More...
std::vector< UInt_t >	fROC_IDs
	Vector of ROC IDs associated with this subsystem. More...
std::vector< std::vector < UInt_t > >	fBank_IDs
	Vector of Bank IDs per ROC ID associated with this subsystem. More...
std::vector< QwSubsystemArray * >	fArrays
	Vector of pointers to subsystem arrays that contain this subsystem. More...
Protected Attributes inherited from MQwHistograms
std::vector< TH1_ptr >	fHistograms
	Histograms associated with this data element. More...

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

Private Attributes
Int_t	fQwMollerErrorCount
UInt_t	fNumberOfEvents

Static Private Attributes
static const Bool_t	bDEBUG = kFALSE

Additional Inherited Members
Static Public Member Functions inherited from VQwSubsystem
static void	DefineOptions ()
	Define options function (note: no virtual static functions in C++) More...
Static Public Member Functions inherited from MQwCloneable< VQwSubsystem, QwMollerDetector >
static VQwSubsystem *	Create (const std::string &name)
	Object creation. More...
static QwMollerDetector *	Cast (QwMollerDetector *type)
	Object dynamic cast. More...
Protected Member Functions inherited from VQwSubsystem
void	UpdatePublishedValue (const TString &name, VQwHardwareChannel *data_channel)
void	ClearAllBankRegistrations ()
	Clear all registration of ROC and Bank IDs for this subsystem. More...
virtual Int_t	RegisterROCNumber (const UInt_t roc_id, const UInt_t bank_id=0)
	Tell the object that it will decode data from this ROC and sub-bank. More...
Int_t	RegisterSubbank (const UInt_t bank_id)
	Tell the object that it will decode data from this sub-bank in the ROC currently open for registration. More...
Int_t	GetSubbankIndex () const
Int_t	GetSubbankIndex (const UInt_t roc_id, const UInt_t bank_id) const
void	SetDataLoaded (Bool_t flag)
Int_t	FindIndex (const std::vector< UInt_t > &myvec, const UInt_t value) const
Bool_t	Compare (VQwSubsystem *source)
Protected Member Functions inherited from MQwHistograms
	MQwHistograms ()
	Default constructor. More...
	MQwHistograms (const MQwHistograms &source)
	Copy constructor. More...
virtual	~MQwHistograms ()
	Virtual destructor. More...
virtual MQwHistograms &	operator= (const MQwHistograms &value)
void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)
void	AddHistogram (TH1 *h)
	Register a histogram. More...

Detailed Description

Implementation of the analysis of Moller data (copied from QwComptonElectronDetector.h)

Author

Andrew Kubera

Date

2010-06-07

Definition at line 62 of file QwMollerDetector.h.

Constructor & Destructor Documentation

QwMollerDetector::QwMollerDetector ( )

private

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

QwMollerDetector::QwMollerDetector ( const TString &  name )

inline

Constructor with name.

Definition at line 73 of file QwMollerDetector.h.

    : VQwSubsystem(name), VQwSubsystemParity(name)
    { };

QwMollerDetector::QwMollerDetector ( const QwMollerDetector &  source )

inline

Copy constructor.

Definition at line 77 of file QwMollerDetector.h.

    : VQwSubsystem(source),VQwSubsystemParity(source),
      fSTR7200_Channel(source.fSTR7200_Channel)
    { }

virtual QwMollerDetector::~QwMollerDetector ( )

inlinevirtual

Virtual destructor.

Definition at line 82 of file QwMollerDetector.h.

{ };

Member Function Documentation

void QwMollerDetector::AccumulateRunningSum ( VQwSubsystem *  value )

virtual

Update the running sums for devices.

Implements VQwSubsystemParity.

Definition at line 317 of file QwMollerDetector.cc.

References Compare(), and fSTR7200_Channel.

                                                              {
  if (Compare(value)) {
    QwMollerDetector* v = dynamic_cast<QwMollerDetector*>(value);

    for (size_t i = 0; i < fSTR7200_Channel.size(); i++){
      for (size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
        fSTR7200_Channel[i][j].AccumulateRunningSum(v->fSTR7200_Channel[i][j]);
      }
    }

  }
}

Here is the call graph for this function:

Bool_t QwMollerDetector::ApplySingleEventCuts ( )

virtual

Apply the single event cuts.

Implements VQwSubsystemParity.

Definition at line 340 of file QwMollerDetector.cc.

References bDEBUG, and fSTR7200_Channel.

                                             {
  std::cout << "QwMoller::ApplySingleEventCuts() ";
  Bool_t test = kTRUE, test_1 = kTRUE;

  for (size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for (size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      test_1 = fSTR7200_Channel[i][j].ApplySingleEventCuts();
      test &= test_1;
      if(!test_1 && bDEBUG){
        std::cout << "***** QwMoller::SingleEventCuts()->Channel[" << i << "][" << j << "]:" << fSTR7200_Channel[i][j].GetElementName() << std::endl;
      }
    }
  }
  return test;
}

void QwMollerDetector::CalculateRunningAverage ( )

virtual

Calculate the average for all good events.

Implements VQwSubsystemParity.

Definition at line 330 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                              {
  for (size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for (size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].CalculateRunningAverage();
    }
  }
}

void QwMollerDetector::ClearEventData ( )

virtual

Implements VQwSubsystem.

Definition at line 137 of file QwMollerDetector.cc.

{}

Bool_t QwMollerDetector::Compare

(

VQwSubsystem *

source

)

Definition at line 418 of file QwMollerDetector.cc.

References fSTR7200_Channel.

Referenced by AccumulateRunningSum(), Difference(), operator+=(), operator-=(), operator=(), Ratio(), and Sum().

                                                    {
  //std::cout << "Beginning QwMollerDetector::Compare" << std::endl;
  
  if (source == 0) return kFALSE;

  Bool_t result = kTRUE;
  if(typeid(*source) != typeid(*this)){
    result = kFALSE;
    std::cout << " Type mismatch! This is bypassed for now but should be fixed eventually." << std::endl;
  } else { //same type, test for # of modules
    QwMollerDetector* input = dynamic_cast<QwMollerDetector*>(source);
    if(input->fSTR7200_Channel.size() != fSTR7200_Channel.size()){
      result = kFALSE;
      std::cout << " Not the same number of Modules" << std::endl;
    }else { //same # modules, loop through and make sure each one has same amount of channels
      for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
        if(input->fSTR7200_Channel[i].size() != fSTR7200_Channel[i].size()){
          result = kFALSE;
          std::cout << " Different number of channels in module " << i << std::endl;
        }
      }
    }
  }
  return result;
}

Here is the caller graph for this function:

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

inlinevirtual

Construct the branch and tree vector.

Implements VQwSubsystem.

Definition at line 125 of file QwMollerDetector.h.

{ };

void QwMollerDetector::ConstructBranch ( TTree *  tree, TString &  prefix, QwParameterFile &  trim_file  )

inlinevirtual

Construct the branch and tree vector based on the trim file.

Implements VQwSubsystem.

Definition at line 126 of file QwMollerDetector.h.

{ };

void QwMollerDetector::ConstructBranchAndVector	(	TTree *	,
		TString &	,
		std::vector< double, std::allocator< double > > &
	)

Definition at line 224 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                                                                                          {
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].ConstructBranchAndVector(tree, prefix, values);
    }
  }
}

void QwMollerDetector::ConstructHistograms ( TDirectory *  folder, TString &  prefix  )

virtual

Construct the histograms for this subsystem in a folder with a prefix.

Implements VQwSubsystem.

Definition at line 208 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                                                             {
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].ConstructHistograms(folder, prefix);
    }
  }
}

void QwMollerDetector::DeaccumulateRunningSum ( VQwSubsystem *  value )

inlinevirtual

remove one entry from the running sums for devices

Implements VQwSubsystemParity.

Definition at line 110 of file QwMollerDetector.h.

                                                    {
    };

void QwMollerDetector::Difference ( VQwSubsystem *  value1, VQwSubsystem *  value2  )

virtual

Implements VQwSubsystemParity.

Definition at line 288 of file QwMollerDetector.cc.

References Compare().

                                                                             {
  if (Compare(value1) && Compare(value2)) {
    *this  = value1;
    *this -= value2;
  }
}

Here is the call graph for this function:

void QwMollerDetector::FillHistograms ( )

virtual

Fill the histograms for this subsystem.

Implements VQwSubsystem.

Definition at line 216 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                     {
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].FillHistograms();
    }
  }
}

void QwMollerDetector::FillTreeVector ( std::vector< Double_t > &  values ) const

virtual

Fill the tree vector.

Implements VQwSubsystem.

Definition at line 232 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                                                       {
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].FillTreeVector(values);
    }
  }
}

Int_t QwMollerDetector::GetChannelIndex	(	TString	channelName,
		UInt_t	module_number
	)

Definition at line 393 of file QwMollerDetector.cc.

References QwMollerChannelID::fChannelName, QwMollerChannelID::fIndex, QwMollerChannelID::fModuleNumber, and fMollerChannelID.

Referenced by GetDataForChannelInModule().

{
  Bool_t ldebug=kFALSE;

  channelName.ToLower();

  if (ldebug){
    std::cout << "QwMollerDetector::GetDetectorIndex" << std::endl;
    std::cout << "module_number: " << module_number << " name=" << channelName << std::endl;
  }

  Int_t result = -1;
  for(size_t i = 0; i < fMollerChannelID.size(); i++){
    QwMollerChannelID *nextChannel = &fMollerChannelID[i];
    //std::cout << ' '<< nextChannel->fChannelName << '=' << channelName << ':' << (nextChannel->fChannelName == channelName) << std::endl;
    if (nextChannel->fChannelName == channelName && nextChannel->fModuleNumber == module_number){
      result = nextChannel->fIndex;
      break;
    }
  }

  return result;
}

Here is the caller graph for this function:

float QwMollerDetector::GetDataForChannelInModule ( Int_t  module_number, Int_t  channel_index  )

inline

Definition at line 135 of file QwMollerDetector.h.

References fSTR7200_Channel.

Referenced by GetDataForChannelInModule().

                                                                             {
      return fSTR7200_Channel[module_number][channel_index].GetValue();
    }

Here is the caller graph for this function:

float QwMollerDetector::GetDataForChannelInModule ( Int_t  module_number, TString  channel_name  )

inline

Definition at line 139 of file QwMollerDetector.h.

References GetChannelIndex(), and GetDataForChannelInModule().

                                                                              {
      return GetDataForChannelInModule(module_number, GetChannelIndex(channel_name,module_number));
    }

Here is the call graph for this function:

UInt_t QwMollerDetector::GetEventcutErrorFlag ( )

virtual

Return the error flag to the top level routines related to stability checks and ErrorFlag updates.

Implements VQwSubsystemParity.

Definition at line 367 of file QwMollerDetector.cc.

                                             {
  return 0;
}

float * QwMollerDetector::GetRawChannelArray

(

)

Definition at line 371 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                           {
  size_t len = 0;
  for (size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for (size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      len++;
    }
  }
  float *result = new float[len];
  
  //float result[96];

  int n = 0;
  for (size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for (size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      result[n+j] = fSTR7200_Channel[i][j].GetValue();
    }
    n=fSTR7200_Channel[i].size();
  }
 
  return result;
}

void QwMollerDetector::IncrementErrorCounters ( )

inlinevirtual

Increment the error counters.

Implements VQwSubsystemParity.

Definition at line 115 of file QwMollerDetector.h.

{};

Int_t QwMollerDetector::LoadChannelMap ( TString  mapfile )

virtual

Load the channel map

Parameters

mapfile

Map file

Returns

Zero if successful

Implements VQwSubsystem.

Definition at line 33 of file QwMollerDetector.cc.

References QwMollerChannelID::fChannelName, QwMollerChannelID::fChannelNumber, VQwSubsystem::fDetectorMapsNames, QwMollerChannelID::fDetectorType, QwMollerChannelID::fIndex, QwMollerChannelID::fModuleNumber, QwMollerChannelID::fModuleType, fMollerChannelID, fPrevious_STR7200_Channel, fSTR7200_Channel, QwMollerChannelID::fWordInSubbank, VQwSubsystem::GetSubbankIndex(), QwParameterFile::GetUInt(), VQwSubsystem::RegisterROCNumber(), and VQwSubsystem::RegisterSubbank().

{
  TString varname, varvalue;
  TString modtype, dettype, name, keyword;
  Int_t modnum, channum;
  Int_t currentrocread = 0;
  Int_t currentbankread = 0;
  Int_t wordsofar = 0;
  Int_t currentsubbankindex = -1;
  
  QwParameterFile mapstr(mapfile.Data());  // Open the file
  fDetectorMapsNames.push_back(mapstr.GetParamFilename());
  while (mapstr.ReadNextLine()) {
    mapstr.TrimComment('!');   // Remove everything after a '!' character.
    mapstr.TrimWhitespace();   // Get rid of leading and trailing whitespace (spaces or tabs).
    if (mapstr.LineIsEmpty())  continue;

    if (mapstr.HasVariablePair("=", varname, varvalue)) {
      // This is a declaration line.  Decode it.
      
      varname.ToLower();
      UInt_t value = QwParameterFile::GetUInt(varvalue);
      if (varname == "roc") {
        currentrocread=value;
        RegisterROCNumber(value,0);
      } else if (varname == "bank") {
        currentbankread=value;
        RegisterSubbank(value);
        if(currentsubbankindex != GetSubbankIndex(currentrocread,currentbankread)){
          currentsubbankindex = GetSubbankIndex(currentrocread,currentbankread);
        }
      }
    } else {
      Bool_t lineok = kTRUE;

      //  Break this line into tokens to process it.
      modtype   = mapstr.GetTypedNextToken<TString>();
      modnum    = mapstr.GetTypedNextToken<Int_t>();
      channum   = mapstr.GetTypedNextToken<Int_t>();
      dettype   = mapstr.GetTypedNextToken<TString>();
      name      = mapstr.GetTypedNextToken<TString>();
      dettype.ToLower();
      name.ToLower();

      QwMollerChannelID newChannelID;
      newChannelID.fModuleType    = modtype;
      newChannelID.fModuleNumber  = modnum;
      newChannelID.fChannelName   = name;
      newChannelID.fDetectorType  = dettype;
      newChannelID.fChannelNumber = channum;
      newChannelID.fWordInSubbank = wordsofar;
       
      if (modtype == "SIS3801"){
        wordsofar += 1;
      }else if(modtype == "SIS7200"){
        wordsofar += 1;
      }else if(modtype == "VQWK"){
        wordsofar += 6;
      }else {
        std::cerr << "Unknown module type: " << modtype << ". Skipping channel "  << name << '.' << std::endl;
        lineok = kFALSE;
      }

      if (name.Length() == 0){
        lineok = kFALSE;
      }

//    add new modules until current number (modnum) is reached 
      std::size_t chan_size;
      chan_size = fSTR7200_Channel.size();
      while ((Int_t) chan_size <= modnum) {
        std::vector<QwSTR7200_Channel> new_module;
        fSTR7200_Channel.push_back(new_module);
      }

      //change this line if names are supposed to be exclusive, as of now, 
      newChannelID.fIndex = -1; // GetChannelIndex(name, modnum);

      if (newChannelID.fIndex == -1 && lineok){
        QwSTR7200_Channel localSTR7200_Channel(newChannelID.fChannelName);
        fSTR7200_Channel[modnum].push_back(localSTR7200_Channel);
//        fSTR7200_Channel[modnum][fMollerChannelID.size() - 1].SetDefaultSampleSize(fSample_size);
        newChannelID.fIndex = fSTR7200_Channel[modnum].size() - 1;
        //std::cout << name << ':' << newChannelID.fIndex << ':' << wordsofar << '\n';
      }

      //  Push a new record into the element array
      if(lineok){
        fMollerChannelID.push_back(newChannelID);
      }
    }
  } // end looping over parameter file

  fPrevious_STR7200_Channel = fSTR7200_Channel;

//  for (size_t i = 0; i < fMollerChannelID.size(); i++){
//     std::cout << i << ": " << fMollerChannelID[i].fChannelName << ' ' << fMollerChannelID[i].fChannelNumber << ' ' << fMollerChannelID[i].fIndex << std::endl;
//   }

  return 0;
}

Here is the call graph for this function:

Int_t QwMollerDetector::LoadEventCuts

(

TString &

filename

)

Int_t QwMollerDetector::LoadEventCuts ( TString  filename )

virtual

Load the event cuts file.

Implements VQwSubsystemParity.

Definition at line 338 of file QwMollerDetector.cc.

{return 0;}

Int_t QwMollerDetector::LoadInputParameters ( TString  mapfile )

virtual

Mandatory parameter file definition.

Implements VQwSubsystem.

Definition at line 136 of file QwMollerDetector.cc.

{ return 0;}

VQwSubsystem & QwMollerDetector::operator+= ( VQwSubsystem *  value )

virtual

Implements VQwSubsystemParity.

Definition at line 254 of file QwMollerDetector.cc.

References Compare(), and fSTR7200_Channel.

                                                              {
  //std::cout << "QwMollerDetector addition assignment (operator+=)" << std::endl;
  if(Compare(value)){
    QwMollerDetector* input = dynamic_cast<QwMollerDetector *> (value);
    for(size_t i = 0; i < input->fSTR7200_Channel.size(); i++){
      for(size_t j = 0; j < input->fSTR7200_Channel[i].size(); j++){
        this->fSTR7200_Channel[i][j] += input->fSTR7200_Channel[i][j];
        //std::cout << "+= " << this->fSTR7200_Channel[i][j].GetValue() << std::endl;
      }
    }
  }
  return *this;
}

Here is the call graph for this function:

VQwSubsystem & QwMollerDetector::operator-= ( VQwSubsystem *  value )

virtual

Implements VQwSubsystemParity.

Definition at line 268 of file QwMollerDetector.cc.

References Compare(), and fSTR7200_Channel.

                                                              {
  //std::cout << "QwMollerDetector subtraction assignment (operator-=)" << std::endl;
  if(Compare(value)){
    QwMollerDetector* input = dynamic_cast<QwMollerDetector *> (value);
    for(size_t i = 0; i < input->fSTR7200_Channel.size(); i++){
      for(size_t j = 0; j < input->fSTR7200_Channel[i].size(); j++){
        this->fSTR7200_Channel[i][j] -= input->fSTR7200_Channel[i][j];
      }
    }
  }
  return *this;
}

Here is the call graph for this function:

VQwSubsystem & QwMollerDetector::operator= ( VQwSubsystem *  value )

virtual

Assignment Note: Must be called at the beginning of all subsystems routine call to operator=(VQwSubsystem *value) by VQwSubsystem::operator=(value)

Implements VQwSubsystemParity.

Definition at line 240 of file QwMollerDetector.cc.

References Compare(), and fSTR7200_Channel.

                                                             {
  // std::cout << "QwMollerDetector assignment (operator=)" << std::endl;
  if(Compare(value)){
    //VQwSubsystem::operator=(value);
    QwMollerDetector* input = dynamic_cast<QwMollerDetector *> (value);
    for(size_t i = 0; i < input->fSTR7200_Channel.size(); i++){
      for(size_t j = 0; j < input->fSTR7200_Channel[i].size(); j++){
        this->fSTR7200_Channel[i][j] = input->fSTR7200_Channel[i][j];
      }
    }
  }
  return *this; 
}

Here is the call graph for this function:

void QwMollerDetector::print

(

)

Definition at line 444 of file QwMollerDetector.cc.

References fMollerChannelID, and fSTR7200_Channel.

                            {
  std::cout << " " << fSTR7200_Channel.size() << std::endl;
  UInt_t max = 0;
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    UInt_t next = fSTR7200_Channel[i].size();
    if (next > max){
      max = next;
    }
  }

  for(size_t i = 0; i < max; i++){
    std::cout << fMollerChannelID[i].fChannelName << ":\t" << (fMollerChannelID[i].fChannelName.Length() < 14 ? "\t" : "");
    for(size_t j = 0; j < fSTR7200_Channel.size(); j++){

      if ( i < fSTR7200_Channel[j].size()){
        std::cout << "0x" << std::hex << (int)fSTR7200_Channel[j][i].GetValue() << std::dec;
      } else {
        std::cout << "            ";
      }
      std::cout << "\t\t";
    }
    std::cout << std::endl;
  }

}

void QwMollerDetector::PrintErrorCounters ( ) const

virtual

Report the number of events failed due to HW and event cut failures.

Implements VQwSubsystemParity.

Definition at line 356 of file QwMollerDetector.cc.

References fQwMollerErrorCount, and fSTR7200_Channel.

                                               {
  std::cout << "***************QwMoller Error Summary****************" << std::endl;
  for (size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for (size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].PrintErrorCounters();
    }
  }
  std::cout << "total failed events: " << fQwMollerErrorCount << std::endl;
  std::cout << "************End QwMoller Error Summary***************" << std::endl;
}

void QwMollerDetector::PrintValue ( ) const

virtual

Print values of all channels.

Reimplemented from VQwSubsystemParity.

Definition at line 470 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                        {
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
         fSTR7200_Channel[i][j].PrintValue();
    }
  }
}

Int_t QwMollerDetector::ProcessConfigurationBuffer ( const UInt_t  roc_id, const UInt_t  bank_id, UInt_t *  buffer, UInt_t  num_words  )

virtual

Implements VQwSubsystem.

Definition at line 140 of file QwMollerDetector.cc.

                                                                                 {
  return 0;
}

Int_t QwMollerDetector::ProcessConfigurationBuffer	(	UInt_t	ev_type,
		const UInt_t	roc_id,
		const UInt_t	bank_id,
		UInt_t *	buffer,
		UInt_t	num_words
	)

Definition at line 144 of file QwMollerDetector.cc.

{
 return 0;
}

Int_t QwMollerDetector::ProcessEvBuffer ( UInt_t  roc_id, UInt_t  bank_id, UInt_t *  buffer, UInt_t  num_words  )

virtual

TODO: The non-event-type-aware ProcessEvBuffer routine should be replaced with the event-type-aware version.

Implements VQwSubsystem.

Definition at line 154 of file QwMollerDetector.cc.

References fMollerChannelID, and fSTR7200_Channel.

Referenced by ProcessEvBuffer().

{
  Int_t index = 0; // GetSubbankIndex(roc_id, bank_id);

  if (index >= 0 && num_words > 0){
    //we want to process this ROC

    //Loop through all the channels as defined in map file.. should be 96 in standard 3 module, 32 channel Moller configuration
    for(size_t i = 0; i < fMollerChannelID.size(); i++){

      Int_t wordsLeft = num_words - fMollerChannelID[i].fWordInSubbank;

//       std::cout << std::dec << fMollerChannelID[i].fChannelName << ' ' << fMollerChannelID[i].fIndex << " in module " << fMollerChannelID[i].fModuleNumber << " has data " << std::hex << buffer[fMollerChannelID[i].fWordInSubbank] << std::dec << std::endl;

      UInt_t mod_index = fMollerChannelID[i].fModuleNumber;
      Int_t chan_index = fMollerChannelID[i].fIndex;
//      char *junk = new char[2];
      if (mod_index < fSTR7200_Channel.size()) {
        //mod_index should be less than 3, next_module is reference to vector where the Channel should be located
     //   std::cout << std::hex << buffer[fMollerChannelID[i].fWordInSubbank] << '\n';
        //std::cin.getline(junk, 1);
//         std::cout << '(' << mod_index << ',' << chan_index << ')' << buffer[fMollerChannelID[i].fWordInSubbank] - fSTR7200_Channel[mod_index][chan_index].GetValue() << ' ' << std::endl;

        fSTR7200_Channel[mod_index][chan_index].ProcessEvBuffer(&(buffer[fMollerChannelID[i].fWordInSubbank]), wordsLeft);
      } else {
        std::cout << "Problem reading buffer, incorrect structure of map file?" << std::endl;
      }

    }
  }
//  std::cout << fSTR7200_Channel[2][0].GetValue() << std::endl;
//  print();
  return 0;
}

Here is the caller graph for this function:

Int_t QwMollerDetector::ProcessEvBuffer ( UInt_t  ev_type, UInt_t  roc_id, UInt_t  bank_id, UInt_t *  buffer, UInt_t  num_words  )

virtual

Reimplemented from VQwSubsystem.

Definition at line 149 of file QwMollerDetector.cc.

References ProcessEvBuffer().

{
  return ProcessEvBuffer(roc_id, bank_id, buffer, num_words);
}

Here is the call graph for this function:

void QwMollerDetector::ProcessEvent ( )

virtual

Implements VQwSubsystem.

Definition at line 191 of file QwMollerDetector.cc.

References fPrevious_STR7200_Channel, and fSTR7200_Channel.

{
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].ProcessEvent();
      QwSTR7200_Channel tempscaler(fSTR7200_Channel[i][j]);
      tempscaler = fSTR7200_Channel[i][j];
      fSTR7200_Channel[i][j] -= fPrevious_STR7200_Channel[i][j];
      fPrevious_STR7200_Channel[i][j] = tempscaler;
      // Store a temproary copy of this channel's raw value as a scaler channel
      // Subtract the corresponding fPrevious_STR7200_Channel from this scaler channel
      // Put the temprary copy into the fPrevious_STR7200_Channel
    }
  }
}

void QwMollerDetector::ProcessOptions ( QwOptions &  options )

virtual

Process the command line options.

Reimplemented from VQwSubsystem.

Definition at line 135 of file QwMollerDetector.cc.

{}

void QwMollerDetector::Ratio ( VQwSubsystem *  value1, VQwSubsystem *  value2  )

virtual

Implements VQwSubsystemParity.

Definition at line 295 of file QwMollerDetector.cc.

References Compare(), and fSTR7200_Channel.

                                                                        {
  if (Compare(value1) && Compare(value2)) {
    QwMollerDetector* v1 = dynamic_cast<QwMollerDetector *> (value1);
    QwMollerDetector* v2 = dynamic_cast<QwMollerDetector *> (value2);

    for(size_t i = 0; i < v1->fSTR7200_Channel.size(); i++){
      for(size_t j = 0; j < v1->fSTR7200_Channel[i].size(); j++){
        fSTR7200_Channel[i][j].Ratio(v1->fSTR7200_Channel[i][j],v2->fSTR7200_Channel[i][j]);
      }
    }
  }
  return;
}

Here is the call graph for this function:

void QwMollerDetector::Scale ( Double_t  factor )

virtual

Implements VQwSubsystemParity.

Definition at line 309 of file QwMollerDetector.cc.

References fSTR7200_Channel.

                                           {
  for(size_t i = 0; i < fSTR7200_Channel.size(); i++){
    for(size_t j = 0; j < fSTR7200_Channel[i].size(); j++){
      fSTR7200_Channel[i][j].Scale(factor);
    }
  }
}

Bool_t QwMollerDetector::SingleEventCuts

(

)

void QwMollerDetector::Sum ( VQwSubsystem *  value1, VQwSubsystem *  value2  )

virtual

Implements VQwSubsystemParity.

Definition at line 281 of file QwMollerDetector.cc.

References Compare().

                                                                      {
  if (Compare(value1) && Compare(value2)) {
    *this  = value1;
    *this += value2;
  }
}

Here is the call graph for this function:

void QwMollerDetector::UpdateErrorFlag ( const VQwSubsystem *  ev_error )

inlinevirtual

update the error flag in the subsystem level from the top level routines related to stability checks. This will uniquely update the errorflag at each channel based on the error flag in the corresponding channel in the ev_error subsystem

Implements VQwSubsystemParity.

Definition at line 120 of file QwMollerDetector.h.

                                                      {
    };

Field Documentation

const Bool_t QwMollerDetector::bDEBUG = kFALSE

staticprivate

Definition at line 153 of file QwMollerDetector.h.

Referenced by ApplySingleEventCuts().

std::vector<QwMollerChannelID> QwMollerDetector::fMollerChannelID

protected

Definition at line 145 of file QwMollerDetector.h.

Referenced by GetChannelIndex(), LoadChannelMap(), print(), and ProcessEvBuffer().

UInt_t QwMollerDetector::fNumberOfEvents

private

Definition at line 154 of file QwMollerDetector.h.

std::vector< std::vector<QwSTR7200_Channel> > QwMollerDetector::fPrevious_STR7200_Channel

protected

Definition at line 149 of file QwMollerDetector.h.

Referenced by LoadChannelMap(), and ProcessEvent().

Int_t QwMollerDetector::fQwMollerErrorCount

private

Definition at line 152 of file QwMollerDetector.h.

Referenced by PrintErrorCounters().

std::vector< std::vector<QwSTR7200_Channel> > QwMollerDetector::fSTR7200_Channel

protected

Definition at line 148 of file QwMollerDetector.h.

Referenced by AccumulateRunningSum(), ApplySingleEventCuts(), CalculateRunningAverage(), Compare(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), GetDataForChannelInModule(), GetRawChannelArray(), LoadChannelMap(), operator+=(), operator-=(), operator=(), print(), PrintErrorCounters(), PrintValue(), ProcessEvBuffer(), ProcessEvent(), Ratio(), and Scale().

---

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

• QwMollerDetector.h
• QwMollerDetector.cc