#include <QwBeamMod.h>

Inheritance diagram for QwBeamMod:

Collaboration diagram for QwBeamMod:

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

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

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

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

void	AccumulateRunningSum (VQwSubsystem *)
	Update the running sums for devices. More...

void	DeaccumulateRunningSum (VQwSubsystem *value)
	remove one entry from the running sums for devices More...

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

Int_t	LoadEventCuts (TString filename)
	Load the event cuts file. More...

Int_t	LoadGeometry (TString mapfile)

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

Bool_t	ApplySingleEventCuts ()
	Apply the single event cuts. More...

void	IncrementErrorCounters ()
	Increment the error counters. More...

void	PrintErrorCounters () const
	Report the number of events failed due to HW and event cut failures. More...

UInt_t	GetEventcutErrorFlag ()
	Return the error flag to the top level routines related to stability checks and ErrorFlag updates. More...

void	UpdateErrorFlag (const VQwSubsystem *ev_error)
	update the error flag in the subsystem level from the top level routines related to stability checks. This will uniquely update the errorflag at each channel based on the error flag in the corresponding channel in the ev_error subsystem More...

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

Int_t	ProcessEvBuffer (const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)
	TODO: The non-event-type-aware ProcessEvBuffer routine should be replaced with the event-type-aware version. More...

void	ClearEventData ()

void	ProcessEvent ()

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

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

VQwSubsystem &	operator= (VQwSubsystem *value)
	Assignment Note: Must be called at the beginning of all subsystems routine call to operator=(VQwSubsystem *value) by VQwSubsystem::operator=(value) More...

VQwSubsystem &	operator+= (VQwSubsystem *value)

VQwSubsystem &	operator-= (VQwSubsystem *value)

void	Sum (VQwSubsystem value1, VQwSubsystem value2)

void	Difference (VQwSubsystem value1, VQwSubsystem value2)

void	Ratio (VQwSubsystem numer, VQwSubsystem denom)

void	Scale (Double_t factor)

void	CalculateRunningAverage ()
	Calculate the average for all good events. More...

void	PrintModChannelID ()

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

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

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

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

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

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

void	FillDB_MPS (QwParityDB *db, TString datatype)
	Fill the database with MPS-based variables Note that most subsystems don't need to do this. More...

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

void	FillErrDB (QwParityDB *db, TString datatype)

void	WritePromptSummary (QwPromptSummary *ps, TString type)

Bool_t	Compare (VQwSubsystem *source)

void	Print ()

void	AtEndOfEventLoop ()
	Perform actions at the end of the event loop. More...

void	AnalyzeOpticsPlots ()

void	ResizeOpticsDataContainers (Int_t)

void	ClearVectors ()

Double_t	GetAmplitudeSign (Double_t, Double_t, Double_t, Double_t)

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

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

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

virtual UInt_t	UpdateErrorFlag ()
	Uses the error flags of contained data elements to update Returns the error flag to the top level routines related to stability checks and ErrorFlag updates. More...

virtual void	Blind (const QwBlinder *blinder)
	Blind the asymmetry of this subsystem. More...

virtual void	Blind (const QwBlinder blinder, const VQwSubsystemParity subsys)
	Blind the difference of this subsystem. More...

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

virtual Bool_t	CheckForEndOfBurst () const

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

	VQwSubsystem (const VQwSubsystem &orig)
	Copy constructor by object. More...

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

TString	GetSubsystemName () const

Bool_t	HasDataLoaded () const

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

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

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

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

virtual Bool_t	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	LoadGeometryDefinition (TString mapfile)
	Optional geometry definition. More...

virtual Int_t	LoadCrosstalkDefinition (TString mapfile)
	Optional crosstalk definition. More...

void	SetEventTypeMask (const UInt_t mask)
	Set event type mask. More...

UInt_t	GetEventTypeMask () const
	Get event type mask. More...

virtual Int_t	ProcessEvBuffer (const UInt_t event_type, const UInt_t roc_id, const UInt_t bank_id, UInt_t *buffer, UInt_t num_words)

virtual void	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 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, QwBeamMod >
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...

Data Fields
std::vector< TString >	fgModTypeNames

Protected Member Functions
Int_t	GetDetectorTypeID (TString name)

Int_t	GetDetectorIndex (Int_t TypeID, TString name)

Protected Member Functions inherited from VQwSubsystem
void	UpdatePublishedValue (const TString &name, VQwHardwareChannel *data_channel)

void	ClearAllBankRegistrations ()
	Clear all registration of ROC and Bank IDs for this subsystem. More...

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

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

Int_t	GetSubbankIndex () const

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

void	SetDataLoaded (Bool_t flag)

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

Bool_t	Compare (VQwSubsystem *source)

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

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

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

virtual MQwHistograms &	operator= (const MQwHistograms &value)

void	Fill_Pointer (TH1_ptr hist_ptr, Double_t value)

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

Protected Attributes
Int_t	fTreeArrayIndex

std::vector< QwVQWK_Channel >	fModChannel

std::vector< QwModChannelID >	fModChannelID

std::vector< QwWord >	fWord

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

Protected Attributes inherited from VQwSubsystem
std::map< TString, VQwHardwareChannel * >	fPublishedInternalValues
	Map of published internal values. More...

std::vector< std::vector < TString > >	fPublishList
	List of parameters to be published (loaded at the channel map) More...

TString	fSystemName
	Name of this subsystem. More...

UInt_t	fEventTypeMask
	Mask of event types. More...

Bool_t	fIsDataLoaded
	Has this subsystem gotten data to be processed? More...

std::vector< TString >	fDetectorMapsNames

std::map< TString, TString >	fDetectorMaps

Int_t	fCurrentROC_ID
	ROC ID that is currently being processed. More...

Int_t	fCurrentBank_ID
	Bank ID that is currently being processed. More...

std::vector< UInt_t >	fROC_IDs
	Vector of ROC IDs associated with this subsystem. More...

std::vector< std::vector < UInt_t > >	fBank_IDs
	Vector of Bank IDs per ROC ID associated with this subsystem. More...

std::vector< QwSubsystemArray * >	fArrays
	Vector of pointers to subsystem arrays that contain this subsystem. More...

Protected Attributes inherited from MQwHistograms
std::vector< TH1_ptr >	fHistograms
	Histograms associated with this data element. More...

Private Member Functions
	QwBeamMod ()

Private Attributes
UInt_t	fFFB_Index

UInt_t	fFFB_holdoff

UInt_t	fFFB_holdoff_Counter

UInt_t	fFFB_ErrorFlag

Bool_t	fFFB_Flag

std::vector< TString >	fMonitorNames

