FairRoot/PandaRoot
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
PndXYResidualTask Class Reference

#include <PndXYResidualTask.h>

Inheritance diagram for PndXYResidualTask:

Public Member Functions

 PndXYResidualTask ()
 
 PndXYResidualTask (const char *name)
 
virtual ~PndXYResidualTask ()
 
virtual InitStatus Init ()
 
virtual InitStatus ReInit ()
 
virtual void Exec (Option_t *opt)
 
virtual void FinishEvent ()
 
virtual void FinishTask ()
 
virtual void SetParContainers ()
 
void SetPersistance (Bool_t p=kTRUE)
 
Bool_t GetPersistance ()
 
TVector3 PropagateToZ (PndTrack *aTrack, Double_t z)
 
void SetBranchName (TString val)
 

Protected Member Functions

 ClassDef (PndXYResidualTask, 2)
 

Protected Attributes

Bool_t fPersistance
 
TClonesArray * fTrackArray
 
TClonesArray * fHitArray
 
TClonesArray * fMissingHitArray
 
TClonesArray * fProjectedHitArray
 
TClonesArray * fGoodHitArray
 
TString fBranchName
 
Int_t fEntryNr
 
TH2D * fHc0c0
 
TH2D * fHc0c1
 
TH2D * fHc0c2
 
TH2D * fHc0c3
 
TH2D * fHc0c0cut
 
TH2D * fHc0c1cut
 
TH2D * fHc0c2cut
 
TH2D * fHc0c3cut
 
TH1D * fHMissingHits
 

Detailed Description

Definition at line 29 of file PndXYResidualTask.h.

Constructor & Destructor Documentation

PndXYResidualTask::PndXYResidualTask ( )
inline

Default constructor

Definition at line 34 of file PndXYResidualTask.h.

References SetVerbose().

34  :
35  FairTask("MvdXYResidualTask"),
36  fPersistance(kTRUE),
37  fTrackArray(0),
38  fHitArray(0),
39  fMissingHitArray(0),
40  fProjectedHitArray(0),
41  fEntryNr(0),
42  fBranchName("MVDHitsPixelCorrectedSorted_event")
43  {
44  SetVerbose(0);
45  }
PndXYResidualTask::fTrackArray
TClonesArray * fTrackArray
Definition: PndXYResidualTask.h:93
PndXYResidualTask::fHitArray
TClonesArray * fHitArray
Definition: PndXYResidualTask.h:94
PndXYResidualTask::fProjectedHitArray
TClonesArray * fProjectedHitArray
Definition: PndXYResidualTask.h:97
PndXYResidualTask::fMissingHitArray
TClonesArray * fMissingHitArray
Definition: PndXYResidualTask.h:96
SetVerbose
drchit SetVerbose(iVerbose)
PndXYResidualTask::PndXYResidualTask ( const char *  name)
inline

Named constructor

Definition at line 48 of file PndXYResidualTask.h.

References SetVerbose().

48  :
49  FairTask(name),
50  fPersistance(kTRUE),
51  fTrackArray(0),
52  fHitArray(0),
53  fMissingHitArray(0),
54  fProjectedHitArray(0),
55  fEntryNr(0),
56  fBranchName("MVDHitsPixelCorrectedSorted_event")
57  {
58  SetVerbose(0);
59  }
PndXYResidualTask::fTrackArray
TClonesArray * fTrackArray
Definition: PndXYResidualTask.h:93
PndXYResidualTask::fHitArray
TClonesArray * fHitArray
Definition: PndXYResidualTask.h:94
PndXYResidualTask::fProjectedHitArray
TClonesArray * fProjectedHitArray
Definition: PndXYResidualTask.h:97
name
TString name
Definition: createRootGeoFile.C:82
PndXYResidualTask::fMissingHitArray
TClonesArray * fMissingHitArray
Definition: PndXYResidualTask.h:96
SetVerbose
drchit SetVerbose(iVerbose)
virtual PndXYResidualTask::~PndXYResidualTask ( )
inlinevirtual

Destructor

Definition at line 63 of file PndXYResidualTask.h.

63  {
64  }

