]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWGHF/hfe/AliHFEVZEROEventPlane.cxx
Merge branch 'master' of https://git.cern.ch/reps/AliRoot
[u/mrichter/AliRoot.git] / PWGHF / hfe / AliHFEVZEROEventPlane.cxx
index 3580da9438ab2ea4e428800c6f7fe6878e698df5..c6212dcc39631724986dd3f6e261a32b9e16d3e7 100644 (file)
@@ -1,13 +1,38 @@
+/**************************************************************************
+* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+*                                                                        *
+* Author: The ALICE Off-line Project.                                    *
+* Contributors are mentioned in the code where appropriate.              *
+*                                                                        *
+* Permission to use, copy, modify and distribute this software and its   *
+* documentation strictly for non-commercial purposes is hereby granted   *
+* without fee, provided that the above copyright notice appears in all   *
+* copies and that both the copyright notice and this permission notice   *
+* appear in the supporting documentation. The authors make no claims     *
+* about the suitability of this software for any purpose. It is          *
+* provided "as is" without express or implied warranty.                  *
+**************************************************************************/
+//
+//
+// VZERO event plane task for 2010
+// Gain equalization + Recentering 
+// Need a root file
+// 
+// Author:
+//   Raphaelle Bailhache <R.Bailhache@gsi.de>
+//
+//
+
 #include "AliHFEVZEROEventPlane.h"
 
 // AliRoot includes
 #include "AliAnalysisManager.h"
 #include "AliInputEventHandler.h"
 #include "AliVEvent.h"
-#include "AliESDEvent.h"
-#include "AliESDtrack.h"
+//#include "AliESDEvent.h"
 #include "AliCentrality.h"
-#include "AliESDVZERO.h"
+//#include "AliESDVZERO.h"
+#include "AliVVZERO.h"
 #include "TFile.h"
 #include "AliOADBContainer.h"
 #include "TH2F.h"