std::vector< QwVQWK_Channel >	fMonitors

std::vector< Double_t >	fBPMPosition

std::vector< std::vector < Double_t > >	fAmplitude

std::vector< std::vector < Double_t > >	fOffset

std::vector< std::vector < Double_t > >	fPhase

std::vector< std::vector < Double_t > >	fAmplitudeError

std::vector< std::vector < Double_t > >	fOffsetError

std::vector< std::vector < Double_t > >	fPhaseError

std::vector< std::vector < Double_t > >	fChisquare

std::vector< std::vector< Int_t > >	fNFitPoints

Int_t	fRampChannelIndex

Int_t	fPatternWordIndex

Static Private Attributes
static const Bool_t	bDEBUG =kFALSE

static const Int_t	fNumberPatterns = 5

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

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

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

Detailed Description

Definition at line 41 of file QwBeamMod.h.

Constructor & Destructor Documentation

QwBeamMod::QwBeamMod ( )

private

QwBeamMod::QwBeamMod ( const TString & name )

inline

Constructor with name.

Definition at line 49 of file QwBeamMod.h.

References fFFB_Flag, fFFB_holdoff_Counter, fgModTypeNames, fPatternWordIndex, and fRampChannelIndex.

   : VQwSubsystem(name),VQwSubsystemParity(name)
     {
       // these declaration need to be coherent with the enum vector EBeamInstrumentType
       fgModTypeNames.push_back(""); // Need to define these wrt to our detector types.
       fgModTypeNames.push_back("");
       fgModTypeNames.push_back("");
       fgModTypeNames.push_back("");
       for(size_t i=0;i<fgModTypeNames.size();i++)
         fgModTypeNames[i].ToLower();
       fFFB_holdoff_Counter=0;
       fFFB_Flag=kTRUE;
       fRampChannelIndex = -1;
       fPatternWordIndex = -1;
     };

QwBeamMod::QwBeamMod ( const QwBeamMod & source )

inline

Copy constructor.

Definition at line 65 of file QwBeamMod.h.

   : VQwSubsystem(source),VQwSubsystemParity(source),
     fModChannel(source.fModChannel),fWord(source.fWord),
     fMonitorNames(source.fMonitorNames),fMonitors(source.fMonitors),
     fBPMPosition(source.fBPMPosition)
   { }

virtual QwBeamMod::~QwBeamMod ( )

inlinevirtual

Virtual destructor.

Definition at line 72 of file QwBeamMod.h.

72 {};

Member Function Documentation

void QwBeamMod::AccumulateRunningSum ( VQwSubsystem * value )

virtual

Update the running sums for devices.

Implements VQwSubsystemParity.

Definition at line 662 of file QwBeamMod.cc.

662 { }

void QwBeamMod::AnalyzeOpticsPlots ( )

Definition at line 827 of file QwBeamMod.cc.

References QwLog::endl(), fAmplitude, fAmplitudeError, fChisquare, MQwHistograms::fHistograms, fMonitorNames, fMonitors, fNFitPoints, fOffset, fOffsetError, fPhase, fPhaseError, GetAmplitudeSign(), QwDebug, and ResizeOpticsDataContainers().

