Fixing bugs (removing a member from AliSpectraAODTrackCuts.h, added Merge method...
authormfloris <mfloris@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 4 May 2012 05:43:29 +0000 (05:43 +0000)
committermfloris <mfloris@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 4 May 2012 05:43:29 +0000 (05:43 +0000)
PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODPID.cxx
PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODPID.h
PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODTrackCuts.h
PWGLF/SPECTRA/PiKaPr/TestAOD/runAODProof.C

index d4ee35b..cc4a721 100644 (file)
 
 ClassImp(AliSpectraAODPID)
 
-AliSpectraAODPID::AliSpectraAODPID() : TObject(), fPIDType(kNSigmaTPCTOF), fNSigmaPID(3), fPIDResponse(0) {
+AliSpectraAODPID::AliSpectraAODPID() : TNamed("PID", "PID object"), fPIDType(kNSigmaTPCTOF), fNSigmaPID(3), fPIDResponse(0) {
 
 }
 
-AliSpectraAODPID::AliSpectraAODPID(AODPIDType_t pidType) : TObject(), fPIDType(pidType), fNSigmaPID(3), fPIDResponse(0) {
+AliSpectraAODPID::AliSpectraAODPID(AODPIDType_t pidType) : TNamed("PID", "PID object"), fPIDType(pidType), fNSigmaPID(3), fPIDResponse(0) {
 
 
 
@@ -168,3 +168,43 @@ Int_t AliSpectraAODPID::GetParticleSpecie(AliAODTrack      * trk, AliSpectraAODT
   return kSpUndefined;
 
 }
+
+Long64_t AliSpectraAODPID::Merge(TCollection* list)
+{
+  // Merging interface.
+  // Returns the number of merged objects (including this).
+
+  Printf("Merging");
+
+  if (!list)
+    return 0;
+
+  if (list->IsEmpty())
+    return 1;
+
+  TIterator* iter = list->MakeIterator();
+  TObject* obj;
+
+  // Actually, we don't do anything here...
+  // collections of all histograms
+  //  TList collections;
+
+  Int_t count = 0;
+
+  while ((obj = iter->Next())) {
+    AliSpectraAODPID* entry = dynamic_cast<AliSpectraAODPID*> (obj);
+    if (entry == 0) 
+      continue;
+
+    // TH1I * histo = entry->GetHistoCuts();      
+    // collections.Add(histo);
+    count++;
+  }
+  
+  //  fHistoCuts->Merge(&collections);
+  
+  delete iter;
+  Printf("OK");
+  return count+1;
+}
+
index 7ed647f..729cf7d 100644 (file)
@@ -40,7 +40,7 @@ namespace AliSpectraNameSpace {
 
 using namespace AliSpectraNameSpace;
 
-class AliSpectraAODPID : public TObject
+class AliSpectraAODPID : public TNamed
 {
 public:
   AliSpectraAODPID() ;
@@ -53,6 +53,8 @@ public:
   Int_t GetParticleSpecie(AliAODTrack      * trk, AliSpectraAODTrackCuts * trackCuts);
   Int_t GetParticleSpecie(AliAODMCParticle * trk);
   
+  Long64_t Merge(TCollection* list);
+
 
 private:
 
index 392469c..68851ac 100644 (file)
@@ -78,7 +78,6 @@ public:
    Float_t          fPtCutTOFMatching;           // TOF Matching
    Float_t          fQvecCutMin;           // Minimum value of Qvec
    Float_t          fQvecCutMax;           // Minimum value of Qvec
-   AliPIDResponse   *fPIDResponse; // PID response  
    TH1I             *fHistoCuts;       // Cuts statistics
    AliAODTrack      *fTrack;           //! Track pointer
    
index 7837f72..cc8a7b4 100644 (file)
@@ -1,4 +1,4 @@
-#include "AliSpectraAODPID.h" 
+
 
 
 void runAODProof(Int_t c=1, const char * proofMode = "full")
@@ -7,10 +7,10 @@ void runAODProof(Int_t c=1, const char * proofMode = "full")
   //3 data AOD086
   //4 MC AOD090
   
-  gEnv->SetValue("XSec.GSI.DelegProxy", "2");
+  //  gEnv->SetValue("XSec.GSI.DelegProxy", "2");
   
   gSystem->Load("libTree.so");
-   gSystem->Load("libGeom.so");
+  gSystem->Load("libGeom.so");
    gSystem->Load("libVMC.so");
    gSystem->Load("libPhysics.so");
    gSystem->Load("libSTEERBase.so");
@@ -21,7 +21,6 @@ void runAODProof(Int_t c=1, const char * proofMode = "full")
    gSystem->Load("libANALYSISalice.so");
    gSystem->AddIncludePath("-I$ALICE_ROOT/include");
 
-
    AliAnalysisAlien * handler = new AliAnalysisAlien("test");
    handler->SetOverwriteMode();
    handler->SetRunMode(proofMode);
@@ -30,11 +29,10 @@ void runAODProof(Int_t c=1, const char * proofMode = "full")
    //handler->SetAliROOTVersion("v5-03-11-AN");
    handler->SetAliROOTVersion("v5-03-17-AN");
    
-   //handler->SetNproofWorkers(2);
-   handler->SetNproofWorkersPerSlave(1);
+   handler->SetNproofWorkers(10);
+   //     handler->SetNproofWorkersPerSlave(1);
    handler->SetProofCluster(Form("%s@alice-caf.cern.ch", gSystem->Getenv("CAFUSER")));
-   //handler->SetProofCluster(Form("%s@skaf.saske.sk",gSystem->Getenv("CAFUSER")));
-   
+   //   handler->SetProofCluster(Form("%s@skaf.saske.sk",gSystem->Getenv("CAFUSER")));
    // Set handler for Real DATA:
    if (c == 1){
      handler->SetProofDataSet("/alice/data/LHC10h_000138653_p2_AOD049#aodTree");
@@ -55,11 +53,11 @@ void runAODProof(Int_t c=1, const char * proofMode = "full")
    gROOT->LoadMacro("AliSpectraAODHistoManager.cxx+g");
    gROOT->LoadMacro("AliSpectraAODPID.cxx+g");
    gROOT->LoadMacro("AliAnalysisTaskSpectraAOD.cxx+g");
-   
+
    handler->SetAliRootMode("default");
    handler->SetAdditionalLibs("AliSpectraAODHistoManager.cxx AliSpectraAODHistoManager.h AliSpectraAODPID.cxx AliSpectraAODPID.h AliSpectraAODEventCuts.cxx AliSpectraAODEventCuts.h AliSpectraAODTrackCuts.cxx AliSpectraAODTrackCuts.h AliAnalysisTaskSpectraAOD.cxx AliAnalysisTaskSpectraAOD.h");
    handler->SetAnalysisSource("Histograms.h HistogramNames.h AliSpectraAODHistoManager.cxx+ AliSpectraAODEventCuts.cxx+ AliSpectraAODTrackCuts.cxx+ AliSpectraAODPID.cxx+ AliAnalysisTaskSpectraAOD.cxx+");
-   handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
+   //   handler->SetFileForTestMode("filelist.txt"); // list of local files for testing
    //  handler->SetAliRootMode("");
    handler->SetClearPackages();
    
@@ -84,7 +82,8 @@ void runAODProof(Int_t c=1, const char * proofMode = "full")
    // Double_t CentCutMax[4]={10,20,30,40};
    // Double_t QvecCutMin[4]={0,0,0,0};
    // Double_t QvecCutMax[4]={100,100,100,100};
-   AliSpectraAODPID * pid = new AliSpectraAODPID(AODPIDType_t::kNSigmaTPCTOF);
+   using namespace AliSpectraNameSpace;
+   AliSpectraAODPID *pid = new AliSpectraAODPID(kNSigmaTPCTOF); 
    pid->SetNSigmaCut(5.);
    
    for(Int_t iCut=1;iCut<2;iCut++){
@@ -138,7 +137,7 @@ void runAODProof(Int_t c=1, const char * proofMode = "full")
         AliAnalysisDataContainer *coutputpt3 = mgr->CreateContainer(Form("ctcutpt%d",iCut), AliSpectraAODTrackCuts::Class(),     AliAnalysisManager::kOutputContainer, 
                                                                     Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
         AliAnalysisDataContainer *coutputpt4 = mgr->CreateContainer(Form("cpidpt%d",iCut),  AliSpectraAODPID::Class(),     AliAnalysisManager::kOutputContainer, 
-                                                                    Form("Pt.AOD.1._data_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
+                                                                    Form("Pt.AOD.1._data_ptcut_Cent%.0fto%.0f_QVec%.1fto%.1f.root",CentCutMin[iCut],CentCutMax[iCut],QvecCutMin[iCut],QvecCutMax[iCut]));
         
        }
      mgr->ConnectInput(task, 0, cinput);