QwComptonPhotonDetector Class Reference

Class for the analysis of Compton photon detector data. More...

#include <QwComptonPhotonDetector.h>

Inheritance diagram for QwComptonPhotonDetector:

Collaboration diagram for QwComptonPhotonDetector:

Public Member Functions
	QwComptonPhotonDetector (TString name)
	Constructor with name. More...
	QwComptonPhotonDetector (const QwComptonPhotonDetector &source)
	Copy constructor. More...
virtual	~QwComptonPhotonDetector ()
	Virtual destructor. More...
void	ProcessOptions (QwOptions &options)
	Process the command line options. More...
Int_t	LoadChannelMap (TString mapfile)
Int_t	LoadInputParameters (TString pedestalfile)
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	ProcessEvBuffer (const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
void	ClearEventData ()
void	ProcessEvent ()
Bool_t	IsGoodEvent ()
void	RandomizeEventData (int helicity=0)
void	EncodeEventData (std::vector< UInt_t > &buffer)
VQwSubsystem &	operator= (VQwSubsystem *value)
VQwSubsystem &	operator+= (VQwSubsystem *value)
VQwSubsystem &	operator-= (VQwSubsystem *value)
void	Sum (VQwSubsystem *value1, VQwSubsystem *value2)
void	Difference (VQwSubsystem *value1, VQwSubsystem *value2)
void	Ratio (VQwSubsystem *numer, VQwSubsystem *denom)
void	Scale (Double_t factor)
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...
Bool_t	CheckRunningAverages (Bool_t)
void	AccumulateRunningSum (VQwSubsystem *value)
void	DeaccumulateRunningSum (VQwSubsystem *value)
	remove one entry from the running sums for devices More...
void	CalculateRunningAverage ()
void	ConstructHistograms (TDirectory *folder, TString &prefix)
void	FillHistograms ()
void	ConstructTree (TDirectory *folder, TString &prefix)
void	FillTree ()
void	DeleteTree ()
void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
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
Bool_t	Compare (VQwSubsystem *source)
Bool_t	CompareADC (VQwSubsystem *source)
Bool_t	CompareTDC (VQwSubsystem *source)
Bool_t	CompareQDC (VQwSubsystem *source)
void	PrintValue () const
void	PrintInfo () const
QwSIS3320_Channel *	GetSIS3320Channel (const TString 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 Int_t	ProcessEvBuffer (const UInt_t event_type, const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
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	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, 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...
Public Member Functions inherited from MQwHistograms
void	ShareHistograms (const MQwHistograms *source)
	Share histogram pointers between objects. More...
Public Member Functions inherited from MQwV775TDC
	MQwV775TDC ()
	~MQwV775TDC ()
void	DecodeTDCWord (UInt_t &word, const UInt_t roc_id=0)
Bool_t	IsValidDataword ()
Bool_t	IsHeaderword ()
UInt_t	GetTDCSlotNumber ()
UInt_t	GetTDCChannelNumber ()
UInt_t	GetTDCData ()
UInt_t	GetTDCMaxChannels ()
UInt_t	GetTDCEventNumber ()
UInt_t	GetTDCTriggerTime ()
void	SetReferenceParameters (Double_t mindiff, Double_t maxdiff, Double_t offset, Double_t shift)
Double_t	SubtractReference (Double_t rawtime, Double_t reftime)
Bool_t	CheckDataIntegrity (const UInt_t roc_id, UInt_t *buffer, UInt_t num_words)
void	PrintTDCHeader (Bool_t flag)
void	PrintTDCData (Bool_t flag)
Public Member Functions inherited from MQwCloneable< VQwSubsystem, QwComptonPhotonDetector >
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...

Static Public Member Functions
static void	DefineOptions (QwOptions &options)
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, QwComptonPhotonDetector >
static VQwSubsystem *	Create (const std::string &name)
	Object creation. More...
static QwComptonPhotonDetector *	Cast (QwComptonPhotonDetector *type)
	Object dynamic cast. More...

Protected Types
enum	ChannelType_t { kUnknown, kSamplingADC, kIntegratingADC, kIntegratingTDC }
	Mapping from ROC/subbank to channel type. More...
typedef std::map< Int_t, std::vector< std::vector < Int_t > > >	SamplingADC_Mapping_t
	List of sampling ADC channels. More...
typedef std::map< Int_t, std::vector< std::vector < Int_t > > >	IntegratingADC_Mapping_t
	List of integrating QDC channels. More...
typedef std::map< Int_t, std::vector< std::vector < Int_t > > >	IntegratingTDC_Mapping_t
	List of integrating TDC channels. More...

Protected Attributes
TTree *	fTree
	Expert tree. More...
Int_t	fTree_fNEvents
	Expert tree fields. More...
std::map< Int_t, ChannelType_t >	fMapping
SamplingADC_Mapping_t	fSamplingADC_Mapping
std::vector< QwSIS3320_Channel >	fSamplingADC
IntegratingADC_Mapping_t	fMultiQDC_Mapping
std::vector< QwPMT_Channel >	fMultiQDC_Channel
std::vector< std::vector < QwPMT_Channel > >	fMultiQDC_Events
IntegratingTDC_Mapping_t	fMultiTDC_Mapping
std::vector< QwPMT_Channel >	fMultiTDC_Channel
std::vector< std::vector < QwPMT_Channel > >	fMultiTDC_Events
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
	QwComptonPhotonDetector ()
	Private default constructor (not implemented, will throw linker error on use) More...

Private Attributes
Int_t	fGoodEventCount

Static Private Attributes
static const Bool_t	kDebug = kTRUE

Additional Inherited Members
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

Class for the analysis of Compton photon detector data.

Author

W. Deconinck

Date

2009-09-04 18:06:23

The QwComptonPhotonDetector class is defined as a parity subsystem that contains all data modules of the photon detector (SIS3320, V775, V792). It reads in a channel map and pedestal file, and defines the histograms and output trees.

Definition at line 32 of file QwComptonPhotonDetector.h.

Member Typedef Documentation

typedef std::map< Int_t, std::vector< std::vector< Int_t > > > QwComptonPhotonDetector::IntegratingADC_Mapping_t

protected

List of integrating QDC channels.

Definition at line 142 of file QwComptonPhotonDetector.h.

typedef std::map< Int_t, std::vector< std::vector< Int_t > > > QwComptonPhotonDetector::IntegratingTDC_Mapping_t

protected

List of integrating TDC channels.

Definition at line 148 of file QwComptonPhotonDetector.h.

typedef std::map< Int_t, std::vector <std::vector <Int_t> > > QwComptonPhotonDetector::SamplingADC_Mapping_t

protected

List of sampling ADC channels.

Definition at line 137 of file QwComptonPhotonDetector.h.

Member Enumeration Documentation

enum QwComptonPhotonDetector::ChannelType_t

protected

Mapping from ROC/subbank to channel type.

Enumerator
kUnknown
kSamplingADC
kIntegratingADC
kIntegratingTDC

Definition at line 133 of file QwComptonPhotonDetector.h.

{ kUnknown, kSamplingADC, kIntegratingADC, kIntegratingTDC };

Constructor & Destructor Documentation

QwComptonPhotonDetector::QwComptonPhotonDetector ( )

private

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

QwComptonPhotonDetector::QwComptonPhotonDetector ( TString  name )

inline

Constructor with name.

Definition at line 44 of file QwComptonPhotonDetector.h.

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

QwComptonPhotonDetector::QwComptonPhotonDetector ( const QwComptonPhotonDetector &  source )

inline

Copy constructor.

Definition at line 46 of file QwComptonPhotonDetector.h.

    : VQwSubsystem(source),VQwSubsystemParity(source),
      fSamplingADC(source.fSamplingADC),
      fMultiQDC_Channel(source.fMultiQDC_Channel),
      fMultiTDC_Channel(source.fMultiTDC_Channel)
    { }

virtual QwComptonPhotonDetector::~QwComptonPhotonDetector ( )

inlinevirtual

Virtual destructor.

Definition at line 53 of file QwComptonPhotonDetector.h.

{ };

Member Function Documentation

void QwComptonPhotonDetector::AccumulateRunningSum ( VQwSubsystem *  value )

virtual

Accumulate the running sum

Implements VQwSubsystemParity.

Definition at line 641 of file QwComptonPhotonDetector.cc.

References Compare(), and fGoodEventCount.

{
  if (Compare(value)) {
    fGoodEventCount++;
    *this  += value;
  }
}

Here is the call graph for this function:

Bool_t QwComptonPhotonDetector::ApplySingleEventCuts ( )

inlinevirtual

Apply the single event cuts.

Implements VQwSubsystemParity.

Definition at line 84 of file QwComptonPhotonDetector.h.

{ return kTRUE; };

void QwComptonPhotonDetector::CalculateRunningAverage ( )

virtual

Normalize the running sum

Implements VQwSubsystemParity.

Definition at line 652 of file QwComptonPhotonDetector.cc.

References fGoodEventCount, and Scale().

{
  if (fGoodEventCount <= 0) {
    Scale(0);
  } else {
    Scale(1.0/fGoodEventCount);
  }
}

Here is the call graph for this function:

Bool_t QwComptonPhotonDetector::CheckRunningAverages ( Bool_t  )

inline

Definition at line 92 of file QwComptonPhotonDetector.h.

{ return kTRUE; };

void QwComptonPhotonDetector::ClearEventData ( )

virtual

Clear the event data in this subsystem

Implements VQwSubsystem.

Definition at line 514 of file QwComptonPhotonDetector.cc.

References fGoodEventCount, fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

{
  // Clear all sampling ADC channels
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].ClearEventData();

  // Clear all TDC and QDC channels
  for (size_t i = 0; i < fMultiTDC_Channel.size(); i++)
    fMultiTDC_Channel[i].ClearEventData();
  for (size_t i = 0; i < fMultiQDC_Channel.size(); i++)
    fMultiQDC_Channel[i].ClearEventData();

  fGoodEventCount = 0;
}

Bool_t QwComptonPhotonDetector::Compare

(

VQwSubsystem *

value

)

Compare two QwComptonPhotonDetector objects

Parameters

value

Object to compare with

Returns

kTRUE if the object is equal

Definition at line 666 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

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

{
  // Immediately fail on null objects
  if (value == 0) return kFALSE;

  // Continue testing on actual object
  Bool_t result = kTRUE;
  if (typeid(*value) != typeid(*this)) {
    result = kFALSE;

  } else {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*> (value);
    if (input->fSamplingADC.size() != fSamplingADC.size()) {
      result = kFALSE;
    } else {
      if (input->fMultiTDC_Channel.size() != fMultiTDC_Channel.size()) {
        result = kFALSE;
      } else {
        if (input->fMultiQDC_Channel.size() != fMultiQDC_Channel.size()) {
          result = kFALSE;
        }
      }
    }
  }
  return result;
}

Here is the caller graph for this function:

Bool_t QwComptonPhotonDetector::CompareADC

(

VQwSubsystem *

value

)

Compare two QwComptonPhotonDetector ADC objects

Parameters

value

Object to compare with

Returns

kTRUE if the object is equal

Definition at line 698 of file QwComptonPhotonDetector.cc.

References fSamplingADC.

Referenced by operator=().

{
  // Immediately fail on null objects
  if (value == 0) return kFALSE;

  // Continue testing on actual object
  Bool_t result = kTRUE;
  if (typeid(*value) != typeid(*this)) {
    result = kFALSE;

  } else {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*> (value);
    if (input->fSamplingADC.size() != fSamplingADC.size()) {
      result = kFALSE;
    }
  }
  return result;
}

Here is the caller graph for this function:

Bool_t QwComptonPhotonDetector::CompareQDC

(

VQwSubsystem *

value

)

Compare two QwComptonPhotonDetector QDC objects

Parameters

value

Object to compare with

Returns

kTRUE if the object is equal

Definition at line 746 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel.

Referenced by operator=().

{
  // Immediately fail on null objects
  if (value == 0) return kFALSE;

  // Continue testing on actual object
  Bool_t result = kTRUE;
  if (typeid(*value) != typeid(*this)) {
    result = kFALSE;

  } else {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*> (value);
    if (input->fMultiQDC_Channel.size() != fMultiQDC_Channel.size()) {
      result = kFALSE;
    }
  }
  return result;
}

Here is the caller graph for this function:

Bool_t QwComptonPhotonDetector::CompareTDC

(

VQwSubsystem *

value

)

Compare two QwComptonPhotonDetector TDC objects

Parameters

value

Object to compare with

Returns

kTRUE if the object is equal

Definition at line 722 of file QwComptonPhotonDetector.cc.

References fMultiTDC_Channel.

Referenced by operator=().

{
  // Immediately fail on null objects
  if (value == 0) return kFALSE;

  // Continue testing on actual object
  Bool_t result = kTRUE;
  if (typeid(*value) != typeid(*this)) {
    result = kFALSE;

  } else {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*> (value);
    if (input->fMultiTDC_Channel.size() != fMultiTDC_Channel.size()) {
      result = kFALSE;
    }
  }
  return result;
}

Here is the caller graph for this function:

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

inlinevirtual

Construct the branch and tree vector.

Implements VQwSubsystem.

Definition at line 111 of file QwComptonPhotonDetector.h.

{ };

void QwComptonPhotonDetector::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 112 of file QwComptonPhotonDetector.h.

{ };

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

virtual

Construct the tree

Parameters

tree	Pointer to the tree in which the branches will be created
prefix	Prefix with information about the type of histogram
values

Implements VQwSubsystem.

Definition at line 830 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

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

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

virtual

Construct the histograms

Parameters

folder	Folder in which the histograms will be created
prefix	Prefix with information about the type of histogram

Implements VQwSubsystem.

Definition at line 770 of file QwComptonPhotonDetector.cc.

References VQwSubsystem::ConstructHistograms(), fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

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

Here is the call graph for this function:

void QwComptonPhotonDetector::ConstructTree ( TDirectory *  folder, TString &  prefix  )

virtual

Construct the tree

Parameters

folder	Folder in which the tree will be created
prefix	Prefix with information about the type of histogram

Reimplemented from VQwSubsystem.

Definition at line 798 of file QwComptonPhotonDetector.cc.

References fTree, and fTree_fNEvents.

{
  folder->cd();
  fTree = new TTree("ComptonPhoton", "Compton Photon Detector");
  fTree->Branch("nevents",&fTree_fNEvents,"nevents/I");
}

void QwComptonPhotonDetector::DeaccumulateRunningSum ( VQwSubsystem *  value )

inlinevirtual

remove one entry from the running sums for devices

Implements VQwSubsystemParity.

Definition at line 96 of file QwComptonPhotonDetector.h.

                                                    {
    };

void QwComptonPhotonDetector::DefineOptions ( QwOptions &  options )

static

Definition at line 31 of file QwComptonPhotonDetector.cc.

{
  // Define command line options
}

void QwComptonPhotonDetector::DeleteTree ( )

virtual

Delete the tree

Reimplemented from VQwSubsystem.

Definition at line 808 of file QwComptonPhotonDetector.cc.

References fTree.

{
  delete fTree;
}

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

virtual

Difference

Parameters

value1	First value
value2	Second value

Implements VQwSubsystemParity.

Definition at line 605 of file QwComptonPhotonDetector.cc.

References Compare().

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

Here is the call graph for this function:

void QwComptonPhotonDetector::EncodeEventData ( std::vector< UInt_t > &  buffer )

virtual

Encode this event into a CODA buffer

Parameters

buffer

Buffer to append data to

Reimplemented from VQwSubsystem.

Definition at line 300 of file QwComptonPhotonDetector.cc.

References VQwSubsystem::fCurrentBank_ID, VQwSubsystem::fCurrentROC_ID, and fSamplingADC.

{
  std::vector<UInt_t> elements;
  elements.clear();

  // Get all buffers in the order they are defined in the map file
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].EncodeEventData(elements);

  // If there is element data, generate the subbank header
  std::vector<UInt_t> subbankheader;
  std::vector<UInt_t> rocheader;
  if (elements.size() > 0) {

    // Form CODA subbank header
    subbankheader.clear();
    subbankheader.push_back(elements.size() + 1); // subbank size
    subbankheader.push_back((fCurrentBank_ID << 16) | (0x01 << 8) | (1 & 0xff));
    // subbank tag | subbank type | event number

    // Form CODA bank/roc header
    rocheader.clear();
    rocheader.push_back(subbankheader.size() + elements.size() + 1);  // bank/roc size
    rocheader.push_back((fCurrentROC_ID << 16) | (0x10 << 8) | (1 & 0xff));
    // bank tag == ROC | bank type | event number

    // Add bank header, subbank header and element data to output buffer
    buffer.insert(buffer.end(), rocheader.begin(), rocheader.end());
    buffer.insert(buffer.end(), subbankheader.begin(), subbankheader.end());
    buffer.insert(buffer.end(), elements.begin(), elements.end());
  }
}

void QwComptonPhotonDetector::FillHistograms ( )

virtual

Fill the histograms with data

Implements VQwSubsystem.

Definition at line 783 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

{
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].FillHistograms();
  for (size_t i = 0; i < fMultiTDC_Channel.size(); i++)
    fMultiTDC_Channel[i].FillHistograms();
  for (size_t i = 0; i < fMultiQDC_Channel.size(); i++)
    fMultiQDC_Channel[i].FillHistograms();
}

