moved definition of a few histograms so that they are created for both real and MC...
authorjklay <jklay@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 3 Mar 2010 23:16:31 +0000 (23:16 +0000)
committerjklay <jklay@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 3 Mar 2010 23:16:31 +0000 (23:16 +0000)
PWG4/PartCorrDep/AliAnaElectron.cxx
PWG4/PartCorrDep/AliAnaElectron.h

index 4b613d3..4f587a4 100755 (executable)
@@ -74,6 +74,9 @@ AliAnaElectron::AliAnaElectron()
   fhPtElectron(0),fhPhiElectron(0),fhEtaElectron(0),\r
   fhPtNPE(0),fhPhiNPE(0),fhEtaNPE(0),\r
   fhPtPE(0),fhPhiPE(0),fhEtaPE(0),\r
+  //for comparisons with tracking detectors\r
+  fhPtTrack(0),fhPtHadron(0),fhPtNPEleTPC(0),fhPtNPEleTPCTRD(0),fhPtNPEleTTE(0),\r
+  fhPtNPEleEMCAL(0),\r
   //DVM B-tagging\r
   fhDVMBtagCut1(0),fhDVMBtagCut2(0),fhDVMBtagCut3(0),fhDVMBtagQA1(0),fhDVMBtagQA2(0),\r
   fhDVMBtagQA3(0),fhDVMBtagQA4(0),fhDVMBtagQA5(0),\r
@@ -89,8 +92,7 @@ AliAnaElectron::AliAnaElectron()
   //reco electrons from various sources\r
   fhPhiConversion(0),fhEtaConversion(0),\r
   //for comparisons with tracking detectors\r
-  fhPtTrack(0),fhPtHadron(0),fhPtNPEleTPC(0),fhPtNPEleTPCTRD(0),fhPtNPEleTTE(0),\r
-  fhPtNPEleEMCAL(0),fhPtNPEBHadron(0),\r
+  fhPtNPEBHadron(0),\r
   //for computing efficiency of B-jet tags\r
   fhBJetPt1x4(0),fhBJetPt2x3(0),fhBJetPt3x2(0),\r
   fhFakeJetPt1x4(0),fhFakeJetPt2x3(0),fhFakeJetPt3x2(0),fhDVMJet(0),\r
@@ -134,6 +136,10 @@ AliAnaElectron::AliAnaElectron(const AliAnaElectron & g)
     fhPtElectron(g.fhPtElectron),fhPhiElectron(g.fhPhiElectron),fhEtaElectron(g.fhEtaElectron),\r
     fhPtNPE(g.fhPtNPE),fhPhiNPE(g.fhPhiNPE),fhEtaNPE(g.fhEtaNPE),\r
     fhPtPE(g.fhPtPE),fhPhiPE(g.fhPhiPE),fhEtaPE(g.fhEtaPE),\r
+    //for comparisons with tracking detectors\r
+    fhPtHadron(g.fhPtHadron),fhPtNPEleTPC(g.fhPtNPEleTPC),\r
+    fhPtNPEleTPCTRD(g.fhPtNPEleTPCTRD),fhPtNPEleTTE(g.fhPtNPEleTTE),\r
+    fhPtNPEleEMCAL(g.fhPtNPEleEMCAL),\r
     //DVM B-tagging\r
     fhDVMBtagCut1(g.fhDVMBtagCut1),fhDVMBtagCut2(g.fhDVMBtagCut2),fhDVMBtagCut3(g.fhDVMBtagCut3),\r
     fhDVMBtagQA1(g.fhDVMBtagQA1),fhDVMBtagQA2(g.fhDVMBtagQA2),\r
@@ -153,9 +159,8 @@ AliAnaElectron::AliAnaElectron(const AliAnaElectron & g)
     //reco electrons from various sources\r
     fhPhiConversion(g.fhPhiConversion),fhEtaConversion(g.fhEtaConversion),\r
     //for comparisons with tracking detectors\r
-    fhPtTrack(g.fhPtTrack),fhPtHadron(g.fhPtHadron),fhPtNPEleTPC(g.fhPtNPEleTPC),\r
-    fhPtNPEleTPCTRD(g.fhPtNPEleTPCTRD),fhPtNPEleTTE(g.fhPtNPEleTTE),\r
-    fhPtNPEleEMCAL(g.fhPtNPEleEMCAL),fhPtNPEBHadron(g.fhPtNPEBHadron),\r
+    fhPtTrack(g.fhPtTrack),\r
+    fhPtNPEBHadron(g.fhPtNPEBHadron),\r
     //for computing efficiency of B-jet tags\r
     fhBJetPt1x4(g.fhBJetPt1x4),fhBJetPt2x3(g.fhBJetPt2x3),\r
     fhBJetPt3x2(g.fhBJetPt3x2),\r
