#include <QwDriftChamberVDC.h>

Inheritance diagram for QwDriftChamberVDC:

Collaboration diagram for QwDriftChamberVDC:

Public Member Functions
	QwDriftChamberVDC (TString name)

virtual	~QwDriftChamberVDC ()

void	SubtractReferenceTimes ()

void	ProcessEvent ()

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

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

void	ReadEvent (TString &)

void	ClearEventData ()

Public Member Functions inherited from QwDriftChamber
	QwDriftChamber (const TString &region_tmp)
	Constructor with name. More...

	QwDriftChamber (const TString &region_tmp, std::vector< QwHit > &fWireHits_TEMP)
	Constructor with name and hit list. More...

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

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

void	FillDriftDistanceToHits ()

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

void	GetTDCHitList (QwHitContainer &grandHitContainer)

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	ConstructBranchAndVector (TTree *tree, TString &prefix, std::vector< Double_t > &values)
	Construct the branch and tree vector. 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...

virtual void	FillTreeVector (std::vector< Double_t > &values) const
	Fill the tree vector. 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, QwDriftChamberVDC >
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 QwDriftChamber
static void	DefineOptions (QwOptions &options)
	Define options. More...

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

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

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

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

Int_t	AddChannelDefinition ()

Int_t	BuildWireDataStructure (const UInt_t chan, const EQwDetectorPackage package, const Int_t octant, const Int_t plane, const Int_t wire)

Double_t	CalculateDriftDistance (Double_t drifttime, QwDetectorID detector)

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

Int_t	LoadTimeWireOffset (TString t0_map)

void	LoadTtoDParameters (TString ttod_map)

void	SubtractWireTimeOffset ()

void	UpdateHits ()

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

Protected Member Functions inherited from QwDriftChamber
Int_t	LinkReferenceChannel (const UInt_t chan, const Int_t plane, const Int_t wire)

void	ClearAllBankRegistrations ()

Int_t	RegisterROCNumber (const UInt_t roc_id, const UInt_t bank_id)
	Tell the object that it will decode data from this ROC and sub-bank. More...

Int_t	RegisterSubbank (const UInt_t bank_id)

Int_t	RegisterSlotNumber (const UInt_t slot_id)

Int_t	GetTDCIndex (size_t bank_index, size_t slot_num) const

Bool_t	IsSlotRegistered (Int_t bank_index, Int_t slot_num) const

void	InitHistogramPointers ()

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

Int_t	RegisterSubbank (const UInt_t bank_id)
	Tell the object that it will decode data from this sub-bank in the ROC currently open for registration. More...

Int_t	GetSubbankIndex () const

Int_t	GetSubbankIndex (const UInt_t roc_id, const UInt_t bank_id) const

void	SetDataLoaded (Bool_t flag)

Int_t	FindIndex (const std::vector< UInt_t > &myvec, const UInt_t value) const

Bool_t	Compare (VQwSubsystem *source)

Protected Member Functions inherited from MQwHistograms
	MQwHistograms ()
	Default constructor. More...

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

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

virtual MQwHistograms &	operator= (const MQwHistograms &value)

void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)

void	AddHistogram (TH1 *h)
	Register a histogram. More...

Protected Attributes
Bool_t	fUseTDCHits

Bool_t	fDisableWireTimeOffset

Int_t	fR3Octant

std::vector< std::vector < QwDelayLine > >	fDelayLineArray

std::vector< std::vector < QwDelayLineID > >	fDelayLinePtrs

std::vector< QwHit >	fWireHitsVDC

std::vector< std::vector < std::vector< Double_t > > >	fTimeWireOffsets

std::vector< Double_t >	fTtoDNumbers

Protected Attributes inherited from QwDriftChamber
TString	fRegion

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

Int_t	fCurrentSlot

Int_t	fCurrentTDCIndex

UInt_t	kMaxNumberOfChannelsPerTDC

Int_t	fNumberOfTDCs

std::vector< std::vector< Int_t > >	fTDC_Index

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

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

std::vector< std::vector < Double_t > >	fReferenceMaster

std::vector< QwHit >	fTDCHits

std::vector< QwHit > &	fWireHits

std::vector< Int_t >	fWiresPerPlane

MQwF1TDC	fF1TDCDecoder

QwF1TDContainer *	fF1TDContainer

F1TDCReferenceContainer *	fF1RefContainer

std::vector< TH1F * >	TotHits

std::vector< TH1F * >	TOFP

std::vector< TH1F * >	TOFP_raw

std::vector< TH1F * >	WiresHit

std::vector< TH2F * >	TOFW

std::vector< TH2F * >	TOFW_raw

std::vector< TH2F * >	HitsWire

std::vector< std::vector < QwDetectorID > >	fTDCPtrs

std::vector< std::vector < QwDetectorInfo > >	fWireData

std::vector< std::vector < UInt_t > >	fDirectionData

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	kBackPlaneNum =16

static const UInt_t	kLineNum =8

Static Protected Attributes inherited from QwDriftChamber
static bool	fPrintF1TDCConfiguration = true

static const UInt_t	kMaxNumberOfSlotsPerROC = 21

static const Int_t	kReferenceChannelPlaneNumber = 99

static const Int_t	kCodaMasterPlaneNumber = 98

Detailed Description

Definition at line 22 of file QwDriftChamberVDC.h.

Constructor & Destructor Documentation

QwDriftChamberVDC::QwDriftChamberVDC ( TString name )

Definition at line 27 of file QwDriftChamberVDC.cc.

References fDelayLineArray, fDelayLinePtrs, kBackPlaneNum, and kLineNum.

                                                  :
   VQwSubsystem(name), QwDriftChamber(name,fWireHitsVDC )
 {
   std::vector<QwDelayLine> temp;
   temp.clear();
   temp.resize ( kLineNum );
   fDelayLineArray.resize ( kBackPlaneNum, temp );
   fDelayLinePtrs.resize ( 21 );
 }

virtual QwDriftChamberVDC::~QwDriftChamberVDC ( )

inlinevirtual

Definition at line 30 of file QwDriftChamberVDC.h.

30 { };

Member Function Documentation

Int_t QwDriftChamberVDC::AddChannelDefinition ( )

protectedvirtual

Implements QwDriftChamber.

Definition at line 443 of file QwDriftChamberVDC.cc.

References QwDriftChamber::fTDC_Index, QwDriftChamber::fTDCPtrs, QwDriftChamber::fWireData, QwDriftChamber::fWiresPerPlane, QwDriftChamber::kCodaMasterPlaneNumber, QwDriftChamber::kReferenceChannelPlaneNumber, and OK.