Member Function Documentation

PndXYResidualTask::ClassDef ( PndXYResidualTask  ,
2   
)
protected
void PndXYResidualTask::Exec ( Option_t *  opt)
virtual

Virtual method Exec

Definition at line 71 of file PndXYResidualTask.cxx.

References fEntryNr, fHc0c0, fHc0c0cut, fHc0c1, fHc0c1cut, fHc0c2, fHc0c2cut, fHc0c3, fHc0c3cut, fHitArray, fHMissingHits, fMissingHitArray, fProjectedHitArray, fTrackArray, fVerbose, PndTrack::GetChi2(), PndTrack::GetNDF(), PndSdsHit::GetPosition(), PndSdsHit::GetSensorID(), PndTrackCand::GetSortedHits(), PndTrack::GetTrackCandPtr(), hits, i, PropagateToZ(), res, and SetVerbose().

72 {
73  SetVerbose(2);
74  if (fVerbose > 1 && ++fEntryNr % 1000 == 0) {
75  std::cout << "-I- PndXYResidualTask: Event " << fEntryNr << std::endl;
76  }
77  std::set<int> sensorIDs;
78  fMissingHitArray->Delete();
79  fProjectedHitArray->Delete();
80 
81  for (int i = 0; i < 4; i++)
82  sensorIDs.insert(i);
83 
84  if (fTrackArray->GetEntriesFast() > 0){
85  for (int i = 0; i < fTrackArray->GetEntriesFast(); i++){
86  PndTrack* myTrack = (PndTrack*)fTrackArray->At(i);
87  std::vector<PndTrackCandHit> hits = myTrack->GetTrackCandPtr()->GetSortedHits();
88  std::set<int> tempIDs = sensorIDs;
89  for (int iHits = 0; iHits < hits.size(); iHits++){
90  PndSdsHit* myHit = (PndSdsHit*)fHitArray->At(hits[iHits].GetIndex());
91  //std::cout << "HitValue: " << myHit->GetPosition().X() << "/" << myHit->GetPosition().Y()<< "/" << myHit->GetPosition().Z() << std::endl;
92 
93  TVector3 predicted = PropagateToZ(myTrack, myHit->GetZ());
94  //std::cout << "PredictedValue: " <<predicted.X() << "/" << predicted.Y()<< "/" << predicted.Z() << std::endl;
95  TVector3 res = predicted - myHit->GetPosition();
96  tempIDs.erase(myHit->GetSensorID());
97  if (myHit->GetSensorID() == 0) fHc0c0->Fill(res.X(), res.Y());
98  if (myHit->GetSensorID() == 1) fHc0c1->Fill(res.X(), res.Y());
99  if (myHit->GetSensorID() == 2) fHc0c2->Fill(res.X(), res.Y());
100  if (myHit->GetSensorID() == 3) fHc0c3->Fill(res.X(), res.Y());
101  if (myTrack->GetChi2() / myTrack->GetNDF() < 2){
102  if (myHit->GetSensorID() == 0) fHc0c0cut->Fill(res.X(), res.Y());
103  if (myHit->GetSensorID() == 1) fHc0c1cut->Fill(res.X(), res.Y());
104  if (myHit->GetSensorID() == 2) fHc0c2cut->Fill(res.X(), res.Y());
105  if (myHit->GetSensorID() == 3) fHc0c3cut->Fill(res.X(), res.Y());
106  }
107 // if (myTrack->GetChi2() / myTrack->GetNDF() < 2){
108 // PndSdsHit* goodHit = new ((*fGoodHitArray)[fGoodHitArray->GetEntriesFast()]) PndSdsHit(*myHit);
109 // TVector3 correction;
110 // TVector3 errorMHit;
111 // Int_t j = goodHit->GetSensorID();
112 // if (j == 0) correction.SetXYZ(-1.348e-3, -4.01e-3, 0);
113 // else if (j == 1) correction.SetXYZ(+4.97e-3, 4.978e-3, 0);
114 // else if (j == 2) correction.SetXYZ(-7.21e-3, 0.42e-3, 0);
115 // else if (j == 3) correction.SetXYZ(+3.556e-3, -7.569e-3, 0);
116 // errorMHit = goodHit->GetPosition() + correction;
117 // goodHit->SetPositionError(errorMHit);
118 // }
119 
120  }
121  if (tempIDs.size() > 0){
122  for (std::set<int>::iterator it = tempIDs.begin(); it != tempIDs.end(); it++){
123  fHMissingHits->Fill(*it);
124  int sensID = *it;
125  TVector3 missingHit = PropagateToZ(myTrack, 10 + 6 * sensID);
126  TVector3 errorMHit(0,0,0);
127  new ((*fMissingHitArray)[fMissingHitArray->GetEntriesFast()]) PndSdsHit(-1, *it, missingHit, errorMHit, -1, -1, -1, -1);
128  }
129  }
130  if (myTrack->GetChi2() / myTrack->GetNDF() < 2){
131  for (int j = 0; j < 4; j++){
132  TVector3 projectedHit = PropagateToZ(myTrack, 10 + 6 * j);
133 
134  TVector3 correction;
135  TVector3 errorMHit;
136  if (j == 0) correction.SetXYZ(-1.348e-3, -4.01e-3, 0);
137  else if (j == 1) correction.SetXYZ(+4.97e-3, 4.978e-3, 0);
138  else if (j == 2) correction.SetXYZ(-7.21e-3, 0.42e-3, 0);
139  else if (j == 3) correction.SetXYZ(+3.556e-3, -7.569e-3, 0);
140  errorMHit = projectedHit + correction;
141 
142 
143  new ((*fProjectedHitArray)[fProjectedHitArray->GetEntriesFast()]) PndSdsHit(-1, j, projectedHit, errorMHit, -1, -1, hits.size(), -1);
144 
145  }
146  }
147  }
148 
149  }
150 
151 }
fVerbose
int fVerbose
Definition: poormantracks.C:24
PndXYResidualTask::fTrackArray
TClonesArray * fTrackArray
Definition: PndXYResidualTask.h:93
res
Int_t res
Definition: anadigi.C:166
i
Int_t i
Definition: run_full.C:25
PndSdsHit::GetPosition
TVector3 GetPosition() const
Definition: PndSdsHit.h:93
PndTrackCand::GetSortedHits
std::vector< PndTrackCandHit > GetSortedHits()
Definition: PndTrackCand.cxx:105
PndXYResidualTask::PropagateToZ
TVector3 PropagateToZ(PndTrack *aTrack, Double_t z)
Definition: PndXYResidualTask.cxx:153
PndXYResidualTask::fHitArray
TClonesArray * fHitArray
Definition: PndXYResidualTask.h:94
PndXYResidualTask::fProjectedHitArray
TClonesArray * fProjectedHitArray
Definition: PndXYResidualTask.h:97
PndTrack::GetNDF
Int_t GetNDF() const
Definition: PndTrack.h:35
PndTrack::GetChi2
Double_t GetChi2() const
Definition: PndTrack.h:34
PndXYResidualTask::fMissingHitArray
TClonesArray * fMissingHitArray
Definition: PndXYResidualTask.h:96
hits
CbmHit * hits[nHits]
Definition: RiemannTest.C:19
SetVerbose
drchit SetVerbose(iVerbose)
PndSdsHit::GetSensorID
Int_t GetSensorID() const
Definition: PndSdsHit.h:90
PndTrack::GetTrackCandPtr
PndTrackCand * GetTrackCandPtr()
Definition: PndTrack.h:48
void PndXYResidualTask::FinishEvent ( )
virtual