@@ -234,6 +239,10 @@ AliAnaElectron & AliAnaElectron::operator = (const AliAnaElectron & g)
   fhPtPE = g.fhPtPE;\r
   fhPhiPE = g.fhPhiPE;\r
   fhEtaPE = g.fhEtaPE; \r
+  //for comparisons with tracking detectors\r
+  fhPtHadron = g.fhPtHadron; fhPtNPEleTPC = g.fhPtNPEleTPC; \r
+  fhPtNPEleTPCTRD = g.fhPtNPEleTPCTRD; fhPtNPEleTTE = g.fhPtNPEleTTE; \r
+  fhPtNPEleEMCAL = g.fhPtNPEleEMCAL; \r
   //DVM B-tagging\r
   fhDVMBtagCut1 = g.fhDVMBtagCut1;\r
   fhDVMBtagCut2 = g.fhDVMBtagCut2; \r
@@ -269,9 +278,7 @@ AliAnaElectron & AliAnaElectron::operator = (const AliAnaElectron & g)
   fhEtaConversion = g.fhEtaConversion;\r
   //for comparisons with tracking detectors\r
   fhPtTrack = g.fhPtTrack;\r
-  fhPtHadron = g.fhPtHadron; fhPtNPEleTPC = g.fhPtNPEleTPC; \r
-  fhPtNPEleTPCTRD = g.fhPtNPEleTPCTRD; fhPtNPEleTTE = g.fhPtNPEleTTE; \r
-  fhPtNPEleEMCAL = g.fhPtNPEleEMCAL; fhPtNPEBHadron = g.fhPtNPEBHadron;\r
+  fhPtNPEBHadron = g.fhPtNPEBHadron;\r
   //for computing efficiency of B-jet tags\r
   fhBJetPt1x4 = g.fhBJetPt1x4; fhBJetPt2x3 = g.fhBJetPt2x3; \r
   fhBJetPt3x2 = g.fhBJetPt3x2;\r
@@ -396,6 +403,26 @@ TList *  AliAnaElectron::GetCreateOutputObjects()
   outputContainer->Add(fhPhiPE) ; \r
   outputContainer->Add(fhEtaPE) ;\r
 \r
