]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ANALYSIS/AliD0toKpi.h
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
Moving lib*.pkg
[u/mrichter/AliRoot.git] / ANALYSIS / AliD0toKpi.h
diff --git a/ANALYSIS/AliD0toKpi.h b/ANALYSIS/AliD0toKpi.h
index fb6800f938b30935585ad66cdb3f1131176b7477..4f5981ea3c4dbc2abd7ac456fa6c77681f0e1c0e 100644 (file)
--- a/ANALYSIS/AliD0toKpi.h
+++ b/ANALYSIS/AliD0toKpi.h
@@ -29,44 +29,44 @@ const Double_t kMPi = 0.13957; // Pi+ mass
 //  *** Pb-Pb dNch/dy=6000 *** 
 //
 //  PIONS
-const Int_t kPiBins_PbPb = 10;
-const Double_t kPiBinWidth_PbPb = 0.250;
-const Double_t kPiTagPi_PbPb[kPiBins_PbPb] = {0.211421,0.652184,0.624421,0.614727,0.610777,0.628015,0.631520,0.630324,0.637551,0.575235};
-const Double_t kPiTagNid_PbPb[kPiBins_PbPb] = {0.788579,0.347816,0.375579,0.385273,0.389223,0.371985,0.368480,0.369676,0.362449,0.424765};
+const Int_t kPiBinsPbPb = 10;
+const Double_t kPiBinWidthPbPb = 0.250;
+const Double_t kPiTagPiPbPb[kPiBinsPbPb] = {0.211421,0.652184,0.624421,0.614727,0.610777,0.628015,0.631520,0.630324,0.637551,0.575235};
+const Double_t kPiTagNidPbPb[kPiBinsPbPb] = {0.788579,0.347816,0.375579,0.385273,0.389223,0.371985,0.368480,0.369676,0.362449,0.424765};
 //  KAONS
-const Int_t kKBins_PbPb = 10;
-const Double_t kKBinWidth_PbPb = 0.250;
-const Double_t kKTagK_PbPb[kKBins_PbPb] = {0.000000,0.101255,0.397662,0.467586,0.517008,0.555023,0.584185,0.519029,0.464117,0.247308};
-const Double_t kKTagPi_PbPb[kKBins_PbPb] = {0.102049,0.289930,0.101930,0.057771,0.040286,0.028567,0.053108,0.094369,0.066302,0.247308};
-const Double_t kKTagNid_PbPb[kKBins_PbPb] = {0.897951,0.608815,0.500408,0.474643,0.442705,0.416410,0.362707,0.386603,0.469580,0.505383};
+const Int_t kKBinsPbPb = 10;
+const Double_t kKBinWidthPbPb = 0.250;
+const Double_t kKTagKPbPb[kKBinsPbPb] = {0.000000,0.101255,0.397662,0.467586,0.517008,0.555023,0.584185,0.519029,0.464117,0.247308};
+const Double_t kKTagPiPbPb[kKBinsPbPb] = {0.102049,0.289930,0.101930,0.057771,0.040286,0.028567,0.053108,0.094369,0.066302,0.247308};
+const Double_t kKTagNidPbPb[kKBinsPbPb] = {0.897951,0.608815,0.500408,0.474643,0.442705,0.416410,0.362707,0.386603,0.469580,0.505383};
 //  PROTONS
-const Int_t kPBins_PbPb = 9;
-const Double_t kPBinWidth_PbPb = 0.500;
-const Double_t kPTagP_PbPb[kPBins_PbPb] = {0.017940,0.350681,0.535286,0.583264,0.562935,0.560524,0.545992,0.598060,0.351245};
-const Double_t kPTagPi_PbPb[kPBins_PbPb] = {0.195955,0.094949,0.039962,0.026039,0.007556,0.016986,0.030333,0.000000,0.000000};
-const Double_t kPTagNid_PbPb[kPBins_PbPb] = {0.786105,0.554370,0.424751,0.390697,0.429508,0.422491,0.423675,0.401940,0.648755};
+const Int_t kPBinsPbPb = 9;
+const Double_t kPBinWidthPbPb = 0.500;
+const Double_t kPTagPPbPb[kPBinsPbPb] = {0.017940,0.350681,0.535286,0.583264,0.562935,0.560524,0.545992,0.598060,0.351245};
+const Double_t kPTagPiPbPb[kPBinsPbPb] = {0.195955,0.094949,0.039962,0.026039,0.007556,0.016986,0.030333,0.000000,0.000000};
+const Double_t kPTagNidPbPb[kPBinsPbPb] = {0.786105,0.554370,0.424751,0.390697,0.429508,0.422491,0.423675,0.401940,0.648755};
 //
 // pp PYTHIA 
 //
 // *** cuts for pp ***
 //
 //  PIONS