Referenced by LoadChannelMap().

 {
   bool temp_local_debug = false;
   //   if (temp_local_debug){
   //     std::cout << " QwDriftChamberVDC::AddChannelDefinition"<<std::endl;
   //     std::cout << "plane " << plane << " wire " << wire << std::endl;
   //   }
 
   std::size_t i =0;
 
   fWireData.resize ( fWiresPerPlane.size() );
   for ( i=1; i<fWiresPerPlane.size(); i++ )
     {
       if ( temp_local_debug )
   {
     std::cout << "wire #" << i
         << " " << fWiresPerPlane.at ( i )
         << std::endl;
   }
       fWireData.at ( i ).resize ( fWiresPerPlane.at ( i ) );
     }
 
   Int_t mytdc = 0;
   for ( i=0; i<fTDC_Index.size(); i++ )
     {
       for ( size_t j=0; j<fTDC_Index.at ( i ).size(); j++ )
   {
     mytdc = fTDC_Index.at ( i ).at ( j );
     if ( mytdc not_eq -1 )
       {
         for ( size_t k=0; k<fTDCPtrs.at ( mytdc ).size(); k++ )
     {
       Int_t package = fTDCPtrs.at ( mytdc ).at ( k ).fPackage;
       Int_t plane   = fTDCPtrs.at ( mytdc ).at ( k ).fPlane;
       if ( ( plane>0 ) and ( plane not_eq ( Int_t ) kReferenceChannelPlaneNumber ) and ( plane not_eq ( Int_t ) kCodaMasterPlaneNumber) )
         {
           Int_t wire  = fTDCPtrs.at ( mytdc ).at ( k ).fElement;
           Int_t window_number=0;
           if ( wire<9 ) window_number=18;
           else if ( wire==152 ) window_number=16;
           else window_number=17;
           for ( Int_t index=0;index<window_number;index++ )
       {
 
         fWireData.at ( 4* ( package-1 ) +plane ).at ( wire-1 ).SetElectronics ( i,j,k );
         wire+=8;
         if ( temp_local_debug )
           {
             std::cout << "mytdc " << mytdc
           << "wire #" << wire
           << " plane  " << plane
           << " (i j k) (" << i  <<" " << j << " "<< k << ")"
           << std::endl;
           }
       }
         }
     }
       }
   }
     }
   // std::cout << " QwDriftChamberVDC::AddChannelDefinition END"<<std::endl;
   return OK;
 }

Here is the caller graph for this function:

Int_t QwDriftChamberVDC::BuildWireDataStructure	(	const UInt_t	chan,
		const EQwDetectorPackage	package,
		const Int_t	octant,
		const Int_t	plane,
		const Int_t	wire
	)

protectedvirtual

Implements QwDriftChamber.

Definition at line 397 of file QwDriftChamberVDC.cc.

References QwDriftChamber::fCurrentTDCIndex, QwDriftChamber::fTDCPtrs, QwDriftChamber::fWiresPerPlane, and OK.

Referenced by LoadChannelMap().

 {
 
   Int_t r3_wire_number_per_plane = 279;
 
 
   // We should not see when plane is equal to "kReferenceChannelPlaneNumber"
   //, because in LocalChannelMap, we assigned "reference channels" already into whatever objects.
   // And Since it is the same structure as HDC, but R3 uses its own method to select the refernece signals, 
   // which I want to undestand how and why soon, all if conditions are invaild for Region 3. 
   // Wednesday, February  8 10:51:20 EST 2012, jhlee
   
   // if ( plane == kReferenceChannelPlaneNumber )
   //   {
   //     printf("Can we See this?t\n");
   //     LinkReferenceChannel ( chan, plane, wire );
   //   }
   // else
   //   {
   fTDCPtrs.at ( fCurrentTDCIndex ).at ( chan ).fOctant  = octant;
   fTDCPtrs.at ( fCurrentTDCIndex ).at ( chan ).fPackage = package;
   fTDCPtrs.at ( fCurrentTDCIndex ).at ( chan ).fPlane   = plane;
   fTDCPtrs.at ( fCurrentTDCIndex ).at ( chan ).fElement = wire;
   
   Int_t real_size = 0;
   real_size = 4 * ( ( Int_t ) package-1 ) +plane;
   if ( real_size >= ( Int_t ) fWiresPerPlane.size() ) // plane is Int_t
     {
       fWiresPerPlane.resize ( real_size + 1 );
       // size() is one more larger than last plane number
       // For VDC, plane    1,2,....,8
       //          vector 0,1,2,....,8
       //          thus, vector.size() returns 9
       // So the magic number "1" is.
       // Wednesday, July 28 21:58:13 EDT 2010, jhlee
     }
   fWiresPerPlane.at ( real_size ) = r3_wire_number_per_plane;
   //    }
   return OK;
 }

Here is the caller graph for this function:

Double_t QwDriftChamberVDC::CalculateDriftDistance	(	Double_t	drifttime,
		QwDetectorID	detector
	)

protectedvirtual

Implements QwDriftChamber.

Definition at line 276 of file QwDriftChamberVDC.cc.

References fTtoDNumbers.

 {
   //    Double_t angle_degree = 45.0;
   Double_t distance_mm = 0.0;
   Double_t distance_cm = 0.0;
   Double_t dt          = drifttime;
   Double_t resolution  = 0.5;     //resolution is 0.5ns
 
   Int_t index = (Int_t) (dt/resolution);
   
   Int_t maximum_allowed_drifttime_by_john = 310;
   Int_t minimum_allowed_drifttime_by_john = 0;
  
   if(dt >= minimum_allowed_drifttime_by_john && dt < maximum_allowed_drifttime_by_john) {
     distance_mm = ( dt-resolution*index ) / resolution * ( fTtoDNumbers.at ( index+1 )-fTtoDNumbers.at ( index ) ) + fTtoDNumbers.at ( index );
   }
   else {
     distance_mm=-50;
   }
   distance_cm = 0.1*distance_mm;
 
   return distance_cm;
 }

void QwDriftChamberVDC::ClearEventData ( )

virtual

Implements QwDriftChamber.

Definition at line 866 of file QwDriftChamberVDC.cc.

References F1TDCReferenceContainer::ClearEventData(), fDelayLineArray, QwDriftChamber::fF1RefContainer, QwDriftChamber::fReferenceData, QwDriftChamber::fReferenceMaster, QwDriftChamber::fTDCHits, QwDriftChamber::fWireHits, and VQwSubsystem::SetDataLoaded().

 {
   SetDataLoaded ( kFALSE );
 
   fTDCHits.clear();
   fWireHits.clear();
 
   std::size_t index = 0;
   std::size_t j = 0;
   std::size_t i = 0;
 
   for ( i=0;i<fDelayLineArray.size();i++ )
     {
       index = fDelayLineArray.at ( i ).size();
 
       for ( j=0;j<index;j++ )
   {
     fDelayLineArray.at ( i ).at ( j ).Processed=kFALSE;
     fDelayLineArray.at ( i ).at ( j ).LeftHits.clear();
     fDelayLineArray.at ( i ).at ( j ).RightHits.clear();
     fDelayLineArray.at ( i ).at ( j ).Hitscount.clear();
     fDelayLineArray.at ( i ).at ( j ).Wire.clear();
   }
     }
   
   fF1RefContainer->ClearEventData();
 
   for ( i=0;i<fReferenceData.size(); i++ )
     {
       fReferenceData.at ( i ).clear();
     }
   for( i=0;i<fReferenceMaster.size();i++)
     fReferenceMaster.at(i).clear();
   return;
 }

