#include <QwMainCerenkovDetector.h>

Inheritance diagram for QwMainCerenkovDetector:

Collaboration diagram for QwMainCerenkovDetector:

Public Member Functions
	QwMainCerenkovDetector (const TString &name)
	Constructor with name. More...

	QwMainCerenkovDetector (const QwMainCerenkovDetector &source)
	Copy constructor. More...

virtual	~QwMainCerenkovDetector ()
	Virtual destructor. More...

void	ProcessOptions (QwOptions &options)

Int_t	LoadChannelMap (TString mapfile)
	Mandatory map file definition. More...

Int_t	LoadInputParameters (TString pedestalfile)
	Mandatory parameter file definition. 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...

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)
	TODO: The non-event-type-aware ProcessEvBuffer routine should be replaced with the event-type-aware version. More...

void	ClearEventData ()

Bool_t	IsGoodEvent ()

void	ProcessEvent ()

void	ExchangeProcessedData ()

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

Bool_t	PublishInternalValues () const

void	SetRandomEventParameters (Double_t mean, Double_t sigma)

void	SetRandomEventAsymmetry (Double_t asymmetry)

void	RandomizeEventData (int helicity=0, Double_t time=0.0)

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

void	ConstructHistograms (TDirectory *folder)
	Construct the histograms for this subsystem in a folder. More...

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

void	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
	Construct the branch and tree vector. More...

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

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

void	FillErrDB (QwParityDB *db, TString datatype)

const QwIntegrationPMT *	GetChannel (const TString name) const

Bool_t	Compare (VQwSubsystem *source)

void	Blind (const QwBlinder *blinder)
	Blind the asymmetry. More...

void	Blind (const QwBlinder blinder, const VQwSubsystemParity subsys)
	Blind the difference using the yield. 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 numer, VQwSubsystem denom)

void	Scale (Double_t factor)

void	Normalize (VQwDataElement *denom)

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

const QwIntegrationPMT *	GetIntegrationPMT (const TString name) const

const QwCombinedPMT *	GetCombinedPMT (const TString name) const

void	DoNormalization (Double_t factor=1.0)

Bool_t	ApplyHWChecks ()

void	PrintValue () const
	Print values of all channels. More...

void	WritePromptSummary (QwPromptSummary *ps, TString type)

void	PrintInfo () const
	Print some information about the subsystem. More...

void	PrintDetectorID () const

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

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, QwMainCerenkovDetector >
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)
	Define options function. More...

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, QwMainCerenkovDetector >
static VQwSubsystem *	Create (const std::string &name)
	Object creation. More...

static QwMainCerenkovDetector *	Cast (QwMainCerenkovDetector *type)
	Object dynamic cast. More...

Protected Member Functions
EQwPMTInstrumentType	GetDetectorTypeID (TString name)

Int_t	GetDetectorIndex (EQwPMTInstrumentType TypeID, TString name)

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

Protected Attributes
Bool_t	fDEBUG

std::vector< QwIntegrationPMT >	fIntegrationPMT

std::vector< QwCombinedPMT >	fCombinedPMT

std::vector < QwMainCerenkovDetectorID >	fMainDetID

QwBeamCharge	fTargetCharge

QwBeamPosition	fTargetX

QwBeamPosition	fTargetY

QwBeamAngle	fTargetXprime

QwBeamAngle	fTargetYprime

QwBeamEnergy	fTargetEnergy

Bool_t	bIsExchangedDataValid

Bool_t	bNormalization

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
	QwMainCerenkovDetector ()
	Private default constructor (not implemented, will throw linker error on use) More...

Private Attributes
Int_t	fMainDetErrorCount

Static Private Attributes
static const Bool_t	bDEBUG =kFALSE

Friends
class	QwCombinedPMT

Detailed Description

Definition at line 28 of file QwMainCerenkovDetector.h.

Constructor & Destructor Documentation

QwMainCerenkovDetector::QwMainCerenkovDetector ( )

private

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

QwMainCerenkovDetector::QwMainCerenkovDetector ( const TString & name )

inline

Constructor with name.

Definition at line 44 of file QwMainCerenkovDetector.h.

References fTargetCharge, fTargetEnergy, fTargetX, fTargetXprime, fTargetY, and fTargetYprime.

   : VQwSubsystem(name),VQwSubsystemParity(name),bNormalization(kFALSE)
   {
     fTargetCharge.InitializeChannel("q_targ","derived");
     fTargetX.InitializeChannel("x_targ","derived");
     fTargetY.InitializeChannel("y_targ","derived");
     fTargetXprime.InitializeChannel("xp_targ","derived");
     fTargetYprime.InitializeChannel("yp_targ","derived");
     fTargetEnergy.InitializeChannel("e_targ","derived");
   };

QwMainCerenkovDetector::QwMainCerenkovDetector ( const QwMainCerenkovDetector & source )

inline

Copy constructor.

Definition at line 55 of file QwMainCerenkovDetector.h.

   : VQwSubsystem(source),VQwSubsystemParity(source),
     fIntegrationPMT(source.fIntegrationPMT),
     fCombinedPMT(source.fCombinedPMT)
   { }

virtual QwMainCerenkovDetector::~QwMainCerenkovDetector ( )

inlinevirtual

Virtual destructor.

Definition at line 61 of file QwMainCerenkovDetector.h.

61 { };

Member Function Documentation

void QwMainCerenkovDetector::AccumulateRunningSum ( VQwSubsystem * value )

virtual

Update the running sums for devices.

Implements VQwSubsystemParity.

Definition at line 1184 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, and fIntegrationPMT.

 {
   if (Compare(value1)) {
     QwMainCerenkovDetector* value = dynamic_cast<QwMainCerenkovDetector*>(value1);
 
     for (size_t i = 0; i < fIntegrationPMT.size(); i++)
       fIntegrationPMT[i].AccumulateRunningSum(value->fIntegrationPMT[i]);
     for (size_t i = 0; i < fCombinedPMT.size(); i++)
       fCombinedPMT[i].AccumulateRunningSum(value->fCombinedPMT[i]);
   }
 }

Here is the call graph for this function:

Bool_t QwMainCerenkovDetector::ApplyHWChecks ( )

inline

Definition at line 150 of file QwMainCerenkovDetector.h.

References fIntegrationPMT.

                         {//Check for harware errors in the devices
     Bool_t status = kTRUE;
     for (size_t i=0; i<fIntegrationPMT.size(); i++){
       status &= fIntegrationPMT.at(i).ApplyHWChecks();
     }
     return status;
   };

Bool_t QwMainCerenkovDetector::ApplySingleEventCuts ( )

virtual

Apply the single event cuts.

Implements VQwSubsystemParity.

Definition at line 712 of file QwMainCerenkovDetector.cc.

References bDEBUG, fCombinedPMT, fIntegrationPMT, and fMainDetErrorCount.

 {
   Bool_t status=kTRUE;
   for(size_t i=0;i<fIntegrationPMT.size();i++){
     status &= fIntegrationPMT[i].ApplySingleEventCuts();
     if(!status && bDEBUG) std::cout<<"******* QwMainCerenkovDetector::SingleEventCuts()->IntegrationPMT[ "<<i<<" , "<<fIntegrationPMT[i].GetElementName()<<" ] ******\n"; 
   }
   for(size_t i=0;i<fCombinedPMT.size();i++){
     status &= fCombinedPMT[i].ApplySingleEventCuts();
     if(!status && bDEBUG) std::cout<<"******* QwMainCerenkovDetector::SingleEventCuts()->CombinedPMT[ "<<i<<" , "<<fCombinedPMT[i].GetElementName()<<" ] ******\n"; 
   }
 
 
   if (!status) 
    fMainDetErrorCount++;//failed  event counter for QwMainCerenkovDetector
 
   return status;
 
 
 }

void QwMainCerenkovDetector::Blind ( const QwBlinder * blinder )

virtual

Blind the asymmetry.

Blind the asymmetry

Parameters

blinder Blinder

Reimplemented from VQwSubsystemParity.

Definition at line 1212 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

Referenced by Blind().

 {
   for (size_t i = 0; i < fIntegrationPMT.size(); i++)
     fIntegrationPMT[i].Blind(blinder);
   for (size_t i = 0; i < fCombinedPMT.size(); i++)
     fCombinedPMT[i].Blind(blinder);
 }

Here is the caller graph for this function:

void QwMainCerenkovDetector::Blind	(	const QwBlinder *	blinder,
		const VQwSubsystemParity *	subsys
	)

virtual

Blind the difference using the yield.

Blind the difference using the yield

Parameters

blinder	Blinder
subsys	Subsystem

