#include <QwTriggerScintillator.h>

Inheritance diagram for QwTriggerScintillator:

Collaboration diagram for QwTriggerScintillator:

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

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

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

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

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

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

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

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

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

void	ReportConfiguration (Bool_t verbose)

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	FillRawTDCWord (Int_t bank_index, Int_t slot_num, Int_t chan, UInt_t data)

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

void	ClearEventData ()

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

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

QwTriggerScintillator &	operator= (const QwTriggerScintillator &value)

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

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

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

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

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

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

Static Public Member Functions
static void	DefineOptions (QwOptions &options)

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

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

static QwTriggerScintillator *	Cast (QwTriggerScintillator *type)
	Object dynamic cast. 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)

Int_t	GetModuleIndex (size_t bank_index, size_t slot_num) const

Bool_t	IsSlotRegistered (Int_t bank_index, Int_t slot_num) const

EQwModuleType	RegisterModuleType (TString moduletype)

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

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

void	SubtractReferenceTimes ()

void	UpdateHits ()

void	AddSoftwareMeantimeToHits (Bool_t option)

Bool_t	IsF1ReferenceChannel (Int_t slot, Int_t chan)

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

Bool_t	fSoftwareMeantimeOption

Double_t	fSoftwareMeantimeTimeWindowNs

TString	fRegion

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

Int_t	fCurrentSlot

Int_t	fCurrentModuleIndex

UInt_t	fBankID [3]

UInt_t	kMaxNumberOfChannelsPerF1TDC

Int_t	fNumberOfModules

MQwV775TDC	fQDCTDC

MQwF1TDC	fF1TDCDecoder

QwF1TDContainer *	fF1TDContainer

F1TDCReferenceContainer *	fF1RefContainer

MeanTimeContainer *	fSoftwareMeantimeContainer [2]

std::vector< QwHit >	fTDCHits

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

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

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

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

std::vector< EQwModuleType >	fModuleTypes

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

Int_t	fRefTime_SlotNum

Int_t	fRefTime_ChanNum

Double_t	fRefTime

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 Int_t	kF1ReferenceChannelNumber = 99

static const Int_t	kMaxNumberOfQwHitPlane = 2

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

Detailed Description

Definition at line 33 of file QwTriggerScintillator.h.

Constructor & Destructor Documentation

QwTriggerScintillator::QwTriggerScintillator ( )

private

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

QwTriggerScintillator::QwTriggerScintillator ( const TString & name )

Constructor with name.

Definition at line 22 of file QwTriggerScintillator.cc.

References ClearAllBankRegistrations(), fBankID, fCurrentBankIndex, fCurrentModuleIndex, fCurrentSlot, fF1RefContainer, fF1TDCDecoder, fF1TDContainer, fNumberOfModules, fRegion, fSoftwareMeantimeContainer, fSoftwareMeantimeOption, fSoftwareMeantimeTimeWindowNs, QwF1TDContainer::GetF1TDCDecoder(), MQwF1TDC::GetTDCMaxChannels(), kMaxNumberOfChannelsPerF1TDC, kMaxNumberOfQwHitPlane, and VQwSubsystem::SetEventTypeMask().

 : VQwSubsystem(name),VQwSubsystemTracking(name)
 {
   SetEventTypeMask(0xffff); 
 
   ClearAllBankRegistrations();
 
   fRegion             = "";
   fCurrentBankIndex   = 0;
   fCurrentSlot        = 0;
   fCurrentModuleIndex = 0;
 
   fBankID[0]          = 0;
   fBankID[1]          = 0;
   fBankID[2]          = 0;
 
   fNumberOfModules    = 0;
 
   fF1TDContainer = new QwF1TDContainer();
   fF1TDCDecoder  = fF1TDContainer->GetF1TDCDecoder();
   kMaxNumberOfChannelsPerF1TDC = fF1TDCDecoder.GetTDCMaxChannels();
   fF1RefContainer = new F1TDCReferenceContainer();
 
 
   for (Int_t plane_idx=0; plane_idx<kMaxNumberOfQwHitPlane; plane_idx++)
     {
       fSoftwareMeantimeContainer[plane_idx] = new MeanTimeContainer();  
     }
 
   fSoftwareMeantimeOption = false;
   fSoftwareMeantimeTimeWindowNs = 0.0;
 
 }

Here is the call graph for this function:

QwTriggerScintillator::~QwTriggerScintillator ( )

virtual

Virtual destructor.

Definition at line 56 of file QwTriggerScintillator.cc.

References fF1RefContainer, fF1TDContainer, fPMTs, fSCAs, fSoftwareMeantimeContainer, and kMaxNumberOfQwHitPlane.

 {
   fPMTs.clear();
   fSCAs.clear();
   delete fF1TDContainer;
   delete fF1RefContainer;
 
   for (Int_t plane_idx=0; plane_idx < kMaxNumberOfQwHitPlane; plane_idx++)
     {
       delete fSoftwareMeantimeContainer[plane_idx];
     }
 
 };

Member Function Documentation

void QwTriggerScintillator::AddSoftwareMeantimeToHits ( Bool_t option )

protected

Definition at line 1355 of file QwTriggerScintillator.cc.

References fSoftwareMeantimeContainer, fSoftwareMeantimeTimeWindowNs, fTDCHits, QwSubsystemArray::GetCodaEventNumber(), MeanTimeContainer::GetMeanTimeObject(), VQwSubsystem::GetParent(), MeanTime::GetPlane(), MeanTimeContainer::GetPlane(), VQwSubsystem::GetSubsystemName(), kDirectionNull, kMaxNumberOfQwHitPlane, kRegionIDTrig, QwHit::SetTimens(), and MeanTimeContainer::SoftwareMTSize().