Referenced by AtEndOfEventLoop().

 {
   //   How to get the run info:
   //   UInt_t runnum = this->GetParent()->GetCodaRunNumber();
   //   UInt_t segnum = this->GetParent()->GetCodaSegmentNumber();
 
   TF1 *sine = new TF1("sine", "[0] + [1]*sin(TMath::DegToRad()*x + [2])", 5, 350);
 
   TCanvas *canvas = new TCanvas("canvas", "canvas", 5);
 
   //Double_t mean; // unused
   Double_t amplitude;
   Double_t phase;
 
   ResizeOpticsDataContainers(fMonitorNames.size());
 
   canvas->cd();
   for(size_t bpm = 0; bpm < fMonitors.size(); bpm++){
 
     for(size_t pattern = 0; pattern < 5; pattern++){
       //  Only do the fits if there are more than three entries.
       if (fHistograms[5*bpm + pattern]->GetEntries()>3){
   sine->SetParameters(fHistograms[5*bpm + pattern]->GetMean(), 0.10, 0);
   sine->SetLineColor(2);
   sine->SetParLimits(2, 0, TMath::Pi()*2 );
   fHistograms[5*bpm + pattern]->Fit("sine","R B");
   
   //mean = sine->GetParameter(0); // unused
   amplitude = TMath::Abs(sine->GetParameter(1));
   phase = sine->GetParameter(2) * TMath::RadToDeg();
   
   amplitude *= GetAmplitudeSign(sine->Derivative(10), 
               sine->Derivative2(10), 
               sine->Derivative3(10), 
               sine->GetParameter(0));
 
   if(phase >= 180){
     phase -= 180;
 //    amplitude = -amplitude;
   } else if (phase < 0){
     phase += 180;
 //    amplitude = -amplitude;
   } 
   fOffset[bpm][pattern] = sine->GetParameter(0);
   fAmplitude[bpm][pattern] = amplitude;
   fPhase[bpm][pattern] = phase;
   fOffsetError[bpm][pattern] = sine->GetParError(0);
   fAmplitudeError[bpm][pattern] = sine->GetParError(1);
   fPhaseError[bpm][pattern] = sine->GetParError(2);
   fChisquare[bpm][pattern] = sine->GetChisquare();
   fNFitPoints[bpm][pattern] = sine->GetNumberFitPoints();
       } else {
   QwDebug << "QwBeamMod can't fit [" << bpm << "][" << pattern
     << "] because there are only " 
     << fHistograms[5*bpm + pattern]->GetEntries()
     << " entries, and we need at least 3."
     << QwLog::endl;
   //  No events in this histogram.  Zero-out the fit results
   fOffset[bpm][pattern]         = 0.0;
   fAmplitude[bpm][pattern]      = 0.0;
   fPhase[bpm][pattern]          = 0.0;
   fOffsetError[bpm][pattern]    = 0.0;
   fAmplitudeError[bpm][pattern] = 0.0;
   fPhaseError[bpm][pattern]     = 0.0;
   fChisquare[bpm][pattern]      = 0.0;
   fNFitPoints[bpm][pattern]     = 0;
       }
     }
   }
   delete canvas;
   delete sine;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Bool_t QwBeamMod::ApplySingleEventCuts ( )

virtual

Apply the single event cuts.

Implements VQwSubsystemParity.

Definition at line 414 of file QwBeamMod.cc.

References bDEBUG, fFFB_ErrorFlag, fFFB_Flag, fFFB_holdoff, fFFB_holdoff_Counter, fFFB_Index, fModChannel, fWord, kBModErrorFlag, kBModFFBErrorFlag, kEventCutMode3, and kGlobalCut.

                                       {
 
   Bool_t test_Mod=kTRUE;
   Bool_t test_BCM1=kTRUE;
 
 
   for(size_t i=0;i<fModChannel.size();i++){
     //std::cout<<"  BCM ["<<i<<"] "<<std::endl;
     test_BCM1 = fModChannel[i].ApplySingleEventCuts();
     test_Mod &= test_BCM1;
     if(!test_BCM1 && bDEBUG) std::cout<<"******* QwBeamMod::SingleEventCuts()->BCM[ "<<i<<" , "<<fModChannel[i].GetElementName()<<" ] ******\n";
   }
   
   if (fWord[fFFB_Index].fValue==8 && fFFB_Flag && fFFB_holdoff_Counter==0){
     fFFB_holdoff_Counter=fFFB_holdoff;
     fFFB_Flag=kFALSE;
     
   } 
   if (!fFFB_Flag &&  (fFFB_holdoff_Counter>0 && fFFB_holdoff_Counter<=fFFB_holdoff) ){
     fFFB_ErrorFlag = (kGlobalCut+kBModErrorFlag+kBModFFBErrorFlag+kEventCutMode3);
     fFFB_holdoff_Counter--;
     if(fFFB_holdoff_Counter==0)
       fFFB_Flag=kTRUE;
   }else
     fFFB_ErrorFlag=0;
 
 
 
   return test_Mod;
 
 }

void QwBeamMod::AtEndOfEventLoop ( )

virtual

Perform actions at the end of the event loop.

Reimplemented from VQwSubsystem.

Definition at line 784 of file QwBeamMod.cc.

References AnalyzeOpticsPlots().

 {
 
   AnalyzeOpticsPlots();
 
 }

Here is the call graph for this function:

void QwBeamMod::CalculateRunningAverage ( )

virtual

Calculate the average for all good events.

Implements VQwSubsystemParity.

Definition at line 660 of file QwBeamMod.cc.

660 { }

void QwBeamMod::ClearEventData ( )

virtual

Implements VQwSubsystem.

Definition at line 521 of file QwBeamMod.cc.

References fModChannel, and fWord.

 {
    for(size_t i=0;i<fModChannel.size();i++)
     fModChannel[i].ClearEventData();
 
    for (size_t i=0;i<fWord.size();i++)
     fWord[i].ClearEventData();
 }

void QwBeamMod::ClearVectors ( )

Definition at line 791 of file QwBeamMod.cc.

References fAmplitude, fAmplitudeError, fChisquare, fNFitPoints, fOffset, fOffsetError, fPhase, and fPhaseError.

 {
     fOffset.clear();
     fAmplitude.clear();
     fPhase.clear();
     fOffsetError.clear();
     fAmplitudeError.clear();
     fPhaseError.clear();
     fChisquare.clear();
     fNFitPoints.clear();
     
 }

Bool_t QwBeamMod::Compare ( VQwSubsystem * source )

Definition at line 664 of file QwBeamMod.cc.

References fModChannel.

Referenced by Difference(), operator+=(), operator-=(), operator=(), Ratio(), and Sum().

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

Here is the caller graph for this function:

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

inlinevirtual

Construct the branch and tree vector.

Implements VQwSubsystem.

Definition at line 129 of file QwBeamMod.h.

129 { };

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

inlinevirtual

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

Implements VQwSubsystem.

Definition at line 130 of file QwBeamMod.h.

130 { };

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

virtual

Construct the branch and tree vector.

Implements VQwSubsystem.

Definition at line 924 of file QwBeamMod.cc.

References fModChannel, fTreeArrayIndex, and fWord.

 {
   TString basename;
   
   for(size_t i = 0; i < fModChannel.size(); i++)
     fModChannel[i].ConstructBranchAndVector(tree, prefix, values);
 //   for (size_t i=0;i<fWord.size();i++)
 //     fWord[i].ConstructBranchAndVector(tree, prefix, values);
   fTreeArrayIndex  = values.size();
   for (size_t i=0; i<fWord.size(); i++)
   {
 //    basename = fWord[i].fWordName;
     basename = prefix+fWord[i].fWordName;
     values.push_back(0.0);
     tree->Branch(basename, &(values.back()), basename+"/D");
   }
 }

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

virtual

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

Implements VQwSubsystem.

Definition at line 689 of file QwBeamMod.cc.

References QwHistogramHelper::Construct1DProf(), MQwHistograms::fHistograms, fModChannel, fMonitorNames, and gQwHists.

 {
   // Go to the specified folder
   if (folder != NULL) folder->cd();
 
   // No histogram creation for asym, yield, diff, etc
   if (prefix != "") return;
 
   TString basename;
   for (size_t bpm = 0; bpm < fMonitorNames.size(); bpm++) {
     // Find histogram with correct name
     basename = TString("bmod_") + prefix + fMonitorNames[bpm];
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_X"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_Y"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_E"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_XP"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_YP"));
   }
 
   // Beam modulation correlations
   for (size_t chan = 0; chan < fModChannel.size(); chan++) {
     // Find histogram with correct name
     basename = TString("bmod_") + prefix + fModChannel[chan].GetElementName();
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_X"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_Y"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_E"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_XP"));
     fHistograms.push_back(gQwHists.Construct1DProf(basename + "_YP"));
   }
 
   // Beam modulation channels
   //for (size_t i = 0; i < fModChannel.size(); i++)
   //  fModChannel[i].ConstructHistograms(folder,prefix);
 }

Here is the call graph for this function:

void QwBeamMod::DeaccumulateRunningSum ( VQwSubsystem * value )

inlinevirtual

remove one entry from the running sums for devices

Implements VQwSubsystemParity.

Definition at line 81 of file QwBeamMod.h.

81 {

82 };

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

virtual

Implements VQwSubsystemParity.

Definition at line 627 of file QwBeamMod.cc.

References Compare().

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

Here is the call graph for this function:

void QwBeamMod::ExchangeProcessedData ( )

virtual

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.

Reimplemented from VQwSubsystem.

Definition at line 494 of file QwBeamMod.cc.

References QwLog::endl(), fMonitorNames, fMonitors, QwError, and VQwSubsystem::RequestExternalValue().

 {
   // Make sure sizes are equal
   if (fMonitorNames.size() != fMonitors.size())
     QwError << "QwBeamMod: Sizes of fBPMnames and fBPMs do not match!" << QwLog::endl;
   // Loop over BPMs
   for (size_t bpm = 0; bpm < fMonitorNames.size(); bpm++) {
     // Get references to external values
     if (! RequestExternalValue(fMonitorNames[bpm],&fMonitors[bpm]))
       QwError << "QwBeamMod: RequestExternalValue for " << fMonitorNames[bpm]
               << " failed!" << QwLog::endl;
   }
 }