Todo:: TODO (wdc) At some point we should introduce const-correctness in the Compare() routine to ensure nothing funny happens. This const_casting is just an ugly stop-gap measure.

Reimplemented from VQwSubsystemParity.

Definition at line 1225 of file QwMainCerenkovDetector.cc.

References Blind(), Compare(), fCombinedPMT, and fIntegrationPMT.

 {
   /// \todo TODO (wdc) At some point we should introduce const-correctness in
   /// the Compare() routine to ensure nothing funny happens.  This const_casting
   /// is just an ugly stop-gap measure.
   if (Compare(const_cast<VQwSubsystemParity*>(subsys))) {
 
     const QwMainCerenkovDetector* yield = dynamic_cast<const QwMainCerenkovDetector*>(subsys);
     if (yield == 0) return;
 
     for (size_t i = 0; i < fIntegrationPMT.size(); i++)
       fIntegrationPMT[i].Blind(blinder, yield->fIntegrationPMT[i]);
     for (size_t i = 0; i < fCombinedPMT.size(); i++)
       fCombinedPMT[i].Blind(blinder, yield->fCombinedPMT[i]);
   }
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::CalculateRunningAverage ( )

virtual

Calculate the average for all good events.

Implements VQwSubsystemParity.

Definition at line 1173 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   for (size_t i=0;i<fIntegrationPMT.size();i++)
     fIntegrationPMT[i].CalculateRunningAverage();
 
   for (size_t i=0;i<fCombinedPMT.size();i++)
     fCombinedPMT[i].CalculateRunningAverage();
 
   return;
 }

void QwMainCerenkovDetector::ClearEventData ( )

virtual

Implements VQwSubsystem.

Definition at line 571 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   for (size_t i=0;i<fIntegrationPMT.size();i++){
     fIntegrationPMT[i].ClearEventData();
   }
   for (size_t i=0;i<fCombinedPMT.size();i++)
     fCombinedPMT[i].ClearEventData();
 
   return;
 }

Bool_t QwMainCerenkovDetector::Compare ( VQwSubsystem * source )

Definition at line 1036 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

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

 {
   //  std::cout<<" Here in QwLumi::Compare \n";
 
   Bool_t res=kTRUE;
   if (typeid(*value)!=typeid(*this))
     {
       res=kFALSE;
       //      std::cout<<" types are not ok \n";
       //      std::cout<<" this is bypassed just for now but should be fixed eventually \n";
     }
   else
     {
       QwMainCerenkovDetector* input = dynamic_cast<QwMainCerenkovDetector*>(value);
       if (input->fIntegrationPMT.size()!=fIntegrationPMT.size() ||
           input->fCombinedPMT.size()!=fCombinedPMT.size() )
         {
           res=kFALSE;
           //    std::cout<<" not the same number of channels \n";
         }
     }
   return res;
 }

Here is the caller graph for this function:

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

virtual

Construct the branch and tree vector.

Implements VQwSubsystem.

Definition at line 982 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

Referenced by ConstructBranch().

 {
   for (size_t i=0;i<fIntegrationPMT.size();i++)
     fIntegrationPMT[i].ConstructBranch(tree, prefix);
 
   for (size_t i=0;i<fCombinedPMT.size();i++)
     fCombinedPMT[i].ConstructBranch(tree, prefix);
 
   return;
 }

Here is the caller graph for this function:

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

virtual

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

Implements VQwSubsystem.

Definition at line 993 of file QwMainCerenkovDetector.cc.

References ConstructBranch(), fCombinedPMT, QwParameterFile::FileHasModuleHeader(), fIntegrationPMT, QwParameterFile::ReadUntilNextModule(), and QwParameterFile::RewindToFileStart().

 {
   TString tmp;
   QwParameterFile* nextmodule;
   trim_file.RewindToFileStart();
   tmp="QwIntegrationPMT";
   trim_file.RewindToFileStart();
   if (trim_file.FileHasModuleHeader(tmp)){
     nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
     for (size_t i=0;i<fIntegrationPMT.size();i++)
 
       fIntegrationPMT[i].ConstructBranch(tree, prefix, *nextmodule);
   }
 
   tmp="QwCombinedPMT";
   trim_file.RewindToFileStart();
    if (trim_file.FileHasModuleHeader(tmp)){
      nextmodule=trim_file.ReadUntilNextModule();//This section contains sub modules and or channels to be included in the tree
      for (size_t i=0;i<fCombinedPMT.size();i++)
        fCombinedPMT[i].ConstructBranch(tree, prefix, *nextmodule );
    }
 
   return;
 }

Here is the call graph for this function:

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

virtual

Construct the branch and tree vector.

Implements VQwSubsystem.

Definition at line 971 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

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

void QwMainCerenkovDetector::ConstructHistograms ( TDirectory * folder )

inlinevirtual

Construct the histograms for this subsystem in a folder.

Reimplemented from VQwSubsystem.

Definition at line 100 of file QwMainCerenkovDetector.h.

References VQwSubsystem::ConstructHistograms().

                                                {
     TString tmpstr("");
     ConstructHistograms(folder,tmpstr);
   };

Here is the call graph for this function:

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

virtual

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

Implements VQwSubsystem.

Definition at line 948 of file QwMainCerenkovDetector.cc.

References VQwSubsystem::ConstructHistograms(), fCombinedPMT, and fIntegrationPMT.

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

Here is the call graph for this function:

void QwMainCerenkovDetector::DeaccumulateRunningSum ( VQwSubsystem * value )

virtual

remove one entry from the running sums for devices

Implements VQwSubsystemParity.

Definition at line 1196 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, and fIntegrationPMT.

                                                                        {
   if (Compare(value1)) {
     QwMainCerenkovDetector* value = dynamic_cast<QwMainCerenkovDetector*>(value1);
 
     for (size_t i = 0; i < fIntegrationPMT.size(); i++)
       fIntegrationPMT[i].DeaccumulateRunningSum(value->fIntegrationPMT[i]);
     for (size_t i = 0; i < fCombinedPMT.size(); i++)
       fCombinedPMT[i].DeaccumulateRunningSum(value->fCombinedPMT[i]);
   }  
 };

Here is the call graph for this function:

void QwMainCerenkovDetector::DefineOptions ( QwOptions & options )

static

Define options function.

Defines configuration options for QwEventBuffer class using QwOptions functionality.

Parameters

options Options object

Definition at line 30 of file QwMainCerenkovDetector.cc.

References QwOptions::AddOptions(), and default_bool_value.

                                                             {
   // Define the execution options
   options.AddOptions()
     ("QwMainCerenkovDetector.normalize",
      po::value<bool>()->default_bool_value(true),
      "Normalize the detectors by beam current");
 }

Here is the call graph for this function:

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

virtual

Implements VQwSubsystemParity.

Definition at line 1124 of file QwMainCerenkovDetector.cc.

References Compare().

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

Here is the call graph for this function:

void QwMainCerenkovDetector::DoNormalization ( Double_t factor = 1.0 )

Definition at line 1312 of file QwMainCerenkovDetector.cc.

References bIsExchangedDataValid, Qw::e, fTargetCharge, and Normalize().