Here is the call graph for this function:

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

protectedvirtual

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

Implements QwDriftChamber.

Definition at line 926 of file QwDriftChamberVDC.cc.

References QwDriftChamber::fWiresPerPlane, VQwSubsystem::GetSubsystemName(), QwDriftChamber::HitsWire, QwDriftChamber::TOFP, QwDriftChamber::TOFP_raw, QwDriftChamber::TOFW, QwDriftChamber::TOFW_raw, QwDriftChamber::TotHits, and QwDriftChamber::WiresHit.

 {
   //  If we have defined a subdirectory in the ROOT file, then change into it.
   if ( folder ) folder->cd();
   //  Now create the histograms...
   TString region = GetSubsystemName();
   //  Loop over the number of planes.
 
   const Short_t buffer_size  = 2000;
   Float_t bin_offset = -0.5;
 
   std::size_t total_plane_number = 0;
   total_plane_number = fWiresPerPlane.size();
 
   TotHits.resize ( total_plane_number );
   TOFP.resize ( total_plane_number );
   TOFP_raw.resize ( total_plane_number );
   WiresHit.resize ( total_plane_number );
   TOFW.resize ( total_plane_number );
   TOFW_raw.resize ( total_plane_number );
   HitsWire.resize ( total_plane_number );
 
   std::size_t iplane = 0;
   // std::cout <<  "QwDriftChamberVDC::ConstructHistograms, "
   // <<  "we are contructing histograms with index from 0 to " <<total_plane_number
   // << "\n"
   // <<  "Thus, fWiresPerPlane.size() returns "
   // << total_plane_number
   // << " and its array definition is ["
   // << total_plane_number
   // << "]."
   // <<  " And hist[i] <-> hist.at(i) <-> fWiresPerplane[i] <-> fWiresPerPlane.at(i)"
   // << std::endl;
 
   // wire_per_plane is the number of wire per plane?
   //
   // we skip the first zero-th plane or wire histogram. thus
   // i starts with '1'. hist[0] is NULL
 
   for ( iplane=1; iplane< total_plane_number; iplane++ )
     {
       //   std::cout << "QwDriftChamberVDC iplane" << iplane << std::endl;
       // push_back can "push" iplane = 1 into TotHits.at(0) ??
       TotHits[iplane] = new TH1F (
           Form ( "%s%sWiresPlane%zu", prefix.Data(), region.Data(), iplane ),
           Form ( "Total hits on all wires in plane %zu", iplane ),
           fWiresPerPlane[iplane], bin_offset, fWiresPerPlane[iplane]+bin_offset
           );
 
       TotHits[iplane]->GetXaxis()->SetTitle ( "Wire #" );
       TotHits[iplane]->GetYaxis()->SetTitle ( "Events" );
 
       WiresHit[iplane] = new TH1F (
            Form ( "%s%sNumberHitsPlane%zu", prefix.Data(), region.Data(), iplane ),
            Form ( "Number of Wires Hit in plane %zu",iplane ),
            20, bin_offset, 20+bin_offset
            );
       WiresHit[iplane]->GetXaxis()->SetTitle ( "Wires Hit per Event" );
       WiresHit[iplane]->GetYaxis()->SetTitle ( "Events" );
 
       HitsWire[iplane] = new TH2F (
            Form ( "%s%sNumberHitsVsWirePlane%zu", prefix.Data(), region.Data(), iplane ),
            Form ( "hits on all wires per event in plane %zu", iplane ),
            fWiresPerPlane[iplane],bin_offset,fWiresPerPlane[iplane]+bin_offset,
            7, -bin_offset, 7-bin_offset
            );
       HitsWire[iplane]->GetXaxis()->SetTitle ( "Wire Number" );
       HitsWire[iplane]->GetYaxis()->SetTitle ( "Hits" );
 
       TOFP[iplane] = new TH1F (
              Form ( "%s%sTimePlane%zu", prefix.Data(), region.Data(), iplane ),
              Form ( "Subtracted time of flight for events in plane %zu", iplane ),
              310,0,0
              );
       TOFP[iplane] -> SetDefaultBufferSize ( buffer_size );
       TOFP[iplane] -> GetXaxis()->SetTitle ( "Time of Flight" );
       TOFP[iplane] -> GetYaxis()->SetTitle ( "Hits" );
 
 
       TOFP_raw[iplane] = new TH1F (
            Form ( "%s%sRawTimePlane%zu", prefix.Data(), region.Data(), iplane ),
            Form ( "Raw time of flight for events in plane %zu", iplane ),
            //          400,-65000,65000);
            310, 0,0
            );
       TOFP_raw[iplane] -> SetDefaultBufferSize ( buffer_size );
       TOFP_raw[iplane]->GetXaxis()->SetTitle ( "Time of Flight" );
       TOFP_raw[iplane]->GetYaxis()->SetTitle ( "Hits" );
 
       TOFW[iplane] = new TH2F (
              Form ( "%s%sTimeVsWirePlane%zu", prefix.Data(), region.Data(), iplane ),
              Form ( "Subtracted time of flight for each wire in plane %zu", iplane ),
              fWiresPerPlane[iplane], bin_offset, fWiresPerPlane[iplane]+bin_offset,
              100,-2000,2000
              );
       // why this range is not -65000 ??
       TOFW[iplane]->GetXaxis()->SetTitle ( "Wire Number" );
       TOFW[iplane]->GetYaxis()->SetTitle ( "Time of Flight" );
 
       TOFW_raw[iplane] = new TH2F (
            Form ( "%s%sRawTimeofFlightperWirePlane%zu", prefix.Data() ,region.Data(),iplane ),
            Form ( "Raw time of flight for each wire in plane %zu",iplane ),
            fWiresPerPlane[iplane], bin_offset, fWiresPerPlane[iplane]+bin_offset,
            100,-40000,65000
            );
       // why this range is not -65000 ??
       TOFW_raw[iplane]->GetXaxis()->SetTitle ( "Wire Number" );
       TOFW_raw[iplane]->GetYaxis()->SetTitle ( "Time of Flight" );
     }
   return;
 }

Here is the call graph for this function:

void QwDriftChamberVDC::DefineOptions ( QwOptions & options )

static

Definition at line 905 of file QwDriftChamberVDC.cc.

