#include <QwRaster.h>

Inheritance diagram for QwRaster:

Collaboration diagram for QwRaster:

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

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

void	ProcessOptions (QwOptions &options)
	Process the command line options. More...

Int_t	LoadEventCuts (TString filename)
	Optional event cut file. More...

Bool_t	ApplySingleEventCuts ()

void	PrintErrorCounters () const

Bool_t	CheckRunningAverages (Bool_t)

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

Int_t	GetDetectorInfo (std::vector< std::vector< QwDetectorInfo > > &detector_info)

Int_t	LoadInputParameters (TString parameterfile)
	Mandatory parameter file definition. More...

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

void	InitializeChannel (TString name, TString datatosave)

void	ClearEventData ()

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

void	ConstructHistograms (TDirectory *folder, TString &prefix)
	Construct the histograms for this subsystem in a folder with a prefix. More...

void	FillHistograms ()
	Fill the histograms for this subsystem. More...

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

void	FillTreeVector (std::vector< Double_t > &values) const
	Fill the tree vector. More...

void	FillDB (QwDatabase *db, TString type)

void	GetHitList (QwHitContainer &grandHitContainer)
	Get the hit list. More...

void	ReportConfiguration ()

void	PrintValue ()

void	PrintInfo ()

Public Member Functions inherited from VQwSubsystemTracking
	VQwSubsystemTracking (TString name)
	Constructor with name. More...

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

Int_t	LoadGeometryDefinition (TString mapfile)
	Load geometry definition for tracking subsystems. More...

Int_t	LoadCrosstalkDefinition (TString mapfile)
	Load crosstalk definition for tracking subsystems. More...

const QwGeometry &	GetDetectorInfo () const
	Get the detector geometry information. More...

virtual void	FillHardwareErrorSummary ()
	Hardware error summary. More...

virtual void	ConstructBranch (TTree *tree, TString &prefix)
	Construct the branch and tree vector. More...

virtual void	ConstructBranch (TTree *tree, TString &prefix, QwParameterFile &trim_file)
	Construct the branch and tree vector based on the trim file. More...

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

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	EncodeEventData (std::vector< UInt_t > &buffer)

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

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

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

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, QwRaster >
virtual	~MQwCloneable ()
	Virtual destructor. More...

virtual VQwSubsystem *	Clone () const
	Concrete clone method. More...

const VQwFactory< VQwSubsystem > *	Factory () const
	Factory getter. More...

Public Member Functions inherited from VQwCloneable< VQwSubsystem >
virtual	~VQwCloneable ()
	Virtual destructor. More...

std::string	GetClassName () const
	Get demangled name of this class. More...

Protected Member Functions
void	FillRawWord (Int_t bank_index, Int_t slot_num, Int_t chan, UInt_t data)

void	ClearAllBankRegistrations ()

Int_t	RegisterROCNumber (const UInt_t roc_id)

Int_t	RegisterSubbank (const UInt_t bank_id)

Int_t	RegisterSlotNumber (const UInt_t slot_id)

EQwModuleType	RegisterModuleType (TString moduletype)

Int_t	GetModuleIndex (size_t bank_index, size_t slot_num) const

Bool_t	IsSlotRegistered (Int_t bank_index, Int_t slot_num) const

Int_t	LinkChannelToSignal (const UInt_t chan, const TString &name)

Int_t	FindSignalIndex (const EQwModuleType modtype, const TString &name) const

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
EQwModuleType	fCurrentType

Bool_t	fDEBUG

MQwV775TDC	fQDCTDC

TString	fRegion

size_t	fCurrentBankIndex
	Name of this subsystem (the region). More...

Int_t	fCurrentSlot

Int_t	fCurrentIndex

UInt_t	fBankID [2]

Int_t	fNumberOfModules

std::vector< std::vector< Int_t > >	fModuleIndex

std::vector< EQwModuleType >	fModuleTypes
	Module index, indexed by bank_index and slot_number. More...

std::vector< std::vector < std::pair< EQwModuleType, Int_t > > >	fModulePtrs

std::vector< std::vector < QwPMT_Channel > >	fPMTs

std::vector< QwSIS3801D24_Channel >	fSCAs

std::map< TString, size_t >	fSCAs_map

std::vector< Int_t >	fSCAs_offset

TH2D *	fRateMap

Protected Attributes inherited from VQwSubsystemTracking
QwGeometry	fDetectorInfo
	Geometry information of this subsystem. More...

size_t	fTreeArrayIndex
	Tree indices. More...

size_t	fTreeArrayNumEntries

Double_t	fF1TDCResolutionNS

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

Static Protected Attributes
static const UInt_t	kMaxNumberOfModulesPerROC = 21

static const UInt_t	kMaxNumberOfChannelsPerModule = 32

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

Private Attributes
Double_t	fPositionOffsetX

Double_t	fPositionOffsetY

Double_t	fChannel_Offset_X

Double_t	fChannel_Offset_Y

Double_t	fCal_Factor_QDC_X

Double_t	fCal_Factor_QDC_Y