Referenced by ProcessEvent_2().

 {
   if (bIsExchangedDataValid)
     {
       try
         {
     this->Normalize(&fTargetCharge);
         }
       catch (std::exception& e)
         {
           std::cerr << e.what() << std::endl;
         }
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Reimplemented from VQwSubsystem.

Definition at line 617 of file QwMainCerenkovDetector.cc.

References VQwSubsystem::fCurrentBank_ID, VQwSubsystem::fCurrentROC_ID, fIntegrationPMT, fMainDetID, and kQwIntegrationPMT.

 {
   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 < fMainDetID.size(); i++)
     {
       // This is a QwIntegrationPMT
       if (fMainDetID.at(i).fTypeID == kQwIntegrationPMT)
         fIntegrationPMT[fMainDetID.at(i).fIndex].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 QwMainCerenkovDetector::ExchangeProcessedData ( )

virtual

Exchange data between subsystems

Reimplemented from VQwSubsystem.

Definition at line 805 of file QwMainCerenkovDetector.cc.

References bDEBUG, bIsExchangedDataValid, bNormalization, QwLog::endl(), fTargetCharge, VQwSubsystem::GetSubsystemName(), QwError, QwWarning, and VQwSubsystem::RequestExternalValue().

 {
   //QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData "<< QwLog::endl;
   bIsExchangedDataValid = kTRUE;
   if (bNormalization){
     // Create a list of all variables that we need
     // TODO This could be a static list to avoid repeated vector initializiations
     /*
     std::vector<VQwDataElement*> variable_list;
     variable_list.push_back(&fTargetCharge);
     variable_list.push_back(&fTargetX);
     variable_list.push_back(&fTargetY);
     variable_list.push_back(&fTargetXprime);
     variable_list.push_back(&fTargetYprime);
     variable_list.push_back(&fTargetEnergy);
     */
 
     // Loop over all variables in the list
     /*
       //IMPORTANT NOTE ABOUT THE COMMENTED LOOP
       //NAMES OF THE variable_list DEVICES WILL BE OVER WRITTEN BY ACTUAL NAMES OF DATA ELEMENTS THAT WE ARE READING FROM THE
       //ROUTINE RequestExternalValue(variable->GetElementName(), variable). DO TRY THIS AT HOME!
       
 
     std::vector<VQwDataElement*>::iterator variable_iter;
     for (variable_iter  = variable_list.begin(); variable_iter != variable_list.end(); variable_iter++){
   VQwDataElement* variable = *variable_iter;
   //QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData for loop "<<variable->GetElementName()<< QwLog::endl;
   if (RequestExternalValue(variable->GetElementName(), variable))
     {
       if (bDEBUG)
         dynamic_cast<QwVQWK_Channel*>(variable)->PrintInfo();
       //QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<variable->GetElementName()<< QwLog::endl;
     }
   else
     {
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << variable->GetElementName() << QwLog::endl;
     }
       } // end of loop over variables
 
     */
 
     if(RequestExternalValue("q_targ", &fTargetCharge)){
       if (bDEBUG){
   QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<fTargetCharge.GetElementName()<< QwLog::endl;
   //QwWarning <<"****QwMainCerenkovDetector****"<< QwLog::endl;
   (dynamic_cast<QwVQWK_Channel*>(&fTargetCharge))->PrintInfo();
       }
     }
     else{
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << fTargetCharge.GetElementName() << QwLog::endl;
     }
     /*
     if(RequestExternalValue("x_targ", &fTargetX)){
       if (bDEBUG){
   dynamic_cast<QwVQWK_Channel*>(&fTargetX)->PrintInfo();
   QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<fTargetX.GetElementName()<< QwLog::endl;  
       }    
     }else{
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << fTargetX.GetElementName() << QwLog::endl;
     }
 
     if(RequestExternalValue("y_targ", &fTargetY)){
       if (bDEBUG){
   dynamic_cast<QwVQWK_Channel*>(&fTargetY)->PrintInfo();
   QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<fTargetY.GetElementName()<< QwLog::endl;
       }
     }else{
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << fTargetY.GetElementName() << QwLog::endl;
     }
 
     if(RequestExternalValue("xp_targ", &fTargetXprime)){
       if (bDEBUG){
   dynamic_cast<QwVQWK_Channel*>(&fTargetXprime)->PrintInfo();
   QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<fTargetXprime.GetElementName()<< QwLog::endl;
       }
     }else{
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << fTargetXprime.GetElementName() << QwLog::endl;
     }
 
     if(RequestExternalValue("yp_targ", &fTargetYprime)){
       if (bDEBUG){
   dynamic_cast<QwVQWK_Channel*>(&fTargetYprime)->PrintInfo();
   QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<fTargetYprime.GetElementName()<< QwLog::endl;
       }
     }else{
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << fTargetYprime.GetElementName() << QwLog::endl;
     }
 
     if(RequestExternalValue("e_targ", &fTargetEnergy)){
       if (bDEBUG){
   dynamic_cast<QwVQWK_Channel*>(&fTargetEnergy)->PrintInfo();
   QwWarning << "QwMainCerenkovDetector::ExchangeProcessedData Found "<<fTargetEnergy.GetElementName()<< QwLog::endl;
       }
     }else{
       bIsExchangedDataValid = kFALSE;
       QwError << GetSubsystemName() << " could not get external value for "
         << fTargetEnergy.GetElementName() << QwLog::endl;
     }
     */
 
     
   }
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::FillDB	(	QwParityDB *	db,
		TString	type
	)

virtual

Fill the database.

Reimplemented from VQwSubsystemParity.

Definition at line 1327 of file QwMainCerenkovDetector.cc.

References QwDatabase::Connect(), QwDBInterface::DetermineMeasurementTypeID(), QwDatabase::Disconnect(), QwLog::endl(), fCombinedPMT, fIntegrationPMT, QwParityDB::GetAnalysisID(), Qw::kBoldRed, Qw::kGreen, Qw::kNormal, QwDatabase::Query(), and QwMessage.

 {
   Bool_t local_print_flag = false;
 
   if(local_print_flag) {
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
     QwMessage << "            QwMainCerenkovDetector::FillDB                       " << QwLog::endl;
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
   }
 
   std::vector<QwDBInterface> interface;
   std::vector<QwParitySSQLS::md_data> entrylist;
 
   UInt_t analysis_id = db->GetAnalysisID();
 
   TString measurement_type;
   measurement_type = QwDBInterface::DetermineMeasurementTypeID(datatype);
 
   UInt_t i,j;
   i = j = 0;
   if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "IntegrationPMT" <<QwLog::endl;
 
   for(i=0; i<fIntegrationPMT.size(); i++) {
     interface.clear();
     interface = fIntegrationPMT[i].GetDBEntry();
     for(j=0; j<interface.size(); j++) {
       interface.at(j).SetAnalysisID( analysis_id );
       interface.at(j).SetMainDetectorID( db );
       interface.at(j).SetMeasurementTypeID( measurement_type );
       interface.at(j).PrintStatus( local_print_flag );
       interface.at(j).AddThisEntryToList( entrylist );
     }
   }
 
   if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Combined PMT" <<QwLog::endl;
 
   for(i=0; i< fCombinedPMT.size(); i++)
     {
       interface.clear();
       interface = fCombinedPMT[i].GetDBEntry();
       for(j=0; j<interface.size(); j++) {
   interface.at(j).SetAnalysisID( analysis_id );
   interface.at(j).SetMainDetectorID( db );
   interface.at(j).SetMeasurementTypeID( measurement_type );
   interface.at(j).PrintStatus( local_print_flag );
   interface.at(j).AddThisEntryToList( entrylist );
       }
     }
   if(local_print_flag) {
     QwMessage << QwColor(Qw::kGreen) << "Entrylist Size : "
         << QwColor(Qw::kBoldRed) << entrylist.size()
               << QwColor(Qw::kNormal) << QwLog::endl;
   }
 
   db->Connect();
   // Check the entrylist size, if it isn't zero, start to query..
   if( entrylist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(entrylist.begin(), entrylist.end());
     query.execute();
   }
   else {
     QwMessage << "QwMainCerenkovDetector::FillDB :: This is the case when the entrlylist contains nothing in "<< datatype.Data() << QwLog::endl;
   }
   db->Disconnect();
   return;
 }

Here is the call graph for this function:

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

virtual

Reimplemented from VQwSubsystemParity.

Definition at line 1433 of file QwMainCerenkovDetector.cc.

References QwDatabase::Connect(), QwDatabase::Disconnect(), QwLog::endl(), fCombinedPMT, fIntegrationPMT, QwParityDB::GetAnalysisID(), Qw::kBoldRed, Qw::kGreen, Qw::kNormal, QwDatabase::Query(), and QwMessage.

 {
 
   Bool_t local_print_flag = false;
   if(local_print_flag){
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
     QwMessage << "            QwMainCerenkovDetectorID::FillErrDB                  " << QwLog::endl;
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
   }
 
 
   std::vector<QwErrDBInterface> interface;
   std::vector<QwParitySSQLS::md_errors> entrylist;
 
   UInt_t analysis_id = db->GetAnalysisID();
 
   UInt_t i,j;
   i = j = 0;
   if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "IntegrationPMT" <<QwLog::endl;
 
   for(i=0; i<fIntegrationPMT.size(); i++) {
     interface.clear();
     interface = fIntegrationPMT[i].GetErrDBEntry();
     for(j=0; j<interface.size(); j++) {
       interface.at(j).SetAnalysisID     ( analysis_id );
       interface.at(j).SetMainDetectorID ( db );
       interface.at(j).PrintStatus       ( local_print_flag );
       interface.at(j).AddThisEntryToList( entrylist );
     }
   }
 
   if(local_print_flag) QwMessage <<  QwColor(Qw::kGreen) << "Combined PMT" <<QwLog::endl;
 
   for(i=0; i< fCombinedPMT.size(); i++)
     {
       interface.clear();
       interface = fCombinedPMT[i].GetErrDBEntry();
       for(j=0; j<interface.size(); j++) {
   interface.at(j).SetAnalysisID     ( analysis_id );
   interface.at(j).SetMainDetectorID ( db );
   interface.at(j).PrintStatus       ( local_print_flag );
   interface.at(j).AddThisEntryToList( entrylist );
       }
     }
   if(local_print_flag) {
     QwMessage << QwColor(Qw::kGreen) << "Entrylist Size : "
         << QwColor(Qw::kBoldRed) << entrylist.size()
               << QwColor(Qw::kNormal) << QwLog::endl;
   }
 
   db->Connect();
   // Check the entrylist size, if it isn't zero, start to query..
   if( entrylist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(entrylist.begin(), entrylist.end());
     query.execute();
   }
   else {
     QwMessage << "QwMainCerenkovDetector::FillErrDB :: This is the case when the entrlylist contains nothing in "<< datatype.Data() << QwLog::endl;
   }
   db->Disconnect();
 
   return;
 };

Here is the call graph for this function:

void QwMainCerenkovDetector::FillHistograms ( )

virtual

Fill the histograms for this subsystem.

Implements VQwSubsystem.

Definition at line 959 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   for (size_t i=0;i<fIntegrationPMT.size();i++)
     fIntegrationPMT[i].FillHistograms();
 
   for (size_t i=0;i<fCombinedPMT.size();i++)
     fCombinedPMT[i].FillHistograms();
 
   return;
 }

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

virtual

Fill the tree vector.

Implements VQwSubsystem.

Definition at line 1018 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   for (size_t i=0;i<fIntegrationPMT.size();i++)
     fIntegrationPMT[i].FillTreeVector(values);
 
   for (size_t i=0;i<fCombinedPMT.size();i++)
     fCombinedPMT[i].FillTreeVector(values);
 
   return;
 }

const QwIntegrationPMT * QwMainCerenkovDetector::GetChannel ( const TString name ) const

Definition at line 1030 of file QwMainCerenkovDetector.cc.

References GetIntegrationPMT().

Referenced by PublishInternalValues().

 {
   return GetIntegrationPMT(name);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

const QwCombinedPMT * QwMainCerenkovDetector::GetCombinedPMT ( const TString name ) const

Definition at line 1293 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), fCombinedPMT, and QwMessage.

Referenced by PublishInternalValues().

 {
   TString tmpname = name;
   tmpname.ToLower();
   if (! fCombinedPMT.empty())
     {
       for (size_t i=0;i<fCombinedPMT.size();i++)
         {
           if (fCombinedPMT.at(i).GetElementName() == tmpname)
             {
               //std::cout<<"Get CombinedPMT "<<tmpname<<std::endl;
               return &(fCombinedPMT.at(i));
             }
         }
     }
   QwMessage << "QwMainCerenkovDetector::GetCombinedPMT: cannot find channel " << tmpname << QwLog::endl;
   return NULL;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwMainCerenkovDetector::GetDetectorIndex	(	EQwPMTInstrumentType	TypeID,
		TString	name
	)

protected

Definition at line 1248 of file QwMainCerenkovDetector.cc.

References fMainDetID.

Referenced by LoadChannelMap(), and LoadEventCuts().

 {
   Bool_t ldebug=kFALSE;
   if (ldebug)
     {
       std::cout<<"QwMainCerenkovDetector::GetDetectorIndex\n";
       std::cout<<"type_id=="<<type_id<<" name="<<name<<"\n";
       std::cout<<fMainDetID.size()<<" already registered detector\n";
     }
 
   Int_t result=-1;
   for (size_t i=0;i<fMainDetID.size();i++)
     {
       if (fMainDetID[i].fTypeID==type_id)
         if (fMainDetID[i].fdetectorname==name)
           {
             result=fMainDetID[i].fIndex;
             if (ldebug)
               std::cout<<"testing against ("<<fMainDetID[i].fTypeID
               <<","<<fMainDetID[i].fdetectorname<<")=>"<<result<<"\n";
           }
     }
 
   return result;
 }

Here is the caller graph for this function:

EQwPMTInstrumentType QwMainCerenkovDetector::GetDetectorTypeID ( TString name )

protected

Definition at line 1242 of file QwMainCerenkovDetector.cc.

References GetQwPMTInstrumentType().

Referenced by LoadChannelMap(), and LoadEventCuts().

 {
   return GetQwPMTInstrumentType(name);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

UInt_t QwMainCerenkovDetector::GetEventcutErrorFlag ( )

virtual

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

Implements VQwSubsystemParity.

Definition at line 734 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   UInt_t ErrorFlag;
   ErrorFlag=0;
   for(size_t i=0;i<fIntegrationPMT.size();i++){
     ErrorFlag |= fIntegrationPMT[i].GetEventcutErrorFlag();
   }
   for(size_t i=0;i<fCombinedPMT.size();i++){
     ErrorFlag |= fCombinedPMT[i].GetEventcutErrorFlag();
   }
   return ErrorFlag;
 }

const QwIntegrationPMT * QwMainCerenkovDetector::GetIntegrationPMT ( const TString name ) const

Definition at line 1274 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), fIntegrationPMT, and QwMessage.

Referenced by GetChannel(), and PublishInternalValues().

 {
   TString tmpname = name;
   tmpname.ToLower();
   if (! fIntegrationPMT.empty())
     {
       for (size_t i=0;i<fIntegrationPMT.size();i++)
         {
           if (fIntegrationPMT.at(i).GetElementName() == tmpname)
             {
               //std::cout<<"Get IntegrationPMT "<<tmpname<<std::endl;
               return &(fIntegrationPMT.at(i));
             }
         }
     }
   QwMessage << "QwMainCerenkovDetector::GetIntegrationPMT: cannot find channel " << tmpname << QwLog::endl;
   return NULL;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwMainCerenkovDetector::IncrementErrorCounters ( )

virtual

Increment the error counters.

Implements VQwSubsystemParity.

Definition at line 747 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   for(size_t i=0;i<fIntegrationPMT.size();i++){
     fIntegrationPMT[i].IncrementErrorCounters();
   }
   for(size_t i=0;i<fCombinedPMT.size();i++){
     fCombinedPMT[i].IncrementErrorCounters();
   }
 }

Bool_t QwMainCerenkovDetector::IsGoodEvent ( )

Definition at line 565 of file QwMainCerenkovDetector.cc.

 {
   Bool_t test=kTRUE;
   return test;
 }

Int_t QwMainCerenkovDetector::LoadChannelMap ( TString mapfile )

virtual

Mandatory map file definition.

Implements VQwSubsystem.

Definition at line 135 of file QwMainCerenkovDetector.cc.

 {
   Bool_t ldebug=kFALSE;
 
   std::vector<TString> combinedchannelnames;
   std::vector<Double_t> weight;
   Int_t currentrocread=0;
   Int_t currentbankread=0;
   Int_t wordsofar=0;
   Int_t currentsubbankindex=-1;
   Int_t sample_size=0;
 
 
 
   // Open the file
   QwParameterFile mapstr(mapfile.Data());
   TString varname, varvalue;
 
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
   while (mapstr.ReadNextLine())
     {
       mapstr.TrimComment('!');   // Remove everything after a '!' character.
       mapstr.TrimWhitespace();   // Get rid of leading and trailing spaces.
       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);
             }
           else if (varname=="sample_size")
             {
               sample_size=value;
             }
         }
       else
         {
           Bool_t  lineok   = kTRUE;
     TString keyword  = "";
     TString keyword2 = "";
           TString modtype  = "";
           TString dettype  = "";
     TString namech   = "";
           Int_t modnum     = 0;
     Int_t channum    = 0;
 
     modtype = mapstr.GetTypedNextToken<TString>();  // module type
 
           if (modtype == "VQWK")
             {
               modnum    = mapstr.GetTypedNextToken<Int_t>();  //slot number
               channum   = mapstr.GetTypedNextToken<Int_t>();  //channel number
               dettype = mapstr.GetTypedNextToken<TString>();  //type-purpose of the detector
               dettype.ToLower();
               namech  = mapstr.GetTypedNextToken<TString>();  //name of the detector
               namech.ToLower();
 
         keyword   = mapstr.GetTypedNextToken<TString>();
         keyword.ToLower();
         keyword2  = mapstr.GetTypedNextToken<TString>();
         keyword2.ToLower();
             }
           else if (modtype == "VPMT")
             {
               channum       = mapstr.GetTypedNextToken<Int_t>();  //channel number
               Int_t combinedchans = mapstr.GetTypedNextToken<Int_t>();  //number of combined channels
               dettype     = mapstr.GetTypedNextToken<TString>();  //type-purpose of the detector
               dettype.ToLower();
               namech      = mapstr.GetTypedNextToken<TString>();  //name of the detector
               namech.ToLower();
               combinedchannelnames.clear();
               for (int i=0; i<combinedchans; i++)
                 {
                   TString nameofcombinedchan = mapstr.GetTypedNextToken<TString>();
                   nameofcombinedchan.ToLower();
                   combinedchannelnames.push_back(nameofcombinedchan);
                 }
               weight.clear();
               for (int i=0; i<combinedchans; i++)
                 {
                   weight.push_back( mapstr.GetTypedNextToken<Double_t>());
                 }
         keyword  = mapstr.GetTypedNextToken<TString>();
         keyword.ToLower();
         keyword2 = mapstr.GetTypedNextToken<TString>();
         keyword2.ToLower();
             }
 
 
           if (currentsubbankindex!=GetSubbankIndex(currentrocread,currentbankread))
             {
               currentsubbankindex=GetSubbankIndex(currentrocread,currentbankread);
               wordsofar=0;
             }
 
           QwMainCerenkovDetectorID localMainDetID;
           localMainDetID.fdetectorname=namech;
           localMainDetID.fmoduletype=modtype;
           localMainDetID.fSubbankIndex=currentsubbankindex;
           localMainDetID.fdetectortype=dettype;
 
     //          localMainDetID.fWordInSubbank=wordsofar;
           if (modtype=="VQWK"){
       Int_t offset = QwVQWK_Channel::GetBufferOffset(modnum, channum);
       if (offset>0){
         localMainDetID.fWordInSubbank = wordsofar + offset;
       }
     }
           else if (modtype=="VPMT")
             {
               localMainDetID.fCombinedChannelNames = combinedchannelnames;
               localMainDetID.fWeight = weight;
               //std::cout<<"Add in a combined channel"<<std::endl;
             }
           else
             {
               QwError << "QwMainCerenkovDetector::LoadChannelMap:  Unknown module type: "
           << modtype <<", the detector "<<namech<<" will not be decoded "
           << QwLog::endl;
               lineok=kFALSE;
               continue;
             }
 
           localMainDetID.fTypeID=GetDetectorTypeID(dettype);
     if (localMainDetID.fTypeID==kQwUnknownPMT) {
       QwError << "QwMainCerenkovDetector::LoadChannelMap:  Unknown detector type: "
         << dettype <<", the detector "<<namech<<" will not be decoded "
         << QwLog::endl;
       lineok=kFALSE;
       continue;
     }
 
           localMainDetID.fIndex= GetDetectorIndex(localMainDetID.fTypeID,
                                                   localMainDetID.fdetectorname);
 
           if (localMainDetID.fIndex==-1)
             {
               if (localMainDetID.fTypeID==kQwIntegrationPMT)
                 {
                   QwIntegrationPMT localIntegrationPMT(GetSubsystemName(),localMainDetID.fdetectorname);
       if (keyword=="not_blindable"
           || keyword2=="not_blindable")
         localIntegrationPMT.SetBlindability(kFALSE);
       else 
         localIntegrationPMT.SetBlindability(kTRUE);
       if (keyword=="not_normalizable"
           || keyword2=="not_normalizable")
         localIntegrationPMT.SetNormalizability(kFALSE);
       else
         localIntegrationPMT.SetNormalizability(kTRUE);
       fIntegrationPMT.push_back(localIntegrationPMT);
                   fIntegrationPMT[fIntegrationPMT.size()-1].SetDefaultSampleSize(sample_size);
       localMainDetID.fIndex=fIntegrationPMT.size()-1;
                 }
 
               else if (localMainDetID.fTypeID==kQwCombinedPMT)
                 {
       QwCombinedPMT localcombinedPMT(GetSubsystemName(),localMainDetID.fdetectorname);
       if (keyword=="not_normalizable" 
           || keyword2=="not_normalizable")
         localcombinedPMT.SetNormalizability(kFALSE);
       else
         localcombinedPMT.SetNormalizability(kTRUE);
       if (keyword=="not_blindable" 
           || keyword2 =="not_blindable") 
         localcombinedPMT.SetBlindability(kFALSE);
       else 
         localcombinedPMT.SetBlindability(kTRUE);
                   fCombinedPMT.push_back(localcombinedPMT);
                   fCombinedPMT[fCombinedPMT.size()-1].SetDefaultSampleSize(sample_size);
                   localMainDetID.fIndex=fCombinedPMT.size()-1;
                 }
             }
 
           if (ldebug)
             {
               localMainDetID.Print();
               std::cout<<"line ok=";
               if (lineok) std::cout<<"TRUE"<<std::endl;
               else
                 std::cout<<"FALSE"<<std::endl;
             }
 
           if (lineok)
             fMainDetID.push_back(localMainDetID);
         }
     }
 
 
   //std::cout<<"linking combined channels"<<std::endl;
 
   for (size_t i=0; i<fMainDetID.size(); i++)
     {
       if (fMainDetID[i].fTypeID==kQwCombinedPMT)
         {
           Int_t ind = fMainDetID[i].fIndex;
 
           //check to see if all required channels are available
           if (ldebug)
             {
               std::cout<<"fMainDetID[i].fCombinedChannelNames.size()="
               <<fMainDetID[i].fCombinedChannelNames.size()<<std::endl<<"name list: ";
               for (size_t n=0; n<fMainDetID[i].fCombinedChannelNames.size(); n++)
                 std::cout<<"  "<<fMainDetID[i].fCombinedChannelNames[n];
               std::cout<<std::endl;
             }
 
           Int_t chanmatched=0;
           for (size_t j=0; j<fMainDetID[i].fCombinedChannelNames.size(); j++)
             {
 
               for (size_t k=0; k<fMainDetID.size(); k++)
                 {
                   if (fMainDetID[i].fCombinedChannelNames[j]==fMainDetID[k].fdetectorname)
                     {
                       if (ldebug)
                         std::cout<<"found a to-be-combined channel candidate"<<std::endl;
                       chanmatched ++;
                       break;
                     }
                 }
             }
 
           if ((Int_t) fMainDetID[i].fCombinedChannelNames.size()==chanmatched)
             {
               for (size_t l=0; l<fMainDetID[i].fCombinedChannelNames.size(); l++)
                 {
                   Int_t ind_pmt = GetDetectorIndex(GetDetectorTypeID("integrationpmt"),
                                                    fMainDetID[i].fCombinedChannelNames[l]);
 
                   fCombinedPMT[ind].Add(&fIntegrationPMT[ind_pmt],fMainDetID[i].fWeight[l]);
                 }
               fCombinedPMT[ind].LinkChannel(fMainDetID[i].fdetectorname);
               if (ldebug)
                 std::cout<<"linked a combined channel"<<std::endl;
             }
           else
             {
               std::cerr<<"cannot combine void channels for "<<fMainDetID[i].fdetectorname<<std::endl;
               fMainDetID[i].fIndex = -1;
               continue;
             }
         }
     }
 
 
   // Now load the variables to publish
   mapstr.RewindToFileStart();
   QwParameterFile *section;
   std::vector<TString> publishinfo;
   while ((section = mapstr.ReadNextSection(varvalue))) {
     if (varvalue == "PUBLISH") {
       fPublishList.clear();
       while (section->ReadNextLine()) {
         section->TrimComment(); // Remove everything after a comment character
         section->TrimWhitespace(); // Get rid of leading and trailing spaces
         for (int ii = 0; ii < 4; ii++) {
           varvalue = section->GetNextToken().c_str();
           if (varvalue.Length()) {
             publishinfo.push_back(varvalue);
           }
         }
         if (publishinfo.size() == 4)
           fPublishList.push_back(publishinfo);
         publishinfo.clear();
       }
     }
   }
   // Print list of variables to publish
   QwMessage << "Variables to publish:" << QwLog::endl;
   for (size_t jj = 0; jj < fPublishList.size(); jj++)
     QwMessage << fPublishList.at(jj).at(0) << " " << fPublishList.at(jj).at(1) << " "
               << fPublishList.at(jj).at(2) << " " << fPublishList.at(jj).at(3) << QwLog::endl;
 
   if (ldebug)
     {
       std::cout<<"Done with Load channel map\n";
       for (size_t i=0;i<fMainDetID.size();i++)
         if (fMainDetID[i].fIndex>=0)
           fMainDetID[i].Print();
     }
   ldebug=kFALSE;
   mapstr.Close(); // Close the file (ifstream)
   return 0;
 }