void QwComptonPhotonDetector::FillTree ( )

virtual

Fill the tree with data

Reimplemented from VQwSubsystem.

Definition at line 816 of file QwComptonPhotonDetector.cc.

References fSamplingADC, fTree, and fTree_fNEvents.

{
  for (size_t i = 0; i < fSamplingADC.size(); i++) {
    fTree_fNEvents = fSamplingADC[i].GetNumberOfEvents();
    fTree->Fill();
  }
}

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

virtual

Fill the tree with data

Parameters

values

Implements VQwSubsystem.

Definition at line 844 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

{
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].FillTreeVector(values);
  for (size_t i = 0; i < fMultiTDC_Channel.size(); i++)
    fMultiTDC_Channel[i].FillTreeVector(values);
  for (size_t i = 0; i < fMultiQDC_Channel.size(); i++)
    fMultiQDC_Channel[i].FillTreeVector(values);
}

UInt_t QwComptonPhotonDetector::GetEventcutErrorFlag ( )

inlinevirtual

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

Implements VQwSubsystemParity.

Definition at line 87 of file QwComptonPhotonDetector.h.

{ return 0; };

QwSIS3320_Channel * QwComptonPhotonDetector::GetSIS3320Channel

(

const TString

name

)

Get the SIS3320 channel for this photon detector

