AlgChopListFar Class Reference

#include <AlgChopListFar.h>

Inheritance diagram for AlgChopListFar:

List of all members.

Public Member Functions
	AlgChopListFar ()
virtual	~AlgChopListFar ()
virtual void	RunAlg (AlgConfig &ac, CandHandle &ch, CandContext &cx)
virtual void	Trace (const char *c) const

---

Detailed Description

Definition at line 13 of file AlgChopListFar.h.

---

Constructor & Destructor Documentation

AlgChopListFar::AlgChopListFar

(

)

Definition at line 46 of file AlgChopListFar.cxx.

{
}

AlgChopListFar::~AlgChopListFar

(

)

[virtual]

Definition at line 51 of file AlgChopListFar.cxx.

{
}

---

Member Function Documentation

void AlgChopListFar::RunAlg	(	AlgConfig &	algConfig,
		CandHandle &	candHandle,
		CandContext &	candContext
	)			[virtual]

Set functions are called to load the CandStrip member variables based on the CandDigit members, and the daughter list is filled.

Implements AlgBase.

Definition at line 60 of file AlgChopListFar.cxx.

References digits(), Form(), CandDigitListHandle::GetAbsTime(), AlgFactory::GetAlgHandle(), CandContext::GetDataIn(), Registry::GetDouble(), AlgFactory::GetInstance(), Registry::GetInt(), CandContext::GetMom(), Msg::kDebug, Msg::kWarning, CandDigitList::MakeCandidate(), MSG, and CandHandle::SetName().

{
  assert(candHandle.InheritsFrom("CandChopListHandle"));
  CandChopListHandle &sliceList = dynamic_cast<CandChopListHandle &>(candHandle);

   assert(candContext.GetDataIn());
   // Check for CandDigitListHandle input
   if (!(candContext.GetDataIn()->InheritsFrom("CandDigitListHandle"))) {
     MSG("Chop",Msg::kWarning ) << "Data into AlgChopListFar is not a digit list." << std::endl;
   }
   
   const CandDigitListHandle *cdlh_ptr = 
     dynamic_cast<const CandDigitListHandle*>(candContext.GetDataIn());
   
   const MomNavigator *mom = candContext.GetMom();
   const RawDigitDataBlock* rddb = DataUtil::GetRawBlock<RawDigitDataBlock>(mom);
   if (!rddb) {
     MSG("Chop", Msg::kWarning) << "No RawDigitDataBlock in RawRecord." << endl;
     return;
   }
   
   // Get setup for the DigitList maker algorithm
   AlgFactory &af = AlgFactory::GetInstance();
   AlgHandle ah = af.GetAlgHandle("AlgChop","default");
   CandContext cxx(this,candContext.GetMom());

   // Configuration:
   double cTimeGap           = algConfig.GetDouble("TimeGap");
   double cExtraShieldWindow = algConfig.GetDouble("ExtraShieldWindow");
   bool   cIgnorePreTrigger  = algConfig.GetInt("IgnorePreTrigger");
 

   // First, make a nice stl vector of the digits.
   DigitVector digits(cdlh_ptr);

   UInt_t ndigits = digits.size();
   if(ndigits==0) return; // No data.

   // Sort the list by time.
   std::sort(digits.begin(), digits.end(), compareDigitTimes());

   // Trigger time.
   double trigTime = cdlh_ptr->GetAbsTime();
   if(!cIgnorePreTrigger) trigTime = -100.0; // Set trigger time infinately early.

   // Find the first digit after the pre-trigger window.
   int first_digit = 0;
   for(UInt_t i=0;i<ndigits;i++) {
     if( (digits[i].GetTime()-trigTime) >0) {
       first_digit = i;
       break;
     }
   }

    
   // Go through, look for gaps. When you find a gap, make a new slice.

   int nslice = 0;
   double slice_start = digits[first_digit].GetTime();    // start time (inclusive)
   double slice_end   = slice_start;                      // stop  time (inclusive)
   double t_last = 99e9;
   bool have_slice = false;

   for(UInt_t i=first_digit; i<ndigits; i++) {     
     double t = digits[i].GetTime();
     bool create_slice = false;

     // Only look at non-veto digits.
     if(!(digits[i].GetPlexSEIdAltL().IsVetoShield()) ) {
       // This digit isn't in veto.. use it.
       double delta_t = t-t_last;
       t_last = t;
       if(!have_slice) {
         slice_start = t;
         have_slice = true;
       } else {
         if( (delta_t >= cTimeGap) ) create_slice = true;
       }
     }

     // Hack for last digit: close off the slice:
     if(i == ndigits-1) {
       if(have_slice) {
         create_slice = true;
         slice_end = 99e9;
       }
     }
     
     if(create_slice) {

       // Set start and end times. Note that the extra window
       // will only catch shield hits, assuming the extra window 
       // is smaller than the gap window.
       double tstart = slice_start - cExtraShieldWindow;
       double tend   = slice_end   + cExtraShieldWindow;
       
       // Make a new slice.
       MSG("Chop",Msg::kDebug) << "Forming new slice, digit " << slice_start
                               << " to " << slice_end << endl;
       DigitVector slice_data;
       for(UInt_t j=0; j<ndigits; j++) {
         double tt = digits[j].GetTime();            
         if((tt>=tstart)&&(tt<=tend))  slice_data.push_back(digits[j]);
       }
         
       cxx.SetDataIn(&(slice_data));
       CandDigitListHandle sliceHandle = CandDigitList::MakeCandidate(ah,cxx);
       sliceHandle.SetName(Form("Chop %d",nslice++));
       sliceList.AddDaughterLink(sliceHandle);
       
       // Prep next slice.
       slice_start = t;
     }

     slice_end = t;


   }      
   
}

void AlgChopListFar::Trace

(

const char *

c

)

const [virtual]

Reimplemented from AlgBase.

Definition at line 185 of file AlgChopListFar.cxx.

{
}

---

The documentation for this class was generated from the following files:

• AlgChopListFar.h
• AlgChopListFar.cxx

---