Here is the call graph for this function:

Int_t QwMainCerenkovDetector::LoadEventCuts ( TString filename )

virtual

Load the event cuts file.

Implements VQwSubsystemParity.

Definition at line 433 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), fCombinedPMT, VQwSubsystem::fDetectorMaps, fIntegrationPMT, fMainDetErrorCount, GetDetectorIndex(), GetDetectorTypeID(), GetGlobalErrorFlag(), GetQwPMTInstrumentTypeName(), QwParameterFile::GetUInt(), kQwCombinedPMT, kQwIntegrationPMT, QwMessage, and QwWarning.

 {
   Int_t eventcut_flag = 1;
 
   // Open the file
   QwParameterFile mapstr(filename.Data());
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
   while (mapstr.ReadNextLine())
     {
       //std::cout<<"********* In the loop  *************"<<std::endl;
       mapstr.TrimComment('!');   // Remove everything after a '!' character.
       mapstr.TrimWhitespace();   // Get rid of leading and trailing spaces.
       if (mapstr.LineIsEmpty())  continue;
       TString varname, varvalue;
       if (mapstr.HasVariablePair("=",varname,varvalue))
         {
           if (varname=="EVENTCUTS")
             {
               //varname="";
               eventcut_flag = QwParameterFile::GetUInt(varvalue);
               //std::cout<<"EVENT CUT FLAG "<<eventcut_flag<<std::endl;
             }
         }
       else
         {
           TString device_type = mapstr.GetTypedNextToken<TString>();
           device_type.ToLower();
           TString device_name = mapstr.GetTypedNextToken<TString>();
           device_name.ToLower();
           Int_t det_index = GetDetectorIndex(GetDetectorTypeID(device_type),device_name);
     if (det_index == -1) {
       QwWarning << " Device not found " << device_name << " of type " << device_type << QwLog::endl;
       continue;
     }
 
     if (device_type == GetQwPMTInstrumentTypeName(kQwIntegrationPMT)){
       Double_t LLX = mapstr.GetTypedNextToken<Double_t>();  //lower limit for IntegrationPMT value
       Double_t ULX = mapstr.GetTypedNextToken<Double_t>();  //upper limit for IntegrationPMT value
       varvalue = mapstr.GetTypedNextToken<TString>();//global/local
             varvalue.ToLower();
       Double_t stabilitycut = mapstr.GetTypedNextToken<Double_t>();
       QwMessage << "QwMainCerenkovDetector Error Code passing to QwIntegrationPMT " << GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut) << QwLog::endl;
 
       //std::cout<<"*****************************"<<std::endl;
       //std::cout<<" Type "<<device_type<<" Name "<<device_name<<" Index ["<<det_index <<"] "<<" device flag "<<check_flag<<std::endl;
 
       fIntegrationPMT[det_index].SetSingleEventCuts(GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut),LLX,ULX,stabilitycut);
       //std::cout<<"*****************************"<<std::endl;
 
     } else if (device_type == GetQwPMTInstrumentTypeName(kQwCombinedPMT)){
       Double_t LLX = mapstr.GetTypedNextToken<Double_t>();  //lower limit for CombinedPMT value
       Double_t ULX = mapstr.GetTypedNextToken<Double_t>();  //upper limit for CombinedPMT value
       varvalue = mapstr.GetTypedNextToken<TString>();//global/local
             varvalue.ToLower();
       Double_t stabilitycut = mapstr.GetTypedNextToken<Double_t>();
       QwMessage << "QwMainCerenkovDetector Error Code passing to QwCombinedPMT " << GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut) << QwLog::endl;
 
       //std::cout<<"*****************************"<<std::endl;
       //std::cout<<" Type "<<device_type<<" Name "<<device_name<<" Index ["<<det_index <<"] "<<" device flag "<<check_flag<<std::endl;
 
       fCombinedPMT[det_index].SetSingleEventCuts(GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut),LLX,ULX,stabilitycut);
       //std::cout<<"*****************************"<<std::endl;
       
     }
 
         }
 
     }
   for (size_t i = 0; i < fIntegrationPMT.size(); i++)
     fIntegrationPMT[i].SetEventCutMode(eventcut_flag);
   for (size_t i = 0; i < fCombinedPMT.size(); i++)
     fCombinedPMT[i].SetEventCutMode(eventcut_flag);
 
   fMainDetErrorCount = 0; //set the error counter to zero
   mapstr.Close(); // Close the file (ifstream)
   return 0;
 }