Referenced by ProcessEvent().

 {
   // reduce if {}, 
 
   // if option is false, end this function
   if( !option )  return;
   // // if fTDCHits size is zero, end this function
   if (fTDCHits.size() ==0)   return;
   
 
   Bool_t           local_debug = false;
 
   Long64_t           ev_num    = 0;
   Int_t              plane     = 0;
   Int_t              octant    = 0;
   Int_t              element   = 0;
   Int_t              hitnumber = 0;
   Double_t           timens    = 0.0;
  
 
   Int_t              bank_index = 0;
   Int_t              slot_num   = 0;
   Int_t              chan_num   = 0;
   EQwDetectorPackage package    = kPackageNull;
   EQwDirectionID     direction  = kDirectionNull;
   EQwRegionID        region     = kRegionIDTrig;
 
   MeanTime*          ts_mt_time = NULL;
 
   Int_t              v_plane_idx = 0;
 
   for(std::vector<QwHit>::iterator iter=fTDCHits.begin(); iter!=fTDCHits.end(); ++iter)
     {
       plane = iter->GetPlane();
       if(plane == 1 || plane == 2) {
   v_plane_idx = plane -1;
   element     = iter->GetElement();
   hitnumber   = iter->GetHitNumber();
   timens      = iter->GetTimeNs();
   fSoftwareMeantimeContainer[v_plane_idx] -> Add(element, hitnumber, timens);
       }
     }
     
   // Reset bank_index, slot_num, and chan_num to -9,-99,-999
   // because we don't have these numbers for software meantime
   // and we don't need these numbers in order to access software meantime from QwHit
   
   ev_num     = VQwSubsystem::GetParent()->GetCodaEventNumber();
   bank_index = -9;
   slot_num   = -99;
   chan_num   = -999;
 
   
   for (v_plane_idx=0; v_plane_idx<kMaxNumberOfQwHitPlane; v_plane_idx++)
     {
    
       if(v_plane_idx == 0) {
   package = kPackage1;
       }
       else if (v_plane_idx ==1){
   package = kPackage2;
       }
       else {
   break;
       }
       
       fSoftwareMeantimeContainer[v_plane_idx] -> SetEventId(ev_num);
       fSoftwareMeantimeContainer[v_plane_idx] -> SetTimeWindow(fSoftwareMeantimeTimeWindowNs);
       fSoftwareMeantimeContainer[v_plane_idx] -> ProcessMeanTime();
       plane = fSoftwareMeantimeContainer[v_plane_idx]->GetPlane();
       
       ts_mt_time = NULL;
       for (Int_t v_smt_idx=0; v_smt_idx < fSoftwareMeantimeContainer[v_plane_idx]->SoftwareMTSize(); v_smt_idx++ )
   {
     ts_mt_time = fSoftwareMeantimeContainer[v_plane_idx]->GetMeanTimeObject(v_smt_idx);
           ts_mt_time -> Print(local_debug);
 
           QwHit software_meantime_hit(bank_index, slot_num, chan_num, v_smt_idx, 
                                       region, package, octant, plane, direction,
                                       ts_mt_time->GetSoftwareMeantimeHitElement()
                                       );
           software_meantime_hit.SetTimens(ts_mt_time->GetMeanTime());
           fTDCHits.push_back(software_meantime_hit);
           
           QwHit software_positive_hit(bank_index, slot_num, chan_num, v_smt_idx, 
                                       region, package, octant, plane, direction,
                                       ts_mt_time->GetSoftwarePositiveHitElement()
                                       );
           software_positive_hit.SetTimens(ts_mt_time->GetPositiveValue());
     software_positive_hit.SetHitNumberR(ts_mt_time->GetPositiveHitId());
           fTDCHits.push_back(software_positive_hit);
           
           QwHit software_negative_hit(bank_index, slot_num, chan_num, v_smt_idx, 
                                       region, package, octant, plane, direction,
                                       ts_mt_time->GetSoftwareNegativeHitElement()
                                       );
           software_negative_hit.SetTimens(ts_mt_time->GetNegativeValue());
     software_negative_hit.SetHitNumberR(ts_mt_time->GetNegativeHitId());
           fTDCHits.push_back(software_negative_hit);
           
           QwHit software_subtract_hit(bank_index, slot_num, chan_num, v_smt_idx, 
                                       region, package, octant, plane, direction,
                                       ts_mt_time->GetSoftwareSubtractHitElement()
                                       );
           software_subtract_hit.SetTimens(ts_mt_time->GetSubtractTime());
     fTDCHits.push_back(software_subtract_hit);
   }
     }
 
 
   if(local_debug) {
 
     TString output    = "";
     
     for(std::vector<QwHit>::iterator iter=fTDCHits.begin(); iter!=fTDCHits.end(); ++iter)
       {
   
   bank_index = iter->GetSubbankID();
   slot_num   = iter->GetModule();
   chan_num   = iter->GetChannel();
   package    = iter->GetPackage();
   //  direction  = iter->GetDirection();
   plane      = iter->GetPlane();
   element    = iter->GetElement();
   hitnumber  = iter->GetHitNumber();
   timens     = iter->GetTimeNs();
   
   
   output += GetSubsystemName();
   output += Form(" Bank ID %+2d", bank_index);
   output += Form(" Slot %+3d", slot_num);
   output += Form(" Chan %+4d", chan_num);
   output += Form(" Package %2d", (Int_t) package);
   output += Form(" Plane %2d", plane);
   output += " Element ";
   output += element;
   output += " Hit ";
   output += hitnumber;
   output += Form(" TimeNs %+10.2f\n", timens);
   
   
       }
     
     std::cout << "Size fTDCHits " << fTDCHits.size()
         << " Event Number " << ev_num << "\n" << output << std::endl;
   }
   
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwTriggerScintillator::ClearAllBankRegistrations ( )

protected

Definition at line 889 of file QwTriggerScintillator.cc.

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

Referenced by QwTriggerScintillator().

 {
   VQwSubsystemTracking::ClearAllBankRegistrations();
   fModuleIndex.clear();
   fModulePtrs.clear();
   fModuleTypes.clear();
 
   fDetectorIDs.clear();
   fTDCHits.clear();
   fNumberOfModules = 0;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwTriggerScintillator::ClearEventData ( )

virtual

Implements VQwSubsystem.

Definition at line 307 of file QwTriggerScintillator.cc.

References fF1RefContainer, fPMTs, fSCAs, fSoftwareMeantimeContainer, fTDCHits, and VQwSubsystem::SetDataLoaded().

 {
   SetDataLoaded(kFALSE);
 
 
   fTDCHits.clear();
   std::size_t i = 0;
   
   fF1RefContainer               -> ClearEventData();
   fSoftwareMeantimeContainer[0] -> ClearEventData();
   fSoftwareMeantimeContainer[1] -> ClearEventData();
 
   // for (i=0; i<fReferenceData.size(); i++) 
   //   {
   //     fReferenceData.at(i).clear();
   //   }
 
   for (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 (i=0; i<fSCAs.size(); i++) 
     {
       fSCAs.at(i).ClearEventData();
     }
 
   return;
 }

Here is the call graph for this function:

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

virtual

Construct the branch and tree vector.

Reimplemented from VQwSubsystemTracking.

Definition at line 804 of file QwTriggerScintillator.cc.

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

 {
   fTreeArrayIndex = values.size();
 
   TString basename;
   if (prefix=="") basename = "trigscint";
   else            basename = prefix;
 
   TString element_name = "";
 
   TString list = "";
   for (size_t i=0; i<fPMTs.size(); i++){
     for (size_t j=0; j<fPMTs.at(i).size(); j++){
       element_name = fPMTs.at(i).at(j).GetElementName();
       if (element_name!="") {
   values.push_back(0.0);
   list += ":" + element_name + "/D";
       }
     }
   }
   
   element_name = "";
 
   for (size_t i=0; i<fSCAs.size(); i++){
     element_name = fSCAs.at(i).GetElementName();
     if (element_name != "") {
       values.push_back(0.0);
       list += ":" + element_name + "/D";
     }
   }
 
   if (list[0]==':') {
     list = list(1,list.Length()-1);
   }
 
   fTreeArrayNumEntries = values.size() - fTreeArrayIndex;
   tree->Branch(basename, &values[fTreeArrayIndex], list);
   return;
 }

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

virtual

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

Implements VQwSubsystem.

Definition at line 776 of file QwTriggerScintillator.cc.

References fPMTs, and fSCAs.

                                                                                    {
 
   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, prefix);
     }
   }
 
   for (size_t i=0; i<fSCAs.size(); i++) {
     fSCAs.at(i).ConstructHistograms(folder, prefix);
   }
 
 }

