DmxDeMuxCosmicsModule.h

Go to the documentation of this file.

// $Id: DmxDeMuxCosmicsModule.h,v 1.55 2003/03/21 17:15:56 rhatcher Exp $
//
// A module used for demultiplexing cosmics ray events in the far detector
//
// Author:  B. Rebel 6/2001
#ifndef DMXDEMUXCOSMICSMODULE_H
#define DMXDEMUXCOSMICSMODULE_H

#include "TObject.h"
#include "DeMux/DmxStatus.h"
#include "TTree.h"
#include "TFile.h"
#include "Registry/Registry.h"

#ifndef JOBCMODULE_H
#include "JobControl/JobCModule.h"
#endif

class MomNavigator;
class JobCommand;
class TObjectItr;
class TFolder;
class CandDeMuxDigitHandleItr;

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

class DmxDeMuxCosmicsModule : public JobCModule, public TObject
{

 public:
  DmxDeMuxCosmicsModule();
  ~DmxDeMuxCosmicsModule();
  
  void BeginJob();
  JobCResult Reco(MomNavigator *mom);
  JobCResult Ana(const MomNavigator *mom);
  void HandleCommand(JobCommand *command);
  const Registry& DefaultConfig() const;
  void Config(const Registry& r);
  void EndJob();
  void Help();
  void FindDigitsChiSq(Double_t &uChiSq, Double_t &vChiSq, CandDeMuxDigitHandleItr &digitItr, 
                       Double_t uIntercept, Double_t uSlope, Double_t vIntercept, 
                       Double_t vSlope);
  void FindPlanesChiSq(Double_t &uChiSq, Double_t &vChiSq, TObjectItr &planeItr, 
                       Double_t uIntercept, Double_t uSlope, Double_t vIntercept, 
                       Double_t vSlope);
private:

  Float_t fAverageTimingOffset;
  Int_t fCrate;
  Int_t fDataType;  //type of data used - real==0 or monte carlo==1
  Int_t fDeltaStrips;  //number of strips a digit can be off of a track
  Int_t fDeMuxedStrip;
  Int_t fDeMuxedStripW;
  Float_t fDigitCharge;
  Float_t fDigitChargeW;
  Int_t fDigitDeMuxedCorrectly;
  Int_t fDigitsInEvent;
  Int_t fDigitSide; //0 = east, 1 = west
  Int_t fDirection;
  bool fDontUseCandDigitMasks;
  Int_t fEndPlane;
  Float_t fEventCharge;
  Int_t fEventDeMuxed;
  bool fEventLength;
  Int_t fEventNumber;
  bool fFilter;
  Float_t fFractionChargeNearTrack;
  const Char_t *fHistoFileName; //name of root file containing histograms from Ana()
  Float_t fHoughSlope;
  Float_t fHoughIntercept;
  Float_t fHoughPeak;
  Float_t fHypothesisCoG;
  Float_t fHypothesisCoGTPos;
  Float_t fHypothesisMatedSignalRatio;
  Int_t fHypothesisNumber;
  Int_t fHypothesisRank;
  Float_t fHypothesisStat;
  Float_t fHypothesisTieBreaker;
  Int_t fHypothesisValid;
  bool fMakeTrees;
  Int_t fMated;
  Double_t fMatedSignalForValid;
  Float_t fMCDigitChargeE;
  Float_t fMCDigitChargeW;  
  Float_t fMCChargeCorrect;
  Float_t fMCChargeCorrect1Hyp;
  Float_t fMCTotalCharge;
  Float_t fMCTotalCharge1Hyp;
  Float_t fMCStripsCorrect;
  Int_t fMCVertexPlane;
  Int_t fMCEndPlane;
  Int_t fMCShowerSpread;
  Int_t fMCStripSpread;
  Int_t fMCStripSpread24Planes;
  Int_t fMultiple;
  Int_t fNonPhysicalFailure;
  Int_t fNonPhysicalPlanes;
  Int_t fNoVertexFailure;
  Int_t fNumberOfPlanesInEvent;
  Int_t fOneSidedUPlanes;
  Int_t fOneSidedVPlanes;
  Int_t fPlane;
  Float_t fPlaneZ;
  Float_t fPlaneDeMuxedCoG;
  Int_t fPlaneStray;
  Int_t fPlanesInWindow;
  Int_t fPlaneType;
  Float_t fSignalFractionLimit;
  bool fSingleModuleDeMux;
  Float_t fSingleModuleCharge;
  Float_t fSingleModuleCorrectCharge;
  Float_t fSingleModuleDigits;
  Float_t fSingleModuleCorrectDigits;
  DmxStatus fStatus;  //object containing info about the event
  Int_t fStrayDeltaStripLimit;
  Int_t fStrip0;
  Int_t fStrip1;
  Int_t fStrip2;
  Int_t fStrip3;
  Int_t fStrip4;
  Int_t fStrip5;
  Int_t fStrip6;
  Int_t fStrip7;
  Double_t fTime;
  Int_t fTruthStripE;
  Int_t fTruthStripW;
  Double_t fUSlope;
  Double_t fUIntercept;
  Double_t fUSlopeRMS;
  Int_t fUStrayPlanes;
  Int_t fUValidPlanes;
  Int_t fUOverlap;
  Double_t fUInterceptRMS;
  Double_t fUChiSq;
  Double_t fVSlope;
  Double_t fVIntercept;
  Double_t fVSlopeRMS;
  Int_t fVStrayPlanes;
  Int_t fVValidPlanes;
  Int_t fVOverlap;
  Double_t fVInterceptRMS;
  Double_t fVChiSq;
  Int_t fXTalkDigit;
  TString fTrackFitter; // type of track fitter to use
  bool fUnsetDigitCheck;
  Int_t fValidPlanesFailure;
  Int_t fVertexPlane;
  Float_t fVertexPlaneZ;
  bool fWriteHistos;
  Double_t fXTalkSignalRatio; //xtalk signal ratio
  bool fUseStrayPlaneCheck;
  Int_t fStrayPlanesLimit;
  
  TFile *fDeMuxFile;

  //TTrees to hold test information
  TTree *fDeMuxTree;  //holds all information that only occurs once per event
                     //ie fraction of correctly demuxed charge, event length, total signal
  TTree *fDeMuxedDigitTree; //tree to hold info about digit - crate, plane, demuxed strip, charge
  TTree *fDmxPlaneTree; //tree to hold info about digit - crate, plane, demuxed strip, charge
  TTree *fSingleModuleTree;  //tree to hold information about digits from the Single Module test
  TTree *fTrackDigitsTree; //tree to hold info about digits from DigitsAlongTrack test
  TTree *fUnusedDigitsTree; //tree to hold info about digits from UnusedDigits test
  TTree *fMCDeMuxTree; //tree to hold mc demuxing info
  TTree *fMCDeMuxDigitTree; //tree to hold mc demuxed digit info
  TFolder *fTFolder;

ClassDef(DmxDeMuxCosmicsModule, 0) //DmxDeMuxCosmicsModule
};

#endif // DMXDEMUXCOSMICSMODULE_H

---