QwFeedback.cc File Reference

main(...) function for the qwanalysis feedback executable to set IA and/or PITA slopes More...

#include <iostream>
#include <fstream>
#include <vector>
#include <new>
#include <boost/shared_ptr.hpp>
#include "QwLog.h"
#include "QwRootFile.h"
#include "QwOptionsParity.h"
#include "QwEventBuffer.h"
#include "QwHistogramHelper.h"
#include "QwSubsystemArrayParity.h"
#include "QwHelicityCorrelatedFeedback.h"
#include "QwEventRing.h"
#include "QwEPICSEvent.h"
#include "QwHelicity.h"
#include "QwFakeHelicity.h"
#include "QwBeamLine.h"
#include "QwVQWK_Channel.h"

Include dependency graph for QwFeedback.cc:

Go to the source code of this file.

Functions
Int_t	main (Int_t argc, Char_t *argv[])

Detailed Description

main(...) function for the qwanalysis feedback executable to set IA and/or PITA slopes

Definition in file QwFeedback.cc.

Function Documentation

main()

Int_t main	(	Int_t	argc,
		Char_t *	argv[] )

First, fill the search paths for the parameter files; this sets a static variable within the QwParameterFile class which will be used by all instances. The "scratch" directory should be first.

Setup screen and file logging

Load the histogram parameter definitions (from parity_hists.txt) into the global histogram helper: QwHistogramHelper

Create the event buffer

Create the database connection

Start loop over all runs

Set the current event number for parameter file lookup

Create an EPICS event

Load the detectors from file

Create the helicity pattern

Create the event ring with the subsysten array

Create the running sum

Definition at line 42 of file QwFeedback.cc.

{
  ///  First, fill the search paths for the parameter files; this sets a
  ///  static variable within the QwParameterFile class which will be used by
  ///  all instances.
  ///  The "scratch" directory should be first.
  //QwParameterFile::AppendToSearchPath(getenv_safe_string("QW_PRMINPUT"));
  QwParameterFile::AppendToSearchPath(getenv_safe_string("QWANALYSIS") + "/Extensions/QwFeedback/prminput");
  //QwParameterFile::AppendToSearchPath(getenv_safe_string("QWANALYSIS") + "/Parity/prminput");
  //QwParameterFile::AppendToSearchPath(getenv_safe_string("QWANALYSIS") + "/Analysis/prminput");
 
 
  // Then set the command line arguments and the configuration filename, 
  // and we define the options that can be used in them (using QwOptions).
  gQwOptions.SetCommandLine(argc, argv);
  gQwOptions.AddConfigFile("qwfeedback.conf");
  //gQwOptions.ListConfigFiles();
  //gQwOptions.AddConfigFile("qweak_mysql.conf");
    // Define the command line options
  DefineOptionsParity(gQwOptions);
  QwHelicityCorrelatedFeedback::DefineOptions(gQwOptions);
  //Load command line options for the histogram/tree helper class
  gQwHists.ProcessOptions(gQwOptions);
   /// Setup screen and file logging
  gQwLog.ProcessOptions(&gQwOptions);
 
  ///  Load the histogram parameter definitions (from parity_hists.txt) into the global
  ///  histogram helper: QwHistogramHelper
  //gQwHists.LoadHistParamsFromFile("qweak_parity_hists.in");
  //gQwHists.LoadTreeParamsFromFile("Qweak_Tree_Trim_List.in");
 
 
  ///  Create the event buffer
  QwEventBuffer eventbuffer;
  eventbuffer.ProcessOptions(gQwOptions);
 
 
 
  ///  Create the database connection
  //  QwDatabase database(gQwOptions);
 
 
  ///  Start loop over all runs
  //  QwRootFile* rootfile = 0;
 
  
  /*
  ///  Create an EPICS control event
  QwEPICSControl fEPICSCtrl;
  GreenMonster   fScanCtrl;
  //  Let's clear the values for the scan data & set "not clean".
  fScanCtrl.Open();
  fScanCtrl.SCNSetStatus(SCN_INT_NOT);
  fScanCtrl.SCNSetValue(1,0);
  fScanCtrl.SCNSetValue(2,0);
  fScanCtrl.CheckScan();
  fScanCtrl.PrintScanInfo();
  fScanCtrl.Close();
  
  fEPICSCtrl.Print_Qasym_Ctrls();
  */
    
  // Loop over all runs
  while (eventbuffer.OpenNextStream() == CODA_OK){
    //  Begin processing for the first run.
 
    ///  Set the current event number for parameter file lookup
    QwParameterFile::SetCurrentRunNumber(eventbuffer.GetRunNumber());
 
 
    ///  Create an EPICS event
    QwEPICSEvent epicsevent;
    epicsevent.LoadChannelMap("EpicsTable.map");
 
    ///  Load the detectors from file
    QwSubsystemArrayParity detectors(gQwOptions);
    detectors.ProcessOptions(gQwOptions);
 
    ///  Create the helicity pattern
    QwHelicityCorrelatedFeedback helicitypattern(detectors);
    helicitypattern.ProcessOptions(gQwOptions);
    helicitypattern.LoadParameterFile("qweak_fb_prm.in");
    /*
      helicitypattern.UpdateGMClean(0);
      helicitypattern.FeedIASetPoint(0);
      helicitypattern.UpdateGMClean(1);
    */
 
    ///  Create the event ring with the subsysten array
    QwEventRing eventring(gQwOptions,detectors);
 
    ///  Create the running sum
    QwSubsystemArrayParity runningsum(detectors);
 
 
    Int_t failed_events_counts = 0; // count failed total events
 
 
    //  Clear the single-event running sum at the beginning of the runlet
    runningsum.ClearEventData();
    helicitypattern.ClearRunningSum();
    //  Clear the running sum of the burst values at the beginning of the runlet
    helicitypattern.ClearBurstSum();
 
 
    
 
    // Loop over events in this CODA file
    while (eventbuffer.GetNextEvent() == CODA_OK) {
 
      if (eventbuffer.GetEventNumber()%1000)
    std::cout<<std::flush;
 
      //  First, do processing of non-physics events...
      if (eventbuffer.IsROCConfigurationEvent()){
    //  Send ROC configuration event data to the subsystem objects.
    eventbuffer.FillSubsystemConfigurationData(detectors);
      }
 
      //  Dump out of the loop when we see the end event.
      if (eventbuffer.GetEndEventCount()>0){
    QwMessage << "Number of events processed at end of run: " << eventbuffer.GetEventNumber() << std::endl;
    break;
      }
      
      //  Now, if this is not a physics event, go back and get a new event.
      if (!eventbuffer.IsPhysicsEvent()) continue;
 
 
      //  Fill the subsystem objects with their respective data for this event.
      eventbuffer.FillSubsystemData(detectors);
 
      //  Process the subsystem data
      detectors.ProcessEvent();
 
       // The event pass the event cut constraints
      if (detectors.ApplySingleEventCuts()) {
 
        // Accumulate the running sum to calculate the event based running average
        runningsum.AccumulateRunningSum(detectors);
 
//         // Fill the histograms
//         rootfile->FillHistograms(detectors);
 
//         // Fill the tree branches
//         rootfile->FillTreeBranches(detectors);
//         rootfile->FillTree("Mps_Tree");
 
        // Add event to the ring
        eventring.push(detectors);
 
        // Check to see ring is ready
        if (eventring.IsReady()) {
 
          // Load the event into the helicity pattern
          helicitypattern.LoadEventData(eventring.pop());
 
          // Calculate helicity pattern asymmetry
          if (helicitypattern.IsCompletePattern()) { 
 
            // Update the blinder if conditions have changed
            helicitypattern.UpdateBlinder(detectors);
 
            // Calculate the asymmetry
            helicitypattern.CalculateAsymmetry();
            if (helicitypattern.IsGoodAsymmetry()) {
          helicitypattern.ApplyFeedbackCorrections();//apply IA feedback
              // Clear the data
              helicitypattern.ClearEventData();       
            }
 
          } // helicitypattern.IsCompletePattern()
 
        } // eventring.IsReady()
 
 
      // Failed single event cuts
      } else {
    failed_events_counts++;
      }
 
 
      
 
    }   // end of loop over events  
 
    QwMessage << "Number of events processed at end of run: "
              << eventbuffer.GetEventNumber() << std::endl;
 
 
 
    //  Close event buffer stream
    eventbuffer.CloseStream();
 
 
 
    //  Print the event cut error summery for each subsystem
    //    detectors.GetEventcutErrorCounters();
 
 
    //  Read from the datebase
    //    database.SetupOneRun(eventbuffer);
 
    // Each sussystem has its own Connect() and Disconnect() functions.
    //    if (database.AllowsWriteAccess()) {
    //      helicitypattern.FillDB(&database);
    //      epicsevent.FillDB(&database);
    //    }
 
 
    QwMessage << "Total events failed " << failed_events_counts << QwLog::endl;
 
    //  Report run summary
    eventbuffer.ReportRunSummary();
    eventbuffer.PrintRunTimes();
 
  } //end of run loop
 
 
  QwMessage << "I have done everything I can do..." << QwLog::endl; 
 
  return 0;
}