Definition at line 167 of file PndXYResidualTask.cxx.

168 {
169  // fOutputArray->Delete();
170 }
void PndXYResidualTask::FinishTask ( )
virtual

Definition at line 172 of file PndXYResidualTask.cxx.

References fHc0c0, fHc0c0cut, fHc0c1, fHc0c1cut, fHc0c2, fHc0c2cut, fHc0c3, and fHc0c3cut.

173 {
174  fHc0c0->Write();
175  fHc0c1->Write();
176  fHc0c2->Write();
177  fHc0c3->Write();
178  fHc0c0cut->Write();
179  fHc0c1cut->Write();
180  fHc0c2cut->Write();
181  fHc0c3cut->Write();
182 }
Bool_t PndXYResidualTask::GetPersistance ( )
inline

Definition at line 79 of file PndXYResidualTask.h.

References fPersistance.

79 {return fPersistance;};
InitStatus PndXYResidualTask::Init ( )
virtual

Virtual method Init

Definition at line 35 of file PndXYResidualTask.cxx.

References fBranchName, fGoodHitArray, fHc0c0, fHc0c0cut, fHc0c1, fHc0c1cut, fHc0c2, fHc0c2cut, fHc0c3, fHc0c3cut, fHitArray, fHMissingHits, fMissingHitArray, fProjectedHitArray, and fTrackArray.