Double_t	fPositionX_ADC

Double_t	fPositionY_ADC

Double_t	fbpm_3h07a_slope

Double_t	fbpm_3h07a_intercept

Double_t	fbpm_3h09b_slope

Double_t	fbpm_3h09b_intercept

Double_t	fbpm_3h07a_xp

Double_t	fbpm_3h07a_xm

Double_t	fbpm_3h07a_yp

Double_t	fbpm_3h07a_ym

Double_t	fbpm_3h09b_xp

Double_t	fbpm_3h09b_xm

Double_t	fbpm_3h09b_yp

Double_t	fbpm_3h09b_ym

Double_t	fbpm_3h07a_pos_x

Double_t	fbpm_3h07a_pos_y

Double_t	fbpm_3h09b_pos_x

Double_t	fbpm_3h09b_pos_y

Additional Inherited Members
Static Public Member Functions inherited from VQwSubsystem
static void	DefineOptions ()
	Define options function (note: no virtual static functions in C++) More...

Static Public Member Functions inherited from MQwCloneable< VQwSubsystem, QwRaster >
static VQwSubsystem *	Create (const std::string &name)
	Object creation. More...

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

Detailed Description

Definition at line 31 of file QwRaster.h.

Constructor & Destructor Documentation

QwRaster::QwRaster ( )

private

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

QwRaster::QwRaster ( const TString & name )

Constructor with name.

Definition at line 24 of file QwRaster.cc.

References ClearAllBankRegistrations().

 : VQwSubsystem(name),
   VQwSubsystemTracking(name)
 {
     ClearAllBankRegistrations();
 }

Here is the call graph for this function:

QwRaster::~QwRaster ( )

virtual

Virtual destructor.

Definition at line 32 of file QwRaster.cc.

References fPMTs, and fSCAs.

 {
   fPMTs.clear();
   fSCAs.clear();
 }

Member Function Documentation

Bool_t QwRaster::ApplySingleEventCuts ( )

inline

Definition at line 47 of file QwRaster.h.

47 { return kTRUE; };

Bool_t QwRaster::CheckRunningAverages ( Bool_t )

inline

Definition at line 49 of file QwRaster.h.

49 { return kTRUE; };

void QwRaster::ClearAllBankRegistrations ( )

protected

Definition at line 603 of file QwRaster.cc.

References VQwSubsystem::ClearAllBankRegistrations(), fModuleIndex, fModulePtrs, fModuleTypes, and fNumberOfModules.

