CDTrackerModule Class Reference

#include <CDTrackerModule.h>

Inheritance diagram for CDTrackerModule:

List of all members.

Public Member Functions
	CDTrackerModule ()
JobCResult	Ana (const MomNavigator *mom)
void	HandleCommand (JobCommand *cmd)
void	EndJob ()
Private Member Functions
void	CleanDataMembers ()
void	FillPID (CandCalDetSIHandle *CDSI)
void	FillPID (const CandCalDetPIDHandle *cdpid)
void	InitialiseTree ()
TFile *	OpenFile (Int_t runNumber, Int_t runNumberSub)
Private Attributes
bool	ffirstrecord
Float_t	fBeamMomentum
Int_t	fNumDeadChips
Float_t	fTemperature
SimFlag::SimFlag_t	fSimFlag
Int_t	frun
Int_t	fsubrun
Int_t	fsnarl
Int_t	fSec
Int_t	fStartSnarl
Int_t	fEndSnarl
Int_t	fStartTime
Int_t	fEndTime
Char_t	ftag [256]
TFile *	fsavefile
CDTrackerOptions *	fTrkOpt
TTree *	ftree
CDTrackInfo *	fTrkInfo
TClonesArray *	fTrkHits
TClonesArray *	fStTrkHits
TClonesArray *	fCCHits
TClonesArray *	fXTalkHits
TClonesArray *	fUnTrkHits
TClonesArray *	fTruthHits
CDPIDInfo *	fPID
Int_t	fJustWriteTracks
Int_t	fOverwriteFile
Int_t	fWritePID
Int_t	fWriteTruth

---

Detailed Description

Definition at line 30 of file CDTrackerModule.h.

---

Constructor & Destructor Documentation

CDTrackerModule::CDTrackerModule

(

)

Definition at line 68 of file CDTrackerModule.cxx.

References fBeamMomentum, fCCHits, fEndSnarl, fEndTime, ffirstrecord, fJustWriteTracks, fNumDeadChips, fOverwriteFile, fPID, frun, fsavefile, fSec, fSimFlag, fsnarl, fStartSnarl, fStartTime, fStTrkHits, fsubrun, ftag, fTemperature, ftree, fTrkHits, fTrkInfo, fTrkOpt, fTruthHits, fUnTrkHits, fWritePID, fWriteTruth, fXTalkHits, and SimFlag::kUnknown.

{
  ffirstrecord=true;

  fBeamMomentum=0;
  fNumDeadChips=-1;
  fTemperature=-1;
  fSimFlag=SimFlag::kUnknown;
  frun=0;
  fsubrun=0;
  fsnarl=0;
  fSec=-1;
  fStartSnarl=0;
  fEndSnarl=0;
  fStartTime=0;
  fEndTime=0;
  fsavefile=0;
  fJustWriteTracks=0;
  fOverwriteFile=0;
  fWritePID=0;
  sprintf(ftag,"SetFileIDToMakeThisShorter");
  fWriteTruth=0;

  //initialise pointers
  fTrkInfo=0;
  ftree=0;
  fPID=0;

  //"beam" tracked hits
  fTrkHits=new TClonesArray("CDTrackedHitInfo",1000);
  //cosmic counter hits
  fCCHits=new TClonesArray("CDTrackedHitInfo",1000); 
  //"cosmic" tracked hits
  fStTrkHits=new TClonesArray("CDTrackedHitInfo",1000);
  //xtalk hits 
  fXTalkHits=new TClonesArray("CDXTalkHitInfo",1000);
  //unidentified hits
  fUnTrkHits=new TClonesArray("CDTrackedHitInfo",1000);
  //the truth of the hits
  fTruthHits=new TClonesArray("CDTruthHitInfo",1000);

  //tracker options
  fTrkOpt=new CDTrackerOptions();
}

---

Member Function Documentation

JobCResult CDTrackerModule::Ana

(

const MomNavigator *

mom

)

[virtual]

Implement this for read only access to the MomNavigator

Reimplemented from JobCModule.

Definition at line 115 of file CDTrackerModule.cxx.

