CondensedNtpModuleAtm.h

Go to the documentation of this file.

// $Id: CondensedNtpModuleAtm.h,v 1.2 2007/02/12 21:25:37 asousa Exp $
//
// A module used for generating analysis ntuples from the SR ntuples
//
// B. Rebel 10/2004
#ifndef CONDENSEDNTPMODULEATM_H
#define CONDENSEDNTPMODULEATM_H

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

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

#include <string>

#include "AnalysisNtuples/ANtpEventInfo.h"
#include "AnalysisNtuples/ANtpTrackInfoAtm.h"
#include "AnalysisNtuples/ANtpTruthInfoAtm.h"
#include "../ANtpHeaderInfo.h"
#include "AnalysisNtuples/ANtpBeamInfo.h"
#include "AnalysisNtuples/Module/ANtpInfoObjectFillerBeam.h"

class MomNavigator;
class JobCommand;
class NtpSRTrack;
class NtpSRShower;
class NtpSREvent;
class NtpMCTruth;
class NtpMCStdHep;
class NtpSRRecord;
class TClonesArray;

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

class CondensedNtpModuleAtm : public JobCModule, public TObject
{

 public:
  CondensedNtpModuleAtm();
  virtual ~CondensedNtpModuleAtm();
  
  void BeginJob();
  JobCResult Ana(const MomNavigator *mom);
  const Registry& DefaultConfig() const;
  void Config(const Registry& r);
  void Help();
  void EndJob();

private:

  std::string fFileName;
  std::string fRockMapFileName;
  std::string fTreeName;
  std::string fBeamTreeName;
  std::string fBeamPath;
  TFile *fNtpFile;
  TTree *fNtuple;
  TTree *fFailTree;
  TTree *fBeamTree;
  Int_t fDetector; //0 --> near, 1 --> far
  Int_t fGoodTrack;
  Int_t fDataType; //0-->data, 1--> mc
  Int_t fEndPlane;
  Double_t fHornCurrent;
  Double_t fCurrentPOT;  
  Double_t fBeamHPos;
  Double_t fBeamVPos;
  Double_t fTargetPos;
  Double_t fBeamTimeSec;
  Int_t fIsMultiple; //0 --> no, 1 --> yes
  Int_t fValidPlaneFail; //0 --> no, 1 --> yes
  Int_t fFailDeMux; //0 --> no, 1 --> yes
  Int_t fMajorRelease;
  Double_t fPathCosZenDeg[1500];
  Double_t fPathAzimuthDeg[1500];
  Double_t fColDenCosZen[1500];
  Double_t fDensity[1500];
  bool fNewMC;

  ANtpEventInfo *fEventInfo;
  ANtpHeaderInfo *fHeaderInfo;
  ANtpBeamInfo *fBeamInfo;
  ANtpTrackInfoAtm *fTrackInfo;
  ANtpTruthInfoAtm *fTruthInfo;
  
  ANtpInfoObjectFillerBeam *fInfoFiller;
 
  Int_t fCtr;

  void FillEventInformation(ANtpRecoNtpManipulator *ntpManipulator,
                            NtpSREvent *ntpEvent);
  void FillMCInformation(NtpMCTruth *ntpMCTruth, NtpMCStdHep *ntpMCStdHep);
  void FillTrackInformation(NtpSRTrack *ntpTrack, TClonesArray *stripArray, double azimuth,
                            double ra, double dec);
  void FillTrackTimingVariables(NtpSRTrack *ntpTrack, TClonesArray *stripArray);
  void FillTrackMagneticBendingVariables(NtpSRTrack *ntpTrack, TClonesArray *stripArray);
  
  void ResetTreeVariables();
  
  Double_t GetTimeWeight(Float_t ph);
  Double_t ArrivalTime_Weight(Double_t npe) const;
  Double_t ArrivalTime_Uncertainty(Double_t npe) const;
  void WeightedAverage(const Int_t n, const Double_t *x, const Double_t *w,
                       Double_t *parm, Double_t *eparm);
  void LinearFit_Weighted(const Int_t n, const Double_t *x, const Double_t *y, 
                          const Double_t *w,Double_t *parm, Double_t *eparm);
  void FitMinimizingResidual(Int_t &n, Double_t *x, Double_t *y, Double_t *w,
                             Double_t *param, Double_t *eparam, bool findSlope = true);
  void FindLinearFit(Double_t *x, Double_t *y, Double_t *weight, 
                     Int_t nPoints, Double_t &a1, Double_t &a2, 
                     Double_t &chiSq);
  void FindChi2(Double_t *x, Double_t *y, Double_t *weight, 
                Int_t nPoints, Double_t &a1, Double_t &a2, 
                Double_t &chiSq);
  void FindStraightLineDeviation(Double_t *x, Double_t *y, Int_t nPoints, 
                                 Double_t &a1, Double_t &a2, Double_t &deviation);
  void FindMeanAndRMS(double *x, double *weight, int nPoints,
                      double &mean, double &rms);
  Float_t CheapBaseline(float cosz);
  

  ClassDef(CondensedNtpModuleAtm, 3) //CondensedNtpModuleAtm
};

#endif // CONDENSEDNTPMODULEATM_H

---