Referenced by QwRaster().

 {
     VQwSubsystem::ClearAllBankRegistrations();
     fModuleIndex.clear();
     fModulePtrs.clear();
     fModuleTypes.clear();
     fNumberOfModules = 0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwRaster::ClearEventData ( )

virtual

Implements VQwSubsystem.

Definition at line 202 of file QwRaster.cc.

References fPMTs, fSCAs, and VQwSubsystem::SetDataLoaded().

 {
     SetDataLoaded(kFALSE);
 
     for (size_t i=0; i<fPMTs.size(); i++)
     {
         for (size_t j=0; j<fPMTs.at(i).size(); j++)
         {
             fPMTs.at(i).at(j).SetValue(0);
         }
     }
 
     /*for (size_t i=0; i<fADC_Data.size(); i++)
     {
         if (fADC_Data.at(i) != NULL)
         {
             fADC_Data.at(i)->ClearEventData();
         }
     }*/
 
     for (size_t i=0; i<fSCAs.size(); i++) {
       fSCAs.at(i).ClearEventData();
     }
 
 }

Here is the call graph for this function:

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

virtual

Construct the branch and tree vector.

Reimplemented from VQwSubsystemTracking.

Definition at line 495 of file QwRaster.cc.

References fPMTs, fSCAs, VQwSubsystemTracking::fTreeArrayIndex, and VQwSubsystemTracking::fTreeArrayNumEntries.

 {
     fTreeArrayIndex = values.size();
 
     TString basename;
     if (prefix=="") basename = "beamline";
     else basename = prefix;
 
     values.push_back(0.0);
     TString list = "PositionX_ADC/D";
     values.push_back(0.0);
     list += ":PositionY_ADC/D";
     list += ":bpm_3h07a_pos_x/D";
     list += ":bpm_3h07a_pos_y/D";
     list += ":bpm_3h09b_pos_x/D";
     list += ":bpm_3h09b_pos_y/D";
 
     for (size_t i=0; i<fPMTs.size(); i++) {
       for (size_t j=0; j<fPMTs.at(i).size(); j++) {
         TString element_name = fPMTs.at(i).at(j).GetElementName();
         if (element_name=="") {
           //  This channel is not used, so skip setting up the tree.
         } else {
             values.push_back(0.0);
             list += ":"+element_name+"_raw/D";
         }
       }
     }
 
     for (size_t i=0; i<fSCAs.size(); i++){
       if (fSCAs.at(i).GetElementName() != "") {
         values.push_back(0.0);
         list += ":" + fSCAs.at(i).GetElementName() + "/D";
       }
     }
 
     if (list[0]==':') {
       list = list(1,list.Length()-1);
     }
 
     fTreeArrayNumEntries = values.size() - fTreeArrayIndex;
     tree->Branch(basename, &values[fTreeArrayIndex], list);
 
     return;
 }

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

virtual

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

Implements VQwSubsystem.

Definition at line 394 of file QwRaster.cc.

References QwHistogramHelper::Construct1DHist(), MQwHistograms::fHistograms, fPMTs, fRateMap, fSCAs, and gQwHists.

 {
     TString basename;
 
     if (prefix == "")  basename = "";
     else  basename = prefix;
 
     if (folder != NULL) folder->cd();
 
     for (size_t i=0; i<fPMTs.size(); i++)
     {
         for (size_t j=0; j<fPMTs.at(i).size(); j++)
             fPMTs.at(i).at(j).ConstructHistograms(folder, basename);
     }
     
     for (size_t i=0; i<fSCAs.size(); i++) {
       fSCAs.at(i).ConstructHistograms(folder, prefix);
     }
 
     fHistograms.push_back( gQwHists.Construct1DHist(TString("raster_position_x")));
     fHistograms.push_back( gQwHists.Construct1DHist(TString("raster_position_y")));
   
     fRateMap  = new TH2D("raster_rate_map","Raster Rate Map",125,0,0,125,0,0);
     
     fRateMap->GetXaxis()->SetTitle(" X [mm]");
     fRateMap->GetYaxis()->SetTitle(" Y [mm]");
 
     gStyle     -> SetPalette(1);
     fRateMap->SetOption("colz");
 
     fHistograms.push_back( gQwHists.Construct1DHist(TString("bpm_3h07a_pos_x")));
     fHistograms.push_back( gQwHists.Construct1DHist(TString("bpm_3h07a_pos_y")));
     fHistograms.push_back( gQwHists.Construct1DHist(TString("bpm_3h09b_pos_x")));
     fHistograms.push_back( gQwHists.Construct1DHist(TString("bpm_3h09b_pos_y")));
 
 }

Here is the call graph for this function:

void QwRaster::FillDB	(	QwDatabase *	db,
		TString	type
	)

inline

Definition at line 70 of file QwRaster.h.

70 {return;};

void QwRaster::FillHistograms ( )

virtual

Fill the histograms for this subsystem.

Implements VQwSubsystem.

Definition at line 431 of file QwRaster.cc.

References fbpm_3h07a_pos_x, fbpm_3h07a_pos_y, fbpm_3h09b_pos_x, fbpm_3h09b_pos_y, MQwHistograms::fHistograms, fPMTs, fPositionX_ADC, fPositionY_ADC, fRateMap, fSCAs, and VQwSubsystem::HasDataLoaded().

 {
 
     if (! HasDataLoaded()) return;
 
     for (size_t i=0; i<fPMTs.size(); i++)
     {
         for (size_t j=0; j<fPMTs.at(i).size(); j++)
         {
             fPMTs.at(i).at(j).FillHistograms();
         }
     }
 
     for (size_t i=0; i<fSCAs.size(); i++) {
       fSCAs.at(i).FillHistograms();
    }
    
     // FR fudge factor is 3.2, by Dave Mack
     // See hclog ttps://hallcweb.jlab.org/hclog/1010_archive/101012134143.html
     // I inverse x axis and multiply the FR fudge factor to x and y values
     // in order to see the unit is mm in the FR plot
     // Tuesday, October 12 17:11:52 EDT 2010, jhlee
 
     Double_t fudge_factor = 3.2;
     
     Double_t raster_x_mm = 0.0;
     Double_t raster_y_mm = 0.0;
 
     raster_x_mm = - fudge_factor*fPositionX_ADC;
     raster_y_mm =   fudge_factor*fPositionY_ADC;
 
     for (size_t j=0; j<fHistograms.size();j++)
     {
       if (fHistograms.at(j)->GetTitle()==TString("raster_position_x"))
         {
     fHistograms.at(j)->Fill(raster_x_mm);
         }
       if (fHistograms.at(j)->GetTitle()==TString("raster_position_y"))
         {
     fHistograms.at(j)->Fill(raster_y_mm);
         }
       if (fHistograms.at(j)->GetTitle()==TString("bpm_3h07a_pos_x"))
         {
     fHistograms.at(j)->Fill(fbpm_3h07a_pos_x);
         }
       if (fHistograms.at(j)->GetTitle()==TString("bpm_3h07a_pos_y"))
         {
     fHistograms.at(j)->Fill(fbpm_3h07a_pos_y);
         }
       if (fHistograms.at(j)->GetTitle()==TString("bpm_3h09b_pos_x"))
         {
     fHistograms.at(j)->Fill(fbpm_3h09b_pos_x);
         }
       if (fHistograms.at(j)->GetTitle()==TString("bpm_3h09b_pos_y"))
         {
     fHistograms.at(j)->Fill(fbpm_3h09b_pos_y);
         }
     }
     
     fRateMap->Fill(raster_x_mm, raster_y_mm);
 
 }

Here is the call graph for this function:

void QwRaster::FillRawWord	(	Int_t	bank_index,
		Int_t	slot_num,
		Int_t	chan,
		UInt_t	data
	)

protected

Definition at line 691 of file QwRaster.cc.

References fModulePtrs, fPMTs, GetModuleIndex(), and kUnknownModuleType.

Referenced by ProcessEvBuffer().

 {
     Int_t modindex = GetModuleIndex(bank_index,slot_num);
 
     // std::cout<<"modtype="<< fModulePtrs.at(modindex).at(chan).first
     // <<"  chanindex="<<fModulePtrs.at(modindex).at(chan).second<<"  data="<<data<<"\n";
 
     if (modindex != -1)
     {
         EQwModuleType modtype = fModulePtrs.at(modindex).at(chan).first;
         Int_t chanindex       = fModulePtrs.at(modindex).at(chan).second;
 
         if (modtype == kUnknownModuleType || chanindex == -1)
         {
             //  This channel is not connected to anything.
             //  Do nothing.
         }
         else
         {
             fPMTs.at(modtype).at(chanindex).SetValue(data);
         }
     }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Fill the tree vector.

Reimplemented from VQwSubsystemTracking.

Definition at line 542 of file QwRaster.cc.

References fbpm_3h07a_pos_x, fbpm_3h07a_pos_y, fbpm_3h09b_pos_x, fbpm_3h09b_pos_y, fPMTs, fPositionX_ADC, fPositionY_ADC, fSCAs, VQwSubsystemTracking::fTreeArrayIndex, and VQwSubsystem::HasDataLoaded().

 {
     if (! HasDataLoaded()) return;
 
     Int_t index = fTreeArrayIndex;
     values[index++] = fPositionX_ADC;
     values[index++] = fPositionY_ADC;
     values[index++] = fbpm_3h07a_pos_x;
     values[index++] = fbpm_3h07a_pos_y;
     values[index++] = fbpm_3h09b_pos_x;
     values[index++] = fbpm_3h09b_pos_y;
 
     //fill trigvalues
     for (size_t i=0; i<fPMTs.size(); i++)
     {
       for (size_t j=0; j<fPMTs.at(i).size(); j++)
       {
         if (fPMTs.at(i).at(j).GetElementName()=="") {}
         else
         {
           values[index++] = fPMTs.at(i).at(j).GetValue();
         }
       }
     }
 
     for (size_t i=0; i<fSCAs.size(); i++) {
       if (fSCAs.at(i).GetElementName()!="") {
         values[index] = fSCAs.at(i).GetValue();
         index++;
       }
     }
 }

Here is the call graph for this function:

Int_t QwRaster::FindSignalIndex	(	const EQwModuleType	modtype,
		const TString &	name
	)		const

protected

Definition at line 733 of file QwRaster.cc.

References fPMTs.

Referenced by LoadChannelMap().

 {
     Int_t chanindex = -1;
     if (modtype < (Int_t) fPMTs.size())
     {
       for (size_t chan = 0; chan < fPMTs.at(modtype).size(); chan++)
       {
         if (name == fPMTs.at(modtype).at(chan).GetElementName())
         {
           chanindex = chan;
           break;
         }
       }
     }
     return chanindex;
 }

Here is the caller graph for this function:

Int_t QwRaster::GetDetectorInfo ( std::vector< std::vector< QwDetectorInfo > > & detector_info )

inline

Definition at line 53 of file QwRaster.h.

53 {return 0;};

void QwRaster::GetHitList ( QwHitContainer & grandHitContainer )

inlinevirtual

Get the hit list.

Implements VQwSubsystemTracking.

Definition at line 72 of file QwRaster.h.

72 {};

Int_t QwRaster::GetModuleIndex	(	size_t	bank_index,
		size_t	slot_num
	)		const

protected

Definition at line 718 of file QwRaster.cc.

References fModuleIndex.

Referenced by FillRawWord(), IsSlotRegistered(), and ProcessEvBuffer().

 {
     Int_t modindex = -1;
     // bank_index and slot_num are unsigned int and always positive
     if (/* bank_index >= 0 && */ bank_index < fModuleIndex.size())
     {
         if (/* slot_num >= 0 && */ slot_num < fModuleIndex.at(bank_index).size())
         {
             modindex = fModuleIndex.at(bank_index).at(slot_num);
         }
     }
     return modindex;
 }

Here is the caller graph for this function:

void QwRaster::InitializeChannel	(	TString	name,
		TString	datatosave
	)

Bool_t QwRaster::IsSlotRegistered	(	Int_t	bank_index,
		Int_t	slot_num
	)		const

inlineprotected

Definition at line 122 of file QwRaster.h.

References GetModuleIndex().

Referenced by ProcessEvBuffer().

                                                                   {
     return (GetModuleIndex(bank_index,slot_num) != -1);
   };

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwRaster::LinkChannelToSignal	(	const UInt_t	chan,
		const TString &	name
	)

protected

Definition at line 682 of file QwRaster.cc.

References fCurrentIndex, fCurrentType, fModulePtrs, and fPMTs.

Referenced by LoadChannelMap().

 {
     fPMTs.at(fCurrentType).push_back(QwPMT_Channel(name));
     fModulePtrs.at(fCurrentIndex).at(chan).first  = fCurrentType;
     fModulePtrs.at(fCurrentIndex).at(chan).second =
     fPMTs.at(fCurrentType).size() -1;
     return 0;
 }

Here is the caller graph for this function:

Int_t QwRaster::LoadChannelMap ( TString mapfile )

virtual

Mandatory map file definition.

Implements VQwSubsystem.

Definition at line 44 of file QwRaster.cc.

References fBankID, fCurrentIndex, fCurrentType, VQwSubsystem::fDetectorMaps, FindSignalIndex(), fModulePtrs, fSCAs, fSCAs_map, fSCAs_offset, QwParameterFile::GetUInt(), kUnknownModuleType, LinkChannelToSignal(), RegisterModuleType(), RegisterROCNumber(), RegisterSlotNumber(), and RegisterSubbank().

 {
     TString varname, varvalue;
     TString modtype, dettype, name;
     Int_t modnum = 0, channum = 0;
 
     QwParameterFile mapstr(mapfile.Data());  //Open the file
     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") {
       RegisterROCNumber(value);
           } else if (varname=="qdc_bank") {
         RegisterSubbank(value);
               fBankID[0] = value;
           } else if (varname=="sca_bank") {
               fBankID[1] = value;
               RegisterSubbank(value);
           } else if (varname=="slot") {
               RegisterSlotNumber(value);
           } else if (varname=="module") {
         RegisterModuleType(varvalue);
           }
         } else {
             //  Break this line into tokens to process it.
             modtype   = mapstr.GetTypedNextToken<TString>();
             modnum    = mapstr.GetTypedNextToken<Int_t>();
             channum   = mapstr.GetTypedNextToken<Int_t>();
             dettype   = mapstr.GetTypedNextToken<TString>();
             name      = mapstr.GetTypedNextToken<TString>();
 
             //  Push a new record into the element array
             if (modtype=="SIS3801") {
               QwSIS3801D24_Channel localchannel(name);
               localchannel.SetNeedsExternalClock(kFALSE);
               fSCAs.push_back(localchannel);
               fSCAs_map[name] = fSCAs.size()-1;
               Int_t offset = QwSIS3801D24_Channel::GetBufferOffset(modnum,channum);
               fSCAs_offset.push_back(offset);
             } else if (modtype=="V792" || modtype=="V775") {
                 RegisterModuleType(modtype);
                 //  Check to see if we've encountered this channel or name yet
                 if (fModulePtrs.at(fCurrentIndex).at(channum).first != kUnknownModuleType) {
                  //  We've seen this channel
                 } else if (FindSignalIndex(fCurrentType, name)>=0) {
                     //  We've seen this signal
                 } else {
                   //  If not, push a new record into the element array
                   LinkChannelToSignal(channum, name);
                 }
             } else {
                 std::cerr << "LoadChannelMap:  Unknown line: " << mapstr.GetLine().c_str() << std::endl;
            }
        }
   }
     mapstr.Close(); // Close the file (ifstream)
   //ReportConfiguration();
   return 0;
 }