void QwTriggerScintillator::DefineOptions ( QwOptions & options )

static

Definition at line 754 of file QwTriggerScintillator.cc.

References QwOptions::AddOptions(), and default_bool_value.

 {
   options.AddOptions()("enable-ts-software-meantime",
            po::value<Bool_t>()->default_bool_value(true),
            "Create Software meantime for TS in QwHits" 
            );
   options.AddOptions()("set-ts-software-meantime-timewindow",
            po::value<Double_t>()->default_value(20.0),
            "TimeWindow (ns) for TS Software meantime"
            );//max window is 2000 ns
   return;
 };

Here is the call graph for this function:

void QwTriggerScintillator::FillHardwareErrorSummary ( )

virtual

Hardware error summary.

Reimplemented from VQwSubsystemTracking.

Definition at line 1536 of file QwTriggerScintillator.cc.

References fF1RefContainer, and fF1TDContainer.

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

void QwTriggerScintillator::FillHistograms ( )

virtual

Fill the histograms for this subsystem.

Implements VQwSubsystem.

Definition at line 790 of file QwTriggerScintillator.cc.

References fPMTs, 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();
   }
 
 }

Here is the call graph for this function:

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

Definition at line 1037 of file QwTriggerScintillator.cc.

References fDetectorIDs, fF1RefContainer, fReferenceData, fTDCHits, GetModuleIndex(), kDirectionNull, kF1ReferenceChannelNumber, kRegionIDTrig, 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) {
     // if ( slot_num == 4 and chan == 4) {
     //   printf("FillRawTDCWord bank_index %2d slot_number %2d chan %2d, data %10d \n", 
     //       bank_index, slot_num, chan, data);
     // }
     Int_t hitcnt  = 0;
 
     EQwDirectionID direction   = kDirectionNull;
 
     TString name         = "";
  
     fF1RefContainer->SetReferenceSignal(bank_index, slot_num, chan, data, local_debug);
 
     // ts2 plane 2, ts1 is plane 1
     Int_t plane   = fDetectorIDs.at(tdcindex).at(chan).fPlane;
     // mt is element 0, p is element 1, and m is element 2
     Int_t element = fDetectorIDs.at(tdcindex).at(chan).fElement;
     EQwDetectorPackage package = fDetectorIDs.at(tdcindex).at(chan).fPackage;
     Int_t 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
       if(local_debug) {
   std::cout << "At QwTriggerScintillator::FillRawTDCWord " 
       << " bank index " << bank_index
       << " slot num " << slot_num
       << " chan num " << chan
       << " hitcnt " << 0
       << " plane  " << plane
       << " wire   " << element
       << " package " << package 
       << " diection " << direction
       << " data " << data
       << " fTDCHits.size() " <<  0
       << std::endl;
       }
     }
     else {
       direction = fDetectorIDs.at(tdcindex).at(chan).fDirection;
 
       hitcnt = std::count_if(fTDCHits.begin(), fTDCHits.end(), 
            boost::bind(
            &QwHit::WireMatches, _1, kRegionIDTrig, boost::ref(package), boost::ref(plane), boost::ref(element)
            ) 
            );
    
       fTDCHits.push_back(
            QwHit(
                  bank_index, 
                  slot_num, 
                  chan, 
                                hitcnt, 
                                kRegionIDTrig, 
                                package, 
                                octant,
                                plane,
                                direction, 
                                element, 
                  data
                  )
            );
 
       if(local_debug) {
   std::cout << "At QwTriggerScintillator::FillRawTDCWord " 
       << " bank index " << bank_index
       << " slot num " << slot_num
       << " chan num " << chan
       << " hitcnt " << hitcnt
       << " plane  " << plane
       << " wire   " << element
       << " package " << package 
       << " diection " << direction
       << " data " << data
       << " 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 QwTriggerScintillator::FillRawWord	(	Int_t	bank_index,
		Int_t	slot_num,
		Int_t	chan,
		UInt_t	data
	)

protected

Definition at line 1011 of file QwTriggerScintillator.cc.

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

Referenced by ProcessEvBuffer().

 {
   Int_t modindex = GetModuleIndex(bank_index,slot_num);
   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 {
       // if (bank_index == 2 and  slot_num == 4 and chan == 4) {
       //  printf("FillRawWord    bank_index %2d slot_number %2d chan %2d, data %10d \n", 
       //         bank_index, slot_num, chan, data);
       // }
       fPMTs.at(modtype).at(chanindex).SetValue(data);
       fPMTs.at(modtype).at(chanindex).SetSubbankID(bank_index);
       fPMTs.at(modtype).at(chanindex).SetModule(slot_num);
     }
   }
   
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Fill the tree vector.

Reimplemented from VQwSubsystemTracking.

Definition at line 844 of file QwTriggerScintillator.cc.

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

 {
   if (! HasDataLoaded()) return;
 
   Int_t index = fTreeArrayIndex;
   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()!="") {
         values[index] = fPMTs.at(i).at(j).GetValue();
         index++;
       }
     }
   }
   
   for (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:

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

protected

Definition at line 1521 of file QwTriggerScintillator.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:

void QwTriggerScintillator::GetHitList ( QwHitContainer & grandHitContainer )

inlinevirtual

Get the hit list.

Implements VQwSubsystemTracking.

Definition at line 78 of file QwTriggerScintillator.h.

References QwHitContainer::Append(), and fTDCHits.

                                                       {
     grandHitContainer.Append(fTDCHits);
   };

Here is the call graph for this function:

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

protected

Definition at line 1508 of file QwTriggerScintillator.cc.

References fModuleIndex.

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

 {
   Int_t modindex = -1;
 
   if ( bank_index < fModuleIndex.size()) {
     if ( 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:

Bool_t QwTriggerScintillator::IsF1ReferenceChannel	(	Int_t	slot,
		Int_t	chan
	)

inlineprotected

Definition at line 180 of file QwTriggerScintillator.h.

References fRefTime_ChanNum, and fRefTime_SlotNum.

Referenced by ProcessEvBuffer().

                                                        { 
     return ( slot == fRefTime_SlotNum &&  chan == fRefTime_ChanNum) ;
   };

Here is the caller graph for this function:

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

inlineprotected

Definition at line 131 of file QwTriggerScintillator.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 QwTriggerScintillator::LinkChannelToSignal	(	const UInt_t	chan,
		const TString &	name
	)

protected

Definition at line 1001 of file QwTriggerScintillator.cc.

References fCurrentModuleIndex, fCurrentType, fModulePtrs, and fPMTs.

Referenced by LoadChannelMap().

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

Here is the caller graph for this function:

Int_t QwTriggerScintillator::LoadChannelMap ( TString mapfile )

virtual

Mandatory map file definition.

Implements VQwSubsystem.

Definition at line 71 of file QwTriggerScintillator.cc.

References fBankID, fCurrentBankIndex, fCurrentModuleIndex, fCurrentType, fDetectorIDs, VQwSubsystem::fDetectorMaps, fF1RefContainer, FindSignalIndex(), fModuleIndex, fModulePtrs, fReferenceChannels, fRefTime_ChanNum, fRefTime_SlotNum, fSCAs, fSCAs_map, fSCAs_offset, QwParameterFile::GetUInt(), kDirectionNull, kF1ReferenceChannelNumber, kRegionIDTrig, kUnknownModuleType, LinkChannelToSignal(), RegisterModuleType(), RegisterROCNumber(), RegisterSlotNumber(), RegisterSubbank(), and ReportConfiguration().

 {
 
   TString varname  = "";
   TString varvalue = "";
   UInt_t  value     = 0;
 
   TString modtype  = "";
   TString dettype  = "";
   TString name     = "";
 
   Int_t   modnum   = 0;
   Int_t   channum  = 0;
   Int_t   slotnum  = 0;
 
   Int_t   ts_chan_num_to_plane    = 0;
   Int_t   ts_chan_type_to_element = 0;
 
   // 
   // u ->? ts2 --> plane 2
   // d ->? ts1 --> plane 1
 
   // mt -> element 0
   // p  -> element 1
   // m  -> element 2
 
 
   Bool_t local_debug = false;
   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 (varname=="roc") {
   RegisterROCNumber(value);
       } else if (varname=="qdc_bank") {
     RegisterSubbank(value);
           fBankID[0] = value;
       } else if (varname=="sca_bank") {
           RegisterSubbank(value);
           fBankID[1] = value;
       } else if (varname=="f1tdc_bank") {
     RegisterSubbank(value);
           fBankID[2] = value;
       } else if (varname=="slot") {
           RegisterSlotNumber(value);
           slotnum = 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>();
 
   if (local_debug) printf("%8s, %d, %d %s, %s\n", modtype.Data(), modnum, channum, dettype.Data(), name.Data());
         //  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;
     fSCAs_offset.push_back(QwSIS3801D24_Channel::GetBufferOffset(modnum,channum));
 
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fRegion    = kRegionIDTrig;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPackage   = package;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPlane     = -1;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fDirection = direction;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fElement   = -1;
     
         } 
   else if  (modtype=="V792") { // || modtype=="F1TDC") {
     RegisterModuleType(modtype);
     //  Check to see if we've encountered this channel or name yet
     if (fModulePtrs.at(fCurrentModuleIndex).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);
     }
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fRegion    = kRegionIDTrig;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPackage   = package;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPlane     = -1;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fDirection = direction;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fElement   = -1;
   }
   else if( modtype=="F1TDC") {
 
     RegisterModuleType(modtype);
     //  Check to see if we've encountered this channel or name yet
     if (fModulePtrs.at(fCurrentModuleIndex).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);
     }
 
 
     if (name=="ts_reftime_f1") {
       
       fF1RefContainer -> AddF1TDCReferenceSignal(new F1TDCReferenceSignal(fCurrentBankIndex, slotnum, channum, "MasterTrigger"));
       reference_counter++;
 
       fRefTime_SlotNum = slotnum;
       fRefTime_ChanNum = channum;
 
       // printf("bank index %d Chan %d reference_counter %d\n", fCurrentBankIndex, channum, reference_counter);
       fReferenceChannels.at ( fCurrentBankIndex ).first  = fCurrentModuleIndex;
       fReferenceChannels.at ( fCurrentBankIndex ).second = channum;
       
       package = kPackageNull;
 
       fDetectorIDs.at(fCurrentModuleIndex).at(channum).fElement = fCurrentBankIndex;
       fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPlane   = kF1ReferenceChannelNumber;
 
     }
     else {
       if (name=="ts2m_f1") {
         ts_chan_num_to_plane    = 2;
         ts_chan_type_to_element = 2;
         package = kPackage2;
       }
       else if (name=="ts2p_f1") {
         ts_chan_num_to_plane    = 2;
         ts_chan_type_to_element = 1;
         package = kPackage2;
       }
       else if (name=="ts2mt_f1") {
         ts_chan_num_to_plane    = 2;
         ts_chan_type_to_element = 0;
         package = kPackage2;
       }
       else if (name=="ts1m_f1") {
         ts_chan_num_to_plane    = 1;
         ts_chan_type_to_element = 2;
         package = kPackage1;
       }
       else if (name=="ts1p_f1") {
         ts_chan_num_to_plane    = 1;
         ts_chan_type_to_element = 1;
         package = kPackage1;
       }
       else if (name=="ts1mt_f1") {
         ts_chan_num_to_plane    = 1;
         ts_chan_type_to_element = 0;        
         package = kPackage1;
 
       }
       else {
         ts_chan_num_to_plane    = -1;
         ts_chan_type_to_element = -1;       
         package                 = kPackageNull;
       }
       fDetectorIDs.at(fCurrentModuleIndex).at(channum).fElement   = ts_chan_type_to_element;
       fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPlane     = ts_chan_num_to_plane;
   
     }
       
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fRegion    = kRegionIDTrig;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fPackage   = package;
     fDetectorIDs.at(fCurrentModuleIndex).at(channum).fDirection = direction;
   
   }
   else {
     std::cerr << "LoadChannelMap:  Unknown line: " << mapstr.GetLine().c_str() << std::endl;
         }
     }
   }
 
   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 " << m_idx
       << 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;
       }
     
     printf("\n------------- TS LoadChannelMap End%s\n\n", mapfile.Data());
 
   }
 
   ReportConfiguration(local_debug);
 
   mapstr.Close(); // Close the file (ifstream)
   return 0;
 }

