Add histograms for cluster with origin a Eta, correct the way the histograms dependen...
authorgconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 14 Sep 2011 16:03:51 +0000 (16:03 +0000)
committergconesab <gconesab@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 14 Sep 2011 16:03:51 +0000 (16:03 +0000)
PWG4/PartCorrDep/AliAnaPi0EbE.cxx
PWG4/PartCorrDep/AliAnaPi0EbE.h

index cb1ebf9..529bb6c 100755 (executable)
@@ -65,7 +65,7 @@ AliAnaPi0EbE::AliAnaPi0EbE() :
 {
   //default ctor
   
-  for(Int_t i = 0; i < 5; i++){
+  for(Int_t i = 0; i < 6; i++){
     fhEMCLambda0[i]     = 0;
     fhEMCLambda1[i]     = 0;
     fhEMCDispersion[i]  = 0;
@@ -223,9 +223,9 @@ TList *  AliAnaPi0EbE::GetCreateOutputObjects()
       outputContainer->Add(fhEtaMCNoPi0) ;
       
       if(fAnaType == kIMCalo){
-        TString ptype[] ={"#gamma","#gamma->e^{#pm}","#pi^{0}","e^{#pm}", "hadron"}; 
-        TString pname[] ={"Photon","Conversion",     "Pi0",    "Electron","Hadron"};
-        for(Int_t i = 0; i < 5; i++){ 
+        TString ptype[] ={"#gamma","#gamma->e^{#pm}","#pi^{0}","#eta","e^{#pm}", "hadron"}; 
+        TString pname[] ={"Photon","Conversion",     "Pi0",    "Eta", "Electron","Hadron"};
+        for(Int_t i = 0; i < 6; i++){ 
           
           fhEMCLambda0[i]  = new TH2F(Form("hELambda0_MC%s",pname[i].Data()),
                                       Form("Selected pair, cluster from %s : E vs #lambda_{0}^{2}",ptype[i].Data()),
@@ -516,61 +516,69 @@ void  AliAnaPi0EbE::MakeInvMassInCalorimeter()
           
           if(IsDataMC()) {
             //Photon1
-            if( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCPhoton) && 
-               !GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCConversion) ){
-              fhEMCLambda0[0]    ->Fill(e1, l01);
-              fhEMCLambda1[0]    ->Fill(e1, l11);
-              fhEMCDispersion[0] ->Fill(e1, disp1);
+            if  ( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCPi0)  ){
+              fhEMCLambda0[mcPi0]    ->Fill(e1, l01);
+              fhEMCLambda1[mcPi0]    ->Fill(e1, l11);
+              fhEMCDispersion[mcPi0] ->Fill(e1, disp1);
+            }//pi0
+            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCEta)  ){
+              fhEMCLambda0[mcEta]    ->Fill(e1, l01);
+              fhEMCLambda1[mcEta]    ->Fill(e1, l11);
+              fhEMCDispersion[mcEta] ->Fill(e1, disp1);
+            }//eta          
+            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCPhoton) &&
+                       GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCConversion) ){
+              fhEMCLambda0[mcConversion]    ->Fill(e1, l01);
+              fhEMCLambda1[mcConversion]    ->Fill(e1, l11);
+              fhEMCDispersion[mcConversion] ->Fill(e1, disp1);
+            }//conversion photon
+            else if( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCPhoton) ){
+              fhEMCLambda0[mcPhoton]    ->Fill(e1, l01);
+              fhEMCLambda1[mcPhoton]    ->Fill(e1, l11);
+              fhEMCDispersion[mcPhoton] ->Fill(e1, disp1);
             }//photon   no conversion
             else if  ( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCElectron)){
-              fhEMCLambda0[3]    ->Fill(e1, l01);
-              fhEMCLambda1[3]    ->Fill(e1, l11);
-              fhEMCDispersion[3] ->Fill(e1, disp1);
+              fhEMCLambda0[mcElectron]    ->Fill(e1, l01);
+              fhEMCLambda1[mcElectron]    ->Fill(e1, l11);
+              fhEMCDispersion[mcElectron] ->Fill(e1, disp1);
             }//electron
-            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCConversion) && 
-                       GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCConversion) ){
-              fhEMCLambda0[1]    ->Fill(e1, l01);
-              fhEMCLambda1[1]    ->Fill(e1, l11);
-              fhEMCDispersion[1] ->Fill(e1, disp1);
-            }//conversion photon
-            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag1,AliMCAnalysisUtils::kMCPi0)  ){
-              fhEMCLambda0[2]    ->Fill(e1, l01);
-              fhEMCLambda1[2]    ->Fill(e1, l11);
-              fhEMCDispersion[2] ->Fill(e1, disp1);
-            }//pi0
             else {
-              fhEMCLambda0[4]    ->Fill(e1, l01);
-              fhEMCLambda1[4]    ->Fill(e1, l11);
-              fhEMCDispersion[4] ->Fill(e1, disp1);
+              fhEMCLambda0[mcHadron]    ->Fill(e1, l01);
+              fhEMCLambda1[mcHadron]    ->Fill(e1, l11);
+              fhEMCDispersion[mcHadron] ->Fill(e1, disp1);
             }//other particles 
             
             //Photon 2