References SimFlag::AsString(), CleanDataMembers(), CDTrackerOptions::Cosmic(), count, fBeamMomentum, fCCHits, fEndSnarl, fEndTime, ffirstrecord, FillPID(), CandRecord::FindCandHandle(), RawRecord::FindRawBlock(), fJustWriteTracks, fNumDeadChips, fPID, frun, fsavefile, fSec, fSimFlag, fsnarl, fStartSnarl, fStartTime, fStTrkHits, fsubrun, fTemperature, ftree, fTrkHits, fTrkInfo, fTrkOpt, fTruthHits, fUnTrkHits, fWritePID, fWriteTruth, fXTalkHits, CandDigitListHandle::GetAbsTime(), CDTrackInfo::GetAngle(), CDTracker::GetCCStripMap(), RawDigitDataBlock::GetDatumIter(), MomNavigator::GetFragment(), RawDeadChipBlock::GetNumberOfEntries(), CandCalDetPIDHandle::GetOLChi2(), CandCalDetPIDHandle::GetPIDType(), RawRecord::GetRawHeader(), CDTrackInfo::GetResult(), RawDaqHeader::GetRun(), VldTimeStamp::GetSec(), VldContext::GetSimFlag(), RawDaqSnarlHeader::GetSnarl(), CDMapMaker::GetStripMap(), RawDaqHeader::GetSubRun(), Calibrator::GetTemperature(), VldContext::GetTimeStamp(), CDTracker::GetTrackedStripMap(), CDTracker::GetTrackInfo(), CDPIDInfo::GetTrueParticleMomentum(), CDTruthifier::GetTruthHitInfo(), RecMinosHdr::GetVldContext(), CandCalDetPIDHandle::InCERTime(), InitialiseTree(), Calibrator::Instance(), RawChannelId::IsSameChannel(), Detector::kCalDet, SimFlag::kData, Msg::kDebug, Dbi::kDefaultTask, Dbi::kDisabled, Msg::kError, Msg::kFatal, Msg::kInfo, JobCResult::kPassed, ElecType::kVA, Msg::kVerbose, Msg::kWarning, PlaneView::kX, PlaneView::kY, MAXMSG, MSG, CandCalDetPIDHandle::NoOverlap(), OpenFile(), CDTruthifier::ParticleMomentum(), CDTruthifier::RecoStdHep(), CalScheme::Reset(), CDTrackInfo::SetAbsTime(), JobCResult::SetFailed(), CDTrackInfo::SetRange(), and CDTrackInfo::SetTriggerTime().