Here is the call graph for this function:

Int_t QwTriggerScintillator::LoadInputParameters ( TString mapfile )

inlinevirtual

Mandatory parameter file definition.

Implements VQwSubsystem.

Definition at line 51 of file QwTriggerScintillator.h.

51 {return 0;};

QwTriggerScintillator & QwTriggerScintillator::operator= ( const QwTriggerScintillator & value )

Definition at line 871 of file QwTriggerScintillator.cc.

References fPMTs.

 {
   if (fPMTs.size() == value.fPMTs.size()){
     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) = value.fPMTs.at(i).at(j);
       }
     }
   } else {
     std::cerr << "QwTriggerScintillator::operator=:  Problems!!!"
         << std::endl;
   }
   return *this;
 }

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

virtual

Implements VQwSubsystem.

Definition at line 340 of file QwTriggerScintillator.cc.

References QwF1TDContainer::AddQwF1TDC(), QwF1TDContainer::DoneF1TDCsConfiguration(), fBankID, fF1RefContainer, fF1TDContainer, VQwSubsystemTracking::fF1TDCResolutionNS, fSoftwareMeantimeContainer, GetModuleIndex(), VQwSubsystem::GetSubbankIndex(), VQwSubsystem::GetSubsystemName(), kMaxNumberOfModulesPerROC, kMaxNumberOfQwHitPlane, QwF1TDContainer::Print(), QwF1TDC::SetF1BankIndex(), QwF1TDC::SetF1SystemName(), QwF1TDC::SetF1TDCBuffer(), QwF1TDC::SetF1TDCIndex(), MeanTimeContainer::SetPlane(), and MeanTimeContainer::SetSystemName().

 {
 
   if( bank_id==fBankID[2] ) {
 
     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);
       fF1RefContainer-> SetSystemName(subsystem_name);
 
       for (Int_t plane_idx=0; plane_idx<kMaxNumberOfQwHitPlane; plane_idx++) 
   {
     fSoftwareMeantimeContainer[plane_idx]->SetSystemName(subsystem_name);
     fSoftwareMeantimeContainer[plane_idx]->SetPlane(plane_idx+1);
   }
 
       if(local_debug) std::cout << "-----------------------------------------------------" << std::endl;
       
       std::cout << "QwTriggerScintillator : " 
     << 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<kMaxNumberOfModulesPerROC; 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 = 0;
 
       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;
       }
     }
   }
 
   return 0;
 }