Here is the call graph for this function:

Int_t QwRaster::LoadEventCuts ( TString mapfile )

inlinevirtual

Optional event cut file.

Reimplemented from VQwSubsystem.

Definition at line 46 of file QwRaster.h.

46 { return 0; };

Int_t QwRaster::LoadInputParameters ( TString mapfile )

virtual

Mandatory parameter file definition.

Implements VQwSubsystem.

Definition at line 113 of file QwRaster.cc.

References fbpm_3h07a_intercept, fbpm_3h07a_slope, fbpm_3h09b_intercept, fbpm_3h09b_slope, fCal_Factor_QDC_X, fCal_Factor_QDC_Y, fChannel_Offset_X, fChannel_Offset_Y, VQwSubsystem::fDetectorMaps, fPositionOffsetX, and fPositionOffsetY.

 {
     Bool_t ldebug=kFALSE;
     TString varname, varvalue;
     Double_t varped;
     Double_t varcal;
     TString localname;
 
     Int_t lineread=0;
 
     QwParameterFile mapstr(parameterfile.Data());  //Open the file
     fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
     if (ldebug) std::cout<<"\nReading raster parameter file: "<<parameterfile<<"\n";
 
     while (mapstr.ReadNextLine())
     {
         lineread+=1;
         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))
         {
             varname.ToLower();
             Double_t value = atof(varvalue.Data());
             if (varname=="position_offset_x")
             {
                 fPositionOffsetX = value;
             }
             else if (varname=="position_offset_y")
             {
                 fPositionOffsetY = value;
             }
             else if (varname=="cal_factor_qdc_x")
             {
                 fCal_Factor_QDC_X = value;
             }
             else if (varname=="cal_factor_qdc_y")
             {
                 fCal_Factor_QDC_Y = value;
             }
             else if (varname=="channel_offset_x")
             {
                 fChannel_Offset_X = value;
             }
             else if (varname=="channel_offset_y")
             {
                 fChannel_Offset_Y = value;
             }
             else if (varname=="bpm_3h07a_slope")
             {
                 fbpm_3h07a_slope = value;
             }
             else if (varname=="bpm_3h07a_intercept")
             {
                 fbpm_3h07a_intercept = value;
             }
             else if (varname=="bpm_3h09b_slope")
             {
                 fbpm_3h09b_slope = value;
             }
             else if (varname=="bpm_3h09b_intercept")
             {
                 fbpm_3h09b_intercept = value;
             }
 
             if (ldebug) std::cout<<"inputs for "<<varname<<": "<<value<<"\n";
         }
 
         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";
         }
     }
     if (ldebug) std::cout<<" line read in the parameter file ="<<lineread<<" \n";
 
     ldebug=kFALSE;
     mapstr.Close(); // Close the file (ifstream)
     return 0;
 }

