add histogram with shower shape of cluster before the split cluster selection
authorgconesab <gustavo.conesa.balbastre@cern.ch>
Fri, 8 Aug 2014 13:44:38 +0000 (15:44 +0200)
committergconesab <gustavo.conesa.balbastre@cern.ch>
Fri, 8 Aug 2014 13:46:49 +0000 (15:46 +0200)
PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.cxx
PWGGA/CaloTrackCorrelations/AliAnaPi0EbE.h

index c9cbeb5..9b6b399 100755 (executable)
@@ -86,8 +86,8 @@ fhSplitPtEta(0),                    fhSplitPtPhi(0),
 fhNLocMaxSplitPt(0),
 fhPtDecay(0),                       fhEDecay(0),
 // Shower shape histos
-fhPtDispersion(0),                  fhPtLambda0(0),                     fhPtLambda1(0),
-fhPtLambda0NoTRD(0),                fhPtLambda0FracMaxCellCut(0),
+fhPtDispersion(0),                  fhPtLambda0(0),                     fhPtLambda0NoSplitCut(0),
+fhPtLambda1(0),                     fhPtLambda0NoTRD(0),                fhPtLambda0FracMaxCellCut(0),
 fhPtFracMaxCell(0),                 fhPtFracMaxCellNoTRD(0),
 fhPtNCells(0),                      fhPtTime(0),                        fhEPairDiffTime(0),
 fhPtDispEta(0),                     fhPtDispPhi(0),
@@ -1096,13 +1096,18 @@ TList *  AliAnaPi0EbE::GetCreateOutputObjects()
   
   if(fAnaType == kSSCalo)
   {
-    
     fhMassPt  = new TH2F
     ("hMassPt","all pairs #it{M}: #it{p}_{T} vs #it{M}",nptbins,ptmin,ptmax, nmassbins,massmin,massmax);
     fhMassPt->SetYTitle("#it{M} (GeV/#it{c}^{2})");
     fhMassPt->SetXTitle("#it{p}_{T} (GeV/#it{c})");
     outputContainer->Add(fhMassPt) ;
-    
+
+    fhPtLambda0NoSplitCut  = new TH2F
+    ("hPtLambda0NoSplitCut","all clusters: #it{p}_{T} vs #lambda_{0}^{2}",nptbins,ptmin,ptmax, ssbins,ssmin,ssmax);
+    fhPtLambda0NoSplitCut->SetYTitle("#lambda_{0}^{2}");
+    fhPtLambda0NoSplitCut->SetXTitle("#it{p}_{T} (GeV/#it{c})");
+    outputContainer->Add(fhPtLambda0NoSplitCut) ;
+
     fhSelectedMassPt  = new TH2F
     ("hSelectedMassPt","Selected #pi^{0} (#eta) pairs #it{M}: #it{p}_{T} vs #it{M}",nptbins,ptmin,ptmax, nmassbins,massmin,massmax);
     fhSelectedMassPt->SetYTitle("#it{M} (GeV/#it{c}^{2})");
@@ -3207,16 +3212,18 @@ void  AliAnaPi0EbE::MakeShowerShapeIdentification()
       continue ;
     }
 
+    Float_t l0 = calo->GetM02();
     Float_t e1 = l1.Energy();
     Float_t e2 = l2.Energy();
     TLorentzVector l12 = l1+l2;
     Float_t ptSplit = l12.Pt();
     Float_t  eSplit = e1+e2;
-    
+
     //mass of all clusters
     fhMass       ->Fill(mom.E() ,mass);
     fhMassPt     ->Fill(mom.Pt(),mass);
     fhMassSplitPt->Fill(ptSplit ,mass);
+    fhPtLambda0NoSplitCut->Fill(mom.Pt(),l0);
     
     // Asymmetry of all clusters
     Float_t asy =-10;
@@ -3335,7 +3342,7 @@ void  AliAnaPi0EbE::MakeShowerShapeIdentification()
     if(nSM < GetCaloUtils()->GetNumberOfSuperModulesUsed() && nSM >=0)
     {
       fhSelectedMassPtLocMaxSM   [indexMax][nSM]->Fill(mom.Pt(),mass);
-      fhSelectedLambda0PtLocMaxSM[indexMax][nSM]->Fill(mom.Pt(),calo->GetM02());
+      fhSelectedLambda0PtLocMaxSM[indexMax][nSM]->Fill(mom.Pt(),l0  );
     }
     
     if(IsDataMC())
@@ -3419,7 +3426,6 @@ void  AliAnaPi0EbE::MakeShowerShapeIdentification()
     if(nMaxima == 2 && fNLMECutMin[1] > mom.E()) continue;
     if(nMaxima >  2 && fNLMECutMin[2] > mom.E()) continue;
 
-    
     //Fill some histograms about shower shape
     if(fFillSelectClHisto && GetReader()->GetDataType()!=AliCaloTrackReader::kMC)
     {
index 74094cb..ca956e3 100755 (executable)
@@ -190,7 +190,7 @@ class AliAnaPi0EbE : public AliAnaCaloTrackCorrBaseClass {
   TH2F         * fhSelectedMass  ;         //! pair mass vs E, for selected pairs
   TH2F         * fhSelectedMassPt  ;       //! pair mass vs pT, for selected pairs
   TH2F         * fhSelectedMassSplitPt  ;  //! pair mass vs pT (split), for selected pairs
-  
+    
   TH2F         * fhMassPtLocMax[3] ;             //! pair mass vs pT, for all pairs, for each NLM case
   TH2F         * fhSelectedMassPtLocMax[3] ;     //! pair mass vs pT, for selected pairs, for each NLM case
   TH2F         * fhSelectedMassPtLocMaxSM[3][22];//! pair mass vs pT, for selected pairs, for each NLM case, for each SM
@@ -247,8 +247,9 @@ class AliAnaPi0EbE : public AliAnaCaloTrackCorrBaseClass {
   TH1F         * fhEDecay   ;              //! Number of identified  pi0/eta decay photons vs E
   
   TH2F         * fhPtDispersion ;           //! pT vs disp of selected cluster
-  TH2F         * fhPtLambda0 ;              //! pT vs lambda0 of selected cluster 
-  TH2F         * fhPtLambda1 ;              //! pT vs lambda1 of selected cluster 
+  TH2F         * fhPtLambda0 ;              //! pT vs lambda0 of selected cluster
+  TH2F         * fhPtLambda0NoSplitCut ;    //! pT vs lambda0 of cluster before the split selection.
+  TH2F         * fhPtLambda1 ;              //! pT vs lambda1 of selected cluster
   TH2F         * fhPtLambda0NoTRD ;         //! pT vs lambda0 of selected cluster, not behind TRD 
   TH2F         * fhPtLambda0FracMaxCellCut ;//! pT vs lambda0 of selected cluster, fraction of cluster energy in max cell cut 
   TH2F         * fhPtFracMaxCell ;          //! pT vs frac max cell of selected cluster 
@@ -406,7 +407,7 @@ class AliAnaPi0EbE : public AliAnaCaloTrackCorrBaseClass {
   AliAnaPi0EbE(              const AliAnaPi0EbE & pi0ebe) ; // cpy ctor
   AliAnaPi0EbE & operator = (const AliAnaPi0EbE & pi0ebe) ; // cpy assignment
   
-  ClassDef(AliAnaPi0EbE,38)
+  ClassDef(AliAnaPi0EbE,39)
 } ;