@@ -76,8 +101,8 @@ AliHFEVZEROEventPlane::AliHFEVZEROEventPlane(const char *name, const Char_t *tit
   // Constructor
   //
 
-  char namelist[100];
-  sprintf(namelist,"QA_hist_%s",name);
+  TString namelist("QA_hist_");
+  namelist += name;
   fOutputList = new TList();
   fOutputList->SetName((const char*)namelist);
   fOutputList->SetOwner(kTRUE);
@@ -91,35 +116,55 @@ AliHFEVZEROEventPlane::AliHFEVZEROEventPlane(const char *name, const Char_t *tit
   fOutputList->Add(fMultV0After);
 
   // Recentering
-  char namecontbefore[100];
-  char namecontafter[100];    
   for(Int_t k = 0; k < fgknCentrBin; k++) {
     for(Int_t iside = 0; iside < 2; iside++) {
       for(Int_t icoord = 0; icoord < 2; icoord++) {
        
-       if(iside==0 && icoord==0)
-         sprintf(namecontbefore,"hQxc2_%i",k);
-       else if(iside==1 && icoord==0)
-         sprintf(namecontbefore,"hQxa2_%i",k);
-       else if(iside==0 && icoord==1)
-         sprintf(namecontbefore,"hQyc2_%i",k);
-       else if(iside==1 && icoord==1)
-         sprintf(namecontbefore,"hQya2_%i",k);
-       //
-       if(iside==0 && icoord==0)
-         sprintf(namecontafter,"hQxc2_%i_after",k);
-       else if(iside==1 && icoord==0)
-         sprintf(namecontafter,"hQxa2_%i_after",k);
-       else if(iside==0 && icoord==1)
-         sprintf(namecontafter,"hQyc2_%i_after",k);
-       else if(iside==1 && icoord==1)
-         sprintf(namecontafter,"hQya2_%i_after",k);
+       TString namecontbefore("");
+       TString namecontafter("");
        
+       if(iside==0 && icoord==0) {
+         namecontbefore += "hQxc2_";
+         namecontbefore += k;
+       }
+       else if(iside==1 && icoord==0){
+         namecontbefore += "hQxa2_";
+         namecontbefore += k;
+       }
+               else if(iside==0 && icoord==1){
+         namecontbefore += "hQyc2_";
+         namecontbefore += k;    
+       }
+       else if(iside==1 && icoord==1) {
+         namecontbefore += "hQya2_";
+         namecontbefore += k;
+       }
+       //
+       if(iside==0 && icoord==0){
+         namecontafter += "hQxc2_";
+         namecontafter += k;
+         namecontafter += "_after";
+       }
+       else if(iside==1 && icoord==0){
+         namecontafter += "hQxa2_";
+         namecontafter += k;
+         namecontafter += "_after";
+       }
+       else if(iside==0 && icoord==1) {
+         namecontafter += "hQyc2_";
+         namecontafter += k;
+         namecontafter += "_after";
+       }
+       else if(iside==1 && icoord==1) {
+         namecontafter += "hQya2_";
+         namecontafter += k;
+         namecontafter += "_after";    
+       }
        //
-       fQBefore[k][iside][icoord] = new TH1F(((const char*)namecontbefore),"",800,-400.0,400.0);
+       fQBefore[k][iside][icoord] = new TH1F(((const char*)namecontbefore),"",1600,-800.0,800.0);
        fQBefore[k][iside][icoord]->Sumw2();
        fOutputList->Add(fQBefore[k][iside][icoord]);
-       fQAfter[k][iside][icoord] = new TH1F(((const char*)namecontafter),"",800,-400.0,400.0);
+       fQAfter[k][iside][icoord] = new TH1F(((const char*)namecontafter),"",1600,-800.0,800.0);
        fQAfter[k][iside][icoord]->Sumw2();
        fOutputList->Add(fQAfter[k][iside][icoord]);
        //
@@ -184,21 +229,10 @@ AliHFEVZEROEventPlane::~AliHFEVZEROEventPlane(){
     delete fOutputList;
   }
   fOutputList = 0x0;
-  if(fMultV0Before) delete fMultV0Before;
-  if(fMultV0After) delete fMultV0After;
-
-  for(Int_t k = 0; k < fgknCentrBin; k++) {
-    for(Int_t iside = 0; iside < 2; iside++) {
-      for(Int_t icoord = 0; icoord < 2; icoord++) {
-       if(fQBefore[k][iside][icoord]) delete fQBefore[k][iside][icoord];
-       if(fQAfter[k][iside][icoord]) delete fQAfter[k][iside][icoord];
-      }
-    }
-  }
-
+  
 }
 //______________________________________________________________________________
-void AliHFEVZEROEventPlane::ProcessEvent(AliESDEvent *event) 
+void AliHFEVZEROEventPlane::ProcessEvent(AliVEvent *event) 
 {
   //
   // Process the event
@@ -225,7 +259,7 @@ void AliHFEVZEROEventPlane::ProcessEvent(AliESDEvent *event)
 }
 
 //________________________________________________________________________
-void AliHFEVZEROEventPlane::Analyze(AliESDEvent* esdEvent)
+void AliHFEVZEROEventPlane::Analyze(AliVEvent* esdEvent)
 {  
   //
   // Do VZERO calibration + centering
@@ -262,7 +296,7 @@ void AliHFEVZEROEventPlane::Analyze(AliESDEvent* esdEvent)
     Double_t qxc2 = 0, qyc2 = 0;
     
     //V0 info    
-    AliESDVZERO* esdV0 = esdEvent->GetVZEROData();
+    AliVVZERO* esdV0 = esdEvent->GetVZEROData();
     
     for (Int_t iv0 = 0; iv0 < 64; iv0++) {
       Double_t phiV0 = TMath::PiOver4()*(0.5 + iv0 % 8);
@@ -371,25 +405,34 @@ Bool_t AliHFEVZEROEventPlane::OpenInfoCalbration(Int_t run)
   for(Int_t iside=0;iside<2;iside++){
     for(Int_t icoord=0;icoord<2;icoord++){
       for(Int_t i=0;i  < fgknCentrBin;i++){
-       char namecont[100];
-       if(iside==0 && icoord==0)
-         sprintf(namecont,"hQxc2_%i",i);
-       else if(iside==1 && icoord==0)
-         sprintf(namecont,"hQxa2_%i",i);
-       else if(iside==0 && icoord==1)
-         sprintf(namecont,"hQyc2_%i",i);
-       else if(iside==1 && icoord==1)
-         sprintf(namecont,"hQya2_%i",i);
+
+       TString namecont("");
+       if(iside==0 && icoord==0) {
+         namecont += "hQxc2_";
+         namecont += i;
+       }
+       else if(iside==1 && icoord==0){
+         namecont += "hQxa2_";
+         namecont += i;
+       }
+       else if(iside==0 && icoord==1){
+         namecont += "hQyc2_";
+         namecont += i;
+       }
+       else if(iside==1 && icoord==1) {
+         namecont += "hQya2_";
+         namecont += i;
+       }
        
        cont = (AliOADBContainer*) foadb->Get(namecont);
        if(!cont){
-         printf("OADB object %s is not available in the file\n",namecont);
+         printf("OADB object %s is not available in the file\n",(const char*)namecont);
          delete fpol0;
          return kFALSE;        
        }
        
        if(!(cont->GetObject(run))){
-         printf("OADB object %s is not available for run %i\n",namecont,run);
+         printf("OADB object %s is not available for run %i\n",(const char*)namecont,run);
          delete fpol0;
          return kFALSE;        
        }