Commit for Simone - Temorarily removal of some histos
authorbdoenigu <benjamin.doenigus@cern.ch>
Tue, 1 Jul 2014 10:24:52 +0000 (12:24 +0200)
committerbdoenigu <benjamin.doenigus@cern.ch>
Tue, 1 Jul 2014 10:24:52 +0000 (12:24 +0200)
PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx
PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.h

index 6c88ef4..2483f30 100644 (file)
@@ -67,9 +67,9 @@ AliAnalysisTaskV0ForRAA::AliAnalysisTaskV0ForRAA()
   fESDTrackCutsLowPt(0),
   fOutputContainer(0),
 //thnf
-  fTHnFK0s(0),
-  fTHnFL(0),
-  fTHnFAL(0),
+//  fTHnFK0s(0),
+// fTHnFL(0),
+//  fTHnFAL(0),
   fTHnFK0sDauEta(0),
   fTHnFLDauEta(0),
   fTHnFALDauEta(0),
@@ -540,9 +540,9 @@ AliAnalysisTaskV0ForRAA::AliAnalysisTaskV0ForRAA(const char *name)
   fESDTrackCutsLowPt(0),
   fOutputContainer(0),
   //thnf
-  fTHnFK0s(0),
-  fTHnFL(0),
-  fTHnFAL(0),
+  //  fTHnFK0s(0),
+  // fTHnFL(0),
+  // fTHnFAL(0),
   fTHnFK0sDauEta(0),
   fTHnFLDauEta(0),
   fTHnFALDauEta(0),
@@ -1111,8 +1111,8 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   // ------------------- add always ---------------------------//
   //THnF
   Double_t piForAx = 2.0*TMath::Pi();
-  Int_t binsTHnV0K0s[4] = {150,100,80,18};
-  Int_t binsTHnV0L[4] = {100,100,80,18};
+  // Int_t binsTHnV0K0s[4] = {150,100,80,18};
+  // Int_t binsTHnV0L[4] = {100,100,80,18};
 
   Int_t binsTHnV0DauEtaK0s[4] = {150,100,40,18};
   Int_t binsTHnV0DauEtaL[4] = {100,100,40,18};
@@ -1120,28 +1120,28 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   Int_t binsTHnV0DauPhiK0s[4] = {150, 100,40,18};
   Int_t binsTHnV0DauPhiL[4] = {100, 100,40,18};
 
-  Double_t minK0s[4] = {0.35,0.0,0.0,0.0};
-  Double_t maxK0s[4] = {0.65,50.0,40.0,piForAx};
+  //Double_t minK0s[4] = {0.35,0.0,0.0,0.0};
+  //Double_t maxK0s[4] = {0.65,50.0,40.0,piForAx};
   Double_t minK0sDauEta[4] = {0.35, 0.0,-0.8,0.0};
   Double_t maxK0sDauEta[4] = {0.65,50.0, 0.8,piForAx};
   Double_t minK0sDauPhi[4] = {0.35,0.0,-0.8,0.0};
   Double_t maxK0sDauPhi[4] = {0.65,50.0,0.8,piForAx};
 
-  Double_t minL[4] = {1.07, 0.0,0.0,0.0};
-  Double_t maxL[4] = {1.17,50.0, 40.0,piForAx};
+  //Double_t minL[4] = {1.07, 0.0,0.0,0.0};
+  // Double_t maxL[4] = {1.17,50.0, 40.0,piForAx};
   Double_t minLDauEta[4] = {1.07, 0.0,-0.8,0.0};
   Double_t maxLDauEta[4] = {1.17,50.0, 0.8, piForAx};
   Double_t minLDauPhi[4] = {1.07,0.0,-0.8,0.0};
   Double_t maxLDauPhi[4] = {1.17,50.0,0.8,piForAx};
   
-  
+  /*
   char histTitK0s[255];
   snprintf(histTitK0s,255,"fTHnFK0s");
   char histTitL[255];
   snprintf(histTitL,255,"fTHnFL");
   char histTitAL[255];
   snprintf(histTitAL,255,"fTHnFAL");
-
+  */
 
   char histTitK0sDauEta[255];
   snprintf(histTitK0sDauEta,255,"fTHnFK0sDauEta");
@@ -1157,14 +1157,14 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   snprintf(histTitLDauPhi,255,"fTHnFLDauPhi");
   char histTitALDauPhi[255];
   snprintf(histTitALDauPhi,255,"fTHnFALDauPhi");
+  /*
   char axTitK0s[255];
   snprintf(axTitK0s,255,"K^{0}_{s};m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
   char axTitL[255];
   snprintf(axTitL,255,"#Lambda;m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
   char axTitAL[255];
   snprintf(axTitAL,255,"#bar{#Lambda};m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
-
+  */
 
   char axTitK0sDauEta[255];
   snprintf(axTitK0sDauEta,255,"K^{0}_{s} daughter;m_{inv} (GeV/c^{2});p_{T} (Gev/c);#eta_{pos};#phi_{pos}");