Here is the call graph for this function:

Int_t QwTriggerScintillator::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 476 of file QwTriggerScintillator.cc.

References QwF1TDContainer::CheckDataIntegrity(), MQwV775TDC::DecodeTDCWord(), MQwF1TDC::DecodeTDCWord(), Qw::e, fBankID, fF1TDCDecoder, fF1TDContainer, FillRawTDCWord(), FillRawWord(), fModulePtrs, fQDCTDC, fRefTime, fRefTime_ChanNum, fRefTime_SlotNum, fSCAs, fSCAs_offset, GetModuleIndex(), VQwSubsystem::GetSubbankIndex(), MQwV775TDC::GetTDCChannelNumber(), MQwF1TDC::GetTDCChannelNumber(), MQwV775TDC::GetTDCData(), MQwF1TDC::GetTDCData(), MQwV775TDC::GetTDCSlotNumber(), MQwF1TDC::GetTDCSlotNumber(), IsF1ReferenceChannel(), IsSlotRegistered(), MQwV775TDC::IsValidDataword(), MQwF1TDC::IsValidDataword(), and VQwSubsystem::SetDataLoaded().

 {
   Int_t index = 0;
 
   index = GetSubbankIndex(roc_id,bank_id);
 
 
   if (bank_id==fBankID[0]) { // V792 or V775
     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 QwTriggerScintillator::FillRawWord: "
           << 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));
       }
     }
   }
 
   else if (bank_id==fBankID[2]) { // F1TDC
     // reset the refrence time 
     fRefTime = 0.0;
 
     Bool_t local_debug_f1 = false;
 
     Int_t  bank_index      = 0;
     UInt_t tdc_data        = 0;
 
 
     Int_t  dummy_slot_number = 0;
     Int_t  dummy_chan_number = 0;
 
     Bool_t data_integrity_flag = false;
     Bool_t temp_print_flag     = false;
 
     UInt_t hit_counter = 0;
 
     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);
   
   
       //
       // 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) {
 
   dummy_slot_number = 0;
   dummy_chan_number = 0;
   hit_counter       = 0;
   if(local_debug_f1) printf("\n");
   for (UInt_t i=0; i<num_words ; i++) {
 
     //  Decode this word as a F1TDC word.
     fF1TDCDecoder.DecodeTDCWord(buffer[i], roc_id); 
       
     // For MQwF1TDC,   roc_id is needed to print out some warning messages.
           
     Int_t tdc_slot_number = fF1TDCDecoder.GetTDCSlotNumber();
     Int_t tdc_chan_number = fF1TDCDecoder.GetTDCChannelNumber();
     //Int_t tdcindex = GetModuleIndex(bank_index, tdc_slot_number);
     
 
     if ( tdc_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, tdc_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 {
     
         tdc_data = fF1TDCDecoder.GetTDCData();
         if (dummy_slot_number == tdc_slot_number && dummy_chan_number == tdc_chan_number) {  
     hit_counter++;
         } 
         else {
     hit_counter = 0;
         }
         
         if(local_debug_f1) {
     printf("TS::ProcessEvBuffer: [%4d] hit counter %d, bank_index %2d slot_number [%2d,%2d] chan [%2d,%2d] data %10d %10.2f\n", i, hit_counter,
            bank_index, tdc_slot_number, fRefTime_SlotNum, tdc_chan_number, fRefTime_ChanNum,tdc_data, fRefTime);
         }
         
         if (hit_counter == 0) {
     FillRawWord   (bank_index, tdc_slot_number, tdc_chan_number, tdc_data);
     if ( IsF1ReferenceChannel(tdc_slot_number,tdc_chan_number) ) {
       fRefTime = (Double_t) tdc_data;
     }
         }
         
         FillRawTDCWord(bank_index, tdc_slot_number, tdc_chan_number, tdc_data);
         
         dummy_slot_number = tdc_slot_number;
         dummy_chan_number = tdc_chan_number;
         //        printf("\n");
   
       }
       catch (std::exception& e) {
         std::cerr << "Standard exception from QwTriggerScintillator::FillRawWord: "
       << e.what() << std::endl;
         std::cerr << "   Parameters:  index==" <<bank_index
       << "; GetF1SlotNumber()=="   <<tdc_slot_number
       << "; GetF1ChannelNumber()=="<<tdc_chan_number
       << "; GetF1Data()=="         <<tdc_data
       << std::endl;
       }
     }//;;
   } // for (UInt_t i=0; i<num_words ; i++) {
       }
     }
   }
 
   return 0;
 }