Parameters

name	Name of the SIS3320 channel

Returns

Pointer to the SIS3320 channel

Definition at line 897 of file QwComptonPhotonDetector.cc.

References fSamplingADC.

{
  if (! fSamplingADC.empty()) {
    for (std::vector<QwSIS3320_Channel>::iterator adc = fSamplingADC.begin(); adc != fSamplingADC.end(); ++adc) {
      if (adc->GetElementName() == name) {
        return &(*adc);
      }
    }
  }
  return 0;
}

void QwComptonPhotonDetector::IncrementErrorCounters ( )

inlinevirtual

Increment the error counters.

Implements VQwSubsystemParity.

Definition at line 85 of file QwComptonPhotonDetector.h.

{};

Bool_t QwComptonPhotonDetector::IsGoodEvent

(

)

Check whether this is a good event

Returns

True if the event is good

Definition at line 495 of file QwComptonPhotonDetector.cc.

References fSamplingADC.

Referenced by SingleEventCuts().

{
  Bool_t test = kTRUE;

  for(size_t i=0;i<fSamplingADC.size();i++)
    test &= fSamplingADC[i].IsGoodEvent();
//  for(size_t i=0;i<fMultiTDC_Channel.size();i++)
//    test &= fMultiTDC_Channel[i].IsGoodEvent();
//  for(size_t i=0;i<fMultiQDC_Channel.size();i++)
//    test &= fMultiQDC_Channel[i].IsGoodEvent();

  if (!test) std::cerr << "This is not a good event!" << std::endl;
  return test;
}