References QwOptions::AddOptions(), and default_bool_value.

 {
   options.AddOptions("Tracking options") ( "use-tdchit",
       po::value<bool>()->default_bool_value(false),
       "create TDC-based tree" );
   options.AddOptions("Tracking options") ( "disable-wireoffset",
       po::value<bool>()->default_bool_value(false),
       "disable subtraction of t0 for every wire" );
   options.AddOptions("Tracking options") ("R3-octant",
       po::value<int>()->default_value(1),
       "MD Package 2 of R3 is in front of" );
 }

Here is the call graph for this function:

void QwDriftChamberVDC::FillHistograms ( )

protectedvirtual

Fill the histograms for this subsystem.

Implements QwDriftChamber.

Definition at line 1040 of file QwDriftChamberVDC.cc.

References QwDetectorInfo::ClearHits(), QwLog::endl(), QwDetectorID::fElement, QwDetectorID::fPlane, QwDriftChamber::fTDCHits, QwDriftChamber::fWireData, QwDriftChamber::fWireHits, QwDriftChamber::fWiresPerPlane, QwDetectorInfo::GetNumHits(), VQwSubsystem::HasDataLoaded(), QwDriftChamber::HitsWire, QwMessage, QwDriftChamber::TOFP, QwDriftChamber::TOFP_raw, QwDriftChamber::TOFW, QwDriftChamber::TOFW_raw, QwDriftChamber::TotHits, and QwDriftChamber::WiresHit.

 {
 
   Bool_t local_debug = false;
   if ( not HasDataLoaded() ) return;
 
   QwDetectorID   this_detid;
   QwDetectorInfo *this_det;
 
   std::vector<Int_t> wireshitperplane ( fWiresPerPlane.size(),0 );
 
   UInt_t raw_time = 0;
   Double_t time   = 0.0;
 
   Int_t plane = 0;
   Int_t element = 0;
   EQwDetectorPackage package = kPackage1; // this is weird..... 
 
   Int_t plane_index = 0;
 
   std::vector<QwHit>::iterator end=fTDCHits.end();
   for ( std::vector<QwHit>::iterator hit=fTDCHits.begin(); hit!=end; hit++ )
     {
 
       this_detid = hit->GetDetectorID();
       plane      = this_detid.fPlane;
       element    = this_detid.fElement;
       //      package    = //kthis_detid.fPackage;
 
       // I guess, plane_index is calculated from kPackageUp instead of kPackageUp and kPackageDown
       // See the FillRawTDCWordtop
       plane_index = 4* ( ( Int_t ) package -1 ) + plane;
 
       if ( plane<=0 or element<=0 )
     {
       if ( local_debug )
         {
           QwMessage << "QwDriftChamberVDC::FillHistograms:  Bad plane or element index:"
         << "  fPlane = "  << plane
         << "  fElement= " << element
         << QwLog::endl;
         }
       continue;
     }
 
       raw_time = hit->GetRawTime();
       time     = hit->GetTimeNs();
 
       //  Fill ToF histograms
       TOFP_raw[plane_index]->Fill ( raw_time );
       TOFW_raw[plane_index]->Fill ( element, raw_time );
 
 
       if ( local_debug )
     {
       QwMessage << " QwDriftChamberVDC::FillHistograms() plane " << plane
           << " element " << element
           << " package " << package
           << " plane_index " << plane_index
           << QwLog::endl;
     }
     }
 
   end=fWireHits.end();
   for ( std::vector<QwHit>::iterator hit1=fWireHits.begin(); hit1!=end; hit1++ )
     {
 
       this_detid = hit1->GetDetectorID();
       plane      = this_detid.fPlane;
       element    = this_detid.fElement;
       package    = this_detid.fPackage;
 
       if ( plane<=0 or element<=0 )
     {
       if ( local_debug )
         {
           QwMessage << "QwDriftChamberVDC::FillHistograms:  Bad plane or element index:"
         << "  fPlane = "  << plane
         << "  fElement= " << element
         << QwLog::endl;
         }
       continue;
     }
 
       plane_index = 4* ( ( Int_t ) package -1 ) + plane;
       this_det   = & ( fWireData.at ( plane_index ).at ( element-1 ) );
 
 
       if ( hit1->IsFirstDetectorHit() )
     {
       //  If this is the first hit for this detector, then let's plot the
       //  total number of hits this wire had.
       //     std::cout << "this_det->GetNumHits: " << this_det->GetNumHits() << std::endl;
       HitsWire[plane_index]->Fill ( element,this_det->GetNumHits() );
 
       // //     //  Also increment the total number of events in whichthis wire was hit.
       TotHits[plane_index]->Fill ( element,1 );
       //  Increment the number of wires hit in this plane
       wireshitperplane.at ( plane_index ) += 1;
     }
       //  Fill ToF histograms
 
       time     = hit1->GetTimeNs();
 
       TOFP[plane_index]->Fill ( time );
       TOFW[plane_index]->Fill ( element,time );
       this_det->ClearHits();
 
 
     }
 
   std::size_t iplane = 0;
 
   for ( iplane=1; iplane<fWiresPerPlane.size(); iplane++ )
     {
       WiresHit[iplane]->Fill ( wireshitperplane[iplane] );
     }
   return;
 }

Here is the call graph for this function:

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

protectedvirtual

Todo:: Should this direction be properly initialized by fTDCPtrs.at(tdcindex).at(chan).fDirection ?

Implements QwDriftChamber.

Definition at line 301 of file QwDriftChamberVDC.cc.

