adding Histos for bck under pi0 vs pt

diff --git a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
index 0986c594806ffe6c86fd7b9f18a09d60ac005cea..94a97dc2b6d2284eb7b69b5de67d1d4210af369d 100644 (file)
--- a/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
+++ b/PWG4/GammaConv/AliAnalysisTaskGammaConversion.cxx
@@ -897,6 +897,9 @@ void AliAnalysisTaskGammaConversion::ProcessMCData(){
         case 331: // Eta'
            fHistograms->FillHistogram("MC_DecayEtapGamma_Pt", particle->Pt());
            break;
+        case 333: // Phi
+           fHistograms->FillHistogram("MC_DecayPhiGamma_Pt", particle->Pt());
+           break;
         }
       }
        
@@ -919,7 +922,10 @@ void AliAnalysisTaskGammaConversion::ProcessMCData(){
       
        fCFManager->GetParticleContainer()->Fill(containerInput,kStepGenerated);                                        // generated gamma
       }                
-      if(particle->GetMother(0) < 0){   // direct gamma
+      if(particle->GetMother(0) < 0 || //Phojet p+p -> Direct Photons have no mother
+        ((particle->GetMother(0) > -1) && 
+         (TMath::Abs(fStack->Particle(particle->GetMother(0))->GetPdgCode()) < 10)) //Pythia p+p -> Direct Photons have quarks as mother
+        ){   // direct gamma
        fHistograms->FillHistogram("MC_allDirectGamma_Energy",particle->Energy());
        fHistograms->FillHistogram("MC_allDirectGamma_Pt", particle->Pt());
        fHistograms->FillHistogram("MC_allDirectGamma_Eta", particle->Eta());
@@ -1539,6 +1545,7 @@ void AliAnalysisTaskGammaConversion::ProcessV0s(){
     if(fV0Reader->GetNegativeTracknTPCFClusters()!=0 && fV0Reader->GetNegativeTracknTPCClusters()!=0 ){
       Double_t eClsToF= (Double_t)fV0Reader->GetNegativeTracknTPCClusters()/(Double_t)fV0Reader->GetNegativeTracknTPCFClusters();
       fHistograms->FillHistogram("ESD_E_nTPCClustersToFP", fV0Reader->GetNegativeTrackP(),eClsToF );
+      fHistograms->FillHistogram("ESD_E_nTPCClustersToFR", fV0Reader->GetXYRadius(),eClsToF );
       fHistograms->FillHistogram("ESD_E_TPCchi2", fV0Reader->GetNegativeTrackTPCchi2()/(Double_t)fV0Reader->GetNegativeTracknTPCClusters());
     }
 
@@ -1553,6 +1560,7 @@ void AliAnalysisTaskGammaConversion::ProcessV0s(){
     if(fV0Reader->GetPositiveTracknTPCFClusters()!=0 && (Double_t)fV0Reader->GetPositiveTracknTPCClusters()!=0 ){
       Double_t pClsToF= (Double_t)fV0Reader->GetPositiveTracknTPCClusters()/(Double_t)fV0Reader->GetPositiveTracknTPCFClusters();
       fHistograms->FillHistogram("ESD_P_nTPCClustersToFP",fV0Reader->GetPositiveTrackP(), pClsToF);
+      fHistograms->FillHistogram("ESD_P_nTPCClustersToFR",fV0Reader->GetXYRadius(), pClsToF);
       fHistograms->FillHistogram("ESD_P_TPCchi2", fV0Reader->GetPositiveTrackTPCchi2()/(Double_t)fV0Reader->GetPositiveTracknTPCClusters());
     }
 
@@ -1783,7 +1791,18 @@ void AliAnalysisTaskGammaConversion::ProcessV0s(){
        fHistograms->FillHistogram("ESD_TrueConvGamma_Pt_Chi2", fV0Reader->GetMotherCandidatePt(), fV0Reader->GetMotherCandidateChi2());
        fHistograms->FillHistogram("ESD_TrueConvGamma_Eta_Chi2", fV0Reader->GetMotherCandidateEta(), fV0Reader->GetMotherCandidateChi2());
                                
-                               
+
+       if(fV0Reader->GetNegativeTracknTPCFClusters()!=0 && fV0Reader->GetNegativeTracknTPCClusters()!=0 ){
+         Double_t eClsToF= (Double_t)fV0Reader->GetNegativeTracknTPCClusters()/
+           (Double_t)fV0Reader->GetNegativeTracknTPCFClusters();
+         fHistograms->FillHistogram("ESD_TrueConversion_E_nTPCClustersToFR", fV0Reader->GetXYRadius(),eClsToF );
+       }
+
+       if(fV0Reader->GetPositiveTracknTPCFClusters()!=0 && fV0Reader->GetPositiveTracknTPCClusters()!=0 ){
+         Double_t pClsToF= (Double_t)fV0Reader->GetPositiveTracknTPCClusters()
+           /(Double_t)fV0Reader->GetPositiveTracknTPCFClusters();
+         fHistograms->FillHistogram("ESD_TrueConversion_P_nTPCClustersToFR",fV0Reader->GetXYRadius(), pClsToF);
+       }
        fHistograms->FillHistogram("ESD_TrueConversion_XY", fV0Reader->GetX(),fV0Reader->GetY());
        fHistograms->FillHistogram("ESD_TrueConversion_R", fV0Reader->GetXYRadius());
        fHistograms->FillHistogram("ESD_TrueConversion_ZR", fV0Reader->GetZ(),fV0Reader->GetXYRadius());
@@ -2375,9 +2394,12 @@ void AliAnalysisTaskGammaConversion::ProcessGammasForNeutralMesonAnalysis(){
 
                    if(!isRealPi0 && !isRealEta){
                      if(gamma1MotherLabel>-1 && gamma2MotherLabel>-1){
-                       fHistograms->FillHistogram("ESD_TrueBckGG_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TrueBckGG_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }else{
-                       fHistograms->FillHistogram("ESD_TrueBckCont_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TrueBckCont_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
+                     }
+                     if(gamma1MotherLabel==-111 || gamma2MotherLabel==-111){
+                       fHistograms->FillHistogram("ESD_TruePi0DalitzCont_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }
                    }
 
@@ -2396,9 +2418,12 @@ void AliAnalysisTaskGammaConversion::ProcessGammasForNeutralMesonAnalysis(){
                    }
                    if(!isRealPi0 && !isRealEta){
                      if(gamma1MotherLabel>-1 && gamma2MotherLabel>-1){
-                       fHistograms->FillHistogram("ESD_TrueBckGG_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TrueBckGG_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }else{
-                       fHistograms->FillHistogram("ESD_TrueBckCont_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TrueBckCont_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
+                     }
+                     if(gamma1MotherLabel==-111 || gamma2MotherLabel==-111){
+                       fHistograms->FillHistogram("ESD_TruePi0DalitzCont_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }
                    }
                  }
@@ -2413,14 +2438,17 @@ void AliAnalysisTaskGammaConversion::ProcessGammasForNeutralMesonAnalysis(){
                      fHistograms->FillHistogram("ESD_TruePi0_InvMass",massTwoGammaCandidate);
                      fHistograms->FillHistogram("ESD_TruePi0_InvMass_vs_Pt_alpha",massTwoGammaCandidate ,momentumVectorTwoGammaCandidate.Pt());
                      if(gamma1MotherLabel > fV0Reader->GetMCStack()->GetNprimary()){
-                       fHistograms->FillHistogram("ESD_TruePi0Sec_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TruePi0Sec_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }
                    }
                    if(!isRealPi0 && !isRealEta){
                      if(gamma1MotherLabel>-1 && gamma2MotherLabel>-1){
-                       fHistograms->FillHistogram("ESD_TrueBckGG_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TrueBckGG_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }else{
-                       fHistograms->FillHistogram("ESD_TrueBckCont_InvMass",massTwoGammaCandidate);
+                       fHistograms->FillHistogram("ESD_TrueBckCont_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
+                     }
+                     if(gamma1MotherLabel==-111 || gamma2MotherLabel==-111){
+                       fHistograms->FillHistogram("ESD_TruePi0DalitzCont_InvMass_vs_Pt",massTwoGammaCandidate,momentumVectorTwoGammaCandidate.Pt());
                      }
                    }
                  }

diff --git a/PWG4/macros/ConfigGammaConversion.C b/PWG4/macros/ConfigGammaConversion.C
index ddcdd5292dd5f152e520f7c4d2eec9e10996bad5..b5475a1fbbc2c5d0d7ee5d1d6fb8e35ad6344739 100644 (file)
--- a/PWG4/macros/ConfigGammaConversion.C
+++ b/PWG4/macros/ConfigGammaConversion.C
@@ -690,9 +690,9 @@ Double_t kGCfirstXBinNITSClusters = -0.5;
 Double_t kGClastXBinNITSClusters = 6.5;
 
 //TPCcluster to Findable-plots
-Int_t kGCnYBinsClsToF = 120;
+Int_t kGCnYBinsClsToF = 200;
 Double_t kGCfirstYBinClsToF = 0.;
-Double_t kGClastYBinClsToF = 1.2;
+Double_t kGClastYBinClsToF = 2.0;
 
 //TPCchi2 -plots
 Int_t kGCnXBinsTPCchi2 = 100;
@@ -1560,6 +1560,11 @@ AliAnalysisTaskGammaConversion* ConfigGammaConversion(TString arguments, AliAnal
   else{
     cout<<"Using mixed event for bg."<<endl;
   }
+  if(kGCUseTrackMultiplicityForBG){
+    cout<<"Using track multiplicity for bck"<<endl;
+  }else{
+    cout<<"Using V0 multiplicity for bck"<<endl;
+  }
   gammaconversion->SetDoRotation(kGCUseRotationMethodInBG);
   cout<<"Using :"<<kGCnumberOfRotationEventsForBG<<" rotations in bg calculation"<<endl;
   gammaconversion->SetNumberOfRotationsBG(kGCnumberOfRotationEventsForBG);
@@ -1880,6 +1885,7 @@ void AddHistograms(AliGammaConversionHistograms *histograms){
      histograms->AddHistogram("MC_DecayOmegaGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
      histograms->AddHistogram("MC_DecayK0sGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
      histograms->AddHistogram("MC_DecayEtapGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
+     histograms->AddHistogram("MC_DecayPhiGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
      histograms->AddHistogram("MC_DecayAllGamma_Pt" ,"" , kGCnXBinsPt, kGCfirstXBinPt, kGClastXBinPt, "", "");
   }
 
@@ -1895,6 +1901,7 @@ void AddHistograms(AliGammaConversionHistograms *histograms){
     if(kGCplotESDENITSClusters == kTRUE){ histograms->AddHistogram("ESD_E_nITSClusters" ,"" , kGCnXBinsNITSClusters, kGCfirstXBinNITSClusters, kGClastXBinNITSClusters, "", "");}
     if(kGCplotESDENTPCClustersToFP== kTRUE){ histograms->AddHistogram("ESD_E_nTPCClustersToFP" ,"" ,kGCnXBinsP, kGCfirstXBinP, kGClastXBinP,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
     if(kGCplotESDENTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_E_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
+    if(kGCplotESDENTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_TrueConversion_E_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
     if(kGCplotESDETPCchi2 == kTRUE){ histograms->AddHistogram("ESD_E_TPCchi2" ,"" , kGCnXBinsTPCchi2, kGCfirstXBinTPCchi2, kGClastXBinTPCchi2, "", "");}
 
                
@@ -1906,6 +1913,7 @@ void AddHistograms(AliGammaConversionHistograms *histograms){
     if(kGCplotESDPNITSClusters == kTRUE){ histograms->AddHistogram("ESD_P_nITSClusters" ,"" , kGCnXBinsNITSClusters, kGCfirstXBinNITSClusters, kGClastXBinNITSClusters, "", "");}
     if(kGCplotESDPNTPCClustersToFP== kTRUE){ histograms->AddHistogram("ESD_P_nTPCClustersToFP" ,"" ,kGCnXBinsP, kGCfirstXBinP, kGClastXBinP,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);} 
     if(kGCplotESDPNTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_P_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
+    if(kGCplotESDPNTPCClustersToFR== kTRUE){ histograms->AddHistogram("ESD_TrueConversion_P_nTPCClustersToFR" ,"" ,kGCnXBinsR, kGCfirstXBinR, kGClastXBinR,kGCnYBinsClsToF, kGCfirstYBinClsToF, kGClastYBinClsToF,"", "",0);}
     if(kGCplotESDPTPCchi2 == kTRUE){ histograms->AddHistogram("ESD_P_TPCchi2" ,"" , kGCnXBinsTPCchi2, kGCfirstXBinTPCchi2, kGClastXBinTPCchi2, "", "");}
 
     if(kGCplotESDConvGammaEnergy == kTRUE){ histograms->AddHistogram("ESD_ConvGamma_Energy" ,"" , kGCnXBinsEnergy, kGCfirstXBinEnergy, kGClastXBinEnergy, "", "");}
@@ -2255,8 +2263,8 @@ void AddHistograms(AliGammaConversionHistograms *histograms){
       }
 
       //      if(kGCdoNeutralMesonV0MCCheck == kTRUE){
-      histograms->AddHistogram("ESD_TrueBckGG_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
-      histograms->AddHistogram("ESD_TrueBckCont_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
+      histograms->AddHistogram("ESD_TrueBckGG_InvMass_vs_Pt","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
+      histograms->AddHistogram("ESD_TrueBckCont_InvMass_vs_Pt","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
       histograms->AddHistogram("ESD_TruePi0Sec_InvMass_vs_Pt","Invariant mass vs Pt",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
       histograms->AddHistogram("ESD_TruePi0DalitzCont_InvMass_vs_Pt","Invariant mass vs Pt",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,kGCnYBinsSpectra, kGCfirstYBinSpectra, kGClastYBinSpectra,"InvMass [GeV]","Counts");
        histograms->AddHistogram("ESD_TruePi0_InvMass","Invariant mass",kGCnXBinsSpectra,kGCfirstXBinSpectra, kGClastXBinSpectra,"InvMass [GeV]","Counts");
@@ -2952,11 +2960,22 @@ Int_t SetAnalysisCutSelection(TString analysisCutSelection){
   switch(BackgroundScheme){
   case 0: //Rotation
     kGCUseRotationMethodInBG=kTRUE;
+    kGCdoBGProbability=kFALSE;
     break;
-  case 1: // mixed event
+  case 1: // mixed event with track multiplicity
     kGCUseRotationMethodInBG=kFALSE;
+    kGCUseTrackMultiplicityForBG=kFALSE;
+    kGCdoBGProbability=kFALSE;
+    break;
+  case 2: // mixed event with V0 multiplicity
+    kGCUseRotationMethodInBG=kFALSE;
+    kGCUseTrackMultiplicityForBG=kTRUE;
+    kGCdoBGProbability=kFALSE;
+    break;
+  case 3: //Rotation
+    kGCUseRotationMethodInBG=kTRUE;
+    kGCdoBGProbability=kTRUE;
     break;
-
   default:
     return iResult;
   }
@@ -2998,6 +3017,12 @@ Int_t SetAnalysisCutSelection(TString analysisCutSelection){
   case 5:
     kGCnumberOfRotationEventsForBG = 50;
     break;
+  case 6:
+    kGCnumberOfRotationEventsForBG = 80;
+    break;
+  case 7:
+    kGCnumberOfRotationEventsForBG = 100;
+    break;
   default:
     return iResult;
   }