Here is the caller graph for this function:

Int_t QwComptonPhotonDetector::LoadChannelMap ( TString  mapfile )

virtual

Load the channel map

Parameters

mapfile

Map file

Returns

Zero if successful

Implements VQwSubsystem.

Definition at line 46 of file QwComptonPhotonDetector.cc.

References QwLog::endl(), VQwSubsystem::fDetectorMaps, fMapping, fMultiQDC_Channel, fMultiQDC_Mapping, fMultiTDC_Channel, fMultiTDC_Mapping, fSamplingADC, fSamplingADC_Mapping, VQwSubsystem::GetSubbankIndex(), QwParameterFile::GetUInt(), kIntegratingADC, kIntegratingTDC, kSamplingADC, QwMessage, VQwSubsystem::RegisterROCNumber(), and VQwSubsystem::RegisterSubbank().

{
  Int_t subbank = -1; // subbank index
  Int_t current_roc_id = -1; // current ROC id
  Int_t current_bank_id = -1; // current bank id

  Int_t current_logical_accum = 0;
  std::vector<TString> name;
  std::vector<std::vector<TString> >  accums;
  std::vector<std::vector<Double_t> > weights;

  // Open the file
  QwParameterFile mapstr(mapfile.Data());
  fDetectorMaps.insert(mapstr.GetParamFileNameContents());

  while (mapstr.ReadNextLine()) {
    mapstr.TrimComment();      // Remove everything after a comment character
    mapstr.TrimWhitespace();   // Get rid of leading and trailing whitespace
    if (mapstr.LineIsEmpty())  continue;

    TString varname, varvalue;
    if (mapstr.HasVariablePair("=",varname,varvalue)) {
      // This is a declaration line.  Decode it.
      varname.ToLower();

      if (varname == "roc") {
        current_roc_id = QwParameterFile::GetUInt(varvalue);
        RegisterROCNumber(current_roc_id,0);

      } else if (varname == "bank") {
        current_bank_id = QwParameterFile::GetUInt(varvalue);
        RegisterSubbank(current_bank_id);
        subbank = GetSubbankIndex(current_roc_id,current_bank_id);

      } else if (varname == "begin") {
        // Add new logical accumulator
        name.resize(current_logical_accum+1);
        accums.resize(current_logical_accum+1);
        weights.resize(current_logical_accum+1);
        // Read to end of this block
        while(mapstr.ReadNextLine()) {
          mapstr.TrimComment();      // Remove everything after a comment character
          mapstr.TrimWhitespace();   // Get rid of leading and trailing whitespace
          if (mapstr.LineIsEmpty())  continue;

          // Loop over block
          if (mapstr.HasVariablePair("=",varname,varvalue)) {
            varname.ToLower();
            if (varname == "end") {
              current_logical_accum++;
              break;
            } else if (varname == "name") {
              name[current_logical_accum] = varvalue;
            }
          } else {
            accums[current_logical_accum].push_back(mapstr.GetTypedNextToken<TString>());
            weights[current_logical_accum].push_back(mapstr.GetTypedNextToken<Double_t>());
          }
        }
      }

    } else {
      //  Break this line into tokens to process it.
      TString modtype = mapstr.GetTypedNextToken<TString>();
      UInt_t modnum   = mapstr.GetTypedNextToken<UInt_t>();
      UInt_t channum  = mapstr.GetTypedNextToken<UInt_t>();
      TString dettype = mapstr.GetTypedNextToken<TString>();
      TString name    = mapstr.GetTypedNextToken<TString>();
      modtype.ToUpper();
      dettype.ToUpper();

      //  Push a new record into the element array
      if (modtype == "SIS3320") {
        // Register data channel type
        fMapping[subbank] = kSamplingADC;
        // Increase module mapping size
        if (modnum >= fSamplingADC_Mapping[subbank].size())
          fSamplingADC_Mapping[subbank].resize(modnum+1);
        // Increase channel mapping size
        if (channum >= fSamplingADC_Mapping[subbank].at(modnum).size())
          fSamplingADC_Mapping[subbank].at(modnum).resize(channum+1,-1);
        // Search whether this module/channel is already registered
        for (SamplingADC_Mapping_t::iterator iter = fSamplingADC_Mapping.begin();
            iter != fSamplingADC_Mapping.end(); iter++) {
          if (iter->first != subbank &&
              iter->second.size() > modnum &&
              iter->second.at(modnum).size() > channum &&
              iter->second.at(modnum).at(channum) >= 0) {
            QwMessage << "Connecting SIS3320 " << name << "/" << dettype
                      << std::hex
                      << " in ROC 0x" << current_roc_id << ", bank 0x" << current_bank_id
                      << std::dec
                      << " to existing SIS3320." << QwLog::endl;
            fSamplingADC_Mapping[subbank].at(modnum).at(channum) =
                iter->second.at(modnum).at(channum);
          }
        }
        // Register new SIS3320
        if (fSamplingADC_Mapping[subbank].at(modnum).at(channum) < 0) {
          QwMessage << "Registering SIS3320 " << name << "/" << dettype
                    << std::hex
                    << " in ROC 0x" << current_roc_id << ", bank 0x" << current_bank_id
                    << std::dec
                    << " at mod " << modnum << ", chan " << channum
                    << QwLog::endl;
          UInt_t index = fSamplingADC.size();
          fSamplingADC_Mapping[subbank].at(modnum).at(channum) = index;
          fSamplingADC.push_back(QwSIS3320_Channel(channum, name));
          fSamplingADC.at(index).SetNumberOfAccumulators(6);
          fSamplingADC.at(index).InitializeChannel(channum, name);
        }

      } else if (modtype == "V792") {
        // Register data channel type
        fMapping[subbank] = kIntegratingADC;
        // Add to mapping
        if (modnum >= fMultiQDC_Mapping[subbank].size())
          fMultiQDC_Mapping[subbank].resize(modnum+1);
        if (channum >= fMultiQDC_Mapping[subbank].at(modnum).size())
          fMultiQDC_Mapping[subbank].at(modnum).resize(channum+1,-1);
        // Add scaler channel
        if (fMultiQDC_Mapping[subbank].at(modnum).at(channum) < 0) {
          QwMessage << "Registering V792 " << name
                    << std::hex
                    << " in ROC 0x" << current_roc_id << ", bank 0x" << current_bank_id
                    << std::dec
                    << " at mod " << modnum << ", chan " << channum
                    << QwLog::endl;
          UInt_t index = fMultiQDC_Channel.size();
          fMultiQDC_Mapping[subbank].at(modnum).at(channum) = index;
          fMultiQDC_Channel.push_back(QwPMT_Channel(name));
          fMultiQDC_Channel.at(index).SetModule(modnum);
          fMultiQDC_Channel.at(index).SetSubbankID(current_bank_id);
        }

      } else if (modtype == "V775") {
        // Register data channel type
        fMapping[subbank] = kIntegratingTDC;
        // Add to mapping
        if (modnum >= fMultiTDC_Mapping[subbank].size())
          fMultiTDC_Mapping[subbank].resize(modnum+1);
        if (channum >= fMultiTDC_Mapping[subbank].at(modnum).size())
          fMultiTDC_Mapping[subbank].at(modnum).resize(channum+1,-1);
        // Add scaler channel
        if (fMultiTDC_Mapping[subbank].at(modnum).at(channum) < 0) {
          QwMessage << "Registering V775 " << name
                    << std::hex
                    << " in ROC 0x" << current_roc_id << ", bank 0x" << current_bank_id
                    << std::dec
                    << " at mod " << modnum << ", chan " << channum
                    << QwLog::endl;
          UInt_t index = fMultiTDC_Channel.size();
          fMultiTDC_Mapping[subbank].at(modnum).at(channum) = index;
          fMultiTDC_Channel.push_back(QwPMT_Channel(name));
          fMultiTDC_Channel.at(index).SetModule(modnum);
          fMultiTDC_Channel.at(index).SetSubbankID(current_bank_id);
        }

      } // end of switch by modtype

    } // end of if for token line
  } // end of while over parameter file

  // Now create the logical accumulators
  for (size_t logical = 0; logical < name.size(); logical++) {
    for (size_t adc = 0; adc < fSamplingADC.size(); adc++) {
      fSamplingADC[adc].AddLogicalAccumulator(name[logical],accums[logical],weights[logical]);
    }
  }

  return 0;
}