{
  //set the default result to passed
  JobCResult result(JobCResult::kPassed);

  CandRecord* candrec=dynamic_cast<CandRecord*>
    (mom->GetFragment("CandRecord", "PrimaryCandidateRecord"));
  if (candrec==0) {
    MSG("CDTrackerModule",Msg::kWarning) 
      <<"No PrimaryCandidateRecord in MOM." << endl;
    result.SetFailed();
    return result;
  }

  RawRecord *rr=dynamic_cast<RawRecord *>
    (mom->GetFragment("RawRecord"));
  if (rr==0) {
    MSG("CDTrackerModule",Msg::kWarning) 
      <<"No RawRecord in MOM." << endl;
    result.SetFailed();
    return result;
  }

  const RawDeadChipBlock* dcb=dynamic_cast<const RawDeadChipBlock*>
    (rr->FindRawBlock("RawDeadChipBlock"));
  if(dcb) {
    fNumDeadChips=dcb->GetNumberOfEntries();
    if (fNumDeadChips>0) {
      MAXMSG("CDTrackerModule",Msg::kDebug,100)
        <<"Found dead chips="<<fNumDeadChips<<endl;
    }
  }
  else {
    MAXMSG("CDTrackerModule",Msg::kWarning,10)
      <<"No RawDeadChipBlock present"<<endl;
  }

  const RawDaqSnarlHeader* snarlHdr =
    dynamic_cast<const RawDaqSnarlHeader*>(rr->GetRawHeader());
  if(snarlHdr){
    fsnarl=snarlHdr->GetSnarl();
    fSec=snarlHdr->GetVldContext().GetTimeStamp().GetSec();
    fEndSnarl=fsnarl;
    fEndTime=fSec;

    if(ffirstrecord) {
      fSimFlag=snarlHdr->GetVldContext().GetSimFlag();
      fStartTime=snarlHdr->GetVldContext().GetTimeStamp().GetSec();
      fStartSnarl=snarlHdr->GetSnarl();
      frun=snarlHdr->GetRun();
      fsubrun=snarlHdr->GetSubRun();//this isn't always right!

      /*
      DbiResultPtr<DbuSubRunSummary> dsrs(snarlHdr->GetVldContext(),
                                          Dbi::kDefaultTask,
                                          Dbi::kDisabled);
      MSG("CDTrackerModule",Msg::kInfo)
        <<"Number of rows in DbuSubRunSummary="<<dsrs.GetNumRows()
        <<endl;
      
      //get the database row
      if(dsrs.GetNumRows()==1) {
        fsubrun=dsrs.GetRow(0)->fSubRun;
        MSG("CDTrackerModule",Msg::kInfo)
          <<"Found subrun in database="<<fsubrun
          <<", run="<<dsrs.GetRow(0)->fRun<<endl
          <<"Start time="<<dsrs.GetRow(0)->fStartTime.GetSec()<<endl
          <<"End time  ="<<dsrs.GetRow(0)->fEndTime.GetSec()<<endl;
      }
      else{
        MSG("CDTrackerModule",Msg::kWarning)
          <<"Can't find subrun in database, #rows="<<dsrs.GetNumRows()
          <<", fsubrun="<<fsubrun<<endl;
      }
      */
      MSG("CDTrackerModule",Msg::kVerbose)<<"Subrun="<<fsubrun<<endl;
    }
  }
  else{
    MSG("CDTrackerModule",Msg::kWarning) 
      <<"No RawDaqSnarlHeader in MOM."<<endl;
    result.SetFailed();
    return result;
  }

  //get the temperature from the calibrator
  Calibrator& cal=Calibrator::Instance();
  cal.Reset(snarlHdr->GetVldContext());
  fTemperature=cal.GetTemperature();
  MAXMSG("CDTrackerModule",Msg::kDebug,20) 
    <<"Used Calibrator to get temperature="<<fTemperature<<" C"<<endl;

  int ExtTrigTime=0;//stores external trigger time of specific channel
  const RawDigitDataBlock *rddb=dynamic_cast<const RawDigitDataBlock *>
    (rr->FindRawBlock("RawDigitDataBlock"));

  //looping to find particular channel
  if (rddb) {
    RawChannelId id(Detector::kCalDet,ElecType::kVA,0,0x1602);
    id.SetModeBits(true,true,true);
    TIter it=rddb->GetDatumIter();
    while (TObject *obj=it.Next()) {
      RawVaDigit *rvd=dynamic_cast<RawVaDigit *>(obj);
      if (rvd) {
        RawChannelId rcid=rvd->GetChannel();
        if (rcid.IsSameChannel(id)) ExtTrigTime=rvd->GetTDC();
      }
    }
  }
    
  //work out the directory to write to
  if(ffirstrecord) {
    MSG("CDTrackerModule",Msg::kDebug)
      <<"Opening file..."<<endl;
    //open the output file
    fsavefile=this->OpenFile(frun,fsubrun);
    
    //create tree and initialise branches
    this->InitialiseTree();
    ffirstrecord=false;
    MSG("CDTrackerModule",Msg::kDebug) 
      <<"... file opened and tree initialised"<<endl;

    if (fSimFlag==SimFlag::kData){
      //note that only beam momentum for data is set here
      //for MC it is the largest stdhep TParticle momentum
      DbiResultPtr<CalDetBeamMomentum> cdbm(snarlHdr->GetVldContext(),
                                            Dbi::kDefaultTask,
                                            Dbi::kDisabled);
      
      MSG("CDTrackerModule",Msg::kDebug)
        <<"Number of rows in CalDetBeamMomentum="<<cdbm.GetNumRows()
        <<endl;
      
      //get the database row
      if(cdbm.GetNumRows()==1) {
        fBeamMomentum=cdbm.GetRow(0)->GetMomentum();
        MSG("CDTrackerModule",Msg::kInfo)
          <<"Found beam momentum in database="<<fBeamMomentum<<endl;
      }
      else{
        MSG("CDTrackerModule",Msg::kInfo)
          <<"Problem with beam momentum in database"
          <<", #rows="<<cdbm.GetNumRows()
          <<", beamMomentum="<<fBeamMomentum<<endl;
      }
    }
  }
  static Bool_t firstSnarl=true;
  if (firstSnarl){
    MSG("CDTrackerModule",Msg::kInfo)
      <<"Started tracking... (first snarl number="<<fsnarl
      <<", second="<<fSec<<")"<<endl;
    firstSnarl=false;
  }

  if (fWritePID){
    //make PID Info object
    fPID=new CDPIDInfo();

    //Get CandCalDetSIHandle and pass it to CDPIDInfo constructor
    CandCalDetSIHandle* cdsi=dynamic_cast<CandCalDetSIHandle*>
      (candrec->FindCandHandle("CandCalDetSIHandle"));
    
    //now fill CDPIDInfo with SI info object
    if(cdsi) {
      this->FillPID(cdsi);
    }
    else{      
      if (fSimFlag==SimFlag::kData){
        MSG("CDTrackerModule",Msg::kError)
          <<"No caldet SI handle; not filling CDPIDInfo!"<<endl;
      }
      //else simflag=MC and you don't expect to find such a handle
    }

    CandCalDetPIDHandle* cdpid=dynamic_cast<CandCalDetPIDHandle*>
      (candrec->FindCandHandle("CandCalDetPIDHandle"));

    //now fill CDPIDInfo object from cdpid
    if (cdpid){
      MSG("CalDetPID",Msg::kDebug) 
        <<"Tracker::particle="<<cdpid->GetPIDType()
        <<", nooverlap="<<cdpid->NoOverlap()
        <<", inckvtime="<<cdpid->InCERTime()
        <<", ol chi2="<<cdpid->GetOLChi2()
        <<endl;
      this->FillPID(cdpid);
    }
    else{
      if (fSimFlag==SimFlag::kData){
      MAXMSG("CDTrackerModule",Msg::kWarning,5) 
        <<"No caldet pid handle; not filling CDPIDInfo!"<<endl;
      }
      //else simflag=MC and you don't expect to find such a handle
    }
  }

  //Get CandDigitList
  CandDigitListHandle *digitlist=dynamic_cast<CandDigitListHandle*>
    (candrec->FindCandHandle("CandDigitListHandle"));
  if(digitlist==0) {
    MSG("CDTrackerModule", Msg::kFatal) 
      << "No CandDigitList in MOM." << endl;
    result.SetFailed();
    return result;
  }

  //Get CandStripList
  CandStripListHandle *striplist=dynamic_cast<CandStripListHandle*>
    (candrec->FindCandHandle("CandStripListHandle"));
  if(striplist==0) {
    MSG("CDTrackerModule", Msg::kFatal) 
      << "No CandStripList in MOM." << endl;
    result.SetFailed();
    return result;
  }

  //Write the truth information
  if (fWriteTruth){
    //check if the input file is not data
    if (fSimFlag!=SimFlag::kData){
      MAXMSG("CDTrackerModule", Msg::kVerbose,5)
        <<"fSimFlag="<<SimFlag::AsString(fSimFlag)
        <<". Accessing truth info..."<<endl; 
      //create a CDTruthifier
      CDTruthifier cdTruthifier;
      cdTruthifier.GetTruthHitInfo(mom,fTruthHits);

      Float_t beamP=0;

      //get the particle momentum
      //note that for kData the beam momentum is obtained 
      //from the caldet database table
      if (fWritePID && fPID){
        cdTruthifier.RecoStdHep(mom,fPID);
        beamP=fPID->GetTrueParticleMomentum();
      }
      else if (!fPID){
        //still want to get the beam momentum
        //even if not writing the whole CDPIDInfo object
        beamP=cdTruthifier.ParticleMomentum(mom);
      }

      //calculate the average beam momentum
      static Float_t count=0;
      MAXMSG("CDTrackerModule",Msg::kVerbose,200)
        <<"beamP="<<beamP<<", fBeamMomentum="<<fBeamMomentum
        <<", count="<<count<<endl;
      fBeamMomentum=fBeamMomentum*(count/(count+1))+beamP*(1/(count+1));
      count++;
    }
  }

  //make a map of the strips
  CDMapMaker *mpmkrStrips = new CDMapMaker(striplist);

  //cross talk removal code:
  CDCrossTalker xtalker(mpmkrStrips->GetStripMap(),0.9,
                        fTrkOpt->Cosmic());
  xtalker.GetXTalkStripInfo(fXTalkHits);

  //create tracker object
  CDTracker *tracker;

  //pass crosstalk removed map to the tracker:
  //if cosmic option set, use cosmic tracker otherwise beam tracker
  if(fTrkOpt->Cosmic()){
    tracker=new CDCosmicTracker(xtalker.GetCleanStripMap(),fTrkOpt);
  }
  else{
    tracker=new CDBeamMuTracker(xtalker.GetCleanStripMap(),fTrkOpt);
  }
  
  //Get track info from the tracker
  fTrkInfo=tracker->GetTrackInfo();
  
  //Set a few things that can't be set in the CDTracker objects:
  fTrkInfo->SetAbsTime(digitlist->GetAbsTime());//leave as digit
  fTrkInfo->SetTriggerTime(ExtTrigTime);
  
  //calculate pathlengths etc. for each hit in a track:
  CDTrackedHits trkHits(tracker->GetTrackedStripMap(PlaneView::kX),
                        tracker->GetTrackedStripMap(PlaneView::kY),
                        tracker->GetCCStripMap());
  
  //Use angles from the track info class to calculate track quantities:
  trkHits.GetTrackedStripInfo(fStTrkHits,fTrkInfo->GetAngle(0),
                              fTrkInfo->GetAngle(1));
  
  //Get Cosmic counter hit info
  trkHits.GetCCStripInfo(fCCHits);
  
  //Set trk range in TrackInfo object
  fTrkInfo->SetRange(trkHits.GetTrkRange());
  
  //alternative way to calculate pathlengths
  //only necessary for beam muons
  if(!fTrkOpt->Cosmic()){
    trkHits.GetTrackedStripInfo(fTrkHits);  
    //overwrite trk range if BeamTracker used
    fTrkInfo->SetRange(trkHits.GetTrkRange());
  }

  //Find unidentified hits in event
  //need to change the cross talker
  CDUnidentifiedHits uihStrips
    (mpmkrStrips->GetStripMap(),
     tracker->GetTrackedStripMap(PlaneView::kX),
     tracker->GetTrackedStripMap(PlaneView::kY),
     tracker->GetCCStripMap(),
     xtalker.GetXTalkStripMap());

  //get the unidentified strips
  uihStrips.GetUnidentifiedStripInfo(fUnTrkHits);
  
  //fill tree with info:
  if(fTrkInfo->GetResult(0)==1 && fTrkInfo->GetResult(1)==1){
    ftree->Fill();
  }
  else {
    result.SetFailed();
    if(!fJustWriteTracks)
      ftree->Fill();
  }
  
  //delete objects pointed to by local pointers
  delete tracker;
  delete mpmkrStrips;

  //clear the TClonesArrays and delete objects
  this->CleanDataMembers();

  return result;
}

