LIPatternFinderFancy.cxx

Go to the documentation of this file.

#include <TH1D.h>
#include <TH2D.h>

#include "CandDigit/CandDigitHandle.h"
#include "CandDigit/CandDigitListHandle.h"
#include "DatabaseInterface/DbiResultPtr.h"
#include "MessageService/MsgService.h"

#include "LIPatternFinderFancy.h"

ClassImp(LIPatternFinderFancy)

//_____________________________________________________________________
CVSID(
 "$Id: LIPatternFinderFancy.cxx,v 1.3 2004/04/07 20:44:30 costas Exp $"
);
//_____________________________________________________________________
LIPatternFinderFancy::LIPatternFinderFancy()
{
  MSG("LIPatternFinder", Msg::kVerbose)
                                << "LIPatternFinderFancy ctor" << endl;
  DefaultConfig();
}
//_____________________________________________________________________
LIPatternFinderFancy::~LIPatternFinderFancy()
{
  MSG("LIPatternFinder", Msg::kVerbose)
                                << "LIPatternFinderFancy dtor" << endl;
}
//_____________________________________________________________________
void LIPatternFinderFancy::DefaultConfig(void)
{
  fFractionalChargeThreshold = 0.75;
}  
//_____________________________________________________________________
void LIPatternFinderFancy::Configure(const Registry & registry)
{
  double tmpd;
  
  if (registry.Get("FractionalChargeThreshold", tmpd))
                                     fFractionalChargeThreshold = tmpd;
}  
//_____________________________________________________________________
bool LIPatternFinderFancy::IsLightInjectionTrash(
                                            CandDigitListHandle * cdlh)
{
  // Use the database interface, read database tables, save results in
  // STL maps for easy lookup
  
  BuildPlexMaps(cdlh); 

  TH2D * hPulserBoxLed = new TH2D("hPulserBoxLed",
        "(pulser box, led) q-weighted distribution", 16,0,16,20,0,20);

  CandDigitHandleItr cdh_iter ( cdlh->GetDaughterIterator() );

  double total_charge = 0;
  
  while ( CandDigitHandle * digit = cdh_iter.Ptr() ) {

    total_charge += digit->GetCharge();

    const PlexSEIdAltL & seid_altl = digit->GetPlexSEIdAltL();    

    seid_altl.SetFirst();

    while( seid_altl.IsValid() ) {

       const PlexSEIdAltLItem & seid_item = seid_altl.GetCurrentItem();

       PlexStripEndId seid = seid_item.GetSEId();

       const PlexLedId * led_id = FindLedIdFromStripEndId(seid);

       if(led_id) {

          MSG("LIPatternFinder", Msg::kDebug)
             << " Pbox: " << led_id->GetPulserBox()
             << " LED:  " << led_id->GetLedInBox()
             << " Q: "    << digit->GetCharge() << endl;

          hPulserBoxLed->Fill(led_id->GetPulserBox(),
                           led_id->GetLedInBox(), digit->GetCharge());           
       } // if led-id
                
       seid_altl.Next();
        
    } // loop over SEIdAltL  
                       
    cdh_iter.Next();
    
  } //loop over CandDigitList

  bool   is_LI = false;
  
  double max_bin_content  = hPulserBoxLed -> GetMaximum();

  if(total_charge > 0) {

     double fractional_pulse_height  =
                                    max_bin_content  / total_charge;
                                
     is_LI = (fractional_pulse_height > fFractionalChargeThreshold);

     MSG("LIPatternFinder", Msg::kDebug)
          << " Fractional Pulse Height: " << fractional_pulse_height
          << endl;
  }  
  delete hPulserBoxLed;         
                           
  return is_LI;
}
//_____________________________________________________________________
void LIPatternFinderFancy::BuildPlexMaps(CandDigitListHandle * cdlh)
{
  MSG("LIPatternFinder", Msg::kVerbose)
                      << "LIPatternFinderFancy::BuildPlexMaps" << endl;

  //-- Get VldContext from input data

  const VldContext * vld = cdlh->GetVldContext();

  //-- Use DatabseInterface to extract Plex information from the dbase

  MSG("LIPatternFinder", Msg::kVerbose)
          << "Extracting Plex information from the data-base" << endl;

  DbiResultPtr<PlexStripEndToLed>       strip_to_led   (*vld);

  MSG("LIPatternFinder", Msg::kVerbose)
       << "STRIP->LED / nrows: " << strip_to_led.GetNumRows() << endl;

  //-- loop over STRIP END -> LED rows
  for(unsigned int s2l_irow = 0;
                  s2l_irow < strip_to_led.GetNumRows(); s2l_irow++) {

     const PlexStripEndToLed * s2l_row =
                                       strip_to_led.GetRow(s2l_irow);

     PlexStripEndId stripend_id = s2l_row->GetPlexStripEndId();
     PlexLedId      led_id      = s2l_row->GetPlexLedId();

     fStrip2LedMap.insert(
           map<PlexStripEndId, PlexLedId>::value_type(
                                              stripend_id, led_id));
  }
  MSG("LIPatternFinder", Msg::kVerbose)
             << "PlexStripEndId -> PlexLedId map filled in" << endl;
}
//_____________________________________________________________________
const PlexLedId * LIPatternFinderFancy::FindLedIdFromStripEndId(
                              const PlexStripEndId & stripend_id) const
{
  map<PlexStripEndId, PlexLedId>::const_iterator strip2led_iter;

  for(strip2led_iter = fStrip2LedMap.begin();
             strip2led_iter != fStrip2LedMap.end(); ++strip2led_iter) {

      PlexStripEndId curr_stripend_id = strip2led_iter->first;

      if( curr_stripend_id.IsSameStripEnd(stripend_id) )
                                      return &(strip2led_iter->second);
  }
  return 0;
}
//_____________________________________________________________________

---