Here is the call graph for this function:

Int_t QwMainCerenkovDetector::LoadInputParameters ( TString mapfile )

virtual

Mandatory parameter file definition.

Implements VQwSubsystem.

Definition at line 513 of file QwMainCerenkovDetector.cc.

References VQwSubsystem::fDetectorMaps, and fIntegrationPMT.

 {
   Bool_t ldebug=kFALSE;
   TString varname;
   Double_t varped;
   Double_t varcal;
   TString localname;
 
   Int_t lineread=0;
 
   QwParameterFile mapstr(pedestalfile.Data());  //Open the file
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
   while (mapstr.ReadNextLine())
     {
       lineread+=1;
       if (ldebug)std::cout<<" line read so far ="<<lineread<<"\n";
       mapstr.TrimComment('!');   // Remove everything after a '!' character.
       mapstr.TrimWhitespace();   // Get rid of leading and trailing spaces.
       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
           if (ldebug) std::cout<<"inputs for channel "<<varname
             <<": ped="<<varped<<": cal="<<varcal<<"\n";
           Bool_t notfound=kTRUE;
 
           if (notfound)
             for (size_t i=0;i<fIntegrationPMT.size();i++)
               if (fIntegrationPMT[i].GetElementName()==varname)
                 {
                   fIntegrationPMT[i].SetPedestal(varped);
                   fIntegrationPMT[i].SetCalibrationFactor(varcal);
                   i=fIntegrationPMT.size()+1;
                   notfound=kFALSE;
                   i=fIntegrationPMT.size()+1;
                 }
         }
 
     }
   if (ldebug) std::cout<<" line read in the pedestal + cal file ="<<lineread<<" \n";
 
   ldebug=kFALSE;
   mapstr.Close(); // Close the file (ifstream)
   return 0;
 }

void QwMainCerenkovDetector::Normalize ( VQwDataElement * denom )

Definition at line 1163 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

Referenced by DoNormalization().

 {
   for (size_t i = 0; i < fIntegrationPMT.size(); i++)
     fIntegrationPMT[i].Normalize(denom);
   for (size_t i = 0; i < fCombinedPMT.size(); i++)
     fCombinedPMT[i].Normalize(denom);
 }

Here is the caller graph for this function:

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

virtual

Implements VQwSubsystemParity.

Definition at line 1079 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, and fIntegrationPMT.

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

Here is the call graph for this function:

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

virtual

Implements VQwSubsystemParity.

Definition at line 1096 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, and fIntegrationPMT.

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

Here is the call graph for this function:

VQwSubsystem & QwMainCerenkovDetector::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 1061 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, and fIntegrationPMT.

 {
   //  std::cout<<" here in QwMainCerenkovDetector::operator= \n";
   if (Compare(value))
     {
 //      VQwSubsystem::operator=(value);
       QwMainCerenkovDetector* input = dynamic_cast<QwMainCerenkovDetector*> (value);
 
       for (size_t i=0;i<input->fIntegrationPMT.size();i++)
         this->fIntegrationPMT[i]=input->fIntegrationPMT[i];
 
       for (size_t i=0;i<input->fCombinedPMT.size();i++)
         (this->fCombinedPMT[i])=(input->fCombinedPMT[i]);
     }
   return *this;
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::PrintDetectorID ( ) const

Definition at line 1419 of file QwMainCerenkovDetector.cc.

References fMainDetID.

 {
   for (size_t i=0;i<fMainDetID.size();i++)
     {
       std::cout<<"============================="<<std::endl;
       std::cout<<" Detector ID="<<i<<std::endl;
       fMainDetID[i].Print();
     }
   return;
 }

void QwMainCerenkovDetector::PrintErrorCounters ( ) const

virtual

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

Implements VQwSubsystemParity.

Definition at line 758 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), fCombinedPMT, fIntegrationPMT, QwVQWK_Channel::PrintErrorCounterHead(), QwVQWK_Channel::PrintErrorCounterTail(), and QwMessage.

 {
   QwMessage<<"*********QwMainCerenkovDetector Error Summary****************"<<QwLog::endl;
   QwVQWK_Channel::PrintErrorCounterHead();
   for(size_t i=0;i<fIntegrationPMT.size();i++){
     //std::cout<<"  IntegrationPMT ["<<i<<"] "<<std::endl;
     fIntegrationPMT[i].PrintErrorCounters();
   }
   for(size_t i=0;i<fCombinedPMT.size();i++){
     //std::cout<<"  CombinedPMT ["<<i<<"] "<<std::endl;
     fCombinedPMT[i].PrintErrorCounters();
   }
   QwVQWK_Channel::PrintErrorCounterTail();
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::PrintInfo ( ) const

virtual

Print some information about the subsystem.

Reimplemented from VQwSubsystem.

Definition at line 1404 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, fIntegrationPMT, and VQwSubsystem::fSystemName.

 {
   std::cout<<"Name of the subsystem ="<<fSystemName<<"\n";
 
   std::cout<<"there are "<<fIntegrationPMT.size()<<" IntegrationPMT \n";
   std::cout<<"          "<<fCombinedPMT.size()<<" CombinedPMT \n";
 
   std::cout<<" Printing Running AVG and other channel info"<<std::endl;
 
   for (size_t i = 0; i < fIntegrationPMT.size(); i++)
     fIntegrationPMT[i].PrintInfo();
   for (size_t i = 0; i < fCombinedPMT.size(); i++)
     fCombinedPMT[i].PrintInfo();
 }

void QwMainCerenkovDetector::PrintValue ( ) const

virtual

Print values of all channels.

Reimplemented from VQwSubsystemParity.

Definition at line 1395 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), fCombinedPMT, fIntegrationPMT, VQwSubsystem::GetSubsystemName(), and QwMessage.

 {
   QwMessage << "=== QwMainCerenkovDetector: " << GetSubsystemName() << " ===" << QwLog::endl;
   for (size_t i = 0; i < fIntegrationPMT.size(); i++)
     fIntegrationPMT[i].PrintValue();
   for (size_t i = 0; i < fCombinedPMT.size(); i++)
     fCombinedPMT[i].PrintValue();
 }