void CDTrackerModule::CleanDataMembers

(

)

[private]

Definition at line 454 of file CDTrackerModule.cxx.

References fCCHits, fPID, fStTrkHits, fTrkHits, fTrkInfo, fTruthHits, fUnTrkHits, and fXTalkHits.

Referenced by Ana().

{
  //Clear TClonesArrays for next event
  //Note: This used to use Delete but that is slower and only 
  //necessary for objects (stored in the TClonesArray) that 
  //allocate their own memory
  fTrkHits->Clear();
  fStTrkHits->Clear();
  fCCHits->Clear();
  fXTalkHits->Clear();
  fUnTrkHits->Clear();
  fTruthHits->Clear();
  
  //delete objects 
  if (fTrkInfo) delete fTrkInfo;
  if (fPID) delete fPID;

  //reinitialise pointers
  fTrkInfo=0;
  fPID=0;

  //Note: Don't delete the fTrkOpt, it is only created once
}

void CDTrackerModule::EndJob

(

)

[virtual]

Implement for notification of end of job

Reimplemented from JobCModule.

Definition at line 583 of file CDTrackerModule.cxx.

References fBeamMomentum, fCCHits, fEndSnarl, fEndTime, frun, fsavefile, fSimFlag, fStartSnarl, fStartTime, fStTrkHits, fsubrun, ftree, fTrkHits, fTrkOpt, fTruthHits, fUnTrkHits, fXTalkHits, Msg::kDebug, Msg::kInfo, Msg::kWarning, and MSG.