Here is the call graph for this function:

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

virtual

Fill the database.

Reimplemented from VQwSubsystemParity.

Definition at line 1053 of file QwBeamMod.cc.

 {
   return;
 }

void QwBeamMod::FillDB_MPS	(	QwParityDB *	db,
		TString	type
	)

virtual

Fill the database with MPS-based variables Note that most subsystems don't need to do this.

Reimplemented from VQwSubsystemParity.

Definition at line 995 of file QwBeamMod.cc.

References QwDatabase::Connect(), QwDatabase::Disconnect(), QwLog::endl(), fAmplitude, fAmplitudeError, fChisquare, fMonitorNames, fMonitors, fNFitPoints, fOffset, fOffsetError, fPhase, fPhaseError, QwParityDB::GetAnalysisID(), QwParityDB::GetMonitorID(), Qw::kBoldRed, Qw::kGreen, Qw::kNormal, QwDatabase::Query(), and QwMessage.

 {
   Bool_t local_print_flag = false;
 
   if(local_print_flag) {
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
     QwMessage << "            QwBeamMod::FillDB_MPS                       " << QwLog::endl;
     QwMessage << " --------------------------------------------------------------- " << QwLog::endl;
   }
 
   std::vector<QwParitySSQLS::beam_optics> entrylist;
 
   QwParitySSQLS::beam_optics row;
 
   UInt_t analysis_id = db->GetAnalysisID();
 
   for(size_t bpm = 0; bpm < fMonitors.size(); bpm++){
     for(size_t pattern = 0; pattern < 5; pattern++){
       //  Explicitly zero the beam optics ID to ensure a non-sensical default
       //  is not picked up.
       row.beam_optics_id = 0;
       row.analysis_id = analysis_id;
       row.monitor_id = db->GetMonitorID(fMonitorNames[bpm].Data());
       row.modulation_type_id = pattern;
       row.n = fNFitPoints[bpm][pattern];
       row.offset = fOffset[bpm][pattern];
       row.amplitude = fAmplitude[bpm][pattern];
       row.phase = fPhase[bpm][pattern];
       row.o_error = fOffsetError[bpm][pattern];
       row.a_error = fAmplitudeError[bpm][pattern];
       row.p_error = fPhaseError[bpm][pattern];
       row.gof_para = fChisquare[bpm][pattern];
 
       entrylist.push_back(row);
     }
   }
 
   if(local_print_flag) {
     QwMessage << QwColor(Qw::kGreen) << "Entrylist Size : "
         << QwColor(Qw::kBoldRed) << entrylist.size()
               << QwColor(Qw::kNormal) << QwLog::endl;
   }
 
   db->Connect();
   // Check the entrylist size, if it isn't zero, start to query..
   if( entrylist.size() ) {
     mysqlpp::Query query= db->Query();
     query.insert(entrylist.begin(), entrylist.end());
     query.execute();
   }
   else {
     QwMessage << "QwBeamMod::FillDB_MPS :: Nothing to insert in database." << QwLog::endl;
   }
   db->Disconnect();
 
   return;
 }

Here is the call graph for this function:

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

virtual

Reimplemented from VQwSubsystemParity.

Definition at line 1058 of file QwBeamMod.cc.

 {
   return;
 }

void QwBeamMod::FillHistograms ( )

virtual

Fill the histograms for this subsystem.

Implements VQwSubsystem.

Definition at line 724 of file QwBeamMod.cc.

References MQwHistograms::fHistograms, fModChannel, fMonitors, fPatternWordIndex, fRampChannelIndex, fWord, and VQwSubsystem::HasDataLoaded().