+  //These histograms are mixed REAL/MC:\r
+  //Bins along y-axis are:  \r
+  //0 - unfiltered (filled for both real and MC data) \r
+  //1 - bottom, 2 - charm, 3 - charm from bottom  (MC only)\r
+  //4 - conversion, 5 - Dalitz, 6 - W and Z, 7 - junk/unknown (MC only)\r
+  //8 - misidentified (MC only)\r
+\r
+  //histograms for comparison to tracking detectors\r
+  fhPtHadron = new TH2F("hPtHadron","Charged hadrons w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
+  fhPtNPEleTPC = new TH2F("hPtNPEleTPC","Non-phot. Electrons identified by TPC w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
+  fhPtNPEleTPCTRD = new TH2F("hPtNPEleTPCTRD","Non-phot. Electrons identified by TPC+TRD w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
+  fhPtNPEleTTE = new TH2F("hPtNPEleTTE","Non-phot. Electrons identified by TPC+TRD+EMCAL w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);    \r
+  fhPtNPEleEMCAL = new TH2F("hPtNPEleEMCAL","Non-phot. Electrons identified by EMCAL w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
+  \r
+  outputContainer->Add(fhPtHadron);\r
+  outputContainer->Add(fhPtNPEleTPC);\r
+  outputContainer->Add(fhPtNPEleTPCTRD);\r
+  outputContainer->Add(fhPtNPEleTTE);\r
+  outputContainer->Add(fhPtNPEleEMCAL);\r
+\r
   //B-tagging\r
   fhDVMBtagCut1 = new TH2F("hdvmbtag_cut1","DVM B-tag result cut1", 10,0,10 ,nptbins,ptmin,ptmax);\r
   fhDVMBtagCut2 = new TH2F("hdvmbtag_cut2","DVM B-tag result cut2", 10,0,10 ,nptbins,ptmin,ptmax);\r
@@ -474,19 +501,9 @@ TList *  AliAnaElectron::GetCreateOutputObjects()
 \r
     //histograms for comparison to tracking detectors\r
     fhPtTrack  = new TH2F("hPtTrack","Track w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
-    fhPtHadron = new TH2F("hPtHadron","Charged hadrons w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
-    fhPtNPEleTPC = new TH2F("hPtNPEleTPC","Non-phot. Electrons identified by TPC w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
-    fhPtNPEleTPCTRD = new TH2F("hPtNPEleTPCTRD","Non-phot. Electrons identified by TPC+TRD w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
-    fhPtNPEleTTE = new TH2F("hPtNPEleTTE","Non-phot. Electrons identified by TPC+TRD+EMCAL w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);    \r
-    fhPtNPEleEMCAL = new TH2F("hPtNPEleEMCAL","Non-phot. Electrons identified by EMCAL w/in EMCAL acceptance",nptbins,ptmin,ptmax,10,0,10);\r
     fhPtNPEBHadron = new TH2F("hPtNPEBHadron","Non-phot. b-electrons (TPC+TRD+EMCAL) vs B-hadron pt w/in EMCAL acceptance",nptbins,ptmin,ptmax,nptbins,ptmin,ptmax);\r
 \r
     outputContainer->Add(fhPtTrack);\r
-    outputContainer->Add(fhPtHadron);\r
-    outputContainer->Add(fhPtNPEleTPC);\r
-    outputContainer->Add(fhPtNPEleTPCTRD);\r
-    outputContainer->Add(fhPtNPEleTTE);\r
-    outputContainer->Add(fhPtNPEleEMCAL);\r
     outputContainer->Add(fhPtNPEBHadron);\r
 \r
     //for computing efficiency of IPSig tag\r
index e1c2b75..8a87d0d 100755 (executable)
@@ -202,6 +202,23 @@ public:
   TH2F * fhPhiPE; //! Azimuthal angle of photonic electron vs transverse momentum \r
   TH2F * fhEtaPE; //! Pseudorapidity of photonic electron vs tranvserse momentum \r
 \r
+  //These next set do use some MC info.  The first bin of the second\r
+  //dimension is filled for both REAL and MC data, other bins filled\r
+  //only if MC\r
+  //Histograms for comparison to tracking detectors\r
+  TH2F* fhPtHadron;        //!Pt distribution of reco charged hadrons\r
+                           //!(pi,k,p) in EMCAL acceptance\r
+  TH2F* fhPtNPEleTPC;      //!Pt distribution of non-photonic reco electrons using\r
+                          //!just TPC dEdx info in EMCAL acceptance\r
+  TH2F* fhPtNPEleTPCTRD;   //!Pt distribution of non-photonic reco electrons using\r
+                          //!pid info from tracking detectors only in EMCAL acceptance\r
+  TH2F* fhPtNPEleTTE;      //!Pt distribution of non-photonic reco\r
+                          //!electrons using pid info from TPC+TRD+EMCAL\r
+                          //!in EMCAL acceptance\r
+  TH2F* fhPtNPEleEMCAL;    //!Pt distribution of non-photonic reco\r
+                          //!electrons using EMCAL only\r
+                          //!in EMCAL acceptance\r
+\r
   //DVM B-tagging\r
   TH2F * fhDVMBtagCut1; //! DVM B-tagging result for cut1 (minv>1.0)\r
   TH2F * fhDVMBtagCut2; //! DVM B-tagging result for cut2 (minv>1.5)\r
@@ -241,18 +258,6 @@ public:
 \r
   //Histograms for comparison to tracking detectors\r
   TH2F* fhPtTrack;         //!Pt distribution of reco tracks with MC-ID\r
-  TH2F* fhPtHadron;        //!Pt distribution of reco charged hadrons\r
-                           //!(pi,k,p) in EMCAL acceptance\r
-  TH2F* fhPtNPEleTPC;      //!Pt distribution of non-photonic reco electrons using\r
-                          //!just TPC dEdx info in EMCAL acceptance\r
-  TH2F* fhPtNPEleTPCTRD;   //!Pt distribution of non-photonic reco electrons using\r
-                          //!pid info from tracking detectors only in EMCAL acceptance\r
-  TH2F* fhPtNPEleTTE;      //!Pt distribution of non-photonic reco\r
-                          //!electrons using pid info from TPC+TRD+EMCAL\r
-                          //!in EMCAL acceptance\r
-  TH2F* fhPtNPEleEMCAL;    //!Pt distribution of non-photonic reco\r
-                          //!electrons using EMCAL only\r
-                          //!in EMCAL acceptance\r
 \r
   TH2F* fhPtNPEBHadron;    //!correlate our best reconstructed\r
                           //b-electrons with the b-hadron momentum\r