CDUnidentifiedHits.cxx

Go to the documentation of this file.

// Program name: CDUnidentifiedHits.cxx
//
// Package: CalDetTracker
//
// Purpose: 
//
// Contact: Chris Smith, Ryan Nichol, Leo Jenner or Jeff Hartnell

#include "TClonesArray.h"

#include "Calibrator/Calibrator.h"
#include "Plex/PlexSEIdAltL.h"
#include "Conventions/CalDigitType.h"
#include "MessageService/MsgService.h"
#include "Validity/VldContext.h"

#include "CalDetTracker/CDUnidentifiedHits.h"
#include "CalDetTracker/CDTrackedHitInfo.h"

using std::cout;
using std::endl;
using std::map;

CVSID("$Id: CDUnidentifiedHits.cxx,v 1.4 2004/07/26 11:02:09 hartnell Exp $");

//......................................................................

CDUnidentifiedHits::CDUnidentifiedHits
(map<int,CandStripHandle> allStrips, 
 map<int,CandStripHandle> evenPlStrips, 
 map<int,CandStripHandle> oddPlStrips, 
 map<int,CandStripHandle> CCStrips, 
 map<int,CandStripHandle> xTalkStrips)
{
  this->SetMap(allStrips,0);
  this->SetMap(evenPlStrips,1);
  this->SetMap(oddPlStrips,2);
  this->SetMap(CCStrips,3);
  this->SetMap(xTalkStrips,4);
}

//......................................................................

int CDUnidentifiedHits::SetMap(map<int,CandStripHandle> hits,int a)
{
  
  if(a==0) {
    map<int,CandStripHandle>::iterator digi_it = hits.begin();
    map<int,CandStripHandle>::iterator end_digi_it = hits.end();
    while(digi_it!=end_digi_it) {
      fAllStrips[digi_it->first] = digi_it->second;
      digi_it++;
    }
    return 1;
  }

  else if(a==1) {
    map<int,CandStripHandle>::iterator even_it = hits.begin();
    map<int,CandStripHandle>::iterator end_even_it = hits.end();
    while(even_it!=end_even_it) {
      fEvenPlStrips[even_it->first] = even_it->second;
      even_it++;
    }
    return 1;
  }

  else if(a==2) {
    map<int,CandStripHandle>::iterator odd_it = hits.begin();
    map<int,CandStripHandle>::iterator end_odd_it = hits.end();
    while(odd_it!=end_odd_it) {
      fOddPlStrips[odd_it->first] = odd_it->second;
      odd_it++;
    }
    return 1;
  }

  else if(a==3) {
    map<int,CandStripHandle>::iterator cc_it = hits.begin();
    map<int,CandStripHandle>::iterator end_cc_it = hits.end();
    while(cc_it!=end_cc_it) {
      fCCStrips[cc_it->first] = cc_it->second;
      cc_it++;
    }
    return 1;
  }

  else if(a==4) {
    map<int,CandStripHandle>::iterator xtalk_it = hits.begin();
    map<int,CandStripHandle>::iterator end_xtalk_it = hits.end();
    while(xtalk_it!=end_xtalk_it) {
      fXTalkStrips[xtalk_it->first] = xtalk_it->second;
      xtalk_it++;
    }
    return 1;
  }
  
  else 
    MSG("CDUnidentifiedHits",Msg::kWarning)
    << "CDUnidentifiedHits Error: " 
                 << "Not a valid map input number " 
                 << "0==all, 1==even planes, 2==odd planes, " 
                 << "3==cosmic counters, 4==X talk." << endl;

  return 0;
}

//......................................................................

void CDUnidentifiedHits::GetUnidentifiedStripInfo
(TClonesArray *infoarray)
{
  map<int,CandStripHandle>::iterator stpIt = fAllStrips.begin();
  map<int,CandStripHandle>::iterator stpEnd = fAllStrips.end();
  
  //get a reference to the calibrator
  Calibrator& calibrator=Calibrator::Instance();

  if(stpIt!=stpEnd) {
    const VldContext *vc = stpIt->second.GetVldContext();
    calibrator.Reset(*vc);
  }
  
  int cnt = 0;
  
  TClonesArray &unident = *infoarray;
  while(stpIt!=stpEnd){
    
    if(fEvenPlStrips.find(stpIt->first)==fEvenPlStrips.end() &&
       fOddPlStrips.find(stpIt->first)==fOddPlStrips.end() &&
       fXTalkStrips.find(stpIt->first)==fXTalkStrips.end() &&
       fCCStrips.find(stpIt->first)==fCCStrips.end()){
      
      CandStripHandle csh = stpIt->second;
      
      if(csh.GetNDigit(StripEnd::kNegative)>0){
        new(unident[cnt]) 
          CDTrackedHitInfo
          (StripEnd::kNegative,csh.GetPlane(),csh.GetStrip(),
           csh.GetCharge(CalDigitType::kNone,StripEnd::kNegative),
           csh.GetCharge(CalDigitType::kSigLin,StripEnd::kNegative),
           csh.GetCharge(CalDigitType::kSigCorr,StripEnd::kNegative),
           csh.GetCharge(CalDigitType::kPE,StripEnd::kNegative),
           calibrator.GetMIP(csh.GetCharge(CalDigitType::kSigCorr,
                                           StripEnd::kNegative)),
           0.,0.,0.,0.,0.,0.,csh.GetTime(StripEnd::kNegative));
        cnt++;
      }
      
      if(csh.GetNDigit(StripEnd::kPositive)>0){
        new(unident[cnt]) 
          CDTrackedHitInfo
          (StripEnd::kPositive,csh.GetPlane(),csh.GetStrip(),
           csh.GetCharge(CalDigitType::kNone,StripEnd::kPositive),
           csh.GetCharge(CalDigitType::kSigLin,StripEnd::kPositive),
           csh.GetCharge(CalDigitType::kSigCorr,StripEnd::kPositive),
           csh.GetCharge(CalDigitType::kPE,StripEnd::kPositive),
           calibrator.GetMIP(csh.GetCharge(CalDigitType::kSigCorr,
                                           StripEnd::kPositive)),
           0.,0.,0.,0.,0.,0.,csh.GetTime(StripEnd::kPositive));
        cnt++;
      }
      
    }
    stpIt++;
  }
  unident.Compress();

}

//......................................................................

---