void QwRaster::PrintErrorCounters ( ) const

inline

Definition at line 48 of file QwRaster.h.

48 {};

void QwRaster::PrintInfo ( )

Definition at line 750 of file QwRaster.cc.

References fPMTs, and VQwSubsystem::fSystemName.

 {
     std::cout << "QwRaster: " << fSystemName << std::endl;
 
     for (size_t i=0; i<fPMTs.size(); i++)
     {
         for (size_t j=0; j<fPMTs.at(i).size(); j++)
         {
             fPMTs.at(i).at(j).PrintInfo();
         }
     }
 
     return;
 }

void QwRaster::PrintValue ( )

inline

Definition at line 75 of file QwRaster.h.

75 { };

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

virtual

Implements VQwSubsystem.

Definition at line 229 of file QwRaster.cc.

 {
     return 0;
 }

Int_t QwRaster::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 235 of file QwRaster.cc.

References MQwV775TDC::DecodeTDCWord(), Qw::e, fBankID, FillRawWord(), fModulePtrs, fQDCTDC, fSCAs, fSCAs_offset, GetModuleIndex(), VQwSubsystem::GetSubbankIndex(), MQwV775TDC::GetTDCChannelNumber(), MQwV775TDC::GetTDCData(), MQwV775TDC::GetTDCSlotNumber(), IsSlotRegistered(), MQwV775TDC::IsValidDataword(), and VQwSubsystem::SetDataLoaded().

 {
     Int_t index = 0;
 
     index = GetSubbankIndex(roc_id,bank_id);
 
     //  This is a QADC/TDC bank
     if (bank_id==fBankID[0])
     {
         if (index>=0 && num_words>0)
         {
             //  We want to process this ROC.  Begin looping through the data.
             SetDataLoaded(kTRUE);
             for (size_t i=0; i<num_words ; i++)
             {
                 //  Decode this word as a V775TDC word.
                 fQDCTDC.DecodeTDCWord(buffer[i]);
 
                 if (! IsSlotRegistered(index, fQDCTDC.GetTDCSlotNumber())) continue;
 
                 if (fQDCTDC.IsValidDataword())
                 {
                     // This is a V775 TDC data word
                     try
                     {
                         FillRawWord(index,fQDCTDC.GetTDCSlotNumber(),fQDCTDC.GetTDCChannelNumber(),
                                           fQDCTDC.GetTDCData());
                     }
                     catch (std::exception& e)
                     {
                         std::cerr << "Standard exception from Raster::FillRawTDCWord: "
                         << e.what() << std::endl;
                         Int_t chan = fQDCTDC.GetTDCChannelNumber();
                         std::cerr << "   Parameters:  index=="<<index
                         << "; GetV775SlotNumber()=="<<fQDCTDC.GetTDCSlotNumber()
                         << "; GetV775ChannelNumber()=="<<chan
                         << "; GetV775Data()=="<<fQDCTDC.GetTDCData()
                         << std::endl;
                         Int_t modindex = GetModuleIndex(index, fQDCTDC.GetTDCSlotNumber());
                         std::cerr << "   GetModuleIndex()=="<<modindex
                         << "; fModulePtrs.at(modindex).size()=="
                         << fModulePtrs.at(modindex).size()
                         << "; fModulePtrs.at(modindex).at(chan).first {module type}=="
                         << fModulePtrs.at(modindex).at(chan).first
                         << "; fModulePtrs.at(modindex).at(chan).second {signal index}=="
                         << fModulePtrs.at(modindex).at(chan).second
                         << std::endl;
                     }
                 }
             }
         }
     }
 
     else if (bank_id==fBankID[1]) { // SIS Scalar
 
       // Check if scaler buffer contains more than one event
       if (buffer[0]/32!=1) return 0;
 
       if (index>=0 && num_words>0) {
         SetDataLoaded(kTRUE);
         UInt_t words_read = 0;
         for (size_t i=0; i<fSCAs.size(); i++) {
           words_read += fSCAs.at(i).ProcessEvBuffer(&(buffer[fSCAs_offset.at(i)]),
                                                     num_words-fSCAs_offset.at(i));
         }
       }
     }
 
     return 0;
 }