Here is the call graph for this function:

Int_t QwComptonPhotonDetector::LoadEventCuts

(

TString &

filename

)

Load the event cuts

Parameters

filename

Event cuts file

Returns

Zero if successful

Definition at line 224 of file QwComptonPhotonDetector.cc.

{
  return 0;
}

Int_t QwComptonPhotonDetector::LoadEventCuts ( TString  filename )

inlinevirtual

Load the event cuts file.

Implements VQwSubsystemParity.

Definition at line 83 of file QwComptonPhotonDetector.h.

{ return 0; };

Int_t QwComptonPhotonDetector::LoadInputParameters ( TString  pedestalfile )

virtual

Load the input parameters

Parameters

pedestalfile

Input parameters file

Returns

Zero if successful

Implements VQwSubsystem.

Definition at line 234 of file QwComptonPhotonDetector.cc.

References QwLog::endl(), VQwSubsystem::fDetectorMaps, fSamplingADC, QwMessage, and QwVerbose.

{
  // Open the file
  QwParameterFile mapstr(pedestalfile.Data());
  fDetectorMaps.insert(mapstr.GetParamFileNameContents());
  while (mapstr.ReadNextLine()) {
    mapstr.TrimComment();    // Remove everything after a comment character
    mapstr.TrimWhitespace(); // Get rid of leading and trailing spaces

    TString varname = "";
    Double_t varped = 0.0;
    Double_t varcal = 1.0;

    if (mapstr.LineIsEmpty())  continue;
    else {
      varname = mapstr.GetTypedNextToken<TString>();  // name of the channel
      varname.ToLower();
      varname.Remove(TString::kBoth,' ');
      varped = mapstr.GetTypedNextToken<Double_t>(); // value of the pedestal
      varcal = mapstr.GetTypedNextToken<Double_t>(); // value of the calibration factor
      QwVerbose << varname << ": " << QwLog::endl;
    }

    Bool_t found = kFALSE;
    for (size_t i = 0; i < fSamplingADC.size() && !found; i++) {
      TString localname = fSamplingADC[i].GetElementName();
      localname.ToLower();
      if (localname == varname) {
        QwMessage << "Setting pedestal and calibration of sampling ADC " << localname
                  << " to " << varped << " and " << varcal << "." << QwLog::endl;
        fSamplingADC[i].SetPedestal(varped);
        fSamplingADC[i].SetCalibrationFactor(varcal);
        found = kTRUE;
      }
    } // end of loop over all sampling ADC channels

  } // end of loop reading all lines of the pedestal file

  return 0;
}