-            if( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCPhoton) && 
-               !GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCConversion) ){
-              fhEMCLambda0[0]    ->Fill(e2, l02);
-              fhEMCLambda1[0]    ->Fill(e2, l12);
-              fhEMCDispersion[0] ->Fill(e2, disp2);
+            if  ( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCPi0)  ){
+                 fhEMCLambda0[mcPi0]    ->Fill(e2, l02);
+                 fhEMCLambda1[mcPi0]    ->Fill(e2, l12);
+                 fhEMCDispersion[mcPi0] ->Fill(e2, disp2);
+            }//pi0
+            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCEta)  ){
+              fhEMCLambda0[mcEta]    ->Fill(e2, l02);
+              fhEMCLambda1[mcEta]    ->Fill(e2, l12);
+              fhEMCDispersion[mcEta] ->Fill(e2, disp2);
+            }//eta            
+            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCConversion) && 
+                       GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCPhoton) ){
+                 fhEMCLambda0[mcConversion]    ->Fill(e2, l02);
+                 fhEMCLambda1[mcConversion]    ->Fill(e2, l12);
+                 fhEMCDispersion[mcConversion] ->Fill(e2, disp2);
+            }//conversion photon
+            else if( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCPhoton) ){              
+                fhEMCLambda0[mcPhoton]    ->Fill(e2, l02);
+                fhEMCLambda1[mcPhoton]    ->Fill(e2, l12);
+                fhEMCDispersion[mcPhoton] ->Fill(e2, disp2);
             }//photon   no conversion
             else if  ( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCElectron)){
-              fhEMCLambda0[3]    ->Fill(e2, l02);
-              fhEMCLambda1[3]    ->Fill(e2, l12);
-              fhEMCDispersion[3] ->Fill(e2, disp2);
+              fhEMCLambda0[mcElectron]    ->Fill(e2, l02);
+              fhEMCLambda1[mcElectron]    ->Fill(e2, l12);
+              fhEMCDispersion[mcElectron] ->Fill(e2, disp2);
             }//electron
-            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCConversion) && 
-                       GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCConversion) ){
-              fhEMCLambda0[1]    ->Fill(e2, l02);
-              fhEMCLambda1[1]    ->Fill(e2, l12);
-              fhEMCDispersion[1] ->Fill(e2, disp2);
-            }//conversion photon
-            else if  ( GetMCAnalysisUtils()->CheckTagBit(tag2,AliMCAnalysisUtils::kMCPi0)  ){
-              fhEMCLambda0[2]    ->Fill(e2, l02);
-              fhEMCLambda1[2]    ->Fill(e2, l12);
-              fhEMCDispersion[2] ->Fill(e2, disp2);
-            }//pi0
             else {
-              fhEMCLambda0[4]    ->Fill(e2, l02);
-              fhEMCLambda1[4]    ->Fill(e2, l12);
-              fhEMCDispersion[4] ->Fill(e2, disp2);
+              fhEMCLambda0[mcHadron]    ->Fill(e2, l02);
+              fhEMCLambda1[mcHadron]    ->Fill(e2, l12);
+              fhEMCDispersion[mcHadron] ->Fill(e2, disp2);
             }//other particles 
           }//is datamc
         }//MC histograms
index caf5002..55758c0 100755 (executable)
@@ -76,7 +76,7 @@ class AliAnaPi0EbE : public AliAnaPartCorrBaseClass {
                   fMinDist = m1; fMinDist2 = m2; fMinDist3 = m3                              ; }
 
   //For histograms
-  enum mcTypes   {mcPhoton = 0, mcConversion = 1, mcPi0 = 2, mcElectron = 3, mcOther = 4};
+  enum mcTypes   {mcPhoton = 0, mcConversion = 1, mcPi0 = 2,  mcEta = 3, mcElectron = 4, mcHadron = 5};
 
  private:
   
@@ -107,9 +107,9 @@ class AliAnaPi0EbE : public AliAnaPartCorrBaseClass {
   
   //MC histograms
   
-  TH2F         * fhEMCLambda0[5] ;         //! E vs lambda0 of pi0 pairs but really from MC particle
-  TH2F         * fhEMCLambda1[5] ;         //! E vs lambda1 of pi0 pairs but really from MC particle
-  TH2F         * fhEMCDispersion[5] ;      //! E vs dispersion of pi0 pairs but really from MC particle
+  TH2F         * fhEMCLambda0[6] ;         //! E vs lambda0 of pi0 pairs but really from MC particle
+  TH2F         * fhEMCLambda1[6] ;         //! E vs lambda1 of pi0 pairs but really from MC particle
+  TH2F         * fhEMCDispersion[6] ;      //! E vs dispersion of pi0 pairs but really from MC particle
   
   TH1F         * fhPtMCNoPi0;              //! Number of identified pi0, not coming from pi0
   TH2F         * fhPhiMCNoPi0;             //! Phi of identified pi0, not coming from pi0