{
  MSG("CDTrackerModule",Msg::kDebug)<<"Running EndJob()"<<endl;  
  
  if (fsavefile!=0 && fsavefile->IsWritable()){
    fsavefile->cd();
    
    MSG("CDTrackerModule",Msg::kInfo)
      <<"Writing to file... "<<endl;
    //print where tree will be written
    fsavefile->pwd();
    ftree->Write();
    
    MSG("CDTrackerModule",Msg::kDebug)
      <<"Creating options tree..."<<endl;
    TTree opt_tree("TrackerOptions","Tracker Options");
    opt_tree.Branch("SimFlag",&fSimFlag,"SimFlag/I",32000);
    opt_tree.Branch("StartTime",&fStartTime,"StartTime/I",32000);
    opt_tree.Branch("EndTime",&fEndTime,"EndTime/I",32000);
    opt_tree.Branch("StartSnarl",&fStartSnarl,"StartSnarl/I",32000);
    opt_tree.Branch("EndSnarl",&fEndSnarl,"EndSnarl/I",32000);
    opt_tree.Branch("BeamMomentum",&fBeamMomentum,
                    "BeamMomentum/F",32000);
    opt_tree.Branch("Run",&frun,"run/I",32000);
    opt_tree.Branch("SubRun",&fsubrun,"subrun/I",32000);
    opt_tree.Branch("TrackerOptions","CDTrackerOptions",
                     &fTrkOpt,8000,1);
    
    MSG("CDTrackerModule",Msg::kInfo)
      <<"Time="<<fStartTime<<"->"<<fEndTime<<endl
      <<"Snarl="<<fStartSnarl<<"->"<<fEndSnarl<<endl;

    MSG("CDTrackerModule",Msg::kDebug)<<"Filling tree..."<<endl;
    opt_tree.Fill();
    opt_tree.Write();
    MSG("CDTrackerModule",Msg::kInfo)<<"Closing file..."<<endl;
    fsavefile->Close();
  }
  else {
    MSG("CDTrackerModule",Msg::kWarning)
      <<"File not writable!"<<endl;
  }

  delete fTrkHits;
  delete fCCHits;
  delete fStTrkHits;
  delete fXTalkHits;
  delete fUnTrkHits;
  delete fTruthHits;

  MSG("CDTrackerModule",Msg::kDebug)
    <<"End of CDTrackerModule::EndJob() reached"<<endl;
}

void CDTrackerModule::FillPID

(

const CandCalDetPIDHandle *

cdpid

)

[private]

Definition at line 727 of file CDTrackerModule.cxx.

References fPID, CandCalDetPIDHandle::GetInCERTimeBits(), CandCalDetPIDHandle::GetNoOverlapBits(), CandCalDetPIDHandle::GetOLChi2(), CandCalDetPIDHandle::GetPIDType(), CandCalDetPIDHandle::InCERTime(), Msg::kVerbose, MSG, CandCalDetPIDHandle::NoOverlap(), CDPIDInfo::SetInCERTime(), CDPIDInfo::SetInCERTimeBits(), CDPIDInfo::SetNoOverlap(), CDPIDInfo::SetNoOverlapBits(), CDPIDInfo::SetOLChi2(), and CDPIDInfo::SetPIDType().

{
  //check handle exists and fill
  if(cdpid) {
    MSG("CDTrackerModule",Msg::kVerbose)
      <<"Filling CDPIDInfo object with CalDetPID info"<<endl;
    fPID->SetPIDType(cdpid->GetPIDType());
    fPID->SetNoOverlap(cdpid->NoOverlap());
    fPID->SetInCERTime(cdpid->InCERTime());
    fPID->SetNoOverlapBits(cdpid->GetNoOverlapBits());
    fPID->SetInCERTimeBits(cdpid->GetInCERTimeBits());
    fPID->SetOLChi2(cdpid->GetOLChi2());
  }
}

void CDTrackerModule::FillPID

(

CandCalDetSIHandle *

CDSI

)

[private]

Definition at line 744 of file CDTrackerModule.cxx.