References fDelayLinePtrs, QwDriftChamber::fF1RefContainer, QwDriftChamber::fReferenceData, QwDriftChamber::fReferenceMaster, QwDriftChamber::fTDCHits, QwDriftChamber::fTDCPtrs, QwDriftChamber::GetTDCIndex(), QwDriftChamber::kCodaMasterPlaneNumber, kDirectionNull, QwDriftChamber::kReferenceChannelPlaneNumber, kRegionID3, F1TDCReferenceContainer::SetReferenceSignal(), and QwHit::WireMatches().

 {
   Bool_t local_debug = false;
   
   if (local_debug) printf("\n");
   Int_t tdcindex = 0;
   tdcindex = GetTDCIndex ( bank_index,slot_num );
   if ( tdcindex not_eq -1 )
     {
       Int_t hitCount = 1;
 
       /// \todo Should this direction be properly initialized by fTDCPtrs.at(tdcindex).at(chan).fDirection ?
       EQwDirectionID direction   = kDirectionNull;
 
       fF1RefContainer->SetReferenceSignal(bank_index, slot_num, chan, data, false);
 
       Int_t octant  = fTDCPtrs.at(tdcindex).at(chan).fOctant;
       Int_t plane   = fTDCPtrs.at(tdcindex).at(chan).fPlane;
       Int_t wire    = fTDCPtrs.at(tdcindex).at(chan).fElement;
 
       // And this below line was missing.
       EQwDetectorPackage package = fTDCPtrs.at(tdcindex).at(chan).fPackage;
       // But I did the correct way to assign the correct package in fTDCPtrs
       //  *** Break *** segmentation violation
       // due to FillHistograms()
       // I guess, plane_index is calculated from kPackageUp instead of kPackageUp and kPackageDown
       // In here I use kPackageNull as an initial value,
       // and save the correct package number in fTDCPtrs in order to use
       // the correct reference time according to package number,
       // and use kpackageUp instead of access fPackage in fTDCPtrs in FillHistograms.
       // Monday, February 13 14:27:50 EST 2012, jhlee
 
 
 
       if ( plane == -1 or wire == -1 )
   {
     //  This channel is not connected to anything.  Do nothing.
   }
       else if ( plane == kReferenceChannelPlaneNumber )
   {
     if (local_debug)  {
       std::cout << " -------------------------"
           << " plane " << plane 
           << " slot "  << slot_num 
           << " wire "  << wire
           << " data " << data 
           << std::endl;
     }
     
     fReferenceData.at ( wire ).push_back ( data );
   }
       else if ( plane == kCodaMasterPlaneNumber )
   {
     if (local_debug)  {
       std::cout << " -------------------------" << " plane " << plane 
           << " slot "  << slot_num 
           << " wire "  << wire
           << " data " << data 
           << std::endl;
     }
     fReferenceMaster.at ( wire ).push_back ( data );
   }
       else
   {
     plane=fDelayLinePtrs.at ( slot_num ).at ( chan ).fBackPlane;
     //std::cout<<"At QwDriftChamberVDC::FillRawTDCWord_1"<<endl;
     hitCount = std::count_if ( fTDCHits.begin(), fTDCHits.end(),
              boost::bind (
               &QwHit::WireMatches,_1, kRegionID3, boost::ref ( package ), boost::ref ( plane ), boost::ref ( wire )
               )
              );
 
     fTDCHits.push_back (
             QwHit (
              bank_index,
              slot_num,
              chan,
                                      hitCount,
                                      kRegionID3,
                                      package,
                                      octant,
                                      plane,
                                      direction,
                                      wire,
              data
              )
             );
     //in order-> bank index, slot num, chan, hitcount, region=3, package, plane,,direction, wire,wire hit time
 
   }
 
     }
 }

Here is the call graph for this function:

void QwDriftChamberVDC::GetHitList ( QwHitContainer & grandHitContainer )

inlineprotectedvirtual

Get the hit list.

Reimplemented from QwDriftChamber.

Definition at line 76 of file QwDriftChamberVDC.h.

References QwHitContainer::Append(), QwDriftChamber::fTDCHits, fUseTDCHits, and QwDriftChamber::fWireHits.

   {
     if(fUseTDCHits) grandHitContainer.Append(fTDCHits);
     else            grandHitContainer.Append(fWireHits);
   };

Here is the call graph for this function:

Int_t QwDriftChamberVDC::LoadChannelMap ( TString mapfile )

virtual

Mandatory map file definition.

Implements VQwSubsystem.

Definition at line 641 of file QwDriftChamberVDC.cc.

References AddChannelDefinition(), BuildWireDataStructure(), QwDriftChamber::fCurrentBankIndex, QwDriftChamber::fCurrentSlot, fDelayLineArray, fDelayLinePtrs, VQwSubsystem::fDetectorMaps, QwDriftChamber::fF1RefContainer, fR3Octant, QwParameterFile::GetUInt(), QwOptions::GetValue(), gQwOptions, QwDriftChamber::kCodaMasterPlaneNumber, kDirectionNull, kDirectionU, kDirectionV, QwDriftChamber::kReferenceChannelPlaneNumber, QwDriftChamber::LinkReferenceChannel(), LoadTimeWireOffset(), LoadTtoDParameters(), OK, QwDriftChamber::RegisterROCNumber(), QwDriftChamber::RegisterSlotNumber(), and QwDriftChamber::RegisterSubbank().

 {
   //some type(like string,Int_t)need to be changed to root type
   LoadTimeWireOffset ( "R3_timeoffset.map" );
   LoadTtoDParameters ( "R3_TtoDTable.map" );
   TString varname,varvalue;
   UInt_t value   = 0;
   UInt_t channum = 0;            //store temporary channel number
   UInt_t lnnum   = 0;           //store line number
   Int_t  bpnum   = 0;           // temp backplane
   UInt_t plnum,firstwire,LR;         //store temp package,plane,firstwire and left or right information
 
   TString name = "";
 
   plnum = firstwire = LR = 0;
 
   TString pknum,dir;
   Bool_t IsFirstChannel = kTRUE;
 
   std::vector<Double_t> tmpWindows;
   QwParameterFile mapstr ( mapfile.Data() );
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
   EQwDirectionID   direction = kDirectionNull;
 
   std::string string_a;
   std::pair<Double_t,Double_t> pair_a;
 
   Int_t extra_reference_channel_number_for_R3 = 97;
 
   while ( mapstr.ReadNextLine() )
     {
       mapstr.TrimComment ( '!' );
       mapstr.TrimWhitespace();
       if ( mapstr.LineIsEmpty() ) continue;
 
       if ( mapstr.HasVariablePair ( "=",varname,varvalue ) )    //to judge whether we find a new slot
   {
     varname.ToLower();
     value = QwParameterFile::GetUInt ( varvalue );
     if ( value ==0 )
       {
         value = atol ( varvalue.Data() );
       }
     if ( varname == "roc" )
       {
         RegisterROCNumber ( value , 0 );
       }
     else if ( varname=="bank" )
       {
         RegisterSubbank ( value );
       }
     else if ( varname == "slot" )
       {
         RegisterSlotNumber ( value );
       }
     continue;        //go to the next line
   }
 
       channum = mapstr.GetTypedNextToken<Int_t>();
       bpnum   = mapstr.GetTypedNextToken<Int_t>();
       name    = mapstr.GetTypedNextToken<TString>();
 
    
 
       if ( bpnum == kReferenceChannelPlaneNumber || bpnum == kCodaMasterPlaneNumber  || bpnum == extra_reference_channel_number_for_R3) {
 
   fF1RefContainer -> AddF1TDCReferenceSignal(new F1TDCReferenceSignal(fCurrentBankIndex, fCurrentSlot, channum, name));
   //
   // old map uses 0, so simply set lnnum = 0
   //
   if ( bpnum != extra_reference_channel_number_for_R3 ) {
     lnnum = 0;
     LinkReferenceChannel ( channum, bpnum, lnnum);
   }
   continue;
       }
       
       //
       // read third one as string, so it must be converted to integer
       // out of Reference channels
       //
 
       lnnum = name.Atoi();
 
       if ( channum ==0 and bpnum ==0 )
   {
     if ( IsFirstChannel == kTRUE ) IsFirstChannel = kFALSE;
     else                         continue;
   }
 
 
       LR       = mapstr.GetTypedNextToken<Int_t>();
 
       fDelayLinePtrs.at ( value ).push_back ( QwDelayLineID ( bpnum,lnnum,LR ) );    //the slot and channel number must be in order
 
       pknum    = mapstr.GetTypedNextToken<TString>();
       plnum    = mapstr.GetTypedNextToken<Int_t>();
       dir      = mapstr.GetTypedNextToken<TString>();
       firstwire= mapstr.GetTypedNextToken<Int_t>();
 
       Int_t octant = 0;
       EQwDetectorPackage package = kPackageNull;
       fR3Octant = gQwOptions.GetValue<Int_t>("R3-octant");
       if ( pknum=="u" )
         {
           package = kPackage1;
           octant  = (fR3Octant + 4) % 8;
         }
       else if ( pknum=="v" )
         {
           package = kPackage2;
           octant  = fR3Octant;
         }
 
       BuildWireDataStructure ( channum,package,octant,plnum,firstwire );
 
       if ( fDelayLineArray.at ( bpnum ).at ( lnnum ).Fill == kFALSE )   //if this delay line has not been Filled in the data
         {
     string_a = mapstr.GetTypedNextToken<std::string>();
     while ( string_a.size() !=0 )
       {
         tmpWindows.push_back ( atof ( string_a.c_str() ) );
               string_a = mapstr.GetTypedNextToken<std::string>();
             }
 
           fDelayLineArray.at ( bpnum ).at ( lnnum ).fOctant  = octant;
           fDelayLineArray.at ( bpnum ).at ( lnnum ).fPackage = package;
           fDelayLineArray.at ( bpnum ).at ( lnnum ).fPlane   = plnum;
 
     if      ( dir == "u" ) direction = kDirectionU;
     else if ( dir == "v" ) direction = kDirectionV;
 
     fDelayLineArray.at ( bpnum ).at ( lnnum ).fDirection = direction;
     fDelayLineArray.at ( bpnum ).at ( lnnum ).fFirstWire = firstwire;
 
     for ( size_t i=0;i<tmpWindows.size() /2;i++ )
       {
         pair_a.first  = tmpWindows.at ( 2*i );
         pair_a.second = tmpWindows.at ( 2*i+1 );
         fDelayLineArray.at ( bpnum ).at ( lnnum ).Windows.push_back ( pair_a );
       }
 
     //            std::cout << "DelayLine: back plane: " << bpnum
     //          << "line number " << lnnum
     //          << " Windows.size: "  << fDelayLineArray.at ( bpnum ).at ( lnnum ).Windows.size()
     //          << std::endl;
     fDelayLineArray.at ( bpnum ).at ( lnnum ).Fill = kTRUE;
     tmpWindows.clear();
   }
     }
 
   AddChannelDefinition();
 
   mapstr.Close(); // Close the file (ifstream)
   return OK;
 }