@@ -1181,14 +1181,14 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   char axTitALDauPhi[255];
   snprintf(axTitALDauPhi,255,"#bar{#Lambda} daughter;m_{inv} (GeV/c^{2});p_{T}(GeV/c);#eta_{neg};#phi_{neg}");
 
-
+  /*
   fTHnFK0s = new       THnF(histTitK0s,axTitK0s,4,binsTHnV0K0s,minK0s,maxK0s);
   // fTHnFK0s->Sumw2();
   fTHnFL   = new       THnF(histTitL  ,axTitL  ,4,binsTHnV0L,minL  ,maxL);
   // fTHnFL->Sumw2();
   fTHnFAL  = new       THnF(histTitAL ,axTitAL ,4,binsTHnV0L,minL  ,maxL);
   //  fTHnFAL->Sumw2();
-
+  */
 
   fTHnFK0sDauEta = new         THnF(histTitK0sDauEta,axTitK0sDauEta,4,binsTHnV0DauEtaK0s,minK0sDauEta,maxK0sDauEta);
   //  fTHnFK0sDauEta->Sumw2();
@@ -1339,7 +1339,7 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   fOutputContainer->Add(fHistPiPiMonitorMCCuts);
   fOutputContainer->Add(fHistPiPiK0sVsLambdaMass);
   fOutputContainer->Add(fHistPiPiK0sVsALambdaMass);
-  fOutputContainer->Add(fTHnFK0s);
+  //  fOutputContainer->Add(fTHnFK0s);
   fOutputContainer->Add(fTHnFK0sDauEta);
   fOutputContainer->Add(fTHnFK0sDauPhi);
   //fOutputContainer->Add(fHistPiPiPhiPosVsPtPosVsMass);//xxx      
@@ -1348,7 +1348,7 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   fOutputContainer->Add(fHistPiPK0sVsLambdaMass);
   fOutputContainer->Add(fHistPiPALambdaVsLambdaMass);
   //  fOutputContainer->Add(fHistPiPPhiPosVsPtPosVsMass);//xxx
-  fOutputContainer->Add(fTHnFL);
+  //fOutputContainer->Add(fTHnFL);
   fOutputContainer->Add(fTHnFLDauEta);
   fOutputContainer->Add(fTHnFLDauPhi);
 
@@ -1356,7 +1356,7 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   fOutputContainer->Add(fHistPiAPK0sVsALambdaMass);
   fOutputContainer->Add(fHistPiAPLambdaVsALambdaMass);
   //  if(fMCTruthMode)  fOutputContainer->Add(fHistPiAPPhiPosVsPtPosVsMass);//xxx
-  fOutputContainer->Add(fTHnFAL);
+  //fOutputContainer->Add(fTHnFAL);
   fOutputContainer->Add(fTHnFALDauEta);
   fOutputContainer->Add(fTHnFALDauPhi);
 
@@ -2657,10 +2657,10 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
       fHistPiAPMonitorMCCuts[isSecd]->Fill(16*fillFlagAL);
 
     
-      Double_t phiMC =   p0->Phi(); 
+      //      Double_t phiMC =   p0->Phi(); 
       //Double_t etaMC =   p0->Eta(); 
 
-      Double_t valTHnMC[4]  = {massV0MC,ptV0MC,0.0,phiMC}; 
+      //      Double_t valTHnMC[4]  = {massV0MC,ptV0MC,0.0,phiMC}; 
       Double_t valTHnMCDauEta[4] = {massV0MC,ptV0MC,etaMCPos,phiMCPos};
       Double_t valTHnMCDauPhi[4] = {massV0MC,ptV0MC,etaMCNeg,phiMCNeg};
 
@@ -2691,8 +2691,8 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
        // fHistPiPiPhiPosVsPtPosVsMass->Fill(massV0MC,ctTK0s,ptV0MC);//,ctK0s);//phiPosMC);//xxx
        fHistPiPiK0sVsLambdaMass->Fill(calcLambdamass,calcK0smass);
        fHistPiPiK0sVsALambdaMass->Fill(calcALambdamass,calcK0smass);
-       valTHnMC[2] = ctTK0s; 
-       fTHnFK0s->Fill(valTHnMC);
+       //      valTHnMC[2] = ctTK0s; 
+       //      fTHnFK0s->Fill(valTHnMC);
        fTHnFK0sDauEta->Fill(valTHnMCDauEta);
        fTHnFK0sDauPhi->Fill(valTHnMCDauPhi);
       }
@@ -2724,8 +2724,8 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
        
        fHistPiPK0sVsLambdaMass->Fill(calcLambdamass,calcK0smass);
 
-              valTHnMC[2] = ctTL;
-fTHnFL->Fill(valTHnMC);
+       //valTHnMC[2] = ctTL;
+       //fTHnFL->Fill(valTHnMC);
        fTHnFLDauEta->Fill(valTHnMCDauEta);
        fTHnFLDauPhi->Fill(valTHnMCDauPhi);
       }
@@ -2758,8 +2758,8 @@ fTHnFL->Fill(valTHnMC);
        // if(isSecd <1) fHistPiPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
        //else fHistPiAPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx       
        //fHistPiAPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