References fPID, CandCalDetSIHandle::GetFafErr(), CandCalDetSIHandle::GetKovADC1(), CandCalDetSIHandle::GetKovADC2(), CandCalDetSIHandle::GetKovADC3(), CandCalDetSIHandle::GetKovTimeStamp1(), CandCalDetSIHandle::GetKovTimeStamp2(), CandCalDetSIHandle::GetKovTimeStamp3(), CandCalDetSIHandle::GetSnarlMaxTimeStamp(), CandCalDetSIHandle::GetSnarlMinTimeStamp(), CandCalDetSIHandle::GetSnarlTimeFrame(), CandCalDetSIHandle::GetSparseErr(), CandCalDetSIHandle::GetTickSinceLast(), CandCalDetSIHandle::GetTofADC0(), CandCalDetSIHandle::GetTofADC1(), CandCalDetSIHandle::GetTofADC2(), CandCalDetSIHandle::GetTofADCTimeStamp0(), CandCalDetSIHandle::GetTofADCTimeStamp1(), CandCalDetSIHandle::GetTofADCTimeStamp2(), CandCalDetSIHandle::GetTofTDC0(), CandCalDetSIHandle::GetTofTDC1(), CandCalDetSIHandle::GetTofTDC2(), CandCalDetSIHandle::GetTofTimeStamp(), CandCalDetSIHandle::GetTriggerPMT(), CandCalDetSIHandle::GetTrigSource(), CDPIDInfo::SetFafErr(), CDPIDInfo::SetIsEmpty(), CDPIDInfo::SetKovADC1(), CDPIDInfo::SetKovADC2(), CDPIDInfo::SetKovADC3(), CDPIDInfo::SetKovTimeStamp1(), CDPIDInfo::SetKovTimeStamp2(), CDPIDInfo::SetKovTimeStamp3(), CDPIDInfo::SetSnarlMaxTimeStamp(), CDPIDInfo::SetSnarlMinTimeStamp(), CDPIDInfo::SetSnarlTimeFrame(), CDPIDInfo::SetSparseErr(), CDPIDInfo::SetTickSinceLast(), CDPIDInfo::SetTofADC0(), CDPIDInfo::SetTofADC1(), CDPIDInfo::SetTofADC2(), CDPIDInfo::SetTofADCTimeStamp0(), CDPIDInfo::SetTofADCTimeStamp1(), CDPIDInfo::SetTofADCTimeStamp2(), CDPIDInfo::SetTofTDC0(), CDPIDInfo::SetTofTDC1(), CDPIDInfo::SetTofTDC2(), CDPIDInfo::SetTofTimeStamp(), CDPIDInfo::SetTriggerPMT(), and CDPIDInfo::SetTrigSource().

Referenced by Ana().

{
  //check handle exists and fill
  if(cdsi) {
    fPID->SetIsEmpty(true);
    fPID->SetTriggerPMT(cdsi->GetTriggerPMT());
    fPID->SetFafErr(cdsi->GetFafErr());
    fPID->SetSparseErr(cdsi->GetSparseErr());
    fPID->SetTrigSource(cdsi->GetTrigSource());
    fPID->SetKovADC1(cdsi->GetKovADC1());
    fPID->SetKovADC2(cdsi->GetKovADC2());
    fPID->SetKovADC3(cdsi->GetKovADC3());
    fPID->SetKovTimeStamp1(cdsi->GetKovTimeStamp1());
    fPID->SetKovTimeStamp2(cdsi->GetKovTimeStamp2());
    fPID->SetKovTimeStamp3(cdsi->GetKovTimeStamp3());
    fPID->SetSnarlTimeFrame(cdsi->GetSnarlTimeFrame());
    fPID->SetSnarlMinTimeStamp(cdsi->GetSnarlMinTimeStamp());
    fPID->SetSnarlMaxTimeStamp(cdsi->GetSnarlMaxTimeStamp());
    fPID->SetTofTDC0(cdsi->GetTofTDC0());
    fPID->SetTofTDC1(cdsi->GetTofTDC1());
    fPID->SetTofTDC2(cdsi->GetTofTDC2());    
    fPID->SetTofADC0(cdsi->GetTofADC0());
    fPID->SetTofADC1(cdsi->GetTofADC1());
    fPID->SetTofADC2(cdsi->GetTofADC2());
    fPID->SetTofADCTimeStamp0(cdsi->GetTofADCTimeStamp0());
    fPID->SetTofADCTimeStamp1(cdsi->GetTofADCTimeStamp1());
    fPID->SetTofADCTimeStamp2(cdsi->GetTofADCTimeStamp2());
    fPID->SetTofTimeStamp(cdsi->GetTofTimeStamp());
    fPID->SetTickSinceLast(cdsi->GetTickSinceLast());
  }
}

void CDTrackerModule::HandleCommand

(

JobCommand *

cmd

)

[virtual]

Implement to handle a JobCommand

Reimplemented from JobCModule.

Definition at line 480 of file CDTrackerModule.cxx.

References fJustWriteTracks, fOverwriteFile, ftag, fTrkOpt, fWritePID, fWriteTruth, JobCommand::HaveCmd(), Msg::kWarning, MSG, JobCommand::PopCmd(), JobCommand::PopFloatOpt(), JobCommand::PopIntOpt(), JobCommand::PopOpt(), CDTrackerOptions::SetAcceptance(), CDTrackerOptions::SetAngleLimit(), CDTrackerOptions::SetHitsPerPlane(), CDTrackerOptions::SetMaxHotPlanes(), CDTrackerOptions::SetMinAngle(), CDTrackerOptions::SetMinHit(), CDTrackerOptions::SetMinVertAngle(), CDTrackerOptions::SetThresh(), CDTrackerOptions::SetTrackerOptions(), and CDTrackerOptions::SetTrkLenCut().