Here is the call graph for this function:

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

virtual

Implements VQwSubsystem.

Definition at line 657 of file QwMainCerenkovDetector.cc.

 {
 
   /*  Int_t index = GetSubbankIndex(roc_id,bank_id);
     if (index>=0 && num_words>0){
       //  We want to process the configuration data for this ROC.
       UInt_t words_read = 0;
       for (size_t i = 0; i < fMainDetID.size(); i++) {
         words_read += fIntegrationPMT[i].ProcessConfigurationBuffer(&(buffer[words_read]),
                    num_words-words_read);
       }
     }*/
   return 0;
 }

Int_t QwMainCerenkovDetector::ProcessEvBuffer	(	const UInt_t	roc_id,
		const 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 673 of file QwMainCerenkovDetector.cc.

References fIntegrationPMT, fMainDetID, VQwSubsystem::GetSubbankIndex(), and kQwIntegrationPMT.

 {
   Bool_t lkDEBUG=kFALSE;
 
   Int_t index = GetSubbankIndex(roc_id,bank_id);
 
   if (index>=0 && num_words>0)
     {
       //  We want to process this ROC.  Begin looping through the data.
       if (lkDEBUG)
         std::cout << "QwMainCerenkovDetector::ProcessEvBuffer:  "
         << "Begin processing ROC" << roc_id
         << " and subbank "<<bank_id
         << " number of words="<<num_words<<std::endl;
 
       for (size_t i=0;i<fMainDetID.size();i++)
         {
           if (fMainDetID[i].fSubbankIndex==index)
             {
 
               if (fMainDetID[i].fTypeID == kQwIntegrationPMT)
                 {
                   if (lkDEBUG)
                     {
                       std::cout<<"found IntegrationPMT data for "<<fMainDetID[i].fdetectorname<<std::endl;
                       std::cout<<"word left to read in this buffer:"<<num_words-fMainDetID[i].fWordInSubbank<<std::endl;
                     }
                   fIntegrationPMT[fMainDetID[i].fIndex].ProcessEvBuffer(&(buffer[fMainDetID[i].fWordInSubbank]),
                       num_words-fMainDetID[i].fWordInSubbank);
                 }
             }
         }
     }
 
   return 0;
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::ProcessEvent ( )

virtual

Implements VQwSubsystem.

Definition at line 787 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

 {
   for (size_t i=0;i<fIntegrationPMT.size();i++)
     fIntegrationPMT[i].ProcessEvent();
 
   for (size_t i=0;i<fCombinedPMT.size();i++)
     {
       //std::cout<<"Process combination "<<i<<std::endl;
       fCombinedPMT[i].ProcessEvent();
 
     }
 
   return;
 }

void QwMainCerenkovDetector::ProcessEvent_2 ( )

virtual

Process the event data again, including data from other subsystems. Not all derived classes will require a second stage of event data processing.

Reimplemented from VQwSubsystem.

Definition at line 923 of file QwMainCerenkovDetector.cc.

References bDEBUG, bIsExchangedDataValid, bNormalization, DoNormalization(), QwLog::endl(), fTargetCharge, and QwWarning.

 {
   if (bIsExchangedDataValid)
     {
       //data is valid, process it
       if (bDEBUG)
         {
           Double_t  pedestal = fTargetCharge.GetPedestal();
           Double_t  calfactor = fTargetCharge.GetCalibrationFactor();
           Double_t  volts = fTargetCharge.GetAverageVolts();
           std::cout<<"QwMainCerenkovDetector::ProcessEvent_2(): processing with exchanged data"<<std::endl;
           std::cout<<"pedestal, calfactor, average volts = "<<pedestal<<", "<<calfactor<<", "<<volts<<std::endl;
         }
 
       if (bNormalization) this->DoNormalization();
     }
   else
     {
       QwWarning<<"QwMainCerenkovDetector::ProcessEvent_2(): could not get all external values."<<QwLog::endl;
     }
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::ProcessOptions ( QwOptions & options )

virtual

Loads the configuration options into this instance of QwMainCerenkovDetector from the QwOptions object.

Parameters

options Options object

Reimplemented from VQwSubsystem.

Definition at line 45 of file QwMainCerenkovDetector.cc.

References bNormalization, QwLog::endl(), QwOptions::GetValue(), and QwWarning.

                                                              {
   bNormalization = options.GetValue<bool>("QwMainCerenkovDetector.normalize");
   if (! bNormalization){
     QwWarning << "QwMainCerenkovDetector::ProcessOptions:  "
         << "Detector yields WILL NOT be normalized."
         << QwLog::endl;
   }
 }

Here is the call graph for this function:

Bool_t QwMainCerenkovDetector::PublishInternalValues ( ) const

virtual

Publish internal values

Returns

Reimplemented from VQwSubsystem.

Definition at line 60 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), VQwSubsystem::fPublishList, QwCombinedPMT::GetChannel(), QwIntegrationPMT::GetChannel(), GetChannel(), GetCombinedPMT(), GetIntegrationPMT(), VQwSubsystem::PublishInternalValue(), QwDebug, and QwError.

 {
   // Publish variables
   Bool_t status = kTRUE;
 
   status = status && PublishInternalValue("qwk_md1neg", "qwk_md1neg", GetIntegrationPMT("qwk_md1neg")->GetChannel("qwk_md1neg"));
   status = status && PublishInternalValue("qwk_md1pos", "qwk_md1pos", GetIntegrationPMT("qwk_md1pos")->GetChannel("qwk_md1pos"));
   status = status && PublishInternalValue("qwk_md2neg", "qwk_md2neg", GetIntegrationPMT("qwk_md2neg")->GetChannel("qwk_md2neg"));
   status = status && PublishInternalValue("qwk_md2pos", "qwk_md2pos", GetIntegrationPMT("qwk_md2pos")->GetChannel("qwk_md2pos"));
   status = status && PublishInternalValue("qwk_md3neg", "qwk_md3neg", GetIntegrationPMT("qwk_md3neg")->GetChannel("qwk_md3neg"));
   status = status && PublishInternalValue("qwk_md3pos", "qwk_md3pos", GetIntegrationPMT("qwk_md3pos")->GetChannel("qwk_md3pos"));
   status = status && PublishInternalValue("qwk_md4neg", "qwk_md4neg", GetIntegrationPMT("qwk_md4neg")->GetChannel("qwk_md4neg"));
   status = status && PublishInternalValue("qwk_md4pos", "qwk_md4pos", GetIntegrationPMT("qwk_md4pos")->GetChannel("qwk_md4pos"));
   status = status && PublishInternalValue("qwk_md5neg", "qwk_md5neg", GetIntegrationPMT("qwk_md5neg")->GetChannel("qwk_md5neg"));
   status = status && PublishInternalValue("qwk_md5pos", "qwk_md5pos", GetIntegrationPMT("qwk_md5pos")->GetChannel("qwk_md5pos"));
   status = status && PublishInternalValue("qwk_md6neg", "qwk_md6neg", GetIntegrationPMT("qwk_md6neg")->GetChannel("qwk_md6neg"));
   status = status && PublishInternalValue("qwk_md6pos", "qwk_md6pos", GetIntegrationPMT("qwk_md6pos")->GetChannel("qwk_md6pos"));
   status = status && PublishInternalValue("qwk_md7neg", "qwk_md7neg", GetIntegrationPMT("qwk_md7neg")->GetChannel("qwk_md7neg"));
   status = status && PublishInternalValue("qwk_md7pos", "qwk_md7pos", GetIntegrationPMT("qwk_md7pos")->GetChannel("qwk_md7pos"));
   status = status && PublishInternalValue("qwk_md8neg", "qwk_md8neg", GetIntegrationPMT("qwk_md8neg")->GetChannel("qwk_md8neg"));
   status = status && PublishInternalValue("qwk_md8pos", "qwk_md8pos", GetIntegrationPMT("qwk_md8pos")->GetChannel("qwk_md8pos"));
   status = status && PublishInternalValue("qwk_md9neg", "qwk_md9neg", GetIntegrationPMT("qwk_md9neg")->GetChannel("qwk_md9neg"));
   status = status && PublishInternalValue("qwk_md9pos", "qwk_md9pos", GetIntegrationPMT("qwk_md9pos")->GetChannel("qwk_md9pos"));
 
   status = status && PublishInternalValue("qwk_md1barsum","qwk_md1barsum", GetCombinedPMT("qwk_md1barsum")->GetChannel("qwk_md1barsum"));
   status = status && PublishInternalValue("qwk_md2barsum","qwk_md2barsum", GetCombinedPMT("qwk_md2barsum")->GetChannel("qwk_md2barsum"));
   status = status && PublishInternalValue("qwk_md3barsum","qwk_md3barsum", GetCombinedPMT("qwk_md3barsum")->GetChannel("qwk_md3barsum"));
   status = status && PublishInternalValue("qwk_md4barsum","qwk_md4barsum", GetCombinedPMT("qwk_md4barsum")->GetChannel("qwk_md4barsum"));
   status = status && PublishInternalValue("qwk_md5barsum","qwk_md5barsum", GetCombinedPMT("qwk_md5barsum")->GetChannel("qwk_md5barsum"));
   status = status && PublishInternalValue("qwk_md6barsum","qwk_md6barsum", GetCombinedPMT("qwk_md6barsum")->GetChannel("qwk_md6barsum"));
   status = status && PublishInternalValue("qwk_md7barsum","qwk_md7barsum", GetCombinedPMT("qwk_md7barsum")->GetChannel("qwk_md7barsum"));
   status = status && PublishInternalValue("qwk_md8barsum","qwk_md8barsum", GetCombinedPMT("qwk_md8barsum")->GetChannel("qwk_md8barsum"));
  
   status = status && PublishInternalValue("qwk_mdallbars","qwk_mdallbars", GetCombinedPMT("qwk_mdallbars")->GetChannel("qwk_mdallbars"));
 
   return status;
 
 
   // TODO:
   // The variables should be published based on the parameter file.
   // See QwBeamLine class for an implementation.
 
   // Publish variables through map file
   // This should work with bcm, bpmstripline, bpmcavity, combo bpm and combo bcm
   for (size_t pp = 0; pp < fPublishList.size(); pp++) {
     TString publish_name = fPublishList.at(pp).at(0);
     TString device_type = fPublishList.at(pp).at(1);
     TString device_name = fPublishList.at(pp).at(2);
     TString device_prop = fPublishList.at(pp).at(3);
     device_type.ToLower();
     device_prop.ToLower();
 
     const VQwHardwareChannel* tmp_channel;
     if (device_type == "integrationpmt") {
       tmp_channel = GetIntegrationPMT(device_name)->GetChannel(device_name);
     } else if (device_type == "combinedpmt") {
       tmp_channel = GetCombinedPMT(device_name)->GetChannel(device_name);
     } else
       QwError << "QwBeamLine::PublishInternalValues() error "<< QwLog::endl;
 
     if (tmp_channel == NULL) {
       QwError << "QwBeamLine::PublishInternalValues(): " << publish_name << " not found" << QwLog::endl;
       status |= kFALSE;
     } else {
       QwDebug << "QwBeamLine::PublishInternalValues(): " << publish_name << " found" << QwLog::endl;
     }
     status = status && PublishInternalValue(publish_name, publish_name, tmp_channel);
 
   }
 
   return status;
 }

Here is the call graph for this function:

void QwMainCerenkovDetector::RandomizeEventData	(	int	helicity = `0`,
		Double_t	time = `0.0`
	)

Definition at line 606 of file QwMainCerenkovDetector.cc.

References fIntegrationPMT, fMainDetID, and kQwIntegrationPMT.

 {
   for (size_t i = 0; i < fMainDetID.size(); i++)
     {
       // This is a QwIntegrationPMT
       if (fMainDetID.at(i).fTypeID == kQwIntegrationPMT)
         fIntegrationPMT[fMainDetID.at(i).fIndex].RandomizeEventData(helicity, time);
     }
 
 }

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

virtual

Implements VQwSubsystemParity.

Definition at line 1133 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, and fIntegrationPMT.

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

Here is the call graph for this function:

void QwMainCerenkovDetector::Scale ( Double_t factor )

virtual

Implements VQwSubsystemParity.

Definition at line 1151 of file QwMainCerenkovDetector.cc.

References fCombinedPMT, and fIntegrationPMT.

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

void QwMainCerenkovDetector::SetRandomEventAsymmetry ( Double_t asymmetry )

Definition at line 595 of file QwMainCerenkovDetector.cc.

References fIntegrationPMT, fMainDetID, and kQwIntegrationPMT.

Referenced by main().

 {
   for (size_t i = 0; i < fMainDetID.size(); i++)
     {
       // This is a QwIntegrationPMT
       if (fMainDetID.at(i).fTypeID == kQwIntegrationPMT)
         fIntegrationPMT[fMainDetID.at(i).fIndex].SetRandomEventAsymmetry(asymmetry);
     }
 
 }

Here is the caller graph for this function:

void QwMainCerenkovDetector::SetRandomEventParameters	(	Double_t	mean,
		Double_t	sigma
	)

Definition at line 584 of file QwMainCerenkovDetector.cc.

References fIntegrationPMT, fMainDetID, and kQwIntegrationPMT.

Referenced by main().

 {
   for (size_t i = 0; i < fMainDetID.size(); i++)
     {
       // This is a QwIntegrationPMT
       if (fMainDetID.at(i).fTypeID == kQwIntegrationPMT)
         fIntegrationPMT[fMainDetID.at(i).fIndex].SetRandomEventParameters(mean, sigma);
     }
 
 }

Here is the caller graph for this function:

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

virtual

Implements VQwSubsystemParity.

Definition at line 1115 of file QwMainCerenkovDetector.cc.

References Compare().

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

Here is the call graph for this function:

void QwMainCerenkovDetector::UpdateErrorFlag ( const VQwSubsystem * ev_error )

virtual

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 773 of file QwMainCerenkovDetector.cc.

References Compare(), fCombinedPMT, fIntegrationPMT, and VQwSubsystemParity::UpdateErrorFlag().

                                                                         {
   VQwSubsystem* tmp = const_cast<VQwSubsystem*>(ev_error);
   if(Compare(tmp)){
     const QwMainCerenkovDetector* input = dynamic_cast<const QwMainCerenkovDetector*> (ev_error);
 
     for (size_t i=0;i<input->fIntegrationPMT.size();i++)
       this->fIntegrationPMT[i].UpdateErrorFlag(&(input->fIntegrationPMT[i]));
     
     for (size_t i=0;i<input->fCombinedPMT.size();i++)
       this->fCombinedPMT[i].UpdateErrorFlag(&(input->fCombinedPMT[i]));
   }
 };

Here is the call graph for this function:

void QwMainCerenkovDetector::WritePromptSummary	(	QwPromptSummary *	ps,
		TString	type
	)

virtual

Reimplemented from VQwSubsystemParity.

Definition at line 1500 of file QwMainCerenkovDetector.cc.

References QwLog::endl(), and QwMessage.

 {
 
   Bool_t local_print_flag = true;
   if(local_print_flag){
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
     QwMessage << "        QwMainCerenkovDetectorID::WritePromptSummary()          " << QwLog::endl;
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
   }
 
   //  ps->PrintCSV();
 
   return;
 }