36 {
37 
38  FairRootManager* ioman = FairRootManager::Instance();
39  if ( ! ioman ) {
40  std::cout << "-E- PndXYResidualTaskT::Init: "
41  << "RootManager not instantiated!" << std::endl;
42  return kFATAL;
43  }
44 
45  // Create and register output array
46  fHitArray = (TClonesArray*)FairRootManager::Instance()->GetObject(fBranchName);
47  fTrackArray = (TClonesArray*)FairRootManager::Instance()->GetObject("MvdTrack");
48 
49  fMissingHitArray = FairRootManager::Instance()->Register("MvdMissingHit","PndSdsHit", "MVD", kTRUE);
50  fProjectedHitArray = FairRootManager::Instance()->Register("MvdProjectedHit","PndSdsHit", "MVD", kTRUE);
51  fGoodHitArray = FairRootManager::Instance()->Register("GoodHit","PndSdsHit", "MVD", kTRUE);
52 
53  //if(fVerbose>1) { Info("Init","Registering this branch: %s/%s",fFolder.Data(),fOutputBranch.Data()); }
54  //fOutputArray = ioman->Register(fOutputBranch, fInputArray->GetClass()->GetName(), fFolder, fPersistance);
55  fHc0c0 = new TH2D("fHc0c0", "fHc0c0", 200, -0.1, 0.1, 200, -0.1, 0.1);
56  fHc0c1 = new TH2D("fHc0c1", "fHc0c1", 200, -0.1, 0.1, 200, -0.1, 0.1);
57  fHc0c2 = new TH2D("fHc0c2", "fHc0c2", 200, -0.1, 0.1, 200, -0.1, 0.1);
58  fHc0c3 = new TH2D("fHc0c3", "fHc0c3", 200, -0.1, 0.1, 200, -0.1, 0.1);
59  fHc0c0cut = new TH2D("fHc0c0cut", "fHc0c0cut", 400, -0.02, 0.02, 400, -0.02, 0.02);
60  fHc0c1cut = new TH2D("fHc0c1cut", "fHc0c1cut", 400, -0.02, 0.02, 400, -0.02, 0.02);
61  fHc0c2cut = new TH2D("fHc0c2cut", "fHc0c2cut", 400, -0.02, 0.02, 400, -0.02, 0.02);
62  fHc0c3cut = new TH2D("fHc0c3cut", "fHc0c3cut", 400, -0.02, 0.02, 400, -0.02, 0.02);
63  fHMissingHits = new TH1D("fHMIssingHits", "MissingHits", 4, -0.5,3.5);
64 
65 
66  return kSUCCESS;
67 }
PndXYResidualTask::fTrackArray
TClonesArray * fTrackArray
Definition: PndXYResidualTask.h:93
PndXYResidualTask::fHitArray
TClonesArray * fHitArray
Definition: PndXYResidualTask.h:94
PndXYResidualTask::fProjectedHitArray
TClonesArray * fProjectedHitArray
Definition: PndXYResidualTask.h:97
PndXYResidualTask::fMissingHitArray
TClonesArray * fMissingHitArray
Definition: PndXYResidualTask.h:96
PndXYResidualTask::fGoodHitArray
TClonesArray * fGoodHitArray
Definition: PndXYResidualTask.h:98
TVector3 PndXYResidualTask::PropagateToZ ( PndTrack *  aTrack,
Double_t  z 
)