-       valTHnMC[2] = ctTAL;
-       fTHnFAL->Fill(valTHnMC);
+       //      valTHnMC[2] = ctTAL;
+       //fTHnFAL->Fill(valTHnMC);
        fTHnFALDauEta->Fill(valTHnMCDauEta);
        fTHnFALDauPhi->Fill(valTHnMCDauPhi);
       }
@@ -3474,9 +3474,9 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
     //    if( ppTrack.Angle(pmTrack)<0.001) continue;  
     //    if( ppTrack.Angle(pmTrack)<0.004) continue;   
       
-    Double_t px = v0K0.Px();
-    Double_t py = v0K0.Py();
-    Double_t phi  = TMath::Pi()+TMath::ATan2(-py, -px);
+    //    Double_t px = v0K0.Px();
+    //    Double_t py = v0K0.Py();
+    //    Double_t phi  = TMath::Pi()+TMath::ATan2(-py, -px);
     //    Double_t eta =  v0K0.Eta();
 
     /*     
@@ -3672,11 +3672,11 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
            fHistPiPiPtVSY->Fill(rapK0s,ptK0s);
            fHistPiPiDecayLengthVsMass->Fill(massK0s,dim2V0Radius);//decayLength);
            // fHistPiPiPhiPosVsPtPosVsMass->Fill(massK0s,ctTK0,ptV0MC);//,ctK0);//posDaughterPhi);//xxx
-           Double_t valTHnK0s[4]= {massK0s,ptV0MC,ctTK0,phi};
+           //Double_t valTHnK0s[4]= {massK0s,ptV0MC,ctTK0,phi};
            Double_t valTHnK0sDauEta[4]= {massK0s,ptV0MC,posDaughterEta,posDaughterPhi};
            Double_t valTHnK0sDauPhi[4]= {massK0s,ptV0MC,negDaughterEta,negDaughterPhi};
            if(massK0s >=0.35 && massK0s <= 0.65){
-             fTHnFK0s->Fill(valTHnK0s);            
+             // fTHnFK0s->Fill(valTHnK0s);         
              fTHnFK0sDauEta->Fill(valTHnK0sDauEta);
              fTHnFK0sDauPhi->Fill(valTHnK0sDauPhi);
            }
@@ -3817,13 +3817,13 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
            //Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,negDaughterEta};
            //      Double_t valTHnLDauPhi[5]= {massLambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
 
-           Double_t valTHnL[4]= {massLambda,ptV0MC,ctTL,phi};
+           //      Double_t valTHnL[4]= {massLambda,ptV0MC,ctTL,phi};
            Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,posDaughterPhi};
            Double_t valTHnLDauPhi[4]= {massLambda,ptV0MC,negDaughterEta,negDaughterPhi};
 
 
            if(massLambda >=1.07 && massLambda <= 1.17){
-             fTHnFL->Fill(valTHnL);
+             //              fTHnFL->Fill(valTHnL);
              fTHnFLDauEta->Fill(valTHnLDauEta);
              fTHnFLDauPhi->Fill(valTHnLDauPhi);
            }
@@ -3978,14 +3978,14 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
            Double_t valTHnALDauEta[4]={massALambda,ptV0MC,posDaughterEta,negDaughterEta};
            Double_t valTHnALDauPhi[5]={massALambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
            */
-           Double_t valTHnAL[4]= {massALambda,ptV0MC,ctTAL,phi};
+           //Double_t valTHnAL[4]= {massALambda,ptV0MC,ctTAL,phi};
            Double_t valTHnALDauEta[4]= {massALambda,ptV0MC,posDaughterEta,posDaughterPhi};
            Double_t valTHnALDauPhi[4]= {massALambda,ptV0MC,negDaughterEta,negDaughterPhi};
 
            if(massALambda >=1.07 && massALambda <= 1.17){
              fTHnFALDauEta->Fill(valTHnALDauEta);
              fTHnFALDauPhi->Fill(valTHnALDauPhi);
-             fTHnFAL->Fill(valTHnAL);
+             //fTHnFAL->Fill(valTHnAL);
            }
            if(fMCMode && !fMCTruthMode) fHistPiAPPDGCode->Fill(pdgBG);
            if( massALambda>1.108 && massALambda<1.123 )  fHistPiAPDecayLengthVsPt[isSecd]->Fill(ptV0MC,dim2V0Radius);//decayLength);
index 7c39339..c8f199b 100644 (file)
@@ -164,9 +164,9 @@ class AliAnalysisTaskV0ForRAA : public AliAnalysisTaskSE {
   TList           *fOutputContainer;    // output data container
    
   //----------------------------histograms --------------------------------------------//
-  THnF *fTHnFK0s; 
-  THnF *fTHnFL; 
-  THnF *fTHnFAL; 
+  // THnF *fTHnFK0s; 
+  //THnF *fTHnFL; 
+  //THnF *fTHnFAL; 
 
   THnF *fTHnFK0sDauEta; 
   THnF *fTHnFLDauEta;