Here is the call graph for this function:

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

virtual

Addition-assignment operator

Parameters

value

Right-hand side

Returns

Left-hand side

Implements VQwSubsystemParity.

Definition at line 562 of file QwComptonPhotonDetector.cc.

References Compare(), and fSamplingADC.

{
  if (Compare(value)) {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*>(value);
    for (size_t i = 0; i < input->fSamplingADC.size(); i++)
      this->fSamplingADC[i] += input->fSamplingADC[i];
  }
  return *this;
}

Here is the call graph for this function:

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

virtual

Subtraction-assignment operator

Parameters

value

Right-hand side

Returns

Left-hand side

Implements VQwSubsystemParity.

Definition at line 577 of file QwComptonPhotonDetector.cc.

References Compare(), and fSamplingADC.

{
  if (Compare(value)) {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*> (value);
    for (size_t i = 0; i < input->fSamplingADC.size(); i++)
      this->fSamplingADC[i] -= input->fSamplingADC[i];
  }
  return *this;
}

Here is the call graph for this function:

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

virtual

Assignment operator

Parameters

value

Right-hand side

Returns

Left-hand side

Implements VQwSubsystemParity.

Definition at line 534 of file QwComptonPhotonDetector.cc.

References CompareADC(), CompareQDC(), CompareTDC(), fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

{
  if (CompareADC(value)) {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*>(value);
    for (size_t i = 0; i < input->fSamplingADC.size(); i++)
      this->fSamplingADC[i] = input->fSamplingADC[i];
  }

  if (CompareTDC(value)) {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*>(value);
    for (size_t i = 0; i < input->fMultiTDC_Channel.size();i++)
      this->fMultiTDC_Channel[i] = input->fMultiTDC_Channel[i];
  }

  if (CompareQDC(value)) {
    QwComptonPhotonDetector* input = dynamic_cast<QwComptonPhotonDetector*>(value);
    for (size_t i = 0; i < input->fMultiQDC_Channel.size();i++)
      this->fMultiQDC_Channel[i] = input->fMultiQDC_Channel[i];
  }

  return *this;
}

Here is the call graph for this function:

void QwComptonPhotonDetector::PrintErrorCounters ( ) const

inlinevirtual

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

Implements VQwSubsystemParity.

Definition at line 86 of file QwComptonPhotonDetector.h.

{};

void QwComptonPhotonDetector::PrintInfo ( ) const

virtual

Print some debugging output for the subcomponents

Reimplemented from VQwSubsystem.

Definition at line 870 of file QwComptonPhotonDetector.cc.

References QwLog::endl(), fMultiQDC_Channel, fMultiTDC_Channel, fSamplingADC, VQwSubsystem::PrintInfo(), and QwOut.

{
  VQwSubsystemParity::PrintInfo();

  QwOut << " there are " << fSamplingADC.size() << " sampling ADCs" << QwLog::endl;
  QwOut << " there are " << fMultiTDC_Channel.size() << " integrating TDCs" << QwLog::endl;
  QwOut << " there are " << fMultiQDC_Channel.size() << " integrating ADCs" << QwLog::endl;

  for (size_t i = 0; i < fSamplingADC.size(); i++) {
    QwOut << " sampling ADC " << i << ":";
    fSamplingADC[i].PrintInfo();
  }
  for (size_t i = 0; i < fMultiTDC_Channel.size(); i++) {
    QwOut << " integrating TDC " << i << ":";
    fMultiTDC_Channel[i].PrintInfo();
  }
  for (size_t i = 0; i < fMultiQDC_Channel.size(); i++) {
    QwOut << " integrating ADC " << i << ":";
    fMultiQDC_Channel[i].PrintInfo();
  }
}

Here is the call graph for this function:

void QwComptonPhotonDetector::PrintValue ( ) const

virtual

Print the value for the subcomponents

Reimplemented from VQwSubsystemParity.

Definition at line 857 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

{
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].PrintValue();
  for (size_t i = 0; i < fMultiTDC_Channel.size(); i++)
    fMultiTDC_Channel[i].PrintValue();
  for (size_t i = 0; i < fMultiQDC_Channel.size(); i++)
    fMultiQDC_Channel[i].PrintValue();
}

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

virtual

Process the configuration buffer for this subsystem

Parameters

roc_id	ROC ID
bank_id	Subbank ID
buffer	Buffer to read from
num_words	Number of words left in buffer

Returns

Number of words read

Implements VQwSubsystem.

Definition at line 485 of file QwComptonPhotonDetector.cc.

{
  return 0;
}

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

virtual

Process the event buffer for this subsystem