References QwSubsystemArrayParity::AccumulateRunningSum(), QwParameterFile::AppendToSearchPath(), QwHelicityCorrelatedFeedback::ApplyFeedbackCorrections(), QwSubsystemArrayParity::ApplySingleEventCuts(), QwHelicityCorrelatedFeedback::CalculateAsymmetry(), QwHelicityPattern::ClearEventData(), QwSubsystemArray::ClearEventData(), QwHelicityCorrelatedFeedback::ClearRunningSum(), QwEventBuffer::CloseStream(), QwHelicityCorrelatedFeedback::DefineOptions(), DefineOptionsParity(), QwLog::endl(), QwEventBuffer::FillSubsystemConfigurationData(), QwEventBuffer::FillSubsystemData(), getenv_safe_string(), QwEventBuffer::GetEventNumber(), QwEventBuffer::GetNextEvent(), QwEventBuffer::GetRunNumber(), gQwHists, gQwLog, gQwOptions, QwHelicityPattern::IsCompletePattern(), QwHelicityPattern::IsGoodAsymmetry(), QwEventBuffer::IsPhysicsEvent(), QwEventRing::IsReady(), QwEventBuffer::IsROCConfigurationEvent(), QwEPICSEvent::LoadChannelMap(), QwHelicityPattern::LoadEventData(), QwHelicityCorrelatedFeedback::LoadParameterFile(), QwEventBuffer::OpenNextStream(), QwEventRing::pop(), QwEventBuffer::PrintRunTimes(), QwSubsystemArray::ProcessEvent(), QwEventBuffer::ProcessOptions(), QwHelicityCorrelatedFeedback::ProcessOptions(), QwSubsystemArray::ProcessOptions(), QwEventRing::push(), QwMessage, QwEventBuffer::ReportRunSummary(), QwParameterFile::SetCurrentRunNumber(), and QwHelicityPattern::UpdateBlinder().

Here is the call graph for this function: