FairRoot/PandaRoot
Enumerations | Functions
track_check_4pi_tpc.C File Reference

Go to the source code of this file.

Enumerations

enum  DetectorId {
  kDCH, kDRC, kDSK, kEMC,
  kGEM, kLUMI, kMDT, kMVD,
  kRPC, kSTT, kTPC, kTOF,
  kFTS, kHYPG, kHYP, kDCH,
  kDRC, kDSK, kEMC, kGEM,
  kLUMI, kMDT, kMVD, kRPC,
  kSTT, kTPC, kTOF, kFTS,
  kHYPG, kHYP, kDCH, kDRC,
  kDSK, kEMC, kGEM, kLUMI,
  kMDT, kMVD, kRPC, kSTT,
  kTPC, kTOF, kFTS, kHYPG,
  kHYP, kDCH, kDRC, kDSK,
  kEMC, kGEM, kLUMI, kMDT,
  kMVD, kRPC, kSTT, kTPC,
  kTOF, kFTS, kHYPG, kHYP,
  kDCH, kDRC, kDSK, kEMC,
  kGEM, kLUMI, kMDT, kMVD,
  kRPC, kSTT, kTPC, kTOF,
  kFTS, kHYPG, kHYP, kDCH,
  kDRC, kDSK, kEMC, kGEM,
  kLUMI, kMDT, kMVD, kRPC,
  kSTT, kTPC, kTOF, kFTS,
  kHYPG, kHYP, kDCH, kDRC,
  kDSK, kEMC, kGEM, kLUMI,
  kMDT, kMVD, kRPC, kSTT,
  kTPC, kTOF, kFTS, kHYPG,
  kHYP, kRICH, kDRC, kDSK,
  kEMC, kGEM, kLUMI, kMDT,
  kMVD, kRPC, kSTT, kFTOF,
  kTOF, kFTS, kHYPG, kHYP
}
 

Functions

int track_check_4pi_tpc (Int_t nEntries=0)
 

Enumeration Type Documentation

enum DetectorId
Enumerator
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kDCH 

kDCH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kTPC 
kTOF 
kFTS 
kHYPG 
kHYP 
kRICH 

kRICH must be the 1st id, and kHYP must be the last one. Please put new detectors in between!!

kDRC 
kDSK 
kEMC 
kGEM 
kLUMI 
kMDT 
kMVD 
kRPC 
kSTT 
kFTOF 
kTOF 
kFTS 
kHYPG 
kHYP 

Definition at line 2 of file track_check_4pi_tpc.C.

2  {
4  kDCH,kDRC,kDSK,kEMC,kGEM,kLUMI,kMDT,kMVD,kRPC,kSTT,kTPC,kTOF,kFTS,kHYPG,kHYP};
kMVD
Definition: production/tdrct/eta_c/track_check.C:4
kTOF
Definition: production/tdrct/eta_c/track_check.C:4
kHYP
Definition: production/tdrct/eta_c/track_check.C:4
kLUMI
Definition: production/tdrct/eta_c/track_check.C:4
kTPC
Definition: production/tdrct/eta_c/track_check.C:4
kSTT
Definition: production/tdrct/eta_c/track_check.C:4
kDSK
Definition: production/tdrct/eta_c/track_check.C:4
kDRC
Definition: production/tdrct/eta_c/track_check.C:4
kMDT
Definition: production/tdrct/eta_c/track_check.C:4
kRPC
Definition: production/tdrct/eta_c/track_check.C:4
kHYPG
Definition: production/tdrct/eta_c/track_check.C:4
kEMC
Definition: production/tdrct/eta_c/track_check.C:4
kGEM
Definition: production/tdrct/eta_c/track_check.C:4
kFTS
Definition: production/tdrct/eta_c/track_check.C:4
kDCH
Definition: production/tdrct/eta_c/track_check.C:4

Function Documentation

int track_check_4pi_tpc ( Int_t  nEntries = 0)

Definition at line 6 of file track_check_4pi_tpc.C.

