FairRoot/PandaRoot
Functions
run_ana_eta_c_stt_v2.C File Reference

Go to the source code of this file.

Functions

int run_ana_eta_c_stt_v2 (int nevts=0, bool usePID=true)
 

Function Documentation

int run_ana_eta_c_stt_v2 ( int  nevts = 0,
bool  usePID = true 
)

Definition at line 6 of file run_ana_eta_c_stt_v2.C.

References ctime, Double_t, PndAnalysis::FillList(), RhoFitterBase::Fit(), Rho4CFitter::FitConserveMasses(), fRun, RhoFitterBase::GetChi2(), PndAnalysis::GetEntries(), PndAnalysis::GetEvent(), h_chi2_4c, h_chi2_prefit, h_chi2_vtx, h_chi2b_4c, h_chi2b_vtx, h_etac_4c, h_etac_nocut, h_etac_phimass_4c, h_etac_phimass_vtx, h_etac_phimass_vtx_2, h_etac_phimassfit, h_etac_pid, h_etac_vtx, h_mphi_4c, h_mphi_final_4c, h_mphi_final_massfit, h_mphi_final_vtx, h_mphi_final_vtx_2, h_mphi_nocuts, h_mphi_pid, h_mphi_vtx, h_mphi_vtx_2, hvpos, hvtxresX, hvtxresY, hvtxresZ, hvzpos, i, kDCH, kDRC, kDSK, kEMC, kFTS, kGEM, kHYP, kHYPG, kLUMI, kMDT, kMVD, kRPC, kSTT, kTOF, kTPC, Mass, n_etac_4c, n_etac_vtx, n_etac_vtx_2, n_events, nc, out, outFile, p1, p2, parFile, printf(), recoFile, rtdb, rtime, simFile, timer, and TString.