-const Int_t kPiBins_pp = 10;
-const Double_t kPiBinWidth_pp = 0.250;
-const Double_t kPiTagPi_pp[kPiBins_pp] = {0.194528,0.447097,0.603364,0.646413,0.647125,0.669157,0.688139,0.682564,0.689910,0.665710};
-const Double_t kPiTagNid_pp[kPiBins_pp] = {0.805472,0.552903,0.396636,0.353587,0.352875,0.330843,0.311861,0.317436,0.310090,0.334290};
+const Int_t kPiBinsPP = 10;
+const Double_t kPiBinWidthPP = 0.250;
+const Double_t kPiTagPiPP[kPiBinsPP] = {0.194528,0.447097,0.603364,0.646413,0.647125,0.669157,0.688139,0.682564,0.689910,0.665710};
+const Double_t kPiTagNidPP[kPiBinsPP] = {0.805472,0.552903,0.396636,0.353587,0.352875,0.330843,0.311861,0.317436,0.310090,0.334290};
 //  KAONS
-const Int_t kKBins_pp = 10;
-const Double_t kKBinWidth_pp = 0.250;
-const Double_t kKTagK_pp[kKBins_pp] = {0.000000,0.173393,0.439690,0.519423,0.587025,0.605372,0.586021,0.650139,0.444444,0.299363};
-const Double_t kKTagPi_pp[kKBins_pp] = {0.000000,0.001495,0.000000,-0.000000,-0.000000,0.000000,0.032258,0.060572,0.101449,0.242038};
-const Double_t kKTagNid_pp[kKBins_pp] = {1.000000,0.825112,0.560310,0.480577,0.412975,0.394628,0.381720,0.289289,0.454106,0.458599};
+const Int_t kKBinsPP = 10;
+const Double_t kKBinWidthPP = 0.250;
+const Double_t kKTagKPP[kKBinsPP] = {0.000000,0.173393,0.439690,0.519423,0.587025,0.605372,0.586021,0.650139,0.444444,0.299363};
+const Double_t kKTagPiPP[kKBinsPP] = {0.000000,0.001495,0.000000,-0.000000,-0.000000,0.000000,0.032258,0.060572,0.101449,0.242038};
+const Double_t kKTagNidPP[kKBinsPP] = {1.000000,0.825112,0.560310,0.480577,0.412975,0.394628,0.381720,0.289289,0.454106,0.458599};
 //  PROTONS
-const Int_t kPBins_pp = 9;
-const Double_t kPBinWidth_pp = 0.500;
-const Double_t kPTagP_pp[kPBins_pp] = {0.029404,0.438640,0.613710,0.665152,0.634961,0.657711,0.703704,0.685714,0.235294};
-const Double_t kPTagPi_pp[kPBins_pp] = {0.000000,0.000000,0.000000,-0.000000,0.000000,0.000000,-0.000000,0.014286,-0.000000};
-const Double_t kPTagNid_pp[kPBins_pp] = {0.970596,0.561360,0.386290,0.334848,0.365039,0.342289,0.296296,0.300000,0.764706};
+const Int_t kPBinsPP = 9;
+const Double_t kPBinWidthPP = 0.500;
+const Double_t kPTagPPP[kPBinsPP] = {0.029404,0.438640,0.613710,0.665152,0.634961,0.657711,0.703704,0.685714,0.235294};
+const Double_t kPTagPiPP[kPBinsPP] = {0.000000,0.000000,0.000000,-0.000000,0.000000,0.000000,-0.000000,0.014286,-0.000000};
+const Double_t kPTagNidPP[kPBinsPP] = {0.970596,0.561360,0.386290,0.334848,0.365039,0.342289,0.296296,0.300000,0.764706};
 
 
 