Here is the call graph for this function:

Int_t QwDriftChamberVDC::LoadTimeWireOffset ( TString t0_map )

protectedvirtual

Implements QwDriftChamber.

Definition at line 1162 of file QwDriftChamberVDC.cc.

References VQwSubsystem::fDetectorMaps, fTimeWireOffsets, and OK.

Referenced by LoadChannelMap().

 {
   //std::cout << "beginning to load t0 file... " << std::endl;
   //
   QwParameterFile mapstr ( t0_map.Data() );
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
   TString varname,varvalue;
   Int_t plane=0,wire=0;
   Double_t t0 = 0.0;
   EQwDetectorPackage package = kPackageNull;
 
   while ( mapstr.ReadNextLine() )
     {
       mapstr.TrimComment ( '!' );
       mapstr.TrimWhitespace();
       if ( mapstr.LineIsEmpty() ) continue;
       if ( mapstr.HasVariablePair ( "=",varname,varvalue ) )
   {
     varname.ToLower();
     if ( varname=="package" )
       {
         package = ( EQwDetectorPackage ) atoi ( varvalue.Data() );
         if ( package> ( Int_t ) fTimeWireOffsets.size() ) fTimeWireOffsets.resize ( package );
       }
     else if ( varname=="plane" )
       {
         //std::cout << "package: "  <<  fTimeWireOffsets.at(package-1).size()<< std::endl;
         plane = atoi ( varvalue.Data() );
         if ( plane> ( Int_t ) fTimeWireOffsets.at ( package-1 ).size() ) fTimeWireOffsets.at ( package-1 ).resize ( plane );
         //std::cout << "plane: "  <<  fTimeWireOffsets.at(package-1).size()<< std::endl;
 
         // To Siyuan, * : can package be obtained before plane in while loop? if plane is before package
         //                we have at(-1), thus, if condition is always "false", I guess.
         //            * : if, else if then can we expect something more?
         // from Han
       }
     continue;
   }
 
       wire = mapstr.GetTypedNextToken<Int_t>();
       t0   = mapstr.GetTypedNextToken<Int_t>();
 
       if ( wire > ( Int_t ) fTimeWireOffsets.at ( package-1 ).at ( plane-1 ).size() )
   {
     fTimeWireOffsets.at ( package-1 ).at ( plane-1 ).resize ( wire );
 
     fTimeWireOffsets.at ( package-1 ).at ( plane-1 ).at ( wire-1 ) = t0;
   }
 
     }
   //
   mapstr.Close(); // Close the file (ifstream)
   return OK;
 }

Here is the caller graph for this function:

void QwDriftChamberVDC::LoadTtoDParameters ( TString ttod_map )

protected

Definition at line 1395 of file QwDriftChamberVDC.cc.

References VQwSubsystem::fDetectorMaps, and fTtoDNumbers.

Referenced by LoadChannelMap().

 {
 
   QwParameterFile mapstr ( ttod_map.Data() );
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
   while ( mapstr.ReadNextLine() )
     {
       mapstr.TrimComment ( '!' );
       mapstr.TrimWhitespace();
       if ( mapstr.LineIsEmpty() ) continue;
 
       Double_t t __attribute__((unused)) // unused but function call still required
                  = mapstr.GetTypedNextToken<Double_t>();
       Double_t d = mapstr.GetTypedNextToken<Double_t>();
       fTtoDNumbers.push_back ( d );
     }
 
   mapstr.Close(); // Close the file (ifstream)
   return;
 }

Here is the caller graph for this function:

void QwDriftChamberVDC::ProcessEvent ( )

virtual

Implements QwDriftChamber.

Definition at line 841 of file QwDriftChamberVDC.cc.

References fDisableWireTimeOffset, QwDriftChamber::FillDriftDistanceToHits(), QwDriftChamber::fWireHits, VQwSubsystem::HasDataLoaded(), invalid(), SubtractReferenceTimes(), SubtractWireTimeOffset(), and UpdateHits().

 {
   if ( not HasDataLoaded() ) return;
 
   SubtractReferenceTimes();  // A.U. unit
   UpdateHits();
   
   if ( fDisableWireTimeOffset == false ){
     
     SubtractWireTimeOffset();
     FillDriftDistanceToHits();
     
     std::vector<QwHit>::iterator it;
     it=std::remove_if(fWireHits.begin(),fWireHits.end(),invalid);
     if(it!=fWireHits.end())
       fWireHits.erase(it,fWireHits.end());    
   }
         
   return;
 }