Parameters

roc_id	ROC ID
bank_id	Subbank ID
buffer	Buffer to read from
num_words	Number of words left in buffer

Returns

Number of words read

Implements VQwSubsystem.

Definition at line 342 of file QwComptonPhotonDetector.cc.

References MQwV775TDC::DecodeTDCWord(), QwLog::endl(), fMapping, fMultiQDC_Channel, fMultiQDC_Mapping, fMultiTDC_Channel, fMultiTDC_Mapping, fSamplingADC, fSamplingADC_Mapping, VQwSubsystem::GetSubbankIndex(), MQwV775TDC::GetTDCChannelNumber(), MQwV775TDC::GetTDCData(), MQwV775TDC::IsValidDataword(), kIntegratingADC, kIntegratingTDC, kSamplingADC, kUnknown, QwDebug, and QwError.

{
  UInt_t words_read = 0;

  // Get the subbank index (or -1 when no match)
  Int_t subbank = GetSubbankIndex(roc_id, bank_id);

  if (subbank >= 0 && num_words > 0) {

    //  We want to process this ROC.  Begin looping through the data.

    switch (fMapping[subbank]) {

      // Sampling ADCs
      case kSamplingADC:
      {
        for (size_t modnum = 0; modnum < fSamplingADC_Mapping[subbank].size(); modnum++) {
          for (size_t channum = 0; channum < fSamplingADC_Mapping[subbank].at(modnum).size(); channum++) {
            Int_t index = fSamplingADC_Mapping[subbank].at(modnum).at(channum);
            if (index >= 0) {
              words_read += fSamplingADC[index].ProcessEvBuffer(&(buffer[words_read]), num_words - words_read);
            }
          }
        }
        break;
      }

      // Integrating ADCs and TDCs
      case kIntegratingADC:
      case kIntegratingTDC:
      {
        // Read header word
        UInt_t header = buffer[words_read++];
        if ((header & 0xbad00dc0) == 0xbad00dc0)
          break; // Bad ADC/TDC

        // Parse number of events
        //UInt_t num_events = 0;
        //if ((header & 0xda000dc0) == 0xda000dc0)
        //  num_events = (buffer[0] & 0x00FF0000) >> 16;
        // TODO Multihit functionality

        // Loop over buffered events
        for (size_t i = words_read; i < num_words; i++) {

          // Decode this word as a V775 TDC / V792 QDC word
          DecodeTDCWord(buffer[i]);

          // Is this a valid V775 TDC / V792 QDC data word?
          if (! IsValidDataword()) continue;

          // Is there an QDC channel registered for this slot and channel?
          for (size_t modnum = 0; modnum < fMultiQDC_Mapping[subbank].size(); modnum++) {
            UInt_t channum = GetTDCChannelNumber();
            if (channum < fMultiQDC_Mapping[subbank].at(modnum).size()) {
              Int_t index = fMultiQDC_Mapping[subbank].at(modnum).at(channum);
              if (index >= 0) {
                fMultiQDC_Channel.at(index).SetValue(GetTDCData());
                QwDebug << "QDC " << std::hex << subbank << " "
                        << modnum << "," << channum << ": " << std::dec
                        << GetTDCData() << QwLog::endl;
              }
            }
          }

          // Is there an TDC channel registered for this slot and channel?
          for (size_t modnum = 0; modnum < fMultiTDC_Mapping[subbank].size(); modnum++) {
            UInt_t channum = GetTDCChannelNumber();
            if (channum < fMultiTDC_Mapping[subbank].at(channum).size()) {
              Int_t index = fMultiTDC_Mapping[subbank].at(modnum).at(channum);
              if (index >= 0) {
                fMultiTDC_Channel.at(index).SetValue(GetTDCData());
                QwDebug << "TDC " << std::hex << subbank << " "
                        << modnum << "," << channum << ": " << std::dec
                        << GetTDCData() << QwLog::endl;
              }
            }
          }

        } // end of loop over buffered events

        words_read = num_words;
        break;
      }

      // Unknown data channel type
      case kUnknown:
      default:
      {
        QwError << "QwComptonPhotonDetector: Unknown data channel type" << QwLog::endl;
        break;
      }
    }

    // Check for leftover words
    if (num_words != words_read) {
      QwError << "QwComptonPhotonDetector: There were "
              << num_words - words_read
              << " leftover words after decoding everything we recognize"
              << std::hex
              << " in ROC " << roc_id << ", bank " << bank_id << "."
              << std::dec
              << QwLog::endl;
    }
  }

  return words_read;
}

Here is the call graph for this function:

void QwComptonPhotonDetector::ProcessEvent ( )

virtual

Process this event

Implements VQwSubsystem.

Definition at line 466 of file QwComptonPhotonDetector.cc.

References fMultiQDC_Channel, fMultiTDC_Channel, and fSamplingADC.

{
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].ProcessEvent();
  for (size_t i = 0; i < fMultiTDC_Channel.size(); i++)
    fMultiTDC_Channel[i].ProcessEvent();
  for (size_t i = 0; i < fMultiQDC_Channel.size(); i++)
    fMultiQDC_Channel[i].ProcessEvent();
}

void QwComptonPhotonDetector::ProcessOptions ( QwOptions &  options )

virtual

Process the command line options.

Reimplemented from VQwSubsystem.

Definition at line 36 of file QwComptonPhotonDetector.cc.

{
  // Handle command line options
}

void QwComptonPhotonDetector::RandomizeEventData

(

int

helicity = 0

)

Randomize this event with the given helicity

Parameters

helicity

Helicity of this event (default is zero)

Definition at line 280 of file QwComptonPhotonDetector.cc.

References fSamplingADC.

{
  // Randomize all gQwSIS3320 buffers
  for (size_t i = 0; i < fSamplingADC.size(); i++)
    fSamplingADC[i].RandomizeEventData(helicity);

  // Randomize all MQwV775TDC buffers
  //for (size_t i = 0; i < fMultiTDC_Channel.size(); i++)
  //  fMultiTDC_Channel[i].RandomizeEventData(helicity);

  // Randomize all MQwV792ADC buffers
  //for (size_t i = 0; i < fMultiQDC_Channel.size(); i++)
  //  fMultiQDC_Channel[i].RandomizeEventData(helicity);
}