Here is the call graph for this function:

void QwTriggerScintillator::ProcessEvent ( )

virtual

Implements VQwSubsystem.

Definition at line 655 of file QwTriggerScintillator.cc.

References AddSoftwareMeantimeToHits(), fF1TDContainer, fPMTs, fRefTime, fSCAs, fSoftwareMeantimeOption, VQwSubsystem::HasDataLoaded(), QwF1TDContainer::ReferenceSignalCorrection(), QwF1TDContainer::ReturnTimeCalibration(), SubtractReferenceTimes(), and UpdateHits().

 {
   if (! HasDataLoaded()) return;
 
   TString  elementname = "";
   Double_t rawtime_arb_unit       = 0.0;
   Double_t corrected_time_arb_unit = 0.0;
   Double_t time_ns                 = 0.0;
 
   Int_t    bank_index = 0;
   Int_t    slot_num   = 0;
 
   for (size_t i=0; i<fPMTs.size(); i++) 
     {//;
       for (size_t j=0; j<fPMTs.at(i).size(); j++)
   {//;;
 
     elementname = fPMTs.at(i).at(j).GetElementName();
     
     if (elementname.EndsWith("f1") ) {
       rawtime_arb_unit = fPMTs.at(i).at(j).GetValue(); // returns Double_t
       
       if ( rawtime_arb_unit not_eq 0.0) {
         
         if ( not elementname.Contains("reftime") ) {
         bank_index              = fPMTs.at(i).at(j).GetSubbankID();
         slot_num                = fPMTs.at(i).at(j).GetModule();
         corrected_time_arb_unit = fF1TDContainer->ReferenceSignalCorrection(rawtime_arb_unit, fRefTime, bank_index, slot_num);
     time_ns                 = fF1TDContainer->ReturnTimeCalibration(corrected_time_arb_unit);
     fPMTs.at(i).at(j).SetValue(time_ns);
     // if (bank_index == 2 and  slot_num == 4) {
     //   printf("TS::ProcessBuffer:  bank_index %2d slot_number %2d raw au %10f  reftime au %10.2f, corrected au %10.2f, corrected ns %10.2f\n", 
     //   bank_index, slot_num, rawtime_arb_unit, reftime, corrected_time_arb_unit, time_ns);
     // }
         } 
         else {
     //  // we save the referennce raw time
           fPMTs.at(i).at(j).SetValue(rawtime_arb_unit);
         }
       }
       else {
         fPMTs.at(i).at(j).SetValue(rawtime_arb_unit); // zero
       }
       
     }
     else {
       fPMTs.at(i).at(j).ProcessEvent();
     }
 
 
     // fPMTs.at(i).at(j).ProcessEvent();
     // elementname = fPMTs.at(i).at(j).GetElementName();
     
     // // Check whether the element is "reftime" 
     // if ( not elementname.Contains("reftime") ) {
     //   rawtime_arb_unit = fPMTs.at(i).at(j).GetValue();
 
     //   if (elementname.EndsWith("f1") && rawtime_arb_unit!=0.0) {
 
     //     bank_index = fPMTs.at(i).at(j).GetSubbankID();
     //     slot_num   = fPMTs.at(i).at(j).GetModule();
     //     // if the reference time signal is recorded by (a) channel(s) of F1TDC(s),
     //     // we correct them. And if not, we set them to zero. (jhlee)
 
     //     if ( reftime!=0.0 ) {
     //  corrected_time_arb_unit = fF1TDContainer->ReferenceSignalCorrection(rawtime_arb_unit, reftime, bank_index, slot_num);
     //     }
     //     else {
     //  corrected_time_arb_unit = 0.0;
     //     }
     //     fPMTs.at(i).at(j).SetValue(corrected_time_arb_unit);
     //   } // if (elementname.EndsWith("f1") && rawtime!=0.0) {
     // } 
     // // if ( not elementname.Contains("reftime") ) {
     // //
     // // we keep the raw reference time information in an output ROOT file. 
     // //
   }//;;
     }//;
 
   
   // F1TDCs of  QwHit
   SubtractReferenceTimes();
   UpdateHits();
   AddSoftwareMeantimeToHits(fSoftwareMeantimeOption);
 
 
  // 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:

void QwTriggerScintillator::ProcessOptions ( QwOptions & options )

virtual

Process the command line options.

Reimplemented from VQwSubsystem.

Definition at line 767 of file QwTriggerScintillator.cc.

References fSoftwareMeantimeOption, fSoftwareMeantimeTimeWindowNs, and QwOptions::GetValue().

 {
   fSoftwareMeantimeOption       = options.GetValue<Bool_t>   ( "enable-ts-software-meantime" );
   fSoftwareMeantimeTimeWindowNs = options.GetValue<Double_t> ( "set-ts-software-meantime-timewindow" );
   return;
 };

Here is the call graph for this function:

EQwModuleType QwTriggerScintillator::RegisterModuleType ( TString moduletype )

protected

Definition at line 976 of file QwTriggerScintillator.cc.

References fCurrentModuleIndex, fCurrentType, fModuleTypes, fPMTs, kF1TDC, kSIS3801, kV775_TDC, 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=="V775"){
     fCurrentType = kV775_TDC;
   } else if (moduletype=="F1TDC") {
     fCurrentType = kF1TDC;
   } else if (moduletype=="SIS3801") {
     fCurrentType = kSIS3801;
   }
   fModuleTypes.at(fCurrentModuleIndex) = fCurrentType;
 
   if ((Int_t) fPMTs.size() <= fCurrentType){
     fPMTs.resize(fCurrentType+1);
   }
   return fCurrentType;
 }