Here is the call graph for this function:

void QwRaster::ProcessEvent ( )

virtual

Implements VQwSubsystem.

Definition at line 307 of file QwRaster.cc.

References fbpm_3h07a_intercept, fbpm_3h07a_pos_x, fbpm_3h07a_pos_y, fbpm_3h07a_slope, fbpm_3h07a_xm, fbpm_3h07a_xp, fbpm_3h07a_ym, fbpm_3h07a_yp, fbpm_3h09b_intercept, fbpm_3h09b_pos_x, fbpm_3h09b_pos_y, fbpm_3h09b_slope, fbpm_3h09b_xm, fbpm_3h09b_xp, fbpm_3h09b_ym, fbpm_3h09b_yp, fCal_Factor_QDC_X, fCal_Factor_QDC_Y, fChannel_Offset_X, fChannel_Offset_Y, fPMTs, fPositionOffsetX, fPositionOffsetY, fPositionX_ADC, fPositionY_ADC, fSCAs, and VQwSubsystem::HasDataLoaded().

 {
     if (! HasDataLoaded()) return;
 
     for (size_t i=0; i<fPMTs.size(); i++) {
       for (size_t j=0; j<fPMTs.at(i).size(); j++) {
         fPMTs.at(i).at(j).ProcessEvent();
       }
     }
 
     /*for (size_t i=0; i<fADC_Data.size(); i++)
     {
         if (fADC_Data.at(i) != NULL)
         {
             fADC_Data.at(i)->ProcessEvent();
         }
     }*/
 
     for (size_t i=0; i<fSCAs.size(); i++) {
       fSCAs.at(i).ProcessEvent();
     }
 
     //Fill trigger data
     for (size_t i=0; i<fPMTs.size(); i++)
     {
         for (size_t j=0; j<fPMTs.at(i).size(); j++)
         {
 
             TString element_name = fPMTs.at(i).at(j).GetElementName();
             if (element_name==TString("raster_posx_adc"))
             {
                 fPositionX_ADC = fPMTs.at(i).at(j).GetValue();
                 fPositionX_ADC = (fPositionX_ADC-fChannel_Offset_X)*fCal_Factor_QDC_X + fPositionOffsetX;
             }
             else if (element_name==TString("raster_posy_adc"))
             {
                 fPositionY_ADC = fPMTs.at(i).at(j).GetValue();
                 fPositionY_ADC = (fPositionY_ADC-fChannel_Offset_Y)*fCal_Factor_QDC_Y + fPositionOffsetY;
             }
 
             else if (element_name==TString("bpm_3h07a_xp_adc"))
             {
                 fbpm_3h07a_xp = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h07a_xm_adc"))
             {
                 fbpm_3h07a_xm = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h07a_yp_adc"))
             {
                 fbpm_3h07a_yp = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h07a_ym_adc"))
             {
                 fbpm_3h07a_ym = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h09b_xp_adc"))
             {
                 fbpm_3h09b_xp = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h09b_xm_adc"))
             {
                 fbpm_3h09b_xm = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h09b_yp_adc"))
             {
                 fbpm_3h09b_yp = fPMTs.at(i).at(j).GetValue();
             }
             else if (element_name==TString("bpm_3h09b_ym_adc"))
             {
                 fbpm_3h09b_ym = fPMTs.at(i).at(j).GetValue();
             }
         }
      }
      
      fbpm_3h07a_pos_x = fbpm_3h07a_slope*(fbpm_3h07a_xp-fbpm_3h07a_xm)/
                         (fbpm_3h07a_xp+fbpm_3h07a_xm)+fbpm_3h07a_intercept;
      fbpm_3h07a_pos_y = fbpm_3h07a_slope*(fbpm_3h07a_yp-fbpm_3h07a_ym)/
                         (fbpm_3h07a_yp+fbpm_3h07a_ym)+fbpm_3h07a_intercept;
      fbpm_3h09b_pos_x = fbpm_3h09b_slope*(fbpm_3h09b_xp-fbpm_3h09b_xm)/
                         (fbpm_3h09b_xp+fbpm_3h09b_xm)+fbpm_3h09b_intercept;
      fbpm_3h09b_pos_y = fbpm_3h09b_slope*(fbpm_3h09b_yp-fbpm_3h09b_ym)/
                         (fbpm_3h09b_yp+fbpm_3h09b_ym)+fbpm_3h09b_intercept;
 
 }

Here is the call graph for this function:

void QwRaster::ProcessOptions ( QwOptions & options )

virtual

Process the command line options.

Reimplemented from VQwSubsystem.

Definition at line 39 of file QwRaster.cc.

 {
     //Handle command line options
 }

EQwModuleType QwRaster::RegisterModuleType ( TString moduletype )

protected

Definition at line 662 of file QwRaster.cc.

References fCurrentIndex, fCurrentType, fModuleTypes, fPMTs, kSIS3801, and kV792_ADC.

Referenced by LoadChannelMap().

 {
     moduletype.ToUpper();
 
     //  Check to see if we've already registered a type for the current slot,
     //  if so, throw an error...
 
     if (moduletype=="V792"){
       fCurrentType = kV792_ADC;
     } else if (moduletype=="SIS3801") {
       fCurrentType = kSIS3801;
     }
     fModuleTypes.at(fCurrentIndex) = fCurrentType;
     if ((Int_t) fPMTs.size() <= fCurrentType){
       fPMTs.resize(fCurrentType+1);
     }
     return fCurrentType;
 }

Here is the caller graph for this function:

Int_t QwRaster::RegisterROCNumber ( const UInt_t roc_id )

protected

Definition at line 612 of file QwRaster.cc.

References fCurrentBankIndex, fModuleIndex, VQwSubsystem::GetSubbankIndex(), kMaxNumberOfModulesPerROC, and VQwSubsystem::RegisterROCNumber().

Referenced by LoadChannelMap().

 {
     VQwSubsystem::RegisterROCNumber(roc_id, 0);
     fCurrentBankIndex = GetSubbankIndex(roc_id, 0);
     std::vector<Int_t> tmpvec(kMaxNumberOfModulesPerROC,-1);
     fModuleIndex.push_back(tmpvec);
     return fCurrentBankIndex;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwRaster::RegisterSlotNumber ( const UInt_t slot_id )

protected

Definition at line 632 of file QwRaster.cc.

References fCurrentBankIndex, fCurrentIndex, fCurrentSlot, fModuleIndex, fModulePtrs, fModuleTypes, fNumberOfModules, kMaxNumberOfChannelsPerModule, kMaxNumberOfModulesPerROC, and kUnknownModuleType.

Referenced by LoadChannelMap().

 {
 
     std::pair<EQwModuleType, Int_t> tmppair;
     tmppair.first  = kUnknownModuleType;
     tmppair.second = -1;
     if (slot_id<kMaxNumberOfModulesPerROC)
     {
         // fCurrentBankIndex is unsigned int and always positive
         if (/* fCurrentBankIndex >= 0 && */ fCurrentBankIndex <= fModuleIndex.size())
         {
             fModuleTypes.resize(fNumberOfModules+1);
             fModulePtrs.resize(fNumberOfModules+1);
             fModulePtrs.at(fNumberOfModules).resize(kMaxNumberOfChannelsPerModule,
                                                     tmppair);
             fNumberOfModules = fModulePtrs.size();
             fModuleIndex.at(fCurrentBankIndex).at(slot_id) = fNumberOfModules-1;
             fCurrentSlot     = slot_id;
             fCurrentIndex = fNumberOfModules-1;
         }
     }
     else
     {
         std::cerr << "QwRaster::RegisterSlotNumber:  Slot number "
         << slot_id << " is larger than the number of slots per ROC, "
         << kMaxNumberOfModulesPerROC << std::endl;
     }
     return fCurrentIndex;
 }

Here is the caller graph for this function:

Int_t QwRaster::RegisterSubbank ( const UInt_t bank_id )

protected

Definition at line 621 of file QwRaster.cc.

References fCurrentBankIndex, fModuleIndex, kMaxNumberOfModulesPerROC, and VQwSubsystem::RegisterSubbank().

Referenced by LoadChannelMap().

 {
     Int_t stat = VQwSubsystem::RegisterSubbank(bank_id);
     fCurrentBankIndex++;
     std::vector<Int_t> tmpvec(kMaxNumberOfModulesPerROC,-1);
     fModuleIndex.push_back(tmpvec);
     //std::cout<<"RASTER Register Subbank "<<bank_id<<" with BankIndex "<<fCurrentBankIndex<<std::endl;
     return stat;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwRaster::ReportConfiguration ( )

Definition at line 576 of file QwRaster.cc.

 {
     /*std::cout << "Configuration of raster:"<< std::endl;
     for (size_t i = 0; i<fROC_IDs.size(); i++)
     {
         for (size_t j=0; j<fBank_IDs.at(i).size(); j++)
         {
 
             Int_t ind = GetSubbankIndex(fROC_IDs.at(i),fBank_IDs.at(i).at(j));
             std::cout << "ROC " << fROC_IDs.at(i)
             << ", subbank 0x" << std::hex<<fBank_IDs.at(i).at(j)<<std::dec
             << ":  subbank index==" << ind << std::endl;
 
                 for (size_t k=0; k<kMaxNumberOfModulesPerROC; k++)
                 {
                     Int_t QadcTdcindex = GetModuleIndex(ind,k);
                     if (QadcTdcindex != -1)
                     {
                         std::cout << "    Slot "<<k<<"  Module#"<<QadcTdcindex<< std::endl;
                     }
                 }
         }
     }*/
 
 }