Here is the call graph for this function:

void QwDriftChamberVDC::ProcessOptions ( QwOptions & options )

virtual

Process the command line options.

Reimplemented from VQwSubsystem.

Definition at line 918 of file QwDriftChamberVDC.cc.

References fDisableWireTimeOffset, fR3Octant, fUseTDCHits, and QwOptions::GetValue().

 {
   fUseTDCHits=options.GetValue<bool> ( "use-tdchit" );
   fDisableWireTimeOffset=options.GetValue<bool> ( "disable-wireoffset" );
   fR3Octant=options.GetValue<int> ( "R3-octant" );
 }

Here is the call graph for this function:

void QwDriftChamberVDC::ReadEvent ( TString & eventfile )

Definition at line 801 of file QwDriftChamberVDC.cc.

References VQwSubsystem::fDetectorMaps, QwDriftChamber::fTDCHits, kDirectionNull, and kRegionID3.

 {
   TString varname,varvalue;
   UInt_t slotnum = 0;
   UInt_t channum = 0;            //store temporary channel number
   //  UInt_t pknum,plnum;         //store temp package,plane,firstwire and left or right information
   UInt_t value = 0;
   Double_t signal = 0.0; // Double_t? unsigned Int_t ? by jhlee
   EQwDetectorPackage package = kPackageNull;
   EQwDirectionID direction = kDirectionNull;
   
   QwParameterFile mapstr ( eventfile.Data() );
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
 
   while ( mapstr.ReadNextLine() )
     {
       mapstr.TrimComment ( '!' );
       mapstr.TrimWhitespace();
       if ( mapstr.LineIsEmpty() ) continue;
       if ( mapstr.HasVariablePair ( "=",varname,varvalue ) )   //to judge whether we find a new crate
   {
     varname.ToLower();
     value = atol ( varvalue.Data() );
     continue;
   }
       
       slotnum = mapstr.GetTypedNextToken<Int_t>();
       channum = mapstr.GetTypedNextToken<Int_t>();
       signal  = mapstr.GetTypedNextToken<Double_t>();
       //std::cout << "signal is: " << signal << endl;
       fTDCHits.push_back ( QwHit ( value,slotnum,channum,0, kRegionID3,package,0,0,direction,0, ( UInt_t ) signal ) );
     }        //only know TDC information and time value
 
   mapstr.Close(); // Close the file (ifstream)
   return;
 }

void QwDriftChamberVDC::SubtractReferenceTimes ( )

virtual

Implements QwDriftChamber.

Definition at line 37 of file QwDriftChamberVDC.cc.