Here is the caller graph for this function:

Int_t QwTriggerScintillator::RegisterROCNumber ( const UInt_t roc_id )

protected

Definition at line 901 of file QwTriggerScintillator.cc.

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

Referenced by LoadChannelMap().

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

Here is the call graph for this function:

Here is the caller graph for this function:

Int_t QwTriggerScintillator::RegisterSlotNumber ( const UInt_t slot_id )

protected

Definition at line 936 of file QwTriggerScintillator.cc.

References fCurrentBankIndex, fCurrentModuleIndex, fCurrentSlot, fDetectorIDs, fModuleIndex, fModulePtrs, fModuleTypes, fNumberOfModules, kMaxNumberOfChannelsPerF1TDC, kMaxNumberOfModulesPerROC, and kUnknownModuleType.

Referenced by LoadChannelMap().

 {
   std::pair<EQwModuleType, Int_t> tmppair;
   tmppair.first  = kUnknownModuleType;
   tmppair.second = -1;
  
  if (slot_id < kMaxNumberOfModulesPerROC){
    
    if ( fCurrentBankIndex <= (Int_t) fModuleIndex.size()) {
      
 
      // fModuleTypes and fModulePtrs could be replaced by fDetectorID,
      // but I leave them here while .....
      // Tuesday, January 24 11:49:20 EST 2012, jhlee
 
       fModuleTypes.resize(fNumberOfModules+1);
       fModulePtrs.resize(fNumberOfModules+1);
       fModulePtrs.at(fNumberOfModules).resize(kMaxNumberOfChannelsPerF1TDC,
                 tmppair);
 
       //   fNumberOfModules = fModulePtrs.size();
      
       fDetectorIDs.resize(fNumberOfModules+1);
       fDetectorIDs.at(fNumberOfModules).resize(kMaxNumberOfChannelsPerF1TDC);
 
       // // reassign kNumberOfModules after resize it
 
       fNumberOfModules = (Int_t) fDetectorIDs.size();
       fCurrentModuleIndex = fNumberOfModules-1;
       fCurrentSlot        = slot_id;
       fModuleIndex.at(fCurrentBankIndex).at(fCurrentSlot) = fCurrentModuleIndex;
     }
   } else {
     std::cerr << "QwTriggerScintillator::RegisterSlotNumber:  Slot number "
         << slot_id << " is larger than the number of slots per ROC, "
         << kMaxNumberOfModulesPerROC << std::endl;
   }
   return fCurrentModuleIndex;
 }

Here is the caller graph for this function:

Int_t QwTriggerScintillator::RegisterSubbank ( const UInt_t bank_id )

protected

Definition at line 910 of file QwTriggerScintillator.cc.

References fCurrentBankIndex, VQwSubsystem::fCurrentROC_ID, fModuleIndex, fReferenceChannels, fReferenceData, VQwSubsystem::GetSubbankIndex(), kMaxNumberOfModulesPerROC, 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();
 
   fCurrentBankIndex = GetSubbankIndex(current_roc_id, bank_id);
 
   if ( (Int_t) reference_size <= fCurrentBankIndex) {
     fReferenceChannels.resize(fCurrentBankIndex+1);
     fReferenceData.resize(fCurrentBankIndex+1);
   }
 
   std::vector<Int_t> tmpvec(kMaxNumberOfModulesPerROC,-1);
   fModuleIndex.push_back(tmpvec);
 
   std::cout<< "QwTriggerScintillator::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 QwTriggerScintillator::ReportConfiguration ( Bool_t verbose )

Definition at line 1548 of file QwTriggerScintillator.cc.

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