{

  if (cmd->HaveCmd()) {                      // If we have a command...
    string sc = cmd->PopCmd();               // Get the command
    
    if (sc == "TrackOpt") {
      string st(cmd->PopOpt());
      fTrkOpt->SetTrackerOptions(st);
      return;
    }

    else if (sc == "ADCThreshold") {
      fTrkOpt->SetThresh(cmd->PopFloatOpt());
      return;
    }

    else if (sc == "MinNHits") {
      fTrkOpt->SetMinHit(cmd->PopIntOpt());
      return;
    }
    
    else if (sc == "MinTrkLen") {
      fTrkOpt->SetTrkLenCut(cmd->PopIntOpt());
      return;
    }
    
    else if (sc == "MaxAngle") {
      fTrkOpt->SetAngleLimit(cmd->PopFloatOpt());
      return;
    }
    
    else if (sc == "MinAngle") {
      fTrkOpt->SetMinAngle(cmd->PopFloatOpt());
      return;
    }

    else if (sc == "MinVertAngle") {
      fTrkOpt->SetMinVertAngle(cmd->PopFloatOpt());
      return;
    }

    else if (sc == "SearchRange") {
      fTrkOpt->SetAcceptance(cmd->PopIntOpt());
      return;
    }

    else if (sc == "HitsPerPlane") {
      fTrkOpt->SetHitsPerPlane(cmd->PopIntOpt());
      return;
    }

    else if (sc == "MaxHotPlanes") {
      fTrkOpt->SetMaxHotPlanes(cmd->PopIntOpt());
      return;
    }

    else if (sc == "FileID") {
      sprintf(ftag,"%s",cmd->PopOpt());
      return;      
    }

    else if (sc == "JustWriteTracks") {
      int tempNum=cmd->PopIntOpt();
      if(tempNum==0)
        fJustWriteTracks=0;
      if(tempNum>0) 
        fJustWriteTracks=1;
      return;      
    }

    else if (sc == "OverwriteFile") {
      int tempNum=cmd->PopIntOpt();
      if(tempNum==0)
        fOverwriteFile=0;
      if(tempNum>0) 
        fOverwriteFile=1;
      return;      
    }

    else if (sc == "WritePID") {
      int tempNum=cmd->PopIntOpt();
      if(tempNum==0)
        fWritePID=0;
      if(tempNum>0) 
        fWritePID=1;
      return;      
    }
    else if (sc == "WriteTruth") {
      int tempNum=cmd->PopIntOpt();
      if(tempNum==0)
        fWriteTruth=0;
      if(tempNum>0) 
        fWriteTruth=1;
      return; 
    }
    else MSG("CDTrackerModule",Msg::kWarning) 
      << "Don't understand '" << sc.c_str() << "'\n";
  }
}

void CDTrackerModule::InitialiseTree

(

)

[private]

Definition at line 703 of file CDTrackerModule.cxx.

References fCCHits, fNumDeadChips, fPID, fSec, fsnarl, fStTrkHits, fTemperature, ftree, fTrkHits, fTrkInfo, fTruthHits, fUnTrkHits, fWritePID, fWriteTruth, and fXTalkHits.

Referenced by Ana().

{
  //create the tree
  ftree=new TTree("TrackerTree","Tracker Tree");
  
  //set tree branches  
  ftree->Branch("Snarl",&fsnarl,"snarl/I",32000);
  ftree->Branch("NumDeadChips",&fNumDeadChips,"NumDeadChips/I",32000);
  ftree->Branch("Sec",&fSec,"Sec/I",32000);
  ftree->Branch("Temperature",&fTemperature,"Temperature/F",32000);
  if(fWritePID==1) ftree->Branch("PIDInfo","CDPIDInfo",&fPID,8000,1);
  ftree->Branch("TrackInfo","CDTrackInfo",&fTrkInfo,8000,1);
  ftree->Branch("TrackedHitInfo","TClonesArray",&fTrkHits,8000,1);
  ftree->Branch("StraightTrackedHitInfo","TClonesArray",
                &fStTrkHits,8000,1);
  ftree->Branch("CCHitInfo","TClonesArray",&fCCHits,8000,1);
  ftree->Branch("XTalkHitInfo","TClonesArray",&fXTalkHits,8000,1);
  ftree->Branch("UnTrackedHitInfo","TClonesArray",&fUnTrkHits,8000,1);
  if (fWriteTruth) ftree->Branch("TruthHitInfo","TClonesArray",
                                 &fTruthHits,8000,1);
}

TFile * CDTrackerModule::OpenFile	(	Int_t	runNumber,
		Int_t	runNumberSub
	)			[private]

Definition at line 639 of file CDTrackerModule.cxx.

References Form(), fOverwriteFile, ftag, Msg::kInfo, Msg::kVerbose, and MSG.

Referenced by Ana().