Definition at line 153 of file PndXYResidualTask.cxx.

References Double_t, PndTrack::GetParamFirst(), and t.

Referenced by Exec().

153  {
154  TVector3 origin = aTrack->GetParamFirst().GetOrigin();
155  TVector3 dir = aTrack->GetParamFirst().GetMomentum();
156 
157  //std::cout << "OriginValue: " <<origin.X() << "/" << origin.Y()<< "/" << origin.Z() << std::endl;
158  //std::cout << "DirValue: " <<dir.X() << "/" << dir.Y()<< "/" << dir.Z() << std::endl;
159 
160  Double_t t = (z - origin.Z())/dir.Z();
161 
162  return origin + t * dir;
163 }
Double_t
Double_t
Definition: SimCompleteLinkDef.h:6
z
Double_t z
Definition: createRootFscGeometryFile.C:146
t
TTree * t
Definition: bump_analys.C:13
PndTrack::GetParamFirst
FairTrackParP GetParamFirst()
Definition: PndTrack.h:49
InitStatus PndXYResidualTask::ReInit ( )
virtual

Definition at line 29 of file PndXYResidualTask.cxx.

30 {
31  return kSUCCESS;
32 }
void PndXYResidualTask::SetBranchName ( TString  val)
inline

Definition at line 83 of file PndXYResidualTask.h.

References fBranchName, and val.

83 { fBranchName = val;}
val
Double_t val[nBoxes][nFEBox]
Definition: createCalib.C:11
virtual void PndXYResidualTask::SetParContainers ( )
inlinevirtual

Definition at line 75 of file PndXYResidualTask.h.

75 {};
void PndXYResidualTask::SetPersistance ( Bool_t  p = kTRUE)
inline

Definition at line 78 of file PndXYResidualTask.h.

References fPersistance, and p.

78 {fPersistance=p;};
p
Double_t p
Definition: anasim.C:58

Member Data Documentation

TString PndXYResidualTask::fBranchName
protected

Definition at line 100 of file PndXYResidualTask.h.

Referenced by Init(), and SetBranchName().

Int_t PndXYResidualTask::fEntryNr
protected

Definition at line 102 of file PndXYResidualTask.h.

Referenced by Exec().

TClonesArray* PndXYResidualTask::fGoodHitArray
protected

Definition at line 98 of file PndXYResidualTask.h.

Referenced by Init().

TH2D* PndXYResidualTask::fHc0c0
protected

Definition at line 104 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c0cut
protected

Definition at line 108 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c1
protected

Definition at line 105 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c1cut
protected

Definition at line 109 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c2
protected

Definition at line 106 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c2cut
protected

Definition at line 110 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c3
protected

Definition at line 107 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TH2D* PndXYResidualTask::fHc0c3cut
protected

Definition at line 111 of file PndXYResidualTask.h.

Referenced by Exec(), FinishTask(), and Init().

TClonesArray* PndXYResidualTask::fHitArray
protected

Definition at line 94 of file PndXYResidualTask.h.

Referenced by Exec(), and Init().

TH1D* PndXYResidualTask::fHMissingHits
protected

Definition at line 113 of file PndXYResidualTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndXYResidualTask::fMissingHitArray
protected

Definition at line 96 of file PndXYResidualTask.h.

Referenced by Exec(), and Init().

Bool_t PndXYResidualTask::fPersistance
protected

switch to turn on/off storing the arrays to a file

Definition at line 91 of file PndXYResidualTask.h.

Referenced by GetPersistance(), and SetPersistance().

TClonesArray* PndXYResidualTask::fProjectedHitArray
protected

Definition at line 97 of file PndXYResidualTask.h.

Referenced by Exec(), and Init().

TClonesArray* PndXYResidualTask::fTrackArray
protected

Input array of PndSdsPixelDigis

Definition at line 93 of file PndXYResidualTask.h.

Referenced by Exec(), and Init().

The documentation for this class was generated from the following files:
Generated on Wed Apr 3 2019 10:03:47 for FairRoot/PandaRoot by   doxygen 1.8.5