References fabs(), PndMCTrack::Get4Momentum(), PndPidCandidate::GetMcIndex(), PndPidCandidate::GetMomentum(), PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), PndPidCandidate::GetMvdHits(), PndMCTrack::GetNPoints(), PndMCTrack::GetPdgCode(), PndPidCandidate::GetSttHits(), inFile, inSimFile, kGEM, kMVD, kSTT, kTPC, mc_array, mctrack, out, reco_mvd(), rootlogon(), tree, and TString.

7 {
8  gROOT->LoadMacro("$VMCWORKDIR/gconfig/rootlogon.C");
9  rootlogon();
10  TString inPidFile = " evt_pid_tpc.root";
11  TString inSimFile = " evt_points_tpc.root";
12 
13  TFile *inFile = TFile::Open(inSimFile,"READ");
14 
15  TTree *tree=(TTree *) inFile->Get("pndsim") ;
16  tree->AddFriend("pndsim",inPidFile);
17 
18  TClonesArray* cand_array=new TClonesArray("PndPidCandidate");
19  tree->SetBranchAddress("PidChargedCand", &cand_array);
20 
21  TClonesArray* mc_array=new TClonesArray("PndMCTrack");
22  tree->SetBranchAddress("MCTrack", &mc_array);
23 
24  TFile *out = TFile::Open("out_test_4pi_tpc.root","RECREATE");
25  TNtuple *nt = new TNtuple("nt","nt","evt:id:mc_p:mc_theta:mc_phi:mc_pid:mc_stt:mc_tpc:mc_mvd:mc_gem:p:theta:phi:mult:stt:tpc:mvd");
26  TNtuple *ntEvt = new TNtuple("ntEvt","ntEvt","evt:acc_stt:acc_tpc:eff:effct:cand:mc_fourpi:reco_fourpi");
27  Float_t mass_k = 0.13957018;
28 
29  if (nEntries==0) nEntries = tree->GetEntriesFast();
30  for (Int_t j=0; j< nEntries; j++){
31  tree->GetEntry(j);
32  //if (cand_array->GetEntriesFast()==0) continue;
33  if ((j%100)==0) cout << "processing event " << j << "\n";
34 
35  Float_t mc_mom = 0, mc_theta = 0, mc_phi = 0;
36  Float_t rec_mom = -1, rec_theta = -1, rec_phi = 0;
37  Int_t stt_mccount = 0, tpc_mccount = 0, reco_stt = 0, reco_tpc = 0,reco_mvd = 0, reco_gem = 0, reco_count = 0, reco_ctcount = 0;
38  TLorentzVector mc_k[4], reco_k[4], mc_fourpi, reco_fourpi;
39  for (Int_t mc = 0; mc < mc_array->GetEntriesFast(); mc++)
40  {
41  PndMCTrack *mctrack = (PndMCTrack*)mc_array->At(mc);
42  if (mctrack->GetMotherID()!=-1) continue;
43  if (mctrack->GetNPoints(kSTT)) stt_mccount++;
44  if (mctrack->GetNPoints(kTPC)) tpc_mccount++;
45  mc_k[mc] = mctrack->Get4Momentum();
46  mc_mom = mctrack->GetMomentum().Mag();
47  mc_theta = mctrack->GetMomentum().Theta()*TMath::RadToDeg();
48  mc_phi = mctrack->GetMomentum().Phi()*TMath::RadToDeg();
49  Int_t mc_pid = mctrack->GetPdgCode();
50  Int_t cand_mult = 0;
51  for (Int_t pp=0; pp<cand_array->GetEntriesFast(); pp++)
52  {
53  PndPidCandidate *pidCand = (PndPidCandidate*)cand_array->At(pp);
54  if (pidCand->GetMcIndex()!=mc) continue;
55  if ( (cand_mult==0) || ((cand_mult>0) && (fabs(rec_mom-mc_mom)> fabs(pidCand->GetMomentum().Mag()-mc_mom))) )
56  {
57  rec_mom = pidCand->GetMomentum().Mag();
58  rec_theta = pidCand->GetMomentum().Theta();
59  rec_phi = pidCand->GetMomentum().Phi();
60  reco_k[mc].SetXYZM(pidCand->GetMomentum().X(),pidCand->GetMomentum().Y(),pidCand->GetMomentum().Z(),mass_k);
61  reco_stt = pidCand->GetSttHits();
62  reco_tpc = pidCand->GetTpcHits();
63  reco_mvd = pidCand->GetMvdHits();
64 
65 }
66  cand_mult++;
67 
68  } // end of candidate loop
69 
70  if (cand_mult>0) reco_count++;
71  if ((cand_mult>0) && ((reco_stt>0) || (reco_tpc>0))) reco_ctcount++;
72 
73  Float_t ntuple_nt[] = {
74  j,mc, mc_mom,mc_theta,mc_phi, mc_pid,
75  mctrack->GetNPoints(kSTT), mctrack->GetNPoints(kTPC), mctrack->GetNPoints(kMVD), mctrack->GetNPoints(kGEM),
76  rec_mom, rec_theta*TMath::RadToDeg(), rec_phi*TMath::RadToDeg(), cand_mult, reco_stt, reco_tpc, reco_mvd
77  };
78  nt->Fill(ntuple_nt);
79 
80  } // end of MC loop
81  mc_fourpi = mc_k[0] + mc_k[1]+mc_k[2]+mc_k[3];
82  reco_fourpi = reco_k[0] + reco_k[1]+reco_k[2]+reco_k[3];
83 
84  Float_t ntuple_evt[] = {
85  j, stt_mccount, tpc_mccount, reco_count, reco_ctcount, cand_array->GetEntriesFast(),
86  mc_fourpi.M(), reco_fourpi.M()
87  };
88  ntEvt->Fill(ntuple_evt);
89 
90 
91  } // end of event loop
92 
93  nt->Draw("mc_theta>>hMcTheta(75,0,150)");
94  nt->Draw("mc_theta>>hMcSttTheta(75,0,150)","mc_stt>0");
95  nt->Draw("mc_theta>>hMcSttMvdTheta(75,0,150)","(mc_stt>0)||(mc_mvd>0)");
96  nt->Draw("mc_theta>>hMcTpcTheta(75,0,150)","mc_tpc>0");
97  nt->Draw("mc_theta>>hMcTpcMvdTheta(75,0,150)","(mc_tpc>0)||(mc_tpc>0)");
98  nt->Draw("mc_theta>>hMcMvdTheta(75,0,150)","mc_mvd>0");
99  nt->Draw("mc_theta>>hMcGemTheta(75,0,150)","mc_gem>0");
100  nt->Draw("mc_theta>>hRecoTheta(75,0,150)","mult>0");
101  nt->Draw("mc_theta>>hRecoCtTheta(75,0,150)","mult>0&&(stt>0||tpc>0)");
102  nt->Draw("mc_theta>>hRecoEffTheta(75,0,150)","mult>0");
103  nt->Draw("mc_theta>>hRecoEffCtTheta(75,0,150)","mult>0&&(stt>0||tpc>0)");
104  hRecoEffTheta->Divide(hMcTheta);
105  hRecoEffCtTheta->Divide(hMcTheta);
106  nt->Draw("mc_theta>>hMcAccSttTheta(75,0,150)","mc_stt>0");
107  hMcAccSttTheta->Divide(hMcTheta);
108  nt->Draw("mc_theta>>hMcAccTpcTheta(75,0,150)","mc_tpc>0");
109  hMcAccTpcTheta->Divide(hMcTheta);
110 
111  nt->Draw("mc_phi>>hMcPhi(50,-200,200)");
112  nt->Draw("mc_phi>>hMcSttPhi(50,-200,200)","mc_stt>0");
113  nt->Draw("mc_phi>>hMcSttMvdPhi(50,-200,200)","(mc_stt>0)||(mc_mvd>0)");
114  nt->Draw("mc_phi>>hMcTpcPhi(50,-200,200)","mc_tpc>0");
115  nt->Draw("mc_phi>>hMcTpcMvdPhi(50,-200,200)","(mc_tpc>0)||(mc_mvd>0)");
116  nt->Draw("mc_phi>>hMcMvdPhi(50,-200,200)","mc_mvd>0");
117  nt->Draw("mc_phi>>hMcGemPhi(50,-200,200)","mc_gem>0");
118  nt->Draw("mc_phi>>hRecoPhi(50,-200,200)","mult>0");
119  nt->Draw("mc_phi>>hRecoCtPhi(50,-200,200)","mult>0&&(stt>0||tpc>0)");
120  nt->Draw("mc_phi>>hRecoEffPhi(50,-200,200)","mult>0");
121  nt->Draw("mc_phi>>hRecoEffCtPhi(50,-200,200)","mult>0&&(stt>0||tpc>0)");
122  hRecoEffPhi->Divide(hMcPhi);
123  hRecoEffCtPhi->Divide(hMcPhi);
124  nt->Draw("mc_phi>>hMcAccSttPhi(50,-200,200)","mc_stt>0");
125  hMcAccSttPhi->Divide(hMcPhi);
126  nt->Draw("mc_phi>>hMcAccTpcPhi(50,-200,200)","mc_tpc>0");
127  hMcAccTpcPhi->Divide(hMcPhi);
128 
129  nt->Draw("mc_p>>hMcP(100,0,5)");
130  nt->Draw("mc_p>>hMcSttP(100,0,5)","mc_stt>0");
131  nt->Draw("mc_p>>hMcSttMvdP(100,0,5)","(mc_stt>0)||(mc_mvd>0)");
132  nt->Draw("mc_p>>hMcTpcP(100,0,5)","mc_tpc>0");
133  nt->Draw("mc_p>>hMcTpcMvdP(100,0,5)","(mc_tpc>0)||(mc_mvd>0)");
134  nt->Draw("mc_p>>hMcMvdP(100,0,5)","mc_mvd>0");
135  nt->Draw("mc_p>>hMcGemP(100,0,5)","mc_gem>0");
136  nt->Draw("mc_p>>hRecoP(100,0,5)","mult>0");
137  nt->Draw("mc_p>>hRecoCtP(100,0,5)","mult>0&&(stt>0||tpc>0)");
138  nt->Draw("mc_p>>hRecoEffP(100,0,5)","mult>0");
139  nt->Draw("mc_p>>hRecoEffCtP(100,0,5)","mult>0&&(stt>0||tpc>0)");
140  hRecoEffP->Divide(hMcP);
141  hRecoEffCtP->Divide(hMcP);
142  nt->Draw("mc_p>>hMcAccSttP(100,0,5)","mc_stt>0");
143  hMcAccSttP->Divide(hMcP);
144  nt->Draw("mc_p>>hMcAccTpcP(100,0,5)","mc_tpc>0");
145  hMcAccTpcP->Divide(hMcP);
146 
147  nt->Draw("(mc_p-p):mc_p>>hresp_p(100,0,5,100,-0.5,0.5)","mult>0");
148  nt->Draw("(mc_p-p):mc_theta>>hresp_theta(100,0,100,100,-0.5,0.5)","mult>0","colz");
149  nt->Draw("(mc_theta-theta):mc_theta>>hrestheta_theta(100,0,100,100,-2,2)","mult>0");
150  nt->Draw("(mc_theta-theta):mc_p>>hrestheta_p(100,0,5,100,-2,2)","mult>0");
151  nt->Draw("(mc_phi-phi):mc_phi>>hresphi_phi(100,-200,200,100,-2,2)","mult>0");
152 
153 
154  ntEvt->Draw("reco_fourpi>>hpi(100,2.00,4.00)","eff==4");
155  out->cd();
156 
157  nt->Write();
158  ntEvt->Write();
159 
160  hMcTheta->Write(); hMcSttTheta->Write(); hMcSttMvdTheta->Write(); hMcMvdTheta->Write(); hMcGemTheta->Write(); hRecoTheta->Write();
161  hMcTpcTheta->Write(); hMcTpcMvdTheta->Write(); hRecoEffTheta->Write(); hMcAccSttTheta->Write(); hMcAccTpcTheta->Write();
162  hMcPhi->Write(); hMcSttPhi->Write(); hMcSttMvdPhi->Write(); hMcMvdPhi->Write(); hMcGemPhi->Write(); hRecoPhi->Write();
163  hMcTpcPhi->Write(); hMcTpcMvdPhi->Write(); hRecoEffPhi->Write(); hMcAccSttPhi->Write(); hMcAccTpcPhi->Write();
164  hMcP->Write(); hMcSttP->Write(); hMcSttMvdP->Write(); hMcMvdP->Write(); hMcGemP->Write(); hRecoP->Write();
165  hMcTpcP->Write(); hMcTpcMvdP->Write(); hRecoEffP->Write();hMcAccSttP->Write(); hMcAccTpcP->Write();
166  hRecoCtTheta->Write();
167  hRecoCtPhi->Write();
168  hRecoCtP->Write();
169  hRecoEffCtTheta->Write();
170  hRecoEffCtPhi->Write();
171  hRecoEffCtP->Write();
172 
173 
174 
175  hresp_p->Write();
176  hresp_theta->Write();
177  hrestheta_theta->Write();
178  hrestheta_p->Write();
179  hresphi_phi->Write();
180 
181 
182  out->Save();
183 
184  return 0;
185 }
mctrack
PndMCTrack * mctrack
Definition: analysis_point_hit_fwendcap.C:250
kMVD
Definition: production/tdrct/eta_c/track_check.C:4
tree
TTree * tree
Definition: plot_dirc.C:12
PndMCTrack::GetNPoints
Int_t GetNPoints(DetectorId detId) const
Definition: PndMCTrack.cxx:120
PndMCTrack::GetPdgCode
Int_t GetPdgCode() const
Definition: PndMCTrack.h:73
kTPC
Definition: production/tdrct/eta_c/track_check.C:4
PndMCTrack::GetMomentum
TVector3 GetMomentum() const
Definition: PndMCTrack.h:78
mc_array
TClonesArray * mc_array
Definition: anaLmdCluster.C:25
PndMCTrack::Get4Momentum
TLorentzVector Get4Momentum() const
Definition: PndMCTrack.cxx:102
PndPidCandidate::GetMcIndex
Int_t GetMcIndex() const
Definition: PndPidCandidate.h:83
kSTT
Definition: production/tdrct/eta_c/track_check.C:4
inFile
TString inFile
Definition: hit_dirc.C:8
reco_mvd
int reco_mvd()
Definition: reco_mvd.C:1
inSimFile
TString inSimFile
Definition: run_reco_llbar_evtgen.C:8
rootlogon
rootlogon()
Definition: outdated/mpiTools/macros/emc/rootlogon.C:1
PndPidCandidate::GetSttHits
Int_t GetSttHits() const
Definition: PndPidCandidate.h:110
fabs
friend F32vec4 fabs(const F32vec4 &a)
Definition: P4_F32vec4.h:47
out
TFile * out
Definition: reco_muo.C:20
PndPidCandidate::GetMvdHits
Int_t GetMvdHits() const
Definition: PndPidCandidate.h:101
kGEM
Definition: production/tdrct/eta_c/track_check.C:4
PndMCTrack::GetMotherID
Int_t GetMotherID() const
Definition: PndMCTrack.h:74
PndPidCandidate::GetMomentum
TVector3 GetMomentum() const
Definition: PndPidCandidate.h:55
Generated on Wed Apr 3 2019 10:03:20 for FairRoot/PandaRoot by   doxygen 1.8.5