7 {
8  TString simFile = "evt_points_stt.root";
9  TString recoFile = "evt_reco_stt.root";
10  TString pidFile = "evt_pid_stt.root";
11  TString parFile = "evt_params_stt.root";
12  TString outFile = "etac_histo_stt.root";
13 
14  enum DetectorId {kDCH,kDRC,kDSK,kEMC,kGEM,kLUMI,kMDT,
15  kMVD,kRPC,kSTT,kTPC,kTOF,kFTS,kHYPG,kHYP};
16  gStyle->SetOptFit(1011);
17 
18  TStopwatch timer;
19  timer.Start();
20 
21  gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
22 
24  TH1F *h_etac_nocut=new TH1F("h_etac_nocut","#eta_{c}m(#phi,#phi), (no cuts);E, GeV",100,2.5,3.5);
25  TH1F *h_mphi_nocuts=new TH1F("h_mphi_nocuts","#phi: m(K+ K-) (no cuts);E, GeV",100,0.95,1.5);
26 
27  TH1F *h_etac_pid=new TH1F("h_etac_pid","#eta_{c}m(#phi,#phi), (MC PID);E, GeV",100,2.5,3.5);
28  TH1F *h_mphi_pid=new TH1F("h_mphi_pid","#phi: m(K+ K-) (MC PID);E, GeV",100,0.95,1.5);
29 
30  TH1F *h_etac_4c=new TH1F("h_etac_4c","#eta_{c}m(#phi,#phi), 4C-fit;E, GeV",100,2.5,3.5);
31  TH1F *h_etac_4c_refit=new TH1F("h_etac_4c_refit","#eta_{c}m(#phi,#phi), 4C-fit;E, GeV",100,2.5,3.5);
32  TH1F *h_mphi_4c=new TH1F("h_mphi_4c","#phi: m(K+ K-) (4C-fit);E, GeV",100,0.95,1.5);
33 
34  TH1F *h_etac_vtx=new TH1F("h_etac_vtx","#eta_{c}m(#phi,#phi), Vertex fit;E, GeV",100,2.5,3.5);
35  TH1F *h_mphi_vtx=new TH1F("h_mphi_vtx","#phi: m(K+ K-) (Vertex fit);E, GeV",100,0.95,1.5);
36 
37  TH1F *h_etac_vtx_2=new TH1F("h_etac_vtx_2","#eta_{c}m(#phi,#phi), Vertex fit;E, GeV",100,2.5,3.5);
38  TH1F *h_mphi_vtx_2=new TH1F("h_mphi_vtx_2","#phi: m(K+ K-) (Vertex fit);E, GeV",100,0.95,1.5);
39 
40  TH1F *h_etac_phimass_4c=new TH1F("h_etac_phimass_4c","#eta_{c}m(#phi,#phi), (cut on #phi mass);E,
41  GeV",100,2.8,3.2);
42  TH1F *h_mphi_final_4c=new TH1F("h_mphi_final_4c","#phi: m(K+ K-);E, GeV",100,0.95,1.1);
43 
44  TH1F *h_etac_phimass_vtx=new TH1F("h_etac_phimass_vtx","#eta_{c}m(#phi,#phi), (cut on #phi mass);E,
45  GeV",100,2.8,3.2);
46  TH1F *h_mphi_final_vtx=new TH1F("h_mphi_final_vtx","#phi: m(K+ K-);E, GeV",100,0.95,1.1);
47 
48  TH1F *h_etac_phimass_vtx_2=new TH1F("h_etac_phimass_vtx_2","#eta_{c}m(#phi,#phi), (cut on #phi mass);E,
49  GeV",100,2.8,3.2);
50  TH1F *h_mphi_final_vtx_2=new TH1F("h_mphi_final_vtx_2","#phi: m(K+ K-);E, GeV",100,0.95,1.1);
51 
52  TH1F *h_etac_phimassfit=new TH1F("h_etac_phimassfit","#eta_{c}m(#phi,#phi);E, GeV",100,2.8,3.2);
53  TH1F *h_mphi_final_massfit=new TH1F("h_mphi_final_massfit","#phi: m(K+ K-);E, GeV",100,0.95,1.1);
54 
55  TH1F *nc=new TH1F("nc","n charged",20,0,20);
56 
57  TH1F *h_chi2_4c=new TH1F("h_chi2_4c","#chi^{2} 4C-fit;#chi^{2}/N_{df}",100,0,100);
58  TH1F *h_chi2b_4c=new TH1F("h_chi2b_4c","#chi^{2} 4C-fit;#chi^{2}/N_{df}",100,0,100);
59  TH1F *h_chi2_vtx=new TH1F("h_chi2_vtx","#chi^{2} vertex;#chi^{2}/N_{df}",100,0,100);
60  TH1F *h_chi2b_vtx=new TH1F("h_chi2b_vtx","#chi^{2} vertex;#chi^{2}/N_{df}",100,0,100);
61 
62  TH1F *h_chi2_prefit=new TH1F("h_chi2_prefit","#chi^{2} prefit;#chi^{2}",100,0,100);
63 
64  TH1F *h_chi2_mass=new TH1F("h_chi2_mass","#chi^{2} Mass constraint fit;#chi^{2}",100,0,100);
65 
66  TH2F *hvpos = new TH2F("hvpos","(x,y) projection of fitted decay
67  vertex",100,-5,5,100,-5,5);
68  TH1F *hvzpos = new TH1F("hvzpos","z position of fitted decay vertex",100,-10,10);
69  TH1F *hvtxresX = new TH1F("hvtxresX","X resolution of fitted decay
70  vertex",100,-0.1,0.1);
71  TH1F *hvtxresY = new TH1F("hvtxresY","Y resolution of fitted decay
72  vertex",100,-0.1,0.1);
73  TH1F *hvtxresZ = new TH1F("hvtxresZ","Z resolution of fitted decay
74  vertex",100,-0.1,0.1);
75 
76  TH2F *h_theta_p = new TH2F("h_theta_p","Theta vs p",100,0,180,100,0.,3);
77  TH2F *h_theta_p_nr = new TH2F("h_theta_p_nr","Theta vs p",100,0,180,100,0.,3); // not reconstructed
78 
79  TH2F *h_dp_p = new TH2F("h_dp_p","Delta p vs p",100,-3.,3,100,0,3);
80  TH2F *h_dp_theta = new TH2F("h_dp_theta","Delta p vs theta",100,-3.,3,100,0,180);
81  TH1F *h_dp=new TH1F("h_dp","Delta p",100,-3.,3.);
82  TH1F *h_dp_low=new TH1F("h_dp_low","Delta p",100,-3.,3);
83  TH1F *h_dp_high=new TH1F("h_dp_high","Delta p",100,-3.,3);
84 
85  // Number of events in file and number of reconstructed eta_c to store in root file
86  TH1F *n_events=new TH1F("n_events","total number of events",1,0,1);
87  TH1F *n_etac_4c=new TH1F("n_etac_4c","number of reconstructed eta_c (4C-fit)",1,0,1);
88  TH1F *n_etac_vtx=new TH1F("n_etac_vtx","number of reconstructed eta_c (Vertex fit)",1,0,1);
89  TH1F *n_etac_vtx_2=new TH1F("n_etac_vtx_2","number of reconstructed eta_c (Vertex fit)",1,0,1);
90 
92 
93  // Set up the Analysis
94  FairLogger::GetLogger()->SetLogToFile(kFALSE);
95  //FairLogger::GetLogger()->SetLogScreenLevel("DEBUG4");
96  //FairLogger::GetLogger()->SetLogVerbosityLevel("MEDIUM");
97  FairRunAna* fRun = new FairRunAna();
98  FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
99 
100  fRun->SetInputFile(simFile);
101  fRun->AddFriend(recoFile);
102  fRun->AddFriend(pidFile);
103  FairParRootFileIo* parIO = new FairParRootFileIo();
104  parIO->open(parFile.Data());
105  rtdb->setFirstInput(parIO);
106  rtdb->setOutput(parIO);
107 
108  fRun->SetOutputFile(outFile.Data());
109  fRun->Init();
110 
111  TFile *out = TFile::Open(outFile,"RECREATE");
112 
113  PndAnalysis* theAnalysis = new PndAnalysis();
114 
115  // the candidates lists we need
116  TCandList p1, p2, phi1, phi1_pid, phi1_massfit, phi1_massfit_cut;
117  TCandList etac, etac_pid, etac_nocut, etac_massfit, etac_massfit_cut;
118  TCandList etac_vtx;
119  TCandList mctracks;
120 
121  int evts = theAnalysis->GetEntries();
122  if (nevts>0 && nevts<evts) evts=nevts;
123  n_events->SetBinContent(1,evts);
124 
125  TPidMassSelector *phiMassSel=new TPidMassSelector("phi",1.02,0.2);
126  TPidPlusSelector *kplusSel=new TPidPlusSelector("kplus");
127  TPidMinusSelector *kminusSel=new TPidMinusSelector("kminus");
128 
129  int n_reco_4c=0, n_reco_vtx=0, n_reco_vtx_2=0;
130 
131  TLorentzVector ini(0,0,3.6772,4.7333);
132 
133  int i=0,j=0, k=0, l=0;
134 
135  Double_t mc_mom, mc_theta, rec_mom, rec_theta, delta_p;
136  int ievt;
137  while (theAnalysis->GetEvent() && ievt++<evts)
138  {
139  etac_vtx.Cleanup();
140  phi1_massfit.Cleanup();
141  phi1_massfit_cut.Cleanup();
142  if ((ievt%100)==0) cout<<"evt " << ievt << "\n";
143 
144  theAnalysis->FillList(mctracks,"McTruth");
145  theAnalysis->FillList(p1,"Charged");
146  theAnalysis->FillList(p2,"Charged");
147 
148  p1.Select(kplusSel);
149  p2.Select(kminusSel);
150 
151  int n_removed=0;
152  int ii=0;
153  for (Int_t l=0;l<p1.GetLength();l++){
154  ii=l-n_removed;
155  if((p1[ii].GetMicroCandidate().GetSttHits())==0){
156  p1.Remove(p1[ii]);
157  n_removed++;
158  }
159  }
160 
161  int n_removed=0;
162  int ii=0;
163  for (Int_t l=0;l<p2.GetLength();l++){
164  ii=l-n_removed;
165  if((p2[ii].GetMicroCandidate().GetSttHits())==0){
166  p2.Remove(p2[ii]);
167  n_removed++;
168  }
169  }
170 
171  int nchrg=p1.GetLength()+p2.GetLength();
172  nc->Fill(nchrg);
173 
174  for (j=0;j<p1.GetLength();++j) {
175  p1[j].SetMass(TRho::Instance()->GetPDG()->GetParticle(321)->Mass());
176  }
177  for (j=0;j<p2.GetLength();++j) {
178  p2[j].SetMass(TRho::Instance()->GetPDG()->GetParticle(321)->Mass());
179  }
180 
181  phi1.Combine(p1,p2);
182 
183  for (j=0;j<phi1.GetLength();++j) h_mphi_nocuts->Fill(phi1[j].M());
184 
185  phi1.Select(phiMassSel);
186  etac_nocut.Combine(phi1,phi1);
187 
188  for (l=0;l<etac_nocut.GetLength();++l) {
189  h_etac_nocut->Fill(etac_nocut[l].M());
190  }
191 
192  TVector3 mcVertex, mcD1Vertex, mcD2vertex;
193  FairMCEventHeader* evthead =
194  (FairMCEventHeader*)(FairRootManager::Instance()->GetObject("EvtHeader"));
195  evthead->GetVertex(mcVertex);
196 
197  // MC PID
198  // Leave only kaons in particle lists
199  if (usePID)
200  {
201  int n_removed=0;
202  int ii=0;
203  for (l=0;l<p1.GetLength();++l) {
204  ii=l-n_removed;
205  Int_t mcIndex = p1[ii].GetMicroCandidate().GetMcIndex();
206  if (mcIndex>-1){
207  if ((mctracks[mcIndex].PdgCode()!=321))
208  {
209  p1.Remove(p1[ii]);
210  n_removed++;
211  }
212  }
213  else
214  {
215  std::cout<<"stt h: " << p1[ii].GetMicroCandidate().GetSttHits() << std::endl;
216  std::cout<<"Kaon list 1, element "<<l<<" has no assosiated mcTRack"<<std::endl;
217  p1.Remove(p1[ii]);
218  n_removed++;
219  }
220  }
221  n_removed=0;
222  ii=0;
223  for (l=0;l<p2.GetLength();++l) {
224  ii=l-n_removed;
225  Int_t mcIndex = p2[ii].GetMicroCandidate().GetMcIndex();
226  if (mcIndex>-1){
227  if ((mctracks[mcIndex].PdgCode()!=-321))
228  {
229  p2.Remove(p2[ii]);
230  n_removed++;
231  }
232  }
233  else
234  {
235  std::cout<<"stt h: " << p2[ii].GetMicroCandidate().GetSttHits() << std::endl;
236  std::cout<<"Kaon list 2, element "<<l<<" has no assosiated mcTRack"<<std::endl;
237  p2.Remove(p2[ii]);
238  n_removed++;
239  }
240  }
241  }
242 
243  phi1_pid.Combine(p1,p2);
244 
245  for (j=0;j<phi1_pid.GetLength();++j) h_mphi_pid->Fill(phi1_pid[j].M());
246  phi1_pid.Select(phiMassSel);
247  etac.Combine(phi1_pid,phi1_pid);
248 
249  for (l=0;l<etac.GetLength();++l) {
250  h_etac_pid->Fill(etac[l].M());
251  }
252 
254  int best_i=0;
255  double best_chi2=10000;
256  TCandidate *ccfit = new TCandidate();
257  double m_phi1, m_phi2;
258  for (l=0;l<etac.GetLength();++l) {
259  Rho4CFitter fitter(etac[l],ini);
260  fitter.FitConserveMasses();
261  double chi2=fitter.GetChi2();
262  if (chi2<best_chi2)
263  {
264  best_chi2=chi2;
265  best_i = l;
266  ccfit = (TCandidate*)fitter.FittedCand(etac[l]);
267  TCandidate *phi1best = fitter.FittedCand(*(etac[l].Daughter(0)));
268  TCandidate *phi2best = fitter.FittedCand(*(etac[l].Daughter(1)));
269  TCandidate *k1best = fitter.FittedCand(*(etac[l].Daughter(0)->Daughter(0)));
270  TCandidate *k2best = fitter.FittedCand(*(etac[l].Daughter(0)->Daughter(1)));
271  TCandidate *k3best = fitter.FittedCand(*(etac[l].Daughter(1)->Daughter(0)));
272  TCandidate *k4best = fitter.FittedCand(*(etac[l].Daughter(1)->Daughter(1)));
273  TLorentzVector tlvk1 = k1best->P4();
274  TLorentzVector tlvk2 = k2best->P4();
275  TLorentzVector tlvk3 = k3best->P4();
276  TLorentzVector tlvk4 = k4best->P4();
277  m_phi1= (tlvk1+tlvk2).M();
278  m_phi2= (tlvk3+tlvk4).M();
279  }
280  h_chi2_4c->Fill(chi2/9); // Ndf=3N-3=9
281  }
282 
283  if(/*(best_chi2<270)&&*/(etac.GetLength()!=0))
284  {
285  h_chi2b_4c->Fill(best_chi2/9); // Ndf=3N-3=9
286  h_etac_4c->Fill(etac[best_i].M());
287  h_etac_4c_refit->Fill(ccfit->M());
288  h_mphi_4c->Fill(m_phi1);
289  h_mphi_4c->Fill(m_phi2);
290  h_mphi_final_4c->Fill(m_phi1);
291  h_mphi_final_4c->Fill(m_phi2);
292  if (((m_phi1>1.02-0.02)&&(m_phi1<1.02+0.02))
293  && ((m_phi2>1.02-0.02)&&(m_phi2<1.02+0.02)))
294  {
295  h_etac_phimass_4c->Fill(etac[best_i].M());
296  if ((etac[best_i].M()>2.9)&&(etac[best_i].M()<3.06))
297  n_reco_4c++;
298  }
299  }
300 
302  TCandidate *k1, *k2, *k3, *k4, *phi1tmp, *phi2tmp, *etac_tmp;
303 
304  //Combine 4 kaons directly to candidates
305  for (j=0;j<etac.GetLength();++j)
306  {
307  phi1tmp=etac[j].Daughter(0);
308  phi2tmp=etac[j].Daughter(1);
309 
310  k1=phi1tmp->Daughter(0);
311  k2=phi1tmp->Daughter(1);
312  k3=phi2tmp->Daughter(0);
313  k4=phi2tmp->Daughter(1);
314 
315  etac_tmp=k1->Combine(*k2,*k3,*k4);
316  etac_vtx.Add(*etac_tmp);
317  }
318 
319  for (l=0;l<etac_vtx.GetLength();++l) {
320  h_etac_pid->Fill(etac_vtx[l].M());
321  }
322 
323  int best_i=0;
324  double best_chi2=10000;
325  TCandidate *etacfit_best=0;
326  TCandidate *k1fit_best=0, *k2fit_best=0, *k3fit_best=0, *k4fit_best=0;
327  TCandidate *phi1fit_best, *phi2fit_best;
328  TVector3 bestPos;
329  Float_t etacvtx_mass;
330  for (j=0;j<etac_vtx.GetLength();++j)
331  {
332  RhoKinVtxFitter vtxfitter(etac_vtx[j]); // instantiate a vertex fitter
333  vtxfitter.Fit(); // do the vertex fit
334 
335  TCandidate *etacfit=vtxfitter.FittedCand(etac_vtx[j]); // request the fitted EtaC candidate
336  TVector3 etacVtx=etacfit->Pos(); // and the decay vertex position
337  double chi2_vtx=vtxfitter.GlobalChi2();
338  h_chi2_vtx->Fill(chi2_vtx/5); // Number degree of freedom 2N-3=5
339  // plot mass and vtx x,y projection after fit
340  hvpos->Fill(etacVtx.X(),etacVtx.Y());
341  hvzpos->Fill(etacVtx.Z());
342  if(chi2_vtx<best_chi2)
343  {
344  best_chi2=chi2_vtx;
345  best_i=l;
346  etacfit_best=etacfit;
347  k1fit_best=vtxfitter.FittedCand(*(etacfit_best->Daughter(0)));
348  k2fit_best=vtxfitter.FittedCand(*(etacfit_best->Daughter(1)));
349  k3fit_best=vtxfitter.FittedCand(*(etacfit_best->Daughter(2)));
350  k4fit_best=vtxfitter.FittedCand(*(etacfit_best->Daughter(3)));
351  etacvtx_mass = etacfit_best->M();
352  bestPos = etacfit->Pos();
353  }
354  }
355 
356  if(/*(best_chi2<150)&&*/(etac.GetLength()!=0)&&(k1fit_best!=0)&&(k3fit_best!=0))
357  {
358  //h_etac_vtx->Fill(etacfit_best->M());
359  h_chi2b_vtx->Fill(best_chi2/5);
360  h_etac_vtx->Fill(etacvtx_mass);
361  phi1fit_best=k1fit_best->Combine(*k2fit_best);
362  phi2fit_best=k3fit_best->Combine(*k4fit_best);
363  double m_phi1=phi1fit_best->M();
364  double m_phi2=phi2fit_best->M();
365  h_mphi_vtx->Fill(m_phi1);
366  h_mphi_vtx->Fill(m_phi2);
367  h_mphi_final_vtx->Fill(m_phi1);
368  h_mphi_final_vtx->Fill(m_phi2);
369  hvtxresX->Fill(mcVertex.X()-bestPos.X());
370  hvtxresY->Fill(mcVertex.Y()-bestPos.Y());
371  hvtxresZ->Fill(mcVertex.Z()-bestPos.Z());
372 
373  if (((m_phi1>1.02-0.02)&&(m_phi1<1.02+0.02))&&
374  ((m_phi2>1.02-0.02)&&(m_phi2<1.02+0.02)))
375  {
376  h_etac_phimass_vtx->Fill(etacvtx_mass);
377  if ((etacvtx_mass>2.9)&&(etacvtx_mass<3.06))
378  n_reco_vtx++;
379  }
380 
381  }
385  int best_i=0;
386  double best_chi2_prefit=1e6;
387  TCandidate *etacprefit_best=0;
388  TCandidate *k1prefit_best=0, *k2prefit_best=0, *k3prefit_best=0, *k4prefit_best=0;
389  TCandidate *phi1prefit_best, *phi2prefit_best;
390  Double_t mphi1, mphi2, metac;
391  Double_t sigma_etac=0.0331, sigma_phi=0.00392;
392  Double_t m_etac_pdg=2.98, m_phi_pdg=1.02;
393 
394  for (j=0;j<etac_vtx.GetLength();++j)
395  {
396  metac=etac_vtx[j].M();
397  k1=etac_vtx[j].Daughter(0);
398  k2=etac_vtx[j].Daughter(1);
399  k3=etac_vtx[j].Daughter(2);
400  k4=etac_vtx[j].Daughter(3);
401  phi1prefit_best=k1->Combine(*k2);
402  phi2prefit_best=k3->Combine(*k4);
403  m_phi1=phi1prefit_best->M();
404  m_phi2=phi2prefit_best->M();
405 
406  double chi2_prefit=pow(metac-m_etac_pdg,2)/pow(sigma_etac,2)+
407  pow(m_phi1-m_phi_pdg,2)/pow(sigma_phi,2)+pow(m_phi1-m_phi_pdg,2)/pow(sigma_phi,2);
408  h_chi2_prefit->Fill(chi2_prefit);
409 
410  if(chi2_prefit<best_chi2_prefit)
411  {
412  best_chi2_prefit=chi2_prefit;
413  best_i=j;
414  etacprefit_best=etac_vtx[j];
415  etacvtx_mass=etacprefit_best->M();
416  }
417  }
418 
419  if (etacprefit_best!=0)
420  {
421  RhoKinVtxFitter vtxfitter(*etacprefit_best); // instantiate a vertex fitter
422  vtxfitter.Fit(); // do the vertex fit
423  TCandidate *etacfit=vtxfitter.FittedCand(*etacprefit_best);
424  k1prefit_best=vtxfitter.FittedCand(*(etacfit->Daughter(0)));
425  k2prefit_best=vtxfitter.FittedCand(*(etacfit->Daughter(1)));
426  k3prefit_best=vtxfitter.FittedCand(*(etacfit->Daughter(2)));
427  k4prefit_best=vtxfitter.FittedCand(*(etacfit->Daughter(3)));
428  phi1prefit_best=k1prefit_best->Combine(*k2prefit_best);
429  phi2prefit_best=k3prefit_best->Combine(*k4prefit_best);
430  double m_phi1=phi1prefit_best->M();
431  double m_phi2=phi2prefit_best->M();
432  h_mphi_vtx_2->Fill(m_phi1);
433  h_mphi_vtx_2->Fill(m_phi2);
434  h_mphi_final_vtx_2->Fill(m_phi1);
435  h_mphi_final_vtx_2->Fill(m_phi2);
436  if (((m_phi1>1.02-0.02)&&(m_phi1<1.02+0.02))&&
437  ((m_phi2>1.02-0.02)&&(m_phi2<1.02+0.02)))
438  {
439  h_etac_phimass_vtx_2->Fill(etacvtx_mass);
440  if ((etacvtx_mass>2.9)&&(etacvtx_mass<3.06))
441  n_reco_vtx_2++;
442  }
443  }
444  }
445 
446  std::cout<<"Number of reconstructed eta_c (4C) = "<<n_reco_4c<<std::endl;
447  std::cout<<"Number of reconstructed eta_c (Vertex fit)= "<<n_reco_vtx<<std::endl;
448  std::cout<<"Number of reconstructed eta_c (Vertex fit, preselection)= "<<n_reco_vtx_2<<std::endl;
449  n_etac_4c->SetBinContent(1,n_reco_4c);
450  n_etac_vtx->SetBinContent(1,n_reco_vtx);
451  n_etac_vtx_2->SetBinContent(1,n_reco_vtx_2);
452 
453  out->cd();
454  n_etac_4c->Write();
455  n_etac_vtx->Write();
456  n_etac_vtx_2->Write();
457  n_events->Write();
458  h_etac_nocut->Write();
459  h_etac_pid->Write();
460  h_etac_vtx->Write();
461  h_etac_vtx_2->Write();
462  h_etac_4c->Write();
463  h_etac_4c_refit->Write();
464  h_etac_phimass_4c->Write();
465  h_etac_phimass_vtx->Write();
466  h_etac_phimass_vtx_2->Write();
467  h_etac_phimassfit->Write();
468 
469  h_mphi_nocuts->Write();
470  h_mphi_pid->Write();
471  h_mphi_vtx->Write();
472  h_mphi_vtx_2->Write();
473  h_mphi_4c->Write();
474  h_mphi_final_4c->Write();
475  h_mphi_final_vtx->Write();
476  h_mphi_final_vtx_2->Write();
477  h_mphi_final_massfit->Write();
478 
479  nc->Write();
480 
481  h_chi2_4c->Write();
482  h_chi2b_4c->Write();
483  h_chi2_vtx->Write();
484  h_chi2b_vtx->Write();
485  h_chi2_mass->Write();
486  h_chi2_prefit->Write();
487  hvzpos->Write();
488  hvpos->Write();
489 
490  hvtxresX->Write();
491  hvtxresY->Write();
492  hvtxresZ->Write();
493 
494  out->Save();
495 
496  timer.Stop();
497  Double_t rtime = timer.RealTime();
498  Double_t ctime = timer.CpuTime();
499  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
500 
501  return 0;
502 }
PndAnalysis::GetEntries
Int_t GetEntries()
Definition: PndAnalysis.cxx:471
h_chi2_vtx
TH1F * h_chi2_vtx
Definition: plot_eta_c_stt.C:33
h_mphi_final_vtx_2
TH1F * h_mphi_final_vtx_2
Definition: plot_eta_c_stt.C:23
kMVD
Definition: production/tdrct/eta_c/track_check.C:4
h_mphi_final_massfit
TH1F * h_mphi_final_massfit
Definition: plot_eta_c_tpc.C:21
printf
printf("RealTime=%f seconds, CpuTime=%f seconds\n", rtime, ctime)
h_etac_nocut
TH1F * h_etac_nocut
Definition: plot_eta_c_stt.C:7
hvtxresX
TH1F * hvtxresX
Definition: plot_eta_c_stt.C:38
kTOF
Definition: production/tdrct/eta_c/track_check.C:4
i
Int_t i
Definition: run_full.C:25
kHYP
Definition: production/tdrct/eta_c/track_check.C:4
outFile
TString outFile
Definition: hit_dirc.C:17
h_etac_4c
TH1F * h_etac_4c
Definition: plot_eta_c_stt.C:10
hvzpos
TH1F * hvzpos
Definition: plot_eta_c_stt.C:35
h_chi2_prefit
TH1F * h_chi2_prefit
Definition: plot_eta_c_stt.C:26
kLUMI
Definition: production/tdrct/eta_c/track_check.C:4
PndAnalysis::FillList
Bool_t FillList(RhoCandList &l, TString listkey="All", TString pidTcaNames="", int trackHypothesis=-1)
Definition: PndAnalysis.cxx:319
kTPC
Definition: production/tdrct/eta_c/track_check.C:4
h_etac_phimass_4c
TH1F * h_etac_phimass_4c
Definition: plot_eta_c_stt.C:11
h_mphi_nocuts
TH1F * h_mphi_nocuts
Definition: plot_eta_c_stt.C:14
h_mphi_final_4c
TH1F * h_mphi_final_4c
Definition: plot_eta_c_stt.C:18
n_etac_vtx
TH1F * n_etac_vtx
Definition: plot_eta_c_stt.C:43
h_mphi_4c
TH1F * h_mphi_4c
Definition: plot_eta_c_stt.C:17
kSTT
Definition: production/tdrct/eta_c/track_check.C:4
fRun
FairRunAna * fRun
Definition: hit_dirc.C:58
simFile
TString simFile
Definition: bump_emc.C:11
h_mphi_final_vtx
TH1F * h_mphi_final_vtx
Definition: plot_eta_c_stt.C:19
kDSK
Definition: production/tdrct/eta_c/track_check.C:4
h_etac_vtx
TH1F * h_etac_vtx
Definition: plot_eta_c_stt.C:9
h_etac_phimassfit
TH1F * h_etac_phimassfit
Definition: plot_eta_c_tpc.C:13
Double_t
Double_t
Definition: SimCompleteLinkDef.h:6
kDRC
Definition: production/tdrct/eta_c/track_check.C:4
parFile
TString parFile
Definition: hit_dirc.C:14
h_chi2b_4c
TH1F * h_chi2b_4c
Definition: plot_eta_c_stt.C:31
DetectorId
DetectorId
Definition: production/tdrct/eta_c/track_check.C:2
h_chi2_4c
TH1F * h_chi2_4c
Definition: plot_eta_c_stt.C:30
h_mphi_vtx
TH1F * h_mphi_vtx
Definition: plot_eta_c_stt.C:16
timer
TStopwatch timer
Definition: hit_dirc.C:51
kMDT
Definition: production/tdrct/eta_c/track_check.C:4
h_etac_phimass_vtx
TH1F * h_etac_phimass_vtx
Definition: plot_eta_c_stt.C:12
rtdb
FairRuntimeDb * rtdb
Definition: hit_dirc.C:66
kRPC
Definition: production/tdrct/eta_c/track_check.C:4
kHYPG
Definition: production/tdrct/eta_c/track_check.C:4
out
TFile * out
Definition: reco_muo.C:20
n_events
TH1F * n_events
Definition: plot_eta_c_stt.C:44
h_etac_phimass_vtx_2
TH1F * h_etac_phimass_vtx_2
Definition: plot_eta_c_stt.C:21
ctime
Double_t ctime
Definition: hit_dirc.C:114
kEMC
Definition: production/tdrct/eta_c/track_check.C:4
h_etac_pid
TH1F * h_etac_pid
Definition: plot_eta_c_stt.C:8
p2
TPad * p2
Definition: hist-t7.C:117
h_mphi_vtx_2
TH1F * h_mphi_vtx_2
Definition: plot_eta_c_stt.C:22
recoFile
std::string recoFile
Definition: run_ana_llbar_evtgen.C:26
n_etac_vtx_2
TH1F * n_etac_vtx_2
Definition: plot_eta_c_stt.C:25
hvtxresZ
TH1F * hvtxresZ
Definition: plot_eta_c_stt.C:40
p1
TPad * p1
Definition: hist-t7.C:116
kGEM
Definition: production/tdrct/eta_c/track_check.C:4
PndAnalysis::GetEvent
Int_t GetEvent(Int_t n=-1)
Definition: PndAnalysis.cxx:257
h_chi2b_vtx
TH1F * h_chi2b_vtx
Definition: plot_eta_c_stt.C:34
rtime
Double_t rtime
Definition: hit_dirc.C:113
h_mphi_pid
TH1F * h_mphi_pid
Definition: plot_eta_c_stt.C:15
kFTS
Definition: production/tdrct/eta_c/track_check.C:4
n_etac_4c
TH1F * n_etac_4c
Definition: plot_eta_c_stt.C:42
hvtxresY
TH1F * hvtxresY
Definition: plot_eta_c_stt.C:39
hvpos
TH2F * hvpos
Definition: plot_eta_c_stt.C:36
nc
TH1F * nc
Definition: plot_eta_c_stt.C:28
kDCH
Definition: production/tdrct/eta_c/track_check.C:4
Generated on Sun Aug 18 2019 03:02:49 for FairRoot/PandaRoot by   doxygen 1.8.8