@@ -83,13 +83,13 @@ class AliD0toKpi : public TObject {
   AliD0toKpi(const AliD0toKpi& d0toKpi);
 
   Double_t Alpha() const { return (Ql(0)-Ql(1))/(Ql(0)+Ql(1)); }
-  void     ApplyPID(TString pidScheme="TOFparam_PbPb");
+  void     ApplyPID(TString pidScheme="TOFparamPbPb");
   Double_t ChildrenRelAngle() const; 
   void     ComputeWgts();
   void     CorrectWgt4BR(Double_t factor);
   Double_t CosPointing() const;
   Double_t CosPointingXY() const;
-  void     CosThetaStar(Double_t&,Double_t&) const;
+  void     CosThetaStar(Double_t &ctsD0,Double_t &ctsD0bar) const;
   Double_t Ct() const {return Length()*kMD0/P();}
   Double_t Energy() const { return TMath::Sqrt(P()*P()+kMD0*kMD0); }
   Double_t Eta() const;
@@ -100,13 +100,13 @@ class AliD0toKpi : public TObject {
   Double_t Getd0Child(Int_t child) const { return fd0[child]; }
   Int_t    GetPdgChild(Int_t child) const { return fPdg[child]; }
   Int_t    GetPdgMum(Int_t child) const {return fMum[child]; }
-  void     GetWgts(Double_t&,Double_t&,TString sample) const;
+  void     GetWgts(Double_t &WgtD0,Double_t &WgtD0bar,TString sample) const;
   void     GetPrimaryVtx(Double_t vtx[3]) const 
     { vtx[0]=fV1x; vtx[1]=fV1y; vtx[2]=fV1z; return; }
   void     GetSecondaryVtx(Double_t vtx[3]) const 
     { vtx[0]=fV2x; vtx[1]=fV2y; vtx[2]=fV2z; return; }
 
-  void     InvMass(Double_t&,Double_t&) const;
+  void     InvMass(Double_t &mD0,Double_t &mD0bar) const;
   Bool_t   IsSignal() const { if(fSignal) return kTRUE; return kFALSE; } 
   Double_t Length() const
     { return TMath::Sqrt((fV1x-fV2x)*(fV1x-fV2x)
@@ -122,7 +122,7 @@ class AliD0toKpi : public TObject {
   Double_t Ql(Int_t child) const;
   Double_t Qt() const;
   Double_t Rapidity() const { return 0.5*TMath::Log((Energy()+Pz())/(Energy()-Pz()+1.e-13)); }
-  Bool_t   Select(const Double_t* cuts,Int_t&,Int_t&) const;
+  Bool_t   Select(const Double_t* cuts,Int_t& okD0,Int_t& okD0bar) const;
   void     SetPrimaryVtx(Double_t vtx[3]) 
     { fV1x=vtx[0]; fV1y=vtx[1]; fV1z=vtx[2]; return; }
   void     SetSignal() { fSignal =  kTRUE; return; }
@@ -132,7 +132,7 @@ class AliD0toKpi : public TObject {
   void     SetPdgCodes(Int_t pdg[2]) {fPdg[0]=pdg[0];fPdg[1]=pdg[1];return;}
   void     SetMumPdgCodes(Int_t mum[2]) {fMum[0]=mum[0];fMum[1]=mum[1];return;}
 
-  void     DrawPIDinTOF(TString pidScheme="TOFparam_PbPb") const;
+  void     DrawPIDinTOF(TString pidScheme="TOFparamPbPb") const;
   Double_t LinearInterpolation(Double_t p,Int_t nBins,Double_t Bin,
                               const Double_t *values) const;
   //  void     SetPtWgts4pp();
@@ -145,15 +145,15 @@ class AliD0toKpi : public TObject {
 
   Int_t fTrkNum[2]; // numbers of the two decay tracks
 
-  Double_t fV1x; //
-  Double_t fV1y; // position of the primary vertex of the event
-  Double_t fV1z; //
-  Double_t fV2x; //
-  Double_t fV2y; // position of the reconstructed secondary vertex
-  Double_t fV2z; //
+  Double_t fV1x; // X-position of the primary vertex of the event
+  Double_t fV1y; // Y-position of the primary vertex of the event
+  Double_t fV1z; // Z-position of the primary vertex of the event
+  Double_t fV2x; // X-position of the reconstructed secondary vertex
+  Double_t fV2y; // Y-position of the reconstructed secondary vertex
+  Double_t fV2z; // Z-position of the reconstructed secondary vertex
   Double_t fDCA; // DCA of the two tracks
 
-  Double_t fPx[2];  // 
+  Double_t fPx[2];  // X,Y,Z-
   Double_t fPy[2];  // momenta of the two tracks
   Double_t fPz[2];  // at the reconstructed vertex  
 
@@ -174,7 +174,7 @@ class AliD0toKpi : public TObject {
   Double_t fPIDrespPr[2]; // det. response to be proton
   Double_t fTOFmass[2]; // mass estimated by the TOF (-1000. if track not reached TOF)
 
-  Double_t fWgtAD0,fWgtAD0bar; //
+  Double_t fWgtAD0,fWgtAD0bar; // weights for the 3 samples
   Double_t fWgtBD0,fWgtBD0bar; // weights for the 3 samples 
   Double_t fWgtCD0,fWgtCD0bar; // A: (K,Pi)+(K,?) B: (?,Pi) C: (?,?)
   Double_t fWgtDD0,fWgtDD0bar; // D: all other pairs