{
  //create the tfile pointer to be returned
  TFile* outputFile=0;
  
  //get the environmental variable
  char *trackerdir=getenv("TRACKER_DIR");
  
  //use a string to hold env instead 
  string sTrackerDir="";
  
  if (trackerdir!=NULL) {
    sTrackerDir=trackerdir;
  }
  else {
    MSG("CDTrackerModule",Msg::kInfo) 
      <<"Environmental variable $TRACKER_DIR not set." 
      <<" Writing file(s) to current directory"<<endl;
    sTrackerDir=".";
  }
  
  //convert varaibles to string
  string sRunNumber=Form("%d",runNumber);
  string sRunNumberSub=Form("%d",runNumberSub);
  MSG("CDTrackerModule",Msg::kVerbose) 
    <<"srunnumbersub="<<sRunNumberSub<<endl;
  //string sDetector="";
  string sTag=ftag;
  string sBase=sTrackerDir+"/Tracker"+sRunNumber+"_"+
    sRunNumberSub+"_"+sTag;
  string sFileName=sBase+".root";
  
  //test if file already exists
  ifstream Test(sFileName.c_str());
  
  //open the appropriate file
  if(!Test || fOverwriteFile){
    outputFile=new TFile(sFileName.c_str(),"RECREATE");
    if(fOverwriteFile) {
       if(Test) Test.close();
    }
  }
  else {
    //Need new filename
    Int_t fred=1;
    while(Test) {
      Test.close();
      string sAppendage=Form("%d",fred);
      sFileName=sBase+"_"+sAppendage+".root";
      Test.open(sFileName.c_str());
      fred++;
    }
    outputFile = new TFile(sFileName.c_str(),"NEW");
    outputFile->SetCompressionLevel(9);
  }
  
  MSG("CDTrackerModule",Msg::kInfo) 
    <<"Output file opened: "<<sFileName<<endl;
  return outputFile;
}

---

Member Data Documentation

Float_t CDTrackerModule::fBeamMomentum [private]

Definition at line 54 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

TClonesArray* CDTrackerModule::fCCHits [private]

Definition at line 78 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), EndJob(), and InitialiseTree().

Int_t CDTrackerModule::fEndSnarl [private]

Definition at line 63 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

Int_t CDTrackerModule::fEndTime [private]

Definition at line 65 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

bool CDTrackerModule::ffirstrecord [private]

Definition at line 51 of file CDTrackerModule.h.

Referenced by Ana(), and CDTrackerModule().

Int_t CDTrackerModule::fJustWriteTracks [private]

Definition at line 87 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and HandleCommand().

Int_t CDTrackerModule::fNumDeadChips [private]

Definition at line 55 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and InitialiseTree().

Int_t CDTrackerModule::fOverwriteFile [private]

Definition at line 88 of file CDTrackerModule.h.

Referenced by CDTrackerModule(), HandleCommand(), and OpenFile().

CDPIDInfo* CDTrackerModule::fPID [private]

Definition at line 85 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), FillPID(), and InitialiseTree().

Int_t CDTrackerModule::frun [private]

Definition at line 58 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

TFile* CDTrackerModule::fsavefile [private]

Definition at line 68 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

Int_t CDTrackerModule::fSec [private]

Definition at line 61 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and InitialiseTree().

SimFlag::SimFlag_t CDTrackerModule::fSimFlag [private]

Definition at line 57 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

Int_t CDTrackerModule::fsnarl [private]

Definition at line 60 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and InitialiseTree().

Int_t CDTrackerModule::fStartSnarl [private]

Definition at line 62 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

Int_t CDTrackerModule::fStartTime [private]

Definition at line 64 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

TClonesArray* CDTrackerModule::fStTrkHits [private]

Definition at line 76 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), EndJob(), and InitialiseTree().

Int_t CDTrackerModule::fsubrun [private]

Definition at line 59 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and EndJob().

Char_t CDTrackerModule::ftag[256] [private]

Definition at line 66 of file CDTrackerModule.h.

Referenced by CDTrackerModule(), HandleCommand(), and OpenFile().

Float_t CDTrackerModule::fTemperature [private]

Definition at line 56 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), and InitialiseTree().

TTree* CDTrackerModule::ftree [private]

Definition at line 72 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), EndJob(), and InitialiseTree().

TClonesArray* CDTrackerModule::fTrkHits [private]

Definition at line 75 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), EndJob(), and InitialiseTree().

CDTrackInfo* CDTrackerModule::fTrkInfo [private]

Definition at line 74 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), and InitialiseTree().

CDTrackerOptions* CDTrackerModule::fTrkOpt [private]

Definition at line 70 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), EndJob(), and HandleCommand().

TClonesArray* CDTrackerModule::fTruthHits [private]

Definition at line 83 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), EndJob(), and InitialiseTree().

TClonesArray* CDTrackerModule::fUnTrkHits [private]

Definition at line 80 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), EndJob(), and InitialiseTree().

Int_t CDTrackerModule::fWritePID [private]

Definition at line 89 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), HandleCommand(), and InitialiseTree().

Int_t CDTrackerModule::fWriteTruth [private]

Definition at line 90 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), HandleCommand(), and InitialiseTree().

TClonesArray* CDTrackerModule::fXTalkHits [private]

Definition at line 79 of file CDTrackerModule.h.

Referenced by Ana(), CDTrackerModule(), CleanDataMembers(), EndJob(), and InitialiseTree().

---

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

• CDTrackerModule.h
• CDTrackerModule.cxx

---