#include <QwSciFiDetector.h>

Inheritance diagram for QwSciFiDetector:

Collaboration diagram for QwSciFiDetector:

Public Member Functions
	QwSciFiDetector (const TString &name)

virtual	~QwSciFiDetector ()

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

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

void	FillListOfHits (QwHitContainer &hitlist)

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

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

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

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

void	PrintConfigurationBuffer (UInt_t *buffer, UInt_t num_words)

void	ReportConfiguration (Bool_t verbose)

void	ProcessEvent ()

void	FillRawTDCWord (Int_t bank_index, Int_t slot_num, Int_t chan, UInt_t data)

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

void	ClearEventData ()

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

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

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

TString	GetSubsystemName () const

Bool_t	HasDataLoaded () const

void	SetParent (QwSubsystemArray *parent)
	Set the parent of this subsystem to the specified array. More...

QwSubsystemArray *	GetParent (const unsigned int parent=0) const
	Get the parent of this subsystem. More...

VQwSubsystem *	GetSibling (const std::string &name) const
	Get the sibling with specified name. More...

Bool_t	PublishInternalValue (const TString &name, const TString &desc, const VQwHardwareChannel *value) const
	Publish a variable name to the parent subsystem array. More...

virtual Bool_t	PublishInternalValues () const
	Publish all variables of the subsystem. More...

virtual Bool_t	PublishByRequest (TString device_name)
	Try to publish an internal variable matching the submitted name. More...

Bool_t	RequestExternalValue (const TString &name, VQwHardwareChannel *value) const
	Request a named value which is owned by an external subsystem; the request will be handled by the parent subsystem array. More...

virtual const VQwHardwareChannel *	ReturnInternalValue (const TString &name) const
	Return a pointer to a varialbe to the parent subsystem array to be delivered to a different subsystem. More...

virtual Bool_t	ReturnInternalValue (const TString &name, VQwHardwareChannel *value) const
	Return a named value to the parent subsystem array to be delivered to a different subsystem. More...

virtual std::vector< TString >	GetParamFileNameList ()

virtual std::map< TString, TString >	GetDetectorMaps ()

virtual Int_t	LoadDetectorMaps (QwParameterFile &file)
	Parse parameter file to find the map files. More...

virtual Int_t	LoadEventCuts (TString mapfile)
	Optional event cut file. 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, QwSciFiDetector >
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	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)

Int_t	GetModuleIndex (size_t bank_index, size_t slot_num) const

Bool_t	IsSlotRegistered (Int_t bank_index, Int_t slot_num) const

void	SubtractReferenceTimes ()

void	UpdateHits ()

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
TString	kRegion

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

Int_t	fCurrentSlot

Int_t	fCurrentModuleIndex

UInt_t	kMaxNumberOfChannelsPerF1TDC

Int_t	kNumberOfVMEModules

MQwF1TDC	fF1TDCDecoder

QwF1TDContainer *	fF1TDContainer

F1TDCReferenceContainer *	fF1RefContainer

std::vector< QwHit >	fTDCHits

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

std::vector< std::vector < QwDetectorID > >	fDetectorIDs

std::vector< std::pair< Int_t, Int_t > >	fReferenceChannels

std::vector< std::vector < Double_t > >	fReferenceData

std::vector< QwSIS3801D24_Channel >	fSCAs

std::map< TString, size_t >	fSCAs_map

std::vector< Int_t >	fSCAs_offset

Int_t	fScalerBankIndex

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	kMaxNumberOfSlotsPerROC = 21

static const Int_t	kF1ReferenceChannelNumber = 99

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

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

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

Detailed Description

Definition at line 33 of file QwSciFiDetector.h.

Constructor & Destructor Documentation

QwSciFiDetector::QwSciFiDetector ( )

private

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

QwSciFiDetector::QwSciFiDetector ( const TString & name )

Definition at line 20 of file QwSciFiDetector.cc.

References ClearAllBankRegistrations(), fCurrentBankIndex, fCurrentModuleIndex, fCurrentSlot, fF1RefContainer, fF1TDCDecoder, fF1TDContainer, fScalerBankIndex, QwF1TDContainer::GetF1TDCDecoder(), MQwF1TDC::GetTDCMaxChannels(), kMaxNumberOfChannelsPerF1TDC, kNumberOfVMEModules, kRegion, and VQwSubsystem::SetEventTypeMask().

 : VQwSubsystem(name),
   VQwSubsystemTracking(name)
 {
   // accept all event types to transfer to ProcessEvBuffer(roc_id .....)
   // ProcessEvBuffer(eventtype...) is defined in VQwSubsystem.h uses ProcessEvBuffer(roc_id,...)
   // 
   SetEventTypeMask(0xffff); 
   
   ClearAllBankRegistrations();
 
   kRegion             = "";
   fCurrentBankIndex   = 0;
   fCurrentSlot        = 0;
   fCurrentModuleIndex = 0;
   kNumberOfVMEModules = 0;
   fScalerBankIndex    = -1;
 
   fF1TDContainer = new QwF1TDContainer();
   fF1TDCDecoder  = fF1TDContainer->GetF1TDCDecoder();
   kMaxNumberOfChannelsPerF1TDC = fF1TDCDecoder.GetTDCMaxChannels(); 
   fF1RefContainer = new F1TDCReferenceContainer();
 }

Here is the call graph for this function:

QwSciFiDetector::~QwSciFiDetector ( )

virtual

Definition at line 44 of file QwSciFiDetector.cc.

References fF1RefContainer, fF1TDContainer, and fSCAs.

 {
   fSCAs.clear();
   delete fF1TDContainer;
   delete fF1RefContainer;
 };

Member Function Documentation

void QwSciFiDetector::ClearAllBankRegistrations ( )

protected

Definition at line 883 of file QwSciFiDetector.cc.

References VQwSubsystem::ClearAllBankRegistrations(), fDetectorIDs, fModuleIndex, fTDCHits, and kNumberOfVMEModules.

Referenced by QwSciFiDetector().

 {
   VQwSubsystemTracking::ClearAllBankRegistrations();
 
   std::size_t i = 0;
 
   for ( i=0; i<fModuleIndex.size(); i++)
     {
       fModuleIndex.at(i).clear();
     }
   
   fModuleIndex.clear();
   fDetectorIDs.clear();
   fTDCHits.clear();
   kNumberOfVMEModules = 0;
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSciFiDetector::ClearEventData ( )

virtual

Implements VQwSubsystem.

Definition at line 236 of file QwSciFiDetector.cc.

References F1TDCReferenceContainer::ClearEventData(), fF1RefContainer, fReferenceData, fSCAs, fTDCHits, and VQwSubsystem::SetDataLoaded().

 { 
   SetDataLoaded(kFALSE);
 
   fTDCHits.clear();
   std::size_t i = 0;
   
   fF1RefContainer->ClearEventData();
 
   for (i=0; i<fReferenceData.size(); i++) {
     fReferenceData.at(i).clear();
   }
   for (i=0; i<fSCAs.size(); i++) {
     fSCAs.at(i).ClearEventData();
   }
   return;
 };

Here is the call graph for this function:

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

virtual

Construct the branch and tree vector.

Reimplemented from VQwSubsystemTracking.

Definition at line 312 of file QwSciFiDetector.cc.

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

 {
   fTreeArrayIndex = values.size();
 
   TString basename;
   if (prefix=="") basename = "scifiber";
   else            basename = prefix;
 
   TString list = "";
 
   for (std::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 QwSciFiDetector::ConstructHistograms	(	TDirectory *	folder,
		TString &	prefix
	)

virtual

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

Implements VQwSubsystem.

Definition at line 282 of file QwSciFiDetector.cc.

References fSCAs.

 {
   
   for (std::size_t i=0; i<fSCAs.size(); i++) {
     fSCAs.at(i).ConstructHistograms(folder, prefix);
   }
   return;
 };

void QwSciFiDetector::DeleteHistograms ( )

Definition at line 303 of file QwSciFiDetector.cc.

 {
 
   return;
 };

void QwSciFiDetector::FillHardwareErrorSummary ( )

virtual

Hardware error summary.

Reimplemented from VQwSubsystemTracking.

Definition at line 870 of file QwSciFiDetector.cc.

References fF1RefContainer, and fF1TDContainer.

 {
   fF1RefContainer-> PrintCounters();
   fF1TDContainer -> PrintErrorSummary();
   fF1TDContainer -> WriteErrorSummary();
   return;
 };

void QwSciFiDetector::FillHistograms ( )

virtual

Fill the histograms for this subsystem.

Implements VQwSubsystem.

Definition at line 292 of file QwSciFiDetector.cc.

References fSCAs, and VQwSubsystem::HasDataLoaded().

 {
   if (! HasDataLoaded()) return;
   for (std::size_t i=0; i<fSCAs.size(); i++) {
     fSCAs.at(i).FillHistograms();
   }
   
   return;
 };

Here is the call graph for this function:

void QwSciFiDetector::FillListOfHits ( QwHitContainer & hitlist )

Definition at line 276 of file QwSciFiDetector.cc.

 {
   return;
 };

void QwSciFiDetector::FillRawTDCWord	(	Int_t	bank_index,
		Int_t	slot_num,
		Int_t	chan,
		UInt_t	data
	)

Definition at line 771 of file QwSciFiDetector.cc.

References fDetectorIDs, fF1RefContainer, fReferenceData, fTDCHits, GetModuleIndex(), kDirectionNull, kF1ReferenceChannelNumber, kRegionID1, F1TDCReferenceContainer::SetReferenceSignal(), and QwHit::WireMatches().

Referenced by ProcessEvBuffer().

 {
   Bool_t local_debug = false;
   Int_t tdcindex     =  0;
 
   tdcindex = GetModuleIndex(bank_index,slot_num);
 
   // I think, it is not necessary, but is just "safe" double check routine.
   if (tdcindex not_eq -1) {
 
     Int_t hitcnt  = 0;
 
     EQwDetectorPackage package = kPackageNull;
     EQwDirectionID direction   = kDirectionNull;
 
     Int_t   octant  = 0;
     Int_t   plane   = 0; // fiber number (1357, 2,4,6,8)  for the fibers,  and  99 for the reference channel
     Int_t   element = 0; // fiber orient type (0 (both), 1(a), 2(b)) for the fibers, and  the module index for the reference channel
     TString name         = "";
  
     //    F1TDCReferenceSignal *f1_ref_signal;
     
     fF1RefContainer->SetReferenceSignal(bank_index, slot_num, chan, data, local_debug);
 
     // if( f1_ref_signal ) {
     //   f1_ref_signal -> SetRefTimeAU (data);
     //   //  std::cout << *f1_ref_signal << std::endl;
     // }
     
     plane   = fDetectorIDs.at(tdcindex).at(chan).fPlane;
     element = fDetectorIDs.at(tdcindex).at(chan).fElement;
     package = fDetectorIDs.at(tdcindex).at(chan).fPackage;
     octant  = fDetectorIDs.at(tdcindex).at(chan).fOctant;
     if(local_debug) {
       printf("bank_idx %d, slot %d, plane %d, element %d, package %d\n",
              bank_index, slot_num, (Int_t) plane, (Int_t) element, (Int_t) package);
     }
 
     if (plane == -1 or element == -1){
       //  This channel is not connected to anything. Do nothing.
     }
     else if (plane == kF1ReferenceChannelNumber){
       fReferenceData.at(element).push_back(data);
       // we assign the reference time into fReferenceData according the module index
       // See LocalChannelMap
  
     }
     else {
       direction = fDetectorIDs.at(tdcindex).at(chan).fDirection;
 
       hitcnt = std::count_if(fTDCHits.begin(), fTDCHits.end(), 
            boost::bind(
            &QwHit::WireMatches, _1, kRegionID1, boost::ref(package), boost::ref(plane), boost::ref(element)
            ) 
            );
    
       fTDCHits.push_back(
            QwHit(
                  bank_index, 
                  slot_num, 
                  chan, 
                                hitcnt, 
                                kRegionID1, 
                                package, 
                                octant,
                                plane,
                                direction, 
                                element, 
                  data
                  )
            );
 
       if(local_debug) {
   std::cout << "At QwSciFiDetector::FillRawTDCWord " 
       << " bank index " << bank_index
       << " slot num " << slot_num
       << " chan num " << chan
       << " hitcnt " << hitcnt
       << " plane  " << plane
       << " wire   " << element
       << " package " << package 
       << " diection " << direction
       << " fTDCHits.size() " <<  fTDCHits.size() 
       << std::endl;
       }
 
     }
 
   } // (tdcindex not_eq -1) {
   
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Fill the tree vector.

Reimplemented from VQwSubsystemTracking.

Definition at line 342 of file QwSciFiDetector.cc.

References fSCAs, VQwSubsystemTracking::fTreeArrayIndex, and VQwSubsystem::HasDataLoaded().

 {
   if (! HasDataLoaded()) return;
 
   Int_t index = fTreeArrayIndex;
 
   for (std::size_t i=0; i<fSCAs.size(); i++) 
     {
       if (fSCAs.at(i).GetElementName() != "") {
   values[index] = fSCAs.at(i).GetValue();
   index++;
       }
     }
   return;
 }

Here is the call graph for this function:

void QwSciFiDetector::GetHitList ( QwHitContainer & grandHitContainer )

inlinevirtual

Get the hit list.

Implements VQwSubsystemTracking.

Definition at line 78 of file QwSciFiDetector.h.

References QwHitContainer::Append(), and fTDCHits.

                                                      {
     grandHitContainer.Append(fTDCHits);
   };

Here is the call graph for this function:

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

protected

Definition at line 978 of file QwSciFiDetector.cc.

References fModuleIndex.

Referenced by FillRawTDCWord(), IsSlotRegistered(), ProcessConfigurationBuffer(), ProcessEvBuffer(), and ReportConfiguration().

 {
   Int_t module_index = -1;
   if ( bank_index < fModuleIndex.size() ) {
     if ( slot_num < fModuleIndex.at(bank_index).size() ) {
       module_index = fModuleIndex.at(bank_index).at(slot_num);
     }
   }
   return module_index;
 }

Here is the caller graph for this function:

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

inlineprotected

Definition at line 112 of file QwSciFiDetector.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 QwSciFiDetector::LoadChannelMap ( TString mapfile )

virtual

Mandatory map file definition.

Implements VQwSubsystem.

Definition at line 53 of file QwSciFiDetector.cc.

References fCurrentBankIndex, fCurrentModuleIndex, fDetectorIDs, VQwSubsystem::fDetectorMaps, fF1RefContainer, fModuleIndex, fReferenceChannels, fScalerBankIndex, fSCAs, fSCAs_map, fSCAs_offset, QwParameterFile::GetUInt(), kDirectionNull, kRegionID1, RegisterROCNumber(), RegisterSlotNumber(), RegisterSubbank(), and ReportConfiguration().

 {
   
   Bool_t local_debug = false;
 
   if(local_debug) printf("\n------------- R1 LoadChannelMap %s\n\n", mapfile.Data());
 
   TString varname   = "";
   TString varvalue  = "";
   UInt_t  value     = 0;
 
   TString modtype      = "";
   Int_t   slot_number  = 0; // for SIS3801 scaler. In case of F1TDC, it is a dummy slot number.
   Int_t   chan_number  = 0; 
   Int_t   fiber_number = 0;
   TString name         = "";
   Int_t reference_counter = 0;
 
   EQwDetectorPackage package = kPackageNull;
   EQwDirectionID   direction = kDirectionNull;
 
 
   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();
       value = QwParameterFile::GetUInt(varvalue);
       if (value ==0){
   value = atol(varvalue.Data());
       }
       if (varname=="roc") {
   RegisterROCNumber(value);
       } 
       else if (varname=="bank"){
         RegisterSubbank(value);
       }
       else if (varname=="pkg") {
         package=(EQwDetectorPackage)value;
       }
       else if (varname=="slot") {
         RegisterSlotNumber(value);
       } 
     }
     else {
 
       modtype      = mapstr.GetTypedNextToken<TString>();
       slot_number  = mapstr.GetTypedNextToken<Int_t>();
       chan_number  = mapstr.GetTypedNextToken<Int_t>();
       fiber_number = mapstr.GetTypedNextToken<Int_t>();
       name         = mapstr.GetTypedNextToken<TString>();
 
       if(local_debug) {
   printf("Modtype  %8s Slot_Number %4d  ChanN %4d, FiberN %4d, FiberName %s\n", modtype.Data(), slot_number, chan_number, fiber_number, name.Data());
       }
       if (modtype=="F1TDC") {
 
   if (fiber_number==99) {
 
     fF1RefContainer -> AddF1TDCReferenceSignal(new F1TDCReferenceSignal(fCurrentBankIndex, slot_number, chan_number, name));
 
     if (name=="MasterTrigger" ) {
       //    fReferenceChannels.push_back(std::make_pair(fCurrentModuleIndex, chan_number));
       //      printf("bank index %d Chan %d reference_counter %d\n", fCurrentBankIndex, chan_number, reference_counter);
       fReferenceChannels.at ( fCurrentBankIndex ).first  = fCurrentModuleIndex;
       fReferenceChannels.at ( fCurrentBankIndex ).second = chan_number;
       
       fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fElement   = fCurrentBankIndex;
       reference_counter++;
     }
   }
   else {
     fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fElement   = name.Atoi();
   }
   
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fRegion    = kRegionID1;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fPackage   = package;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fPlane     = fiber_number;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fDirection = direction;
 
       }
       else if (modtype == "SIS3801") {
   
   // We have *one* SIS3801 module, and "8" SIS3801 channels,
   // So, fSCAs size is 8
 
   fScalerBankIndex = fCurrentBankIndex;
   //  printf("Scaler bank index %d\n", fScalerBankIndex);
   QwSIS3801D24_Channel localchannel(name);
   localchannel.SetNeedsExternalClock(kFALSE);
   fSCAs.push_back(localchannel);
   fSCAs_map[name] = fSCAs.size()-1;
   fSCAs_offset.push_back( QwSIS3801D24_Channel::GetBufferOffset(slot_number,chan_number) );
 
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fRegion    = kRegionID1;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fPackage   = package;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fPlane     = fiber_number;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fDirection = direction;
   fDetectorIDs.at(fCurrentModuleIndex).at(chan_number).fElement   = -1;
 
       }
       
    
       
     } 
 
   }
   
   if(local_debug) {
   
   
     Int_t unused_size_counter = 0;
     std::size_t i = 0;
     
     
 
     for( i = 0; i < fModuleIndex.size(); i++) 
       {
   for(size_t slot_size =0; slot_size < fModuleIndex.at(i).size(); slot_size++)
     {
       Int_t m_idx = 0;
       m_idx = fModuleIndex.at(i).at(slot_size);
       if(m_idx != -1 ) {
         std::cout << "[" << i <<","<< slot_size << "] "
       << " module index " << fModuleIndex.at(i).at(slot_size)
       << std::endl;
       }
       else {
         unused_size_counter++;
       }
     }
       }
     // why the bank index is always odd number?
     // 
 
     printf("Total unused size of fModuleIndex vector %6d\n", unused_size_counter);
     
     for(i = 0; i < fReferenceChannels.size(); i++) 
       {
   std::cout << "[" << i <<"] "
       << " fRerenceChannel " << fReferenceChannels.at(i).first
       << " " << fReferenceChannels.at(i).second
       << std::endl;
       }
 
 
     // for(i = 0; i < fReferenceData.size(); i++) 
     //   {
     //  for(std::size_t j=0; j< fReferenceData.at(i).size(); j++) 
     //    {
     //      std::cout << "[" << i << "," << j << "]"
     //          << " fReferenceData " << fReferenceData.at(i).at(j)
     //          << std::endl;
     //    }
     //   }
     //  fReferenceChannels.resize(fCurrentBankIndex+1);
     //    fReferenceData.resize(fCurrentBankIndex+1);
 
     printf("\n------------- R1 LoadChannelMap End%s\n\n", mapfile.Data());
 
   }
 
   mapstr.Close(); // Close the file (ifstream)
 
   ReportConfiguration(local_debug);
   return 0;
 };

Here is the call graph for this function:

Int_t QwSciFiDetector::LoadInputParameters ( TString mapfile )

virtual

Mandatory parameter file definition.

Implements VQwSubsystem.

Definition at line 230 of file QwSciFiDetector.cc.

 {
   return 0;
 };

void QwSciFiDetector::PrintConfigurationBuffer	(	UInt_t *	buffer,
		UInt_t	num_words
	)

Definition at line 465 of file QwSciFiDetector.cc.

 {
   UInt_t ipt = 0;
   UInt_t j   = 0;
   UInt_t k   = 0;
   
   for ( j=0; j<(num_words/5); j++ ) 
     {
       printf ( "buffer[%5d] = 0x:", ipt );
       for ( k=j; k<j+5; k++ ) 
   {
     printf ( "%12x", buffer[ipt++] );
   }
       printf ( "\n" );
     }
   
   if ( ipt<num_words ) {
     printf ( "buffer[%5d] = 0x:", ipt );
     for ( k=ipt; k<num_words; k++ ) 
       {
   printf ( "%12x", buffer[ipt++] );
       }
     printf ( "\n" );
   }
   printf ( "\n" );
   
   return;
 }

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

virtual

Implements VQwSubsystem.

Definition at line 498 of file QwSciFiDetector.cc.

References QwF1TDContainer::AddQwF1TDC(), QwF1TDContainer::DoneF1TDCsConfiguration(), fF1TDContainer, VQwSubsystemTracking::fF1TDCResolutionNS, GetModuleIndex(), VQwSubsystem::GetSubbankIndex(), VQwSubsystem::GetSubsystemName(), kMaxNumberOfSlotsPerROC, QwF1TDContainer::Print(), QwF1TDC::SetF1BankIndex(), QwF1TDC::SetF1SystemName(), QwF1TDC::SetF1TDCBuffer(), and QwF1TDC::SetF1TDCIndex().

 {
 
   TString subsystem_name = "";
 
   Int_t bank_index     = 0;
   Int_t tdc_index      = 0;
   UInt_t slot_id       = 0;
   UInt_t vme_slot_num  = 0;
 
   Bool_t local_debug   = false;
   QwF1TDC *local_f1tdc = NULL;
    
   bank_index = GetSubbankIndex(roc_id, bank_id);
   
   if(bank_index >=0) {
     if(local_debug) {
       std::cout << "fF1TDContainer " << fF1TDContainer
     <<" local_f1tdc    " << local_f1tdc << "\n";
     }
     subsystem_name = this->GetSubsystemName();
     fF1TDContainer -> SetSystemName(subsystem_name);
 
     if(local_debug) {
       std::cout << "-----------------------------------------------------" << std::endl;
   
       std::cout << "\nQwSciFiDetector : " 
     << subsystem_name
     << ", "
     << "ProcessConfigurationBuffer"
     << std::endl;
       std::cout << "ROC " 
     << std::setw(2) << roc_id
     << " Bank [index,id]["
     <<  bank_index
     << ","
     << bank_id
     << "]"
     << std::endl;
     }
     for ( slot_id=0; slot_id<kMaxNumberOfSlotsPerROC; slot_id++ ) 
       { 
   // slot id starts from 2, because 0 is one offset (1) difference between QwAnalyzer and VME definition, 
   // and 1 and 2 are used for CPU and TI. Tuesday, August 31 10:57:07 EDT 2010, jhlee
   
   tdc_index    = GetModuleIndex(bank_index, slot_id);
   vme_slot_num = slot_id;
   if(local_debug) {
     std::cout << "    "
         << "Slot [id, VME num] [" 
         << std::setw(2) << slot_id
         << ","
         << std::setw(2) << vme_slot_num
         << "]";
     std::cout << "    ";
   }
   local_f1tdc = NULL;
 
   if(slot_id > 2) { // save time
   
     if (tdc_index not_eq -1) {
 
       if(local_f1tdc) delete local_f1tdc; local_f1tdc = NULL;
 
       local_f1tdc = new QwF1TDC(roc_id, vme_slot_num);
 
       local_f1tdc->SetF1BankIndex(bank_index);
       local_f1tdc->SetF1TDCIndex(tdc_index);
       local_f1tdc->SetF1TDCBuffer(buffer, num_words);
       local_f1tdc->SetF1SystemName(subsystem_name);
 
       fF1TDContainer->AddQwF1TDC(local_f1tdc);
       if(local_debug) { 
         std::cout << "F1TDC index " 
       << std::setw(2) 
       << tdc_index
       << std::setw(16) 
       << " local_f1tdc " 
       << *local_f1tdc
       << " at " 
       << local_f1tdc;
       }
     }
     else {
       if(local_debug){
         std::cout << "Unused in "  
       << std::setw(4) 
       << subsystem_name 
       << std::setw(16) 
       << " local_f1tdc  at " 
       << local_f1tdc;
       }
     }
     
   }
   else { // slot_id == only 0, 1, & 2
     if(local_debug) {
       if      (slot_id == 0) std::cout << "         ";
       else if (slot_id == 1) std::cout << "MVME CPU ";
       else                   std::cout << "Trigger Interface"; // slot_id == 2;
     }
   }
   if(local_debug) {
     std::cout << std::endl;
   }
       }
   
     fF1TDCResolutionNS = fF1TDContainer->DoneF1TDCsConfiguration();
 
     if(local_debug) fF1TDContainer->Print();
     std::cout << "-----------------------------------------------------" << std::endl;
 
     
     // kMaxNumberOfChannelsPerTDC = fF1TDContainer->GetF1TDCChannelNumber();
 
     return 0;
   }
   else {
    if(local_debug) {
      printf("\nQwSciFiDetector::ProcessConfigurationBuffer: \n");
      printf("Bank index is not defined with ROC id %d and Bank id %d. Thus, this is not F1TDC data. Skip this bank\n\n", roc_id, bank_id);
     }
  
     return -1;
   }
 }

Here is the call graph for this function:

Int_t QwSciFiDetector::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 643 of file QwSciFiDetector.cc.

References QwF1TDContainer::CheckDataIntegrity(), MQwF1TDC::DecodeTDCWord(), Qw::e, fDetectorIDs, fF1TDCDecoder, fF1TDContainer, FillRawTDCWord(), fScalerBankIndex, fSCAs, fSCAs_offset, GetModuleIndex(), VQwSubsystem::GetSubbankIndex(), VQwSubsystem::GetSubsystemName(), MQwF1TDC::GetTDCChannelNumber(), MQwF1TDC::GetTDCData(), MQwF1TDC::GetTDCSlotNumber(), IsSlotRegistered(), MQwF1TDC::IsValidDataword(), and VQwSubsystem::SetDataLoaded().

 {
 
   Int_t  bank_index             = 0;
   Int_t  vme_module_slot_number = 0;
   Int_t  vme_module_chan_number = 0;
   UInt_t vme_module_data        = 0;
   Int_t  module_index           = 0;
 
   Bool_t data_integrity_flag    = false;
   Bool_t temp_print_flag        = false;
 
   bank_index = GetSubbankIndex(roc_id, bank_id);
   
   if (bank_index>=0 && num_words>0) {
     //  We want to process this ROC.  Begin looping through the data.
     SetDataLoaded(kTRUE);
 
     if ( temp_print_flag ) {
       std::cout << "\nQwSciFiDetector::ProcessEvBuffer:  "
       << "Begin processing ROC" 
       << std::setw(2)
       << roc_id 
       << " bank id " 
       << bank_id 
       << " Subbbank Index "
       << bank_index
       << " Region "
       << GetSubsystemName()
       << std::endl;
     }
     
     if( bank_index != fScalerBankIndex) {
     //
     // CheckDataIntegrity() do "counter" whatever errors in each F1TDC 
     // and check whether data is OK or not.
 
       data_integrity_flag = fF1TDContainer->CheckDataIntegrity(roc_id, buffer, num_words);
     // if it is false (TFO, EMM, and SYN), the whole buffer is excluded for
     // the further process, because of multiblock data transfer.
 
       if (data_integrity_flag) {
   
   //  printf("\n");
   for (UInt_t i=0; i<num_words ; i++) {
     
     //  Decode this word as a F1TDC word.
     fF1TDCDecoder.DecodeTDCWord(buffer[i], roc_id); 
     
     vme_module_slot_number = fF1TDCDecoder.GetTDCSlotNumber();
     vme_module_chan_number = fF1TDCDecoder.GetTDCChannelNumber();
     module_index           = GetModuleIndex(bank_index, vme_module_slot_number);
     
     if ( vme_module_slot_number == 31) {
       //  This is a custom word which is not defined in
       //  the F1TDC, so we can use it as a marker for
       //  other data; it may be useful for something.
     }
     
     // Each subsystem has its own interesting slot(s), thus
     // here, if this slot isn't in its slot(s) (subsystem map file)
     // we skip this buffer to do the further process
     
     if (not IsSlotRegistered(bank_index, vme_module_slot_number) ) continue;
     
     if(temp_print_flag) std::cout << fF1TDCDecoder << std::endl;
     
     if ( fF1TDCDecoder.IsValidDataword() ) {//;;
       // if decoded F1TDC data has a valid slot, resolution locked, data word, no overflow (0xFFFF), and no fake data
       // try get time and fill it into raw QwHit.. FillRawTDCWord function is in each subsystem.
       
       try {
         vme_module_data = fF1TDCDecoder.GetTDCData();
         FillRawTDCWord(bank_index, vme_module_slot_number, vme_module_chan_number, vme_module_data);
 
       }
       catch (std::exception& e) {
         std::cerr << "Standard exception from QwSciFiDetector::FillRawTDCWord: "
       << e.what() << std::endl;
         std::cerr << "   Parameters:  index==" <<bank_index
       << "; Slot    Number=="   <<vme_module_slot_number
       << "; Channel Number=="<<vme_module_chan_number
       << "; Data=="         <<vme_module_data
       << std::endl;
         
         std::cerr << " Module Index=="<<module_index
       << "; fDetectorIDs.at(module_index).size()=="
       << fDetectorIDs.at(module_index).size()
       << "; fDetectorIDs.at(module_index).at(chan).fPlane=="
       << fDetectorIDs.at(module_index).at(vme_module_chan_number).fPlane
       << "; fDetectorIDs.at(module_index).at(chan).fElement=="
       << fDetectorIDs.at(module_index).at(vme_module_chan_number).fElement
       << std::endl;
       }
     }//;;
   } // for (UInt_t i=0; i<num_words ; i++) {
       }
  
     }
     else {
  
       // printf("Scaler %d\n", (Int_t)fSCAs.size());
       // Check if scaler buffer contains more than one event, and if so, reject it to decode.
       // 
       if (buffer[0]/32!=1) {
   if(temp_print_flag) printf("This buffer contains more than one event, thus, we reject them\n"); // Why? jhlee
   return 0;
       }
 
       UInt_t words_read = 0;
       for (std::size_t k=0; k<fSCAs.size(); k++) {
   words_read += fSCAs.at(k).ProcessEvBuffer(&(buffer[fSCAs_offset.at(k)]), num_words-fSCAs_offset.at(k));
   if(temp_print_flag) printf("k %d words_read %d\n", (Int_t)k, words_read);
       }
       
     }
   }
   
   // // fF1TDContainer-> PrintErrorSummary();
   return 0;
 }

Here is the call graph for this function:

void QwSciFiDetector::ProcessEvent ( )

virtual

Implements VQwSubsystem.

Definition at line 255 of file QwSciFiDetector.cc.

References fSCAs, VQwSubsystem::HasDataLoaded(), SubtractReferenceTimes(), and UpdateHits().

 {
   if (not HasDataLoaded()) return;
 
   // F1TDCs
   SubtractReferenceTimes();
   UpdateHits();
 
 
   // SIS3801 scaler
   std::size_t i = 0;
   for (i=0; i<fSCAs.size(); i++)
     {
       fSCAs.at(i).ProcessEvent();
     }
   return;
 };

Here is the call graph for this function:

Int_t QwSciFiDetector::RegisterROCNumber ( const UInt_t roc_id )

protected

Definition at line 902 of file QwSciFiDetector.cc.

References fModuleIndex, kMaxNumberOfSlotsPerROC, and VQwSubsystem::RegisterROCNumber().

Referenced by LoadChannelMap().

 {
   Int_t status = 0;
   status = VQwSubsystemTracking::RegisterROCNumber(roc_id, 0);
   std::vector<Int_t> tmpvec(kMaxNumberOfSlotsPerROC,-1);
   fModuleIndex.push_back(tmpvec);
   //std::cout<<"Registering ROC "<<roc_id<<std::endl;
 
   return status;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwSciFiDetector::RegisterSlotNumber ( const UInt_t slot_id )

protected

Definition at line 951 of file QwSciFiDetector.cc.

References fCurrentBankIndex, fCurrentModuleIndex, fCurrentSlot, fDetectorIDs, fModuleIndex, kMaxNumberOfChannelsPerF1TDC, kMaxNumberOfSlotsPerROC, and kNumberOfVMEModules.

Referenced by LoadChannelMap().

 {
   if (slot_id<kMaxNumberOfSlotsPerROC) {
     
     if (fCurrentBankIndex <= (Int_t) fModuleIndex.size()) {
  
       fDetectorIDs.resize(kNumberOfVMEModules+1);
       fDetectorIDs.at(kNumberOfVMEModules).resize(kMaxNumberOfChannelsPerF1TDC);
 
       // reassign kNumberOfVMEModules after resize it
 
       kNumberOfVMEModules = (Int_t) fDetectorIDs.size();
       fCurrentModuleIndex = kNumberOfVMEModules-1;
       fCurrentSlot        = slot_id;
       fModuleIndex.at(fCurrentBankIndex).at(fCurrentSlot) = fCurrentModuleIndex;
     }
   } 
   else {
     std::cout << "QwSciFiDetector::RegisterSlotNumber:  Slot number "
        << slot_id << " is larger than the number of slots per ROC, "
        << kMaxNumberOfSlotsPerROC << std::endl;
   }
   return fCurrentModuleIndex;
     
 }

Here is the caller graph for this function:

Int_t QwSciFiDetector::RegisterSubbank ( const UInt_t bank_id )

protected

Definition at line 915 of file QwSciFiDetector.cc.

References fCurrentBankIndex, VQwSubsystem::fCurrentROC_ID, fModuleIndex, fReferenceChannels, fReferenceData, VQwSubsystem::GetSubbankIndex(), kMaxNumberOfSlotsPerROC, and VQwSubsystem::RegisterSubbank().

Referenced by LoadChannelMap().

 {
 
   Int_t status         = VQwSubsystem::RegisterSubbank(bank_id);
   Int_t current_roc_id = VQwSubsystem::fCurrentROC_ID;
 
   std::size_t reference_size = fReferenceChannels.size();
  
   // from here, just expand some vectors in QwSciFiDetector...
   //
   fCurrentBankIndex = GetSubbankIndex(current_roc_id, bank_id);
   //subbank id is directly related to the ROC
 
   if ( (Int_t) reference_size <= fCurrentBankIndex) {
     fReferenceChannels.resize(fCurrentBankIndex+1);
     fReferenceData.resize(fCurrentBankIndex+1);
   }
 
   // Why do we need to expand this vector size
   // even if most of them are empty? --- jhlee
 
   std::vector<Int_t> tmpvec(kMaxNumberOfSlotsPerROC,-1);
 
   fModuleIndex.push_back(tmpvec);
 
   std::cout<< "QwSciFiDetector::RegisterSubbank()" 
      <<" ROC "            << current_roc_id
      <<" Subbank "        << bank_id
      <<" with BankIndex " << fCurrentBankIndex
      << std::endl;
 
   return status;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSciFiDetector::ReportConfiguration ( Bool_t verbose )

Definition at line 362 of file QwSciFiDetector.cc.

References VQwSubsystem::fBank_IDs, fDetectorIDs, VQwSubsystem::fROC_IDs, GetModuleIndex(), VQwSubsystem::GetSubbankIndex(), VQwSubsystem::GetSubsystemName(), and kMaxNumberOfSlotsPerROC.

Referenced by LoadChannelMap().

 {
 
   // Bank index 1 : ROC4  F1TDC
   // Bank index 3 : ROC9  F1TDCs
   // Bank index 5 : ROC11 SIS3801 
 
 
   if(verbose) {
     std::size_t i       = 0;
     std::size_t j       = 0;
 
     Int_t roc_num       = 0;
     Int_t bank_flag     = 0;
     Int_t bank_index    = 0;
     Int_t module_index  = 0;
 
     UInt_t slot_id      = 0;
     UInt_t vme_slot_num = 0;
     
     std::cout << "QwSciFiDetector Region : " 
         << this->GetSubsystemName()
         << "::ReportConfiguration fDetectorIDs.size() " 
         << fDetectorIDs.size() << std::endl;
 
 
     for ( i=0; i<fROC_IDs.size(); i++ ) 
       {
 
   roc_num = fROC_IDs.at(i);
 
   for ( j=0; j<fBank_IDs.at(i).size(); j++ ) 
     {
       bank_flag = fBank_IDs.at(i).at(j);
       if(bank_flag == 0) continue; 
       // must be uncommented if one doesn't define "bank_flag" in a CRL file
       // but, now we use "bank_flag" in our crl files, thus skip to print
       // unnecessary things on a screen
       // Monday, August 30 14:45:34 EDT 2010, jhlee
 
       bank_index = GetSubbankIndex(roc_num, bank_flag);
   
       std::cout << "ROC [index, Num][" 
           << i
           << ","
           << std::setw(2) << roc_num
           << "]"
           << " Bank [index,id]["
           <<  bank_index
           << ","
           << bank_flag
           << "]"
           << std::endl;
   
       for ( slot_id=2; slot_id<kMaxNumberOfSlotsPerROC; slot_id++ ) 
         { 
     // slot id starts from 2, because 0 and 1 are used for CPU and TI.
     // Tuesday, August 31 10:57:07 EDT 2010, jhlee
 
     module_index = GetModuleIndex(bank_index, slot_id);
     
     vme_slot_num = slot_id;
     
     std::cout << "    "
         << "Slot [id, VME num] [" 
         << std::setw(2) << slot_id
         << ","
         << std::setw(2) << vme_slot_num
         << "]";
     if ( module_index == -1 ) {
       std::cout << "    "
           << "Unused in R1 SciFiDetector" 
           << std::endl;
     }
     else {
       std::cout << "    "
           << "Module index " 
           << module_index << std::endl;
     }
         }
     }
       }
     
     for( size_t midx = 0; midx < fDetectorIDs.size(); midx++ )
       {
   for (size_t chan = 0 ; chan< fDetectorIDs.at(midx).size(); chan++)
     {
       std::cout << "[" << midx <<","<< chan << "] "
           << " detectorID " << fDetectorIDs.at(midx).at(chan)
           << std::endl;
     }
       }
   }
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSciFiDetector::SubtractReferenceTimes ( )

protected

Definition at line 993 of file QwSciFiDetector.cc.

References fF1RefContainer, fF1TDContainer, fTDCHits, F1TDCReferenceContainer::GetReferenceTimeAU(), VQwSubsystem::GetSubsystemName(), QwMessage, and QwF1TDContainer::ReferenceSignalCorrection().

Referenced by ProcessEvent().

 {
 
   UInt_t   bank_index        = 0;
   Double_t raw_time_arb_unit = 0.0;
   Double_t ref_time_arb_unit = 0.0;
   Double_t time_arb_unit     = 0.0;
  // EQwDetectorPackage package = kPackageNull;  // Region 1 has only Package 2, so kPackageDown
 
   Bool_t local_debug = false;
   Int_t slot_num = 0;
 
   TString reference_name1 = "MasterTrigger";
   TString reference_name2 = "CopyMasterTrigger";
 
 
   for ( std::vector<QwHit>::iterator hit=fTDCHits.begin(); hit!=fTDCHits.end(); hit++ ) 
     {
 
       bank_index        = hit -> GetSubbankID();
       slot_num          = hit -> GetModule();
       raw_time_arb_unit = (Double_t) hit -> GetRawTime();
       ref_time_arb_unit = fF1RefContainer -> GetReferenceTimeAU(bank_index, reference_name1);
       //
       // if there is no reference time due to a channel error, try to use a copy of mater trigger
       // 
       if(ref_time_arb_unit==0.0) {
   ref_time_arb_unit =  fF1RefContainer->GetReferenceTimeAU(bank_index, reference_name2);
       }
       // second time, it returns 0.0, we simply ignore this event .... 
       // set time zero. ReferenceSignalCorrection() will return zero, and increase RFM counter...
       //
       time_arb_unit = fF1TDContainer->ReferenceSignalCorrection(raw_time_arb_unit, ref_time_arb_unit, bank_index, slot_num);
 
       hit -> SetTime(time_arb_unit); 
       hit -> SetRawRefTime((UInt_t) ref_time_arb_unit);
 
       if(local_debug) {
   QwMessage << this->GetSubsystemName()
       << " BankIndex " << std::setw(2) << bank_index
       << " Slot "      << std::setw(2) << slot_num
       << " RawTime : " << std::setw(6) << raw_time_arb_unit
       << " RefTime : " << std::setw(6) << ref_time_arb_unit
       << " time : "    << std::setw(6) << time_arb_unit
       << std::endl;
   
       }
     }
 
   return;
   
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwSciFiDetector::UpdateHits ( )

protected

Definition at line 1167 of file QwSciFiDetector.cc.

References VQwSubsystemTracking::fF1TDCResolutionNS, and fTDCHits.

Referenced by ProcessEvent().

 {
  
   for(std::vector<QwHit>::iterator iter=fTDCHits.begin(); iter!=fTDCHits.end(); ++iter)
     {
   
       iter->ApplyTimeCalibration(fF1TDCResolutionNS); // Fill fTimeRes and fTimeNs in QwHit
 
     }
 
   return;
 }