References QwDriftChamber::fF1RefContainer, QwDriftChamber::fF1TDContainer, QwDriftChamber::fTDCHits, F1TDCReferenceContainer::GetReferenceTimeAU(), VQwSubsystem::GetSubsystemName(), kPackage1, kPackage2, 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";
 
   std::vector<QwHit>::iterator end=fTDCHits.end();
   for ( std::vector<QwHit>::iterator hit=fTDCHits.begin(); hit!=end; hit++ ) 
     {
 
       bank_index        =            hit  -> GetSubbankID();
       slot_num          =            hit  -> GetModule();
       raw_time_arb_unit = (Double_t) hit  -> GetRawTime();
       package           =            hit  -> GetPackage();
 
 
       // John L. introduced this package match reference time subtraction for only Region 3,
       // I was told it would be better to extend this concept to Region2, so I did, and I 
       // got the worst result.....  Q2 = +2.735199e-02  its RMS +1.192504e-02 
       // if I don't use this method in Region2, Q2 = +2.292274e-02  and RMS = +7.133603e-03
       // if I don't use this method in Region3, I couldn't get "peak" structure 
       // in time distributions of four VDC chambers. 
       // 
       // Why does this method give us a good result for only Region 3? MUX? 
       // I don't know ..
       // Friday, February 17 14:01:12 EST 2012, jhlee
 
       if (package == kPackage1) { // Up is 1
   reference_name1 = "TrigScintPkg1";
       }
       else if (package == kPackage2) {// Down is 2
   reference_name1 = "TrigScintPkg2";
       }
       else {
   reference_name1 = "MasterTrigger";
       }
       
       ref_time_arb_unit =  fF1RefContainer->GetReferenceTimeAU(bank_index, reference_name1);
     
       //   printf("%f\n", ref_time_arb_unit);
 
       if(ref_time_arb_unit== 0.0 ) {
   ref_time_arb_unit =  fF1RefContainer->GetReferenceTimeAU(bank_index, reference_name2);
       }
       //      printf("%f\n", ref_time_arb_unit);
       // 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()
       << " Package "   << std::setw(2) << package
       << " 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 QwDriftChamberVDC::SubtractWireTimeOffset ( )

protectedvirtual

Implements QwDriftChamber.

Definition at line 1219 of file QwDriftChamberVDC.cc.

References VQwSubsystemTracking::fF1TDCResolutionNS, fTimeWireOffsets, and QwDriftChamber::fWireHits.

Referenced by ProcessEvent().

 {
   Int_t plane = 0;
   Int_t wire  = 0;
   EQwDetectorPackage package = kPackageNull;
 
   Double_t t0_NS = 0.0;
 
   Double_t educated_guess_t0_correction_NS = 91; 
   // Double_t educated_guess_t0_correction_AU = educated_guess_t0_correction_NS/fF1TDCResolutionNS;
 
   std::size_t nhits = fWireHits.size();
 
   for(std::size_t i=0;i<nhits;i++)
     {
       package = fWireHits.at(i).GetPackage();
       plane   = fWireHits.at(i).GetPlane();
       wire    = fWireHits.at(i).GetElement();
       t0_NS   = fTimeWireOffsets.at ( package-1 ).at ( plane-1 ).at ( wire-1 );
       t0_NS  += educated_guess_t0_correction_NS;
 
       // real_time=fWireHits.at(i).GetTime()-t0-91;
       fWireHits.at(i).SubtractTimeNsOffset(t0_NS);                     // time unit is ns Replace fTimeNs
       fWireHits.at(i).SubtractTimeAuOffset(t0_NS/fF1TDCResolutionNS);  // time unit is a.u. Replace fTime
     }
   return;
 }

Here is the caller graph for this function:

void QwDriftChamberVDC::UpdateHits ( )

protected

Definition at line 1249 of file QwDriftChamberVDC.cc.

References QwHit::ApplyTimeCalibration(), QwElectronicsID::fChannel, fDelayLineArray, fDelayLinePtrs, VQwSubsystemTracking::fDetectorInfo, VQwSubsystemTracking::fF1TDCResolutionNS, QwElectronicsID::fModule, QwDriftChamber::fTDCHits, QwDriftChamber::fWireData, QwDriftChamber::fWireHits, QwGeometry::in(), kRegionID3, QwHit::SetAmbiguityID(), VQwTrackingElement::SetDetectorInfo(), QwHit::SetHitNumberR(), and QwHit::SetTime().

Referenced by ProcessEvent().

 {
 
   Double_t real_time_au  = 0.0;
   Double_t tmpTime_au    = 0.0;
   Double_t left_time_au  = 0.0;
   Double_t right_time_au = 0.0;
 
   Int_t    tmpCrate   = 0;
   Int_t    tmpModule  = 0;
   Int_t    tmpChan    = 0;
   Int_t    tmpbp      = 0;
   Int_t    tmpln      = 0;
   Int_t    wire_hit   = 0;
   Int_t    mycount    = 0;
   
   Bool_t  kDir        = kTRUE;
   Bool_t  tmpAM       = kFALSE;
 
   std::vector<Int_t> wire_array;
   wire_array.clear();
 
   Int_t channel_offset = 64;
 
   QwElectronicsID tmpElectronicsID;
 
   // processing the delay line starts....
   std::vector<QwHit>::iterator end=fTDCHits.end();
   for ( std::vector<QwHit>::iterator iter=fTDCHits.begin();iter!=end;iter++ )
     {
       //this for loop will Fill in the tdc hits data Int_to the corresponding delay line
       tmpElectronicsID = iter->GetElectronicsID();
       tmpCrate         = iter->GetSubbankID();
       tmpTime_au       = iter->GetTime();  // a.u.
       tmpModule        = tmpElectronicsID.fModule;
       tmpChan          = tmpElectronicsID.fChannel;
       
       if ( tmpCrate==3 ) tmpChan+=channel_offset; // ROC10
 
       if ( fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fSide == 0 ) {
   fDelayLineArray
     .at( fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fBackPlane )
     .at( fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fLineNumber)
     .LeftHits.push_back ( tmpTime_au );
       }
       else {
   fDelayLineArray
     .at ( fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fBackPlane )
     .at ( fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fLineNumber )
     .RightHits.push_back ( tmpTime_au );
       }
     }
 
   for ( std::vector<QwHit>::iterator iter=fTDCHits.begin();iter!=end;iter++ )
     {
       tmpElectronicsID = iter->GetElectronicsID();
       tmpCrate         = iter->GetSubbankID();
       // if(tmpCrate == 0) break;
       //    tmpTime_au        = iter->GetTime();
       tmpModule        = tmpElectronicsID.fModule;
       tmpChan          = tmpElectronicsID.fChannel;
       // Int_t tmpROC     = tmpElectronicsID.fROC;
       if ( tmpCrate==3 ) tmpChan+=channel_offset; // ROC10
 
       tmpbp  = fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fBackPlane;
       tmpln  = fDelayLinePtrs.at ( tmpModule ).at ( tmpChan ).fLineNumber;
 
       Int_t plane     = fDelayLineArray.at ( tmpbp ).at ( tmpln ).fPlane ;
       Int_t octant    = fDelayLineArray.at ( tmpbp ).at ( tmpln ).fOctant;
       EQwDetectorPackage package = fDelayLineArray.at ( tmpbp ).at ( tmpln ).fPackage;
       EQwDirectionID direction = fDelayLineArray.at ( tmpbp ).at ( tmpln ).fDirection;
 
       // QwMessage << "QwDriftChamberVDC::ProcessEvent() :"
       //          << " plane = "  << plane
       //        << " direction = " << direction
       //        << " package = " << package
       //    << " roc = " << tmpROC
       //    << " slot = " << tmpModule
       //    << " chan = " << tmpChan
 
       //        << QwLog::endl;
 
       if ( fDelayLineArray.at ( tmpbp ).at ( tmpln ).Processed == kFALSE )   //if this delay line has been Processed
   {
     wire_array.clear();
     
     if ( tmpbp==0 || tmpbp ==3 || tmpbp==4 || tmpbp==7 || tmpbp==8 || tmpbp==11 || tmpbp==12 || tmpbp==15 ) {
       kDir = kTRUE;         //true means left-right
     }
     else { 
       kDir = kFALSE;
     }
     
     fDelayLineArray.at ( tmpbp ).at ( tmpln ).ProcessHits ( kDir );
 
     std::size_t Wirecount = fDelayLineArray.at ( tmpbp ).at ( tmpln ).Wire.size();
     //      std::cout << "Wirecount " << Wirecount << std::endl;
     for ( std::size_t i=0;i<Wirecount;i++ )
       {
         std::size_t  Ambiguitycount = fDelayLineArray.at ( tmpbp ).at ( tmpln ).Wire.at ( i ).size();   //if there's a ambiguity, it's 2; if not, this is 1
         if ( Ambiguitycount==1 ) tmpAM = kFALSE;
         else                     tmpAM = kTRUE;
         
         Int_t order_L = fDelayLineArray.at ( tmpbp ).at ( tmpln ).Hitscount.at ( i ).first;
         Int_t order_R = fDelayLineArray.at ( tmpbp ).at ( tmpln ).Hitscount.at ( i ).second;
 
         left_time_au  = fDelayLineArray.at ( tmpbp ).at ( tmpln ).LeftHits.at  ( order_L );
         right_time_au = fDelayLineArray.at ( tmpbp ).at ( tmpln ).RightHits.at ( order_R );
 
         for ( std::size_t j=0;j<Ambiguitycount;++j )
     {
       real_time_au = ( left_time_au+right_time_au ) /2.0; // a.u.
 
       wire_hit     = fDelayLineArray.at ( tmpbp ).at ( tmpln ).Wire.at ( i ).at ( j );
 
       wire_array.push_back ( wire_hit );
 
       mycount = count ( wire_array.begin(), wire_array.end(), wire_hit ) - 1;
 
       Int_t temp_Chan=0;
       if ( tmpCrate==3 ) temp_Chan=tmpChan-channel_offset;
       else               temp_Chan=tmpChan;
 
                   const QwDetectorInfo* local_info = fDetectorInfo.in(package).at(plane-1);
 
                   QwHit NewQwHit ( tmpCrate, tmpModule, temp_Chan, mycount, kRegionID3, package, octant, plane, direction, wire_hit );
 
                   NewQwHit.SetHitNumberR       ( order_R );
                   NewQwHit.SetDetectorInfo     ( local_info );
       NewQwHit.SetAmbiguityID      ( tmpAM, j );
 
       NewQwHit.SetTime             ( real_time_au );
       NewQwHit.ApplyTimeCalibration( fF1TDCResolutionNS );  // Fill fTimeRes and fTimeNs in QwHit
 
       fWireData.at ( 4* ( package-1 ) +plane ).at ( wire_hit-1 ).PushHit ( ( Int_t ) real_time_au );
 
       fWireHits.push_back ( NewQwHit );
     }
             }
         }
     }
 }