Referenced by ProcessEvent_2().

 {
   // No data loaded
   if (! HasDataLoaded()) return;
 
   // No histograms defined
   if (fHistograms.size() == 0) return;
 
   // Beam modulation channels
   //for (size_t i = 0; i < fModChannel.size(); i++)
   //  fModChannel[i].FillHistograms();
 
   // Do we have a ramp channel?
   if (fRampChannelIndex < 0 || fRampChannelIndex > Int_t(fModChannel.size())) return;
 
   // Do we have a pattern word?
   if (fPatternWordIndex < 0 || fPatternWordIndex > Int_t(fWord.size())) return;
 
   // Determine the ramp/phase
   Double_t ramp = fModChannel[fRampChannelIndex].GetValue();
   if (ramp < 0) return;
 
   // Determine the pattern number -- the pattern number for single coil is
   // between [0, 4] so we need to check for this.
   Int_t pattern = -1;
   
   if(fWord[fPatternWordIndex].fValue < 11 && fWord[fPatternWordIndex].fValue > 0)
     pattern = fWord[fPatternWordIndex].fValue;
   else
     pattern = fWord[fPatternWordIndex].fValue - 11;
   
   if (pattern < 0 || pattern > 4) return;
 
   // Fill histograms for all BPMs and each of the modulation patterns
   //
   // Due to the the way the ADC averages the ramp signal we want to filter
   // out events at the edged of the signal.
   //
   // Seperated the ramp cut here because it is ridiculously long... 
   //
 
   Double_t ramp_block_41 = fModChannel[fRampChannelIndex].GetValue(4) + fModChannel[fRampChannelIndex].GetValue(1); 
   Double_t ramp_block_32 = fModChannel[fRampChannelIndex].GetValue(3) + fModChannel[fRampChannelIndex].GetValue(2);
   Double_t ramp_block    = ramp_block_41 - ramp_block_32;  
 
   //  Require the difference between the block1/4 average and 
   //  the block2/3 average to be within 5 degrees.
   //  In run 9636, the distribution of ramp_block has a sigma of 0.5 degrees.
   if( ramp_block < 5.0 && ramp_block > -5.0 ){
     for (size_t bpm = 0; bpm < fMonitors.size(); bpm++){
       fHistograms[5 * bpm + pattern]->Fill(ramp,fMonitors[bpm].GetValue());
     }
   
     // Beam modulation correlations
     for (size_t chan = 0; chan < fModChannel.size(); chan++){
       fHistograms[5 * (fMonitors.size() + chan) + pattern]->Fill(ramp,fModChannel[chan].GetValue());
     }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

virtual

Fill the tree vector.

Implements VQwSubsystem.

Definition at line 942 of file QwBeamMod.cc.

References fModChannel, fTreeArrayIndex, and fWord.

 {
   size_t index = fTreeArrayIndex;
 
   for (size_t i = 0; i < fModChannel.size(); i++)
     fModChannel[i].FillTreeVector(values);
   for (size_t i = 0; i < fWord.size(); i++)
     values[index++] = fWord[i].fValue;
 }

Double_t QwBeamMod::GetAmplitudeSign	(	Double_t	d1,
		Double_t	d2,
		Double_t	d3,
		Double_t	fmean
	)

Definition at line 900 of file QwBeamMod.cc.

Referenced by AnalyzeOpticsPlots().

 {
 
   Double_t sign = 0.0;
 
   if(d1 > 0.0 && d2 < 0.0)          sign =  1.0;
   else if(d1 == 0.0 && fmean > 0.0) sign =  1.0;
   else if(d1 < 0.0 && d2 < 0.0)     sign =  1.0;
 
   else if(d1 < 0.0 && d3 < 0)       sign =  1.0;
 
   else if(d1 < 0.0 && d2 > 0.0)     sign = -1.0;
   else if(d1 == 0.0 && fmean < 0.0) sign = -1.0;
   else if(d1 > 0.0 && d2 > 0.0)     sign = -1.0;
 
   else if(d1 < 0.0 && d3 > 0)       sign = -1.0;
 
   else
     sign = 1.0;
 
   return(sign);
 
 }

Here is the caller graph for this function:

Int_t QwBeamMod::GetDetectorIndex	(	Int_t	TypeID,
		TString	name
	)

protected

Definition at line 541 of file QwBeamMod.cc.

References fModChannelID.

Referenced by LoadEventCuts().

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

Here is the caller graph for this function:

Int_t QwBeamMod::GetDetectorTypeID ( TString name )

protected

Definition at line 531 of file QwBeamMod.cc.

References fgModTypeNames.

Referenced by LoadEventCuts().

 {
   Int_t result=-1;
   for(size_t i=0;i<fgModTypeNames.size();i++)
     if(name==fgModTypeNames[i])
       {result=i;i=fgModTypeNames.size()+1;}
   return result;
 }

Here is the caller graph for this function:

UInt_t QwBeamMod::GetEventcutErrorFlag ( )

virtual

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

Implements VQwSubsystemParity.

Definition at line 472 of file QwBeamMod.cc.

References fFFB_ErrorFlag, and fModChannel.

                                       {//return the error flag
   UInt_t ErrorFlag;
   ErrorFlag=0;
   for(size_t i=0;i<fModChannel.size();i++){
     ErrorFlag |= fModChannel[i].GetEventcutErrorFlag();
   }
   ErrorFlag |= fFFB_ErrorFlag;
   return ErrorFlag;
 
 }

void QwBeamMod::IncrementErrorCounters ( )

virtual

Increment the error counters.

Implements VQwSubsystemParity.

Definition at line 446 of file QwBeamMod.cc.

References fModChannel.

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

Int_t QwBeamMod::LoadChannelMap ( TString mapfile )

virtual

Mandatory map file definition.

Implements VQwSubsystem.

Definition at line 36 of file QwBeamMod.cc.

References QwLog::endl(), VQwSubsystem::fDetectorMaps, fFFB_Index, QwModChannelID::fIndex, fModChannel, fModChannelID, QwModChannelID::fmodulename, QwWord::fModuleType, fMonitorNames, fMonitors, fPatternWordIndex, fRampChannelIndex, QwWord::fSubbankIndex, fWord, QwWord::fWordInSubbank, QwWord::fWordName, fWordsPerSubbank, VQwSubsystem::GetSubbankIndex(), VQwSubsystem::GetSubsystemName(), QwVQWK_Channel::InitializeChannel(), Print(), QwModChannelID::Print(), QwDebug, QwMessage, QwWarning, QwParameterFile::ReadNextSection(), VQwSubsystem::RegisterROCNumber(), VQwSubsystem::RegisterSubbank(), and ResizeOpticsDataContainers().

 {
   Bool_t ldebug=kFALSE;
 
   Int_t currentrocread=0;
   Int_t currentbankread=0;
   Int_t wordsofar=0;
   Int_t currentsubbankindex=-1;
   Int_t sample_size=0;
 
   TString varname, varvalue;
 
   // Open the file
   QwParameterFile mapstr(mapfile.Data());
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
   mapstr.EnableGreediness();
   mapstr.SetCommentChars("!");
   while (mapstr.ReadNextLine() && mapstr.SkipSection("Monitors"))
   {
     UInt_t value = 0;
     if (mapstr.PopValue("roc",value)) {
       currentrocread = value;
       RegisterROCNumber(value,0);
       fWordsPerSubbank.push_back( std::pair<Int_t, Int_t>(fWord.size(),fWord.size()));
     }
     if (mapstr.PopValue("bank",value)) {
       currentbankread = value;
       RegisterSubbank(value);
       fWordsPerSubbank.push_back( std::pair<Int_t, Int_t>(fWord.size(),fWord.size()));
       //    std::cout<<"bank " <<  currentbankread <<std::endl;
     }
     if (mapstr.PopValue("sample_size",value)) {
       sample_size = value;
     }
 
     if (mapstr.HasVariablePair("=",varname,varvalue))
     {
       //  This is a declaration line.  Decode it.
       varname.ToLower();
       QwWarning << "QwBeamMod::LoadChannelMap: Unrecognized declaration "
           << varname << " = " << varvalue << QwLog::endl;
 
     } else {
       Bool_t lineok=kTRUE;
       //  Break this line into tokens to process it.
       TString modtype = mapstr.GetTypedNextToken<TString>();  // module type
       /* Int_t modnum    = */ mapstr.GetTypedNextToken<Int_t>();  //slot number
       /* Int_t channum   = */ mapstr.GetTypedNextToken<Int_t>();  //channel number
       //  TString dettype   = mapstr.GetTypedNextToken<TString>();  //type-purpose of the detector
       //  dettype.ToLower();
       TString namech  = mapstr.GetTypedNextToken<TString>();  //name of the detector
       namech.ToLower();
       TString keyword = mapstr.GetTypedNextToken<TString>();
       keyword.ToLower();
 
 
       if (currentsubbankindex != GetSubbankIndex(currentrocread,currentbankread))
       {
         currentsubbankindex = GetSubbankIndex(currentrocread,currentbankread);
         wordsofar = 0;
       }
 
 
       QwModChannelID localModChannelID(currentsubbankindex, wordsofar,namech, modtype, this);
 
 
       if(modtype=="VQWK")
       {
         wordsofar+=6;
 
         if (lineok){
           QwVQWK_Channel localchan;
           localchan.InitializeChannel(GetSubsystemName(),"QwBeamMod",localModChannelID.fmodulename,"raw");
           fModChannel.push_back(localchan);
           fModChannel[fModChannel.size()-1].SetDefaultSampleSize(sample_size);
           localModChannelID.fIndex=fModChannel.size()-1;
           fModChannelID.push_back(localModChannelID);
 
           // Store reference to ramp channel
           if (localModChannelID.fmodulename == "ramp") {
             fRampChannelIndex = fModChannel.size() - 1;
           }
         }
 
         if(ldebug)
         {
           localModChannelID.Print();
           std::cout<<"line ok=";
           if(lineok) std::cout<<"TRUE"<<std::endl;
           else std::cout<<"FALSE"<<std::endl;
         }
       }
 
 
       if(modtype == "WORD")
       {
         //  std::cout << "Decoding QwWord :: " << namech << std::endl;
 
 
         QwWord localword;
         localword.fSubbankIndex=currentsubbankindex;
         localword.fWordInSubbank=wordsofar;
         wordsofar+=1;
         // I assume that one data = one word here. But it is not always the case, for
         // example the triumf adc gives 6 words per channel
         localword.fModuleType=modtype;
         localword.fWordName=namech;
         //localword.fWordType=dettype; // FIXME dettype is undefined so commented this out
         fWord.push_back(localword);
 
         // Store reference to pattern number
         if (localword.fWordName == "bm_pattern_number") {
           fPatternWordIndex = fWord.size() - 1;
         }
 
         if(namech=="ffb_status")//save the location of this word to access this later
           fFFB_Index=fWord.size()-1;
 
 
         fWordsPerSubbank[currentsubbankindex].second = fWord.size();
         QwDebug << "--" << namech << "--" << fWord.size()-1 << QwLog::endl;
 
       }
     }
   }
   if(ldebug)
   {
     std::cout<<"Done with Load map channel \n";
     for(size_t i=0;i<fModChannelID.size();i++)
       fModChannelID[i].Print();
 
     for(size_t i=0;i<fWord.size();i++)
       fWord[i].PrintID();
   }
 
   // Now load the variables to monitor
   mapstr.RewindToFileStart();
   while (QwParameterFile *section = mapstr.ReadNextSection(varvalue)) {
     if (varvalue == "Monitors") {
       fMonitorNames.clear();
       while (section->ReadNextLine()) {
         section->TrimComment();         // Remove everything after a comment character
         section->TrimWhitespace();      // Get rid of leading and trailing spaces
         varvalue = section->GetTypedNextToken<TString>();
   if (varvalue.Length() > 0) {
           // Add names of monitor channels for each degree of freedom
     //
     //
     // Consider whether or not the input monitor is a bpm or not.  This should
     // allow for usage of monitors other than bpms, .eg bcm{1,2}, lumi, pmtonl,..ect.
     //
 
     if(varvalue.Contains("qwk_bpm", TString::kExact)){
       fMonitorNames.push_back(Form("%sX",varvalue.Data()));
       fMonitorNames.push_back(Form("%sY",varvalue.Data()));
     }
     else if(varvalue.Contains("target", TString::kExact)){
       fMonitorNames.push_back(Form("%sX",varvalue.Data()));
       fMonitorNames.push_back(Form("%sY",varvalue.Data()));
     }
     else{
       fMonitorNames.push_back(varvalue);
     }
         }
       }
     }
     delete section;
   }
   // Resize local version of the BPMs
   QwVQWK_Channel dummy("dummy");
   fMonitors.resize(fMonitorNames.size(),dummy);
   ResizeOpticsDataContainers(fMonitorNames.size());
 
   // Debug output
   if (ldebug) {
     QwMessage << "Done with loading monitor channels:" << QwLog::endl;
     for (size_t i = 0; i < fMonitorNames.size(); i++)
       QwMessage << fMonitorNames[i] << QwLog::endl;
   }
 
   return 0;
 }

Here is the call graph for this function:

Int_t QwBeamMod::LoadEventCuts ( TString filename )

virtual

Load the event cuts file.

Implements VQwSubsystemParity.

Definition at line 244 of file QwBeamMod.cc.

References QwLog::endl(), VQwSubsystem::fDetectorMaps, fFFB_holdoff, fModChannel, GetDetectorIndex(), GetDetectorTypeID(), GetGlobalErrorFlag(), QwParameterFile::GetUInt(), kBModErrorFlag, kQwUnknownDeviceType, and QwMessage.

 {
   Int_t eventcut_flag = 1;
 
   // Open the file
   QwParameterFile mapstr(filename.Data());
   fDetectorMaps.insert(mapstr.GetParamFileNameContents());
   fFFB_holdoff=0;//Default holdoff for the FFB pause
   while (mapstr.ReadNextLine()){
     //std::cout<<"********* In the loop  *************"<<std::endl;
     mapstr.TrimComment('!');   // Remove everything after a '!' character.
     mapstr.TrimWhitespace();   // Get rid of leading and trailing spaces.
     if (mapstr.LineIsEmpty())  continue;
 
     TString varname, varvalue;
     if (mapstr.HasVariablePair("=",varname,varvalue)){
       if (varname == "EVENTCUTS"){
   eventcut_flag = QwParameterFile::GetUInt(varvalue);
   //std::cout<<"EVENT CUT FLAG "<<eventcut_flag<<std::endl;
       }
     }
     else{
       TString device_type = mapstr.GetTypedNextToken<TString>();
       device_type.ToLower();
       TString device_name = mapstr.GetTypedNextToken<TString>();
       device_name.ToLower();
 
 
       if (device_type == "vqwk"){
 
   Double_t LLX = mapstr.GetTypedNextToken<Double_t>();  //lower limit for BCM value
   Double_t ULX = mapstr.GetTypedNextToken<Double_t>();  //upper limit for BCM value
   varvalue = mapstr.GetTypedNextToken<TString>();//global/loacal
   varvalue.ToLower();
   Double_t stabilitycut = mapstr.GetTypedNextToken<Double_t>();
   QwMessage<<"QwBeamMod Error Code  "<<GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)<<QwLog::endl;
   Int_t det_index=GetDetectorIndex(GetDetectorTypeID(kQwUnknownDeviceType),device_name);
   QwMessage << "*****************************" << QwLog::endl;
   QwMessage << " Type " << device_type << " Name " << device_name << " Index [" << det_index << "] "
             << " device flag " << eventcut_flag << QwLog::endl;
   fModChannel[det_index].SetSingleEventCuts((GetGlobalErrorFlag(varvalue,eventcut_flag,stabilitycut)|kBModErrorFlag),LLX,ULX,stabilitycut);
   QwMessage << "*****************************" << QwLog::endl;
 
       }
       else if (device_type == "word" && device_name== "ffb_status"){
   fFFB_holdoff=mapstr.GetTypedNextToken<UInt_t>();//Read the FFB OFF interval
       }
 
     }
 
   }
   //update the event cut ON/OFF for all the devices
   //std::cout<<"EVENT CUT FLAG"<<eventcut_flag<<std::endl;
   for (size_t i=0;i<fModChannel.size();i++)
     fModChannel[i].SetEventCutMode(eventcut_flag);
 
   return 0;
 }

Here is the call graph for this function:

Int_t QwBeamMod::LoadGeometry ( TString mapfile )

Int_t QwBeamMod::LoadInputParameters ( TString mapfile )

virtual

Mandatory parameter file definition.

Implements VQwSubsystem.

Definition at line 305 of file QwBeamMod.cc.

References fModChannel, and QwParameterFile::TrimComment().

 {
   Bool_t ldebug=kFALSE;
   TString varname;
 
   Int_t lineread=1;
 
   if(ldebug)std::cout<<"QwBeamMod::LoadInputParameters("<< pedestalfile<<")\n";
 
   QwParameterFile mapstr(pedestalfile.Data());  //Open the file
   while (mapstr.ReadNextLine())
     {
       lineread+=1;
       if(ldebug)std::cout<<" line read so far ="<<lineread<<"\n";
       mapstr.TrimComment('!');   // Remove everything after a '!' character.
       mapstr.TrimWhitespace();   // Get rid of leading and trailing spaces.
       if (mapstr.LineIsEmpty())  continue;
       else
   {
     TString varname = mapstr.GetTypedNextToken<TString>();  //name of the channel
     varname.ToLower();
     varname.Remove(TString::kBoth,' ');
     Double_t varped = mapstr.GetTypedNextToken<Double_t>(); // value of the pedestal
     Double_t varcal = mapstr.GetTypedNextToken<Double_t>(); // value of the calibration factor
     /* Double_t varweight = */ mapstr.GetTypedNextToken<Double_t>(); // value of the statistical weight
 
     //if(ldebug) std::cout<<"inputs for channel "<<varname
     //        <<": ped="<<varped<<": cal="<<varcal<<": weight="<<varweight<<"\n";
     for(size_t i=0;i<fModChannel.size();i++)
       if(fModChannel[i].GetElementName()==varname)
     {
       fModChannel[i].SetPedestal(varped);
       fModChannel[i].SetCalibrationFactor(varcal);
       break;
     }
   }
 
     }
   if(ldebug) std::cout<<" line read in the pedestal + cal file ="<<lineread<<" \n";
 
   ldebug=kFALSE;
   return 0;
 }

Here is the call graph for this function:

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

virtual

Implements VQwSubsystemParity.

Definition at line 584 of file QwBeamMod.cc.

References Compare(), fFFB_ErrorFlag, and fModChannel.

 {
   if(Compare(value))
     {
       //QwBeamMod* input= (QwBeamMod*)value ;
       QwBeamMod* input = dynamic_cast<QwBeamMod*>(value);
 
       for(size_t i=0;i<input->fModChannel.size();i++)
   this->fModChannel[i]+=input->fModChannel[i];
 //         for(size_t i=0;i<input->fWord.size();i++)
 //      this->fWord[i]+=input->fWord[i];
       this->fFFB_ErrorFlag |= input->fFFB_ErrorFlag;
     }
   return *this;
 }

Here is the call graph for this function:

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

virtual

Implements VQwSubsystemParity.

Definition at line 600 of file QwBeamMod.cc.

References Compare(), fFFB_ErrorFlag, and fModChannel.

 {
 
   if(Compare(value))
     {
       //QwBeamMod* input= (QwBeamMod*)value;
       QwBeamMod* input = dynamic_cast<QwBeamMod*>(value);
 
       for(size_t i=0;i<input->fModChannel.size();i++)
   this->fModChannel[i]-=input->fModChannel[i];
 //       for(size_t i=0;i<input->fWord.size();i++)
 //         this->fWord[i]-=input->fWord[i];
       this->fFFB_ErrorFlag |= input->fFFB_ErrorFlag;
 
     }
   return *this;
 }

Here is the call graph for this function:

VQwSubsystem & QwBeamMod::operator= ( VQwSubsystem * value )

virtual

Assignment Note: Must be called at the beginning of all subsystems routine call to operator=(VQwSubsystem *value) by VQwSubsystem::operator=(value)

Implements VQwSubsystemParity.

Definition at line 567 of file QwBeamMod.cc.

References Compare(), fFFB_ErrorFlag, fModChannel, and fWord.

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

Here is the call graph for this function:

void QwBeamMod::Print ( )

Definition at line 954 of file QwBeamMod.cc.

References QwLog::endl(), fModChannel, VQwSubsystem::fSystemName, fWord, VQwSubsystemParity::PrintValue(), and QwMessage.

Referenced by LoadChannelMap().

 {
   QwMessage << "Name of the subsystem =" << fSystemName << QwLog::endl;
 
   QwMessage << "there are " << fModChannel.size() << " mods" << QwLog::endl;
 
   QwMessage << " Printing Running AVG and other channel info for fModChannel" << QwLog::endl;
   for(size_t i=0;i<fModChannel.size();i++)
     fModChannel[i].PrintValue();
   for(size_t i=0;i<fWord.size();i++)
     fWord[i].Print();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void QwBeamMod::PrintErrorCounters ( ) const

virtual

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

Implements VQwSubsystemParity.

Definition at line 453 of file QwBeamMod.cc.

References fModChannel, QwVQWK_Channel::PrintErrorCounterHead(), and QwVQWK_Channel::PrintErrorCounterTail().

                                         {//inherited from the VQwSubsystemParity; this will display the error summary
 
   std::cout<<"*********QwBeamMod Error Summary****************"<<std::endl;
   QwVQWK_Channel::PrintErrorCounterHead();
   for(size_t i=0;i<fModChannel.size();i++){
     fModChannel[i].PrintErrorCounters();
   }
   QwVQWK_Channel::PrintErrorCounterTail();
 }

Here is the call graph for this function:

void QwBeamMod::PrintModChannelID ( )

Definition at line 967 of file QwBeamMod.cc.

References fModChannelID.

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

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

virtual

Implements VQwSubsystem.

Definition at line 514 of file QwBeamMod.cc.

 {
   return 0;
 }

Int_t QwBeamMod::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 350 of file QwBeamMod.cc.

References QwLog::endl(), fModChannel, fModChannelID, fWord, fWordsPerSubbank, VQwSubsystem::GetSubbankIndex(), QwDebug, QwMessage, QwWarning, and VQwSubsystem::SetDataLoaded().

 {
   Bool_t lkDEBUG=kFALSE;
 
   Int_t index = GetSubbankIndex(roc_id,bank_id);
 
   if (index>=0 && num_words>0){
     //  We want to process this ROC.  Begin looping through the data.
     if (lkDEBUG)
       QwMessage << "QwBeamMod::ProcessEvBuffer:  "
     << "Begin processing ROC" << roc_id
     << " and subbank "<<bank_id
     << " number of words="<<num_words<<QwLog::endl;
     if (buffer[0]==0xf0f0f0f0 && num_words%2==1){
       buffer++;
       if (lkDEBUG)
   QwMessage << "QwBeamMod::ProcessEvBuffer:  "
       << "Skipped padding word 0xf0f0f0f0 at beginning of buffer."
       << QwLog::endl;
     }
 
     for(size_t i=0;i<fModChannelID.size();i++)
       {
   if(fModChannelID[i].fSubbankIndex==index)
     {
 
     if (lkDEBUG)
       {
         std::cout<<"found modulation data for "<<fModChannelID[i].fmodulename<<std::endl;
         std::cout<<"word left to read in this buffer:"<<num_words-fModChannelID[i].fWordInSubbank<<std::endl;
       }
     fModChannel[fModChannelID[i].fIndex].  // Instead of BCM this will become our data word.If we use detectorID we need to add a specific ID
       ProcessEvBuffer(&(buffer[fModChannelID[i].fWordInSubbank]),
           num_words-fModChannelID[i].fWordInSubbank);
 
     }
       }
 
      for(Int_t i=fWordsPerSubbank[index].first; i<fWordsPerSubbank[index].second; i++) {
       if(fWord[i].fWordInSubbank+1<= (Int_t) num_words) {
   fWord[i].fValue=buffer[fWord[i].fWordInSubbank];
       } else {
   QwWarning << "QwBeamMod::ProcessEvBuffer:  There is not enough word in the buffer to read data for "
       << fWord[i].fWordName << QwLog::endl;
   QwWarning << "QwBeamMod::ProcessEvBuffer:  Words in this buffer:" << num_words
       << " trying to read word number =" << fWord[i].fWordInSubbank << QwLog::endl;
       }
     }
     if(lkDEBUG) {
       QwDebug << "QwBeamMod::ProcessEvBuffer:  Done with Processing this event" << QwLog::endl;
       for(size_t i=0;i<fWord.size();i++) {
   std::cout << " word number = " << i << " ";
   fWord[i].Print();
       }
     }
 
   }
   lkDEBUG=kFALSE;
 
   SetDataLoaded(kTRUE);
 
   return 0;
 }

Here is the call graph for this function:

void QwBeamMod::ProcessEvent ( )

virtual

Implements VQwSubsystem.

Definition at line 485 of file QwBeamMod.cc.

References fModChannel.

 {
   for (size_t i = 0; i < fModChannel.size(); i++) {
     // First apply HW checks and update HW  error flags.
     fModChannel[i].ApplyHWChecks();
     fModChannel[i].ProcessEvent();
   }
 }

void QwBeamMod::ProcessEvent_2 ( )

virtual

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

Reimplemented from VQwSubsystem.

Definition at line 508 of file QwBeamMod.cc.

References FillHistograms().

 {
   // Fill histograms here to bypass event cuts
   FillHistograms();
 }

Here is the call graph for this function:

void QwBeamMod::ProcessOptions ( QwOptions & options )

virtual

Process the command line options.

Reimplemented from VQwSubsystem.

Definition at line 32 of file QwBeamMod.cc.

                                                 {
       //Handle command line options
 }

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

virtual

Implements VQwSubsystemParity.

Definition at line 636 of file QwBeamMod.cc.

References Compare(), fModChannel, and fWord.

 {
   if(Compare(numer)&&Compare(denom))
     {
       //QwBeamMod* innumer= (QwBeamMod*)numer ;
       QwBeamMod* innumer = dynamic_cast<QwBeamMod*>(numer);
       //QwBeamMod* indenom= (QwBeamMod*)denom ;
       QwBeamMod* indenom = dynamic_cast<QwBeamMod*>(denom);
 
       for(size_t i=0;i<innumer->fModChannel.size();i++)
   this->fModChannel[i].Ratio(innumer->fModChannel[i],indenom->fModChannel[i]);
       for(size_t i=0;i<innumer->fWord.size();i++)
   this->fWord[i].fValue=innumer->fWord[i].fValue;
 
     }
 }

Here is the call graph for this function:

void QwBeamMod::ResizeOpticsDataContainers ( Int_t size )

Definition at line 804 of file QwBeamMod.cc.

References fAmplitude, fAmplitudeError, fChisquare, fNFitPoints, fNumberPatterns, fOffset, fOffsetError, fPhase, and fPhaseError.

Referenced by AnalyzeOpticsPlots(), and LoadChannelMap().

 {
   fOffset.resize(size);
   fAmplitude.resize(size);
   fPhase.resize(size);
   fOffsetError.resize(size);
   fAmplitudeError.resize(size);
   fPhaseError.resize(size);
   fChisquare.resize(size);
   fNFitPoints.resize(size);
 
   for(Int_t i = 0; i < size; i++){
     fOffset[i].resize(fNumberPatterns);
     fAmplitude[i].resize(fNumberPatterns);
     fPhase[i].resize(fNumberPatterns);
     fOffsetError[i].resize(fNumberPatterns);
     fAmplitudeError[i].resize(fNumberPatterns);
     fPhaseError[i].resize(fNumberPatterns);
     fChisquare[i].resize(fNumberPatterns);
     fNFitPoints[i].resize(fNumberPatterns);
   }
 }

Here is the caller graph for this function:

void QwBeamMod::Scale ( Double_t factor )

virtual

Implements VQwSubsystemParity.

Definition at line 654 of file QwBeamMod.cc.

References fModChannel.

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

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

virtual

Implements VQwSubsystemParity.

Definition at line 618 of file QwBeamMod.cc.

References Compare().

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

Here is the call graph for this function:

void QwBeamMod::UpdateErrorFlag ( const VQwSubsystem * ev_error )

virtual

update the error flag in the subsystem level from the top level routines related to stability checks. This will uniquely update the errorflag at each channel based on the error flag in the corresponding channel in the ev_error subsystem

Implements VQwSubsystemParity.

Definition at line 463 of file QwBeamMod.cc.

References fModChannel.

 {
   const QwBeamMod* input = dynamic_cast<const QwBeamMod*> (ev_error);
 
   for(size_t i=0;i<fModChannel.size();i++){
     fModChannel[i].UpdateErrorFlag(input->fModChannel[i]);
   }
 }