TMap object for scaling VZERO A response in 7 TeV pp added
authorlmilano <lmilano@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 17 Jul 2013 10:33:49 +0000 (10:33 +0000)
committerlmilano <lmilano@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 17 Jul 2013 10:33:49 +0000 (10:33 +0000)
PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx
PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.h

index a3bae53..76cd906 100644 (file)
@@ -24,6 +24,7 @@
 #include <TH2F.h>
 #include <TH3F.h>
 #include <TRandom.h>
+#include <TParameter.h>
 
 #include "AliAnalysisTaskPhiCorrelations.h"
 #include "AliAnalyseLeadingTrackUE.h"
@@ -56,6 +57,7 @@
 
 #include "AliHelperPID.h"
 #include "AliAnalysisUtils.h"
+#include "TMap.h"
 
 ////////////////////////////////////////////////////////////////////////
 //
@@ -96,9 +98,10 @@ fUseVtxAxis(kFALSE),
 fCourseCentralityBinning(kFALSE),
 fSkipTrigger(kFALSE),
 fInjectedSignals(kFALSE),
-// pointers to UE classes
 fHelperPID(0x0),
 fAnalysisUtils(0x0),
+fMap(0x0),
+// pointers to UE classes
 fAnalyseUE(0x0),
 fHistos(0x0),
 fHistosMixed(0),
@@ -302,6 +305,9 @@ void  AliAnalysisTaskPhiCorrelations::CreateOutputObjects()
   // add HelperPID to list
   if (fHelperPID)
     fListOfHistos->Add(fHelperPID);
+  // add TMap to list
+  if (fMap)
+    fListOfHistos->Add(fMap);
   
   fListOfHistos->Add(new TH2F("trackletsVsV0Cent", ";L1 clusters;v0 centrality", 100, -0.5, 9999.5, 101, 0, 101));
   fListOfHistos->Add(new TH2F("processIDs", ";#Delta#phi;process id", 100, -0.5 * TMath::Pi(), 1.5 * TMath::Pi(), kPNoProcess + 1, -0.5, kPNoProcess + 0.5));
@@ -309,7 +315,7 @@ void  AliAnalysisTaskPhiCorrelations::CreateOutputObjects()
   fListOfHistos->Add(new TH2F("mixedDist", ";centrality;tracks;events", 101, 0, 101, 200, 0, fMixingTracks * 1.5));
   fListOfHistos->Add(new TH1F("pids", ";pdg;tracks", 2001, -1000.5, 1000.5));
   fListOfHistos->Add(new TH2F("referenceMultiplicity", ";centrality;tracks;events", 101, 0, 101, 200, 0, 200));
-  fListOfHistos->Add(new TH1F("V0AMult", ";;V0A multiplicity", 2000, -.5, 1999.5));
+  fListOfHistos->Add(new TH2F("V0AMult", "V0A multiplicity;V0A multiplicity;V0A multiplicity (scaled)", 1000, -.5, 999.5, 1000, -.5, 999.5));
   
   PostData(0,fListOfHistos);
   
@@ -929,6 +935,15 @@ void  AliAnalysisTaskPhiCorrelations::AnalyseDataMode()
     }
     else if (fCentralityMethod == "V0A_MANUAL")
     {
+      //Total multiplicity in the VZERO A detector
+      Float_t MV0A=inputEvent->GetVZEROData()->GetMTotV0A();
+      Float_t MV0AScaled=0.;
+      if (fMap){
+       TParameter<float>* sf=(TParameter<float>*)fMap->GetValue(Form("%d",inputEvent->GetRunNumber()));
+       if(sf)MV0AScaled=MV0A*sf->GetVal();
+      }
+      ((TH2F*) fListOfHistos->FindObject("V0AMult"))->Fill(MV0A,MV0AScaled);
+      
       // for pp, to be implemented
       centrality = 1;
     }
@@ -1037,9 +1052,6 @@ void  AliAnalysisTaskPhiCorrelations::AnalyseDataMode()
     return;
   }
 
-  //Total multiplicity in the VZERO A detector
-  ((TH1F*) fListOfHistos->FindObject("V0AMult"))->Fill(inputEvent->GetVZEROData()->GetMTotV0A());
-  
   // correlate particles with...
   TObjArray* tracksCorrelate = 0;
   if (fParticleSpeciesAssociated != fParticleSpeciesTrigger)
index aace72f..9654ff2 100644 (file)
@@ -47,6 +47,7 @@ class AliESDEvent;
 class AliHelperPID;
 class AliAnalysisUtils;
 class TFormula;
+class TMap;
 
 
 class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
@@ -120,6 +121,9 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
    
     AliAnalysisUtils* GetAnalysisUtils() { return fAnalysisUtils; }
     void   SetAnalysisUtils(AliAnalysisUtils* utils){ fAnalysisUtils = utils; }
+   
+    TMap* GetMap() { return fMap; }
+    void   SetMap(TMap* map){ fMap = map; }
     
   private:
     AliAnalysisTaskPhiCorrelations(const  AliAnalysisTaskPhiCorrelations &det);
@@ -148,9 +152,11 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     Bool_t             fSkipTrigger;             // skip trigger selection
     Bool_t             fInjectedSignals;         // check header to skip injected signals in MC
     
-    // Pointers to external UE classes
-    AliHelperPID*     fHelperPID;      // points to class containing common analysis algorithms
-    AliAnalysisUtils*     fAnalysisUtils;      // points to class containing common analysis algorithms
+    AliHelperPID*     fHelperPID;      // points to class for PID
+    AliAnalysisUtils*     fAnalysisUtils;      // points to class with common analysis utilities
+    TMap*     fMap;                   // points to TMap class containing scaling factors for VZERO A signal
+   // Pointers to external UE classes
     AliAnalyseLeadingTrackUE*     fAnalyseUE;      //! points to class containing common analysis algorithms
     AliUEHistograms*  fHistos;       //! points to class to handle histograms/containers  
     AliUEHistograms*  fHistosMixed;       //! points to class to handle mixed histograms/containers  
@@ -210,7 +216,7 @@ class  AliAnalysisTaskPhiCorrelations : public AliAnalysisTask
     
     Bool_t fFillpT;                // fill sum pT instead of number density
     
-    ClassDef( AliAnalysisTaskPhiCorrelations, 35); // Analysis task for delta phi correlations
+    ClassDef( AliAnalysisTaskPhiCorrelations, 36); // Analysis task for delta phi correlations
   };
 
 class AliDPhiBasicParticle : public AliVParticle