Referenced by LoadChannelMap().

 {
 
   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 << "QwTriggerScintillator 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<kMaxNumberOfModulesPerROC; 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 "
           << GetSubsystemName()
           << 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 QwTriggerScintillator::SubtractReferenceTimes ( )

protected

Definition at line 1145 of file QwTriggerScintillator.cc.

References fF1RefContainer, fF1TDContainer, fTDCHits, 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;  
 
   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;
   
       }
     }
   // std::vector<Double_t> reftimes;
   // std::vector<Bool_t>   refchecked;
   // std::vector<Bool_t>   refokay;
   // Bool_t allrefsokay;
 
 
   // std::size_t ref_size = 0;
   // std::size_t i = 0;
   // std::size_t j = 0;
 
   // ref_size = fReferenceData.size();
 
   // reftimes.resize  ( ref_size );
   // refchecked.resize( ref_size );
   // refokay.resize   ( ref_size );
 
   // for ( i=0; i<ref_size; i++ ) {
   //   reftimes.at(i)   = 0.0;
   //   refchecked.at(i) = kFALSE;
   //   refokay.at(i)    = kFALSE;
   // }
 
   // allrefsokay = kTRUE;
 
 
 
   // 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;
 
   // Bool_t local_debug = false;
 
   // for ( std::vector<QwHit>::iterator hit=fTDCHits.begin(); hit!=fTDCHits.end(); hit++ ) 
   //   {
   //     //  Only try to check the reference time for a bank if there is at least one
   //     //  non-reference hit in the bank.
   //     bank_index = hit->GetSubbankID();
 
   //     //     if(local_debug) printf("QwHit :: bank index %d\n", bank_index);
 
   //     if ( not refchecked.at(bank_index) ) {
   
   //  if ( fReferenceData.at(bank_index).empty() ) {
   //    std::cout << "QwTriggerScintillator::SubtractReferenceTimes:  Subbank ID "
   //            << bank_index << " is missing a reference time." << std::endl;
   //    refokay.at(bank_index) = kFALSE;
   //    allrefsokay            = kFALSE;
   //  }
   //  else {
   //    if(fReferenceData.at(bank_index).size() not_eq 1) {
   //      std::cout << "Multiple hits are recorded in the reference channel, we use the first hit signal as the refererence signal." << std::endl;
   //    }
   //    reftimes.at(bank_index) = fReferenceData.at(bank_index).at(0);
   //    refokay.at(bank_index)  = kTRUE;
   //  }
 
   //  if (refokay.at(bank_index)){
   //    for ( j=0; j<fReferenceData.at(bank_index).size(); j++ ) 
   //      {
   //        // printf("Reference time %f fReferenceData.at(%d).at(%d) %f\n", 
   //        //       reftimes.at(bank_index), (Int_t) bank_index, (Int_t) j, fReferenceData.at(bank_index).at(j));
   //        fReferenceData.at(bank_index).at(j) -= reftimes.at(bank_index);
   //        // printf("Reference time %f fReferenceData.at(%d).at(%d) %f\n", 
   //        //       reftimes.at(bank_index), (Int_t) bank_index, (Int_t) j, fReferenceData.at(bank_index).at(j));
   //      }
   //  }
 
   //      refchecked.at(bank_index) = kTRUE;
   //     }
       
   //     if ( refokay.at(bank_index) ) {
   //    Int_t slot_num    = hit -> GetModule();
   //      raw_time_arb_unit = (Double_t) hit -> GetRawTime();
   //      ref_time_arb_unit = (Double_t) reftimes.at(bank_index);
 
   //      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;
     
   //      }
   //     }
   //   }
   
   // // bank_index = 0;
   
   // if ( not allrefsokay ) {
   //   std::vector<QwHit> tmp_hits;
   //   tmp_hits.clear();
   //   for ( std::vector<QwHit>::iterator hit=fTDCHits.begin(); hit!=fTDCHits.end(); hit++ ) 
   //     {
   //  bank_index = hit->GetSubbankID();
   //  if ( refokay.at(bank_index) ) tmp_hits.push_back(*hit);
   //     }
   //   //    std::cout << "FTDC size " << fTDCHits.size() << "tmp hit size " << tmp_hits.size() << std::endl;
   //   fTDCHits.clear();
   //   fTDCHits = tmp_hits;
   //   //     std::cout << "FTDC size " << fTDCHits.size() << "tmp hit size " << tmp_hits.size() << std::endl;
   // }
   
   return;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwTriggerScintillator::UpdateHits ( )

protected

Definition at line 1309 of file QwTriggerScintillator.cc.

References VQwSubsystemTracking::fF1TDCResolutionNS, and fTDCHits.

Referenced by ProcessEvent().

 {
  
   // std::vector... is weird
   // I use a plane number to assign a f1tdc channel. Channel 99 is the F1TDC reference channel
   // , becaue I see plan is 0 in qweak_new.geo file and I think, nowhere plane is used in analyzer
   // But, I am wrong. This double std::vector fDetectorInfo uses plane number as second index,
   // So I had the notorious following error
   //  Terminate called after throwing an instance of 'std::out_of_range'
   //   what():  vector::_M_range_check
   // Abort
   //
   //
   // I don't care about QwDetectorInfo for other tracking subsystems except VDC/HDC,
   // I disable to fill QwDetectorInfo into QwHit
 
   // Thursday, January 26 10:30:09 EST 2012, jhlee
 
 
   // EQwDetectorPackage package = kPackageNull;
   // //  Int_t              plane   = 0;
 
   // QwDetectorID       local_id;
   // QwDetectorInfo  *local_info;
 
   for(std::vector<QwHit>::iterator iter=fTDCHits.begin(); iter!=fTDCHits.end(); ++iter)
     {
       // local_id   = iter->GetDetectorID();
       // package    = local_id.fPackage;
       // //      plane      = local_id.fPlane;
       // // For TS, we have only 1 plane in each package according to its geometry, but I use a plane
       // // number to assign a f1tdc channel. Channel 99 is not assigned into this vector
       // // 
       // local_info = fDetectorInfo.in(package).at(0); 
       // //      local_info = fDetectorInfo.in(package).at(plane); 
 
       // iter->SetDetectorInfo(local_info);
       iter->ApplyTimeCalibration(fF1TDCResolutionNS); // Fill fTimeRes and fTimeNs in QwHit
 
     }
  
   return;
 }