void QwComptonPhotonDetector::Ratio ( VQwSubsystem *  numer, VQwSubsystem *  denom  )

virtual

Determine the ratio of two photon detectors

Parameters

numer	Numerator
denom	Denominator

Implements VQwSubsystemParity.

Definition at line 618 of file QwComptonPhotonDetector.cc.

References Compare(), and fSamplingADC.

{
  if (Compare(numer) && Compare(denom)) {
    QwComptonPhotonDetector* innumer = dynamic_cast<QwComptonPhotonDetector*> (numer);
    QwComptonPhotonDetector* indenom = dynamic_cast<QwComptonPhotonDetector*> (denom);
    for (size_t i = 0; i < innumer->fSamplingADC.size(); i++)
      this->fSamplingADC[i].Ratio(innumer->fSamplingADC[i], indenom->fSamplingADC[i]);
  }
}

Here is the call graph for this function:

void QwComptonPhotonDetector::Scale ( Double_t  factor )

virtual

Scale the photon detector

Parameters

factor

Scale factor

Implements VQwSubsystemParity.

Definition at line 632 of file QwComptonPhotonDetector.cc.

References fSamplingADC.

Referenced by CalculateRunningAverage().

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

Here is the caller graph for this function:

Bool_t QwComptonPhotonDetector::SingleEventCuts

(

)

Process the single event cuts

Returns

Definition at line 456 of file QwComptonPhotonDetector.cc.

References QwLog::endl(), IsGoodEvent(), and QwDebug.

{
  QwDebug << "QwComptonPhotonDetector::SingleEventCuts()" << QwLog::endl;
  return IsGoodEvent();
}

Here is the call graph for this function:

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

virtual

Summation

Parameters

value1	First value
value2	Second value

Implements VQwSubsystemParity.

Definition at line 592 of file QwComptonPhotonDetector.cc.

References Compare().

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

Here is the call graph for this function:

void QwComptonPhotonDetector::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 90 of file QwComptonPhotonDetector.h.

                                                      {
    };

Field Documentation

Int_t QwComptonPhotonDetector::fGoodEventCount

private

Definition at line 156 of file QwComptonPhotonDetector.h.

Referenced by AccumulateRunningSum(), CalculateRunningAverage(), and ClearEventData().

std::map< Int_t, ChannelType_t > QwComptonPhotonDetector::fMapping

protected

Definition at line 134 of file QwComptonPhotonDetector.h.

Referenced by LoadChannelMap(), and ProcessEvBuffer().

std::vector< QwPMT_Channel > QwComptonPhotonDetector::fMultiQDC_Channel

protected

Definition at line 144 of file QwComptonPhotonDetector.h.

Referenced by ClearEventData(), Compare(), CompareQDC(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), LoadChannelMap(), operator=(), PrintInfo(), PrintValue(), ProcessEvBuffer(), and ProcessEvent().

std::vector< std::vector< QwPMT_Channel > > QwComptonPhotonDetector::fMultiQDC_Events

protected

Definition at line 145 of file QwComptonPhotonDetector.h.

IntegratingADC_Mapping_t QwComptonPhotonDetector::fMultiQDC_Mapping

protected

Definition at line 143 of file QwComptonPhotonDetector.h.

Referenced by LoadChannelMap(), and ProcessEvBuffer().

std::vector< QwPMT_Channel > QwComptonPhotonDetector::fMultiTDC_Channel

protected

Definition at line 150 of file QwComptonPhotonDetector.h.

Referenced by ClearEventData(), Compare(), CompareTDC(), ConstructBranchAndVector(), ConstructHistograms(), FillHistograms(), FillTreeVector(), LoadChannelMap(), operator=(), PrintInfo(), PrintValue(), ProcessEvBuffer(), and ProcessEvent().

std::vector< std::vector< QwPMT_Channel > > QwComptonPhotonDetector::fMultiTDC_Events

protected

Definition at line 151 of file QwComptonPhotonDetector.h.

IntegratingTDC_Mapping_t QwComptonPhotonDetector::fMultiTDC_Mapping

protected

Definition at line 149 of file QwComptonPhotonDetector.h.

Referenced by LoadChannelMap(), and ProcessEvBuffer().

std::vector< QwSIS3320_Channel > QwComptonPhotonDetector::fSamplingADC

protected

Definition at line 139 of file QwComptonPhotonDetector.h.

Referenced by ClearEventData(), Compare(), CompareADC(), ConstructBranchAndVector(), ConstructHistograms(), EncodeEventData(), FillHistograms(), FillTree(), FillTreeVector(), GetSIS3320Channel(), IsGoodEvent(), LoadChannelMap(), LoadInputParameters(), operator+=(), operator-=(), operator=(), PrintInfo(), PrintValue(), ProcessEvBuffer(), ProcessEvent(), RandomizeEventData(), Ratio(), and Scale().

SamplingADC_Mapping_t QwComptonPhotonDetector::fSamplingADC_Mapping

protected

Definition at line 138 of file QwComptonPhotonDetector.h.

Referenced by LoadChannelMap(), and ProcessEvBuffer().

TTree* QwComptonPhotonDetector::fTree

protected

Expert tree.

Definition at line 128 of file QwComptonPhotonDetector.h.

Referenced by ConstructTree(), DeleteTree(), and FillTree().

Int_t QwComptonPhotonDetector::fTree_fNEvents

protected

Expert tree fields.

Definition at line 130 of file QwComptonPhotonDetector.h.

Referenced by ConstructTree(), and FillTree().

const Bool_t QwComptonPhotonDetector::kDebug = kTRUE

staticprivate

Definition at line 158 of file QwComptonPhotonDetector.h.

---

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

• QwComptonPhotonDetector.h
• QwComptonPhotonDetector.cc