adding TPC cosmics selector based on ESD

diff --git a/PWG0/PWG0selectorsLinkDef.h b/PWG0/PWG0selectorsLinkDef.h
index 3b70f65f3d0fbf55985e95e2d5b067e006737ab2..52a9de4d9ffd47a79aced8f33ff71750bb56e462 100644 (file)
--- a/PWG0/PWG0selectorsLinkDef.h
+++ b/PWG0/PWG0selectorsLinkDef.h
@@ -15,5 +15,6 @@
 #pragma link C++ class AliMultiplicityMCSelector+;
 #pragma link C++ class AlidNdEtaSystematicsSelector+;
 #pragma link C++ class AliTestESDtrackCutsSelector+;
+#pragma link C++ class AliROCESDAnalysisSelector+;
 
 #endif

diff --git a/PWG0/TPC/AliROCESDAnalysisSelector.cxx b/PWG0/TPC/AliROCESDAnalysisSelector.cxx
new file mode 100644 (file)
index 0000000..31d608d
--- /dev/null
+++ b/PWG0/TPC/AliROCESDAnalysisSelector.cxx
@@ -0,0 +1,199 @@
+/**************************************************************************
+ * 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.                  *
+ **************************************************************************/
+
+/* $Id$ */
+
+// The ESD is available as member fESD
+//
+// The Process function is nearly empty. Implement your analysis there and look at the other listed below functions you
+// might need.
+//
+// The following methods can be overrriden. Please do not forgot to call the base class function.
+//
+//    Begin():        called everytime a loop on the tree starts,
+//                    a convenient place to create your histograms.
+//    SlaveBegin():   called after Begin(), when on PROOF called only on the
+//                    slave servers.
+//    Init():         called for each new tree. Enable/Disable branches here.
+//    Process():      called for each event, in this function you decide what
+//                    to read and fill your histograms.
+//    SlaveTerminate: called at the end of the loop on the tree, when on PROOF
+//                    called only on the slave servers.
+//    Terminate():    called at the end of the loop on the tree,
+//                    a convenient place to draw/fit your histograms.
+//
+//  Author: Jan.Fiete.Grosse-Oetringhaus@cern.ch
+
+#include "AliROCESDAnalysisSelector.h"
+
+#include <AliLog.h>
+#include <AliESD.h>
+#include <AliESDfriend.h>
+#include <../TPC/AliTPCclusterMI.h>
+#include <../TPC/AliTPCseed.h>
+
+#include <TH2F.h>
+#include <TFile.h>
+#include <TTree.h>
+#include <TCanvas.h>
+#include <TProfile.h>
+
+ClassImp(AliROCESDAnalysisSelector)
+
+AliROCESDAnalysisSelector::AliROCESDAnalysisSelector() :
+  AliSelector(),
+  fESDfriend(0),
+  fhQmaxVsRow(0)
+{
+  //
+  // Constructor. Initialization of pointers
+  //
+}
+
+AliROCESDAnalysisSelector::~AliROCESDAnalysisSelector()
+{
+  //
+  // Destructor
+  //
+}
+
+void AliROCESDAnalysisSelector::SlaveBegin(TTree* tree)
+{
+  //
+  
+  AliSelector::SlaveBegin(tree);
+  
+  fhQmaxVsRow = new TH2F("QmaxVsRow", "Qmax vs. pad row;Pad row;Qmax", 91, -0.5, 90, 201, -0.5, 200.5);
+} 
+
+void AliROCESDAnalysisSelector::Init(TTree *tree)
+{
+  // The Init() function is called when the selector needs to initialize
+  // a new tree or chain. Typically here the branch addresses of the tree
+  // will be set. It is normaly not necessary to make changes to the
+  // generated code, but the routine can be extended by the user if needed.
+  // Init() will be called many times when running with PROOF.
+
+  AliSelector::Init(tree);
+
+  // Set branch address
+  if (tree)
+    tree->SetBranchAddress("ESDfriend", &fESDfriend);
+    
+  if (fESDfriend != 0)
+    AliDebug(AliLog::kInfo, "INFO: Found ESDfriend branch in chain.");
+}
+
+Bool_t AliROCESDAnalysisSelector::Process(Long64_t entry)
+{
+  //
+  // Implement your analysis here. Do not forget to call the parent class Process by
+  // if (AliSelector::Process(entry) == kFALSE)
+  //   return kFALSE;
+  //
+
+  if (AliSelector::Process(entry) == kFALSE)
+    return kFALSE;
+
+  // Check prerequisites
+  if (!fESD)
+  {
+    AliDebug(AliLog::kError, "ESD branch not available");
+    return kFALSE;
+  }
+
+  // Check prerequisites
+  if (!fESDfriend)
+  {
+    AliDebug(AliLog::kError, "ESDfriend branch not available");
+    return kFALSE;
+  }
+
+  fESD->SetESDfriend(fESDfriend);
+
+  Int_t nTracks = fESD->GetNumberOfTracks();
+  
+  // loop over esd tracks
+  for (Int_t t=0; t<nTracks; t++)
+  {
+    AliESDtrack* esdTrack = dynamic_cast<AliESDtrack*> (fESD->GetTrack(t));
+    if (!esdTrack)
+    {
+      AliDebug(AliLog::kError, Form("ERROR: Could not retrieve track %d.", t));
+      continue;
+    }
+    
+    AliESDfriendTrack* friendtrack = const_cast<AliESDfriendTrack*> (dynamic_cast<const AliESDfriendTrack*> (esdTrack->GetFriendTrack()));
+    if (!friendtrack)
+    {
+      AliDebug(AliLog::kError, Form("ERROR: Could not retrieve friend of track %d.", t));
+      continue;
+    }
+    
+    const AliTPCseed* seed = dynamic_cast<const AliTPCseed*> (friendtrack->GetCalibObject(0));
+    if (!seed)
+    {
+      AliDebug(AliLog::kError, Form("ERROR: Could not retrieve seed of track %d.", t));
+      continue;
+    }
+    
+    for (Int_t clusterID = 0; clusterID < 160; clusterID++)
+    {
+      AliTPCclusterMI* cluster = seed->GetClusterPointer(clusterID);
+      if (!cluster)
+      {
+        //AliDebug(AliLog::kError, Form("ERROR: Could not retrieve cluster %d of track %d.", clusterID, t));
+        continue;
+      }
+      
+      //AliDebug(AliLog::kDebug, Form("We found a cluster from sector %d", cluster->GetDetector()));
+
+      if (cluster->GetDetector() != 5)
+        continue;
+      
+      fhQmaxVsRow->Fill(cluster->GetRow(), cluster->GetMax());
+    }
+  }
+   
+  return kTRUE;
+}
+
+void AliROCESDAnalysisSelector::SlaveTerminate()
+{
+  //
+  
+  fOutput->Add(fhQmaxVsRow);
+} 
+
+void AliROCESDAnalysisSelector::Terminate()
+{
+  //
+  
+  fhQmaxVsRow = dynamic_cast<TH2F*> (fOutput->FindObject("QmaxVsRow"));
+  if (!fhQmaxVsRow)
+    return;
+
+  new TCanvas;
+  fhQmaxVsRow->DrawCopy("COLZ");
+    
+  new TCanvas;
+  TProfile* profile = fhQmaxVsRow->ProfileX("fhQmaxVsRowProfile");
+  profile->DrawCopy();
+    
+  TFile* file = TFile::Open("rocESD.root", "RECREATE");
+  fhQmaxVsRow->Write();
+  profile->Write();
+  file->Close();
+} 

diff --git a/PWG0/TPC/AliROCESDAnalysisSelector.h b/PWG0/TPC/AliROCESDAnalysisSelector.h
new file mode 100644 (file)
index 0000000..9c59020
--- /dev/null
+++ b/PWG0/TPC/AliROCESDAnalysisSelector.h
@@ -0,0 +1,36 @@
+/* $Id$ */
+
+#ifndef AliROCESDAnalysisSelector_H
+#define AliROCESDAnalysisSelector_H
+
+#include "AliSelector.h"
+
+class AliESDfriend;
+class TH2F;
+
+// this is an empty selector that can be used to create an analysis
+
+class AliROCESDAnalysisSelector : public AliSelector {
+  public:
+    AliROCESDAnalysisSelector();
+    virtual ~AliROCESDAnalysisSelector();
+
+    virtual void    SlaveBegin(TTree* tree);
+    virtual void    Init(TTree *tree);
+    virtual Bool_t  Process(Long64_t entry);
+    virtual void    SlaveTerminate();
+    virtual void    Terminate();
+
+ protected:
+    AliESDfriend* fESDfriend;
+    
+    TH2F* fhQmaxVsRow;
+
+ private:
+    AliROCESDAnalysisSelector(const AliROCESDAnalysisSelector&);
+    AliROCESDAnalysisSelector& operator=(const AliROCESDAnalysisSelector&);
+
+  ClassDef(AliROCESDAnalysisSelector, 0);
+};
+
+#endif

diff --git a/PWG0/TPC/runROCESDAnalysis.C b/PWG0/TPC/runROCESDAnalysis.C
new file mode 100644 (file)
index 0000000..e02da46
--- /dev/null
+++ b/PWG0/TPC/runROCESDAnalysis.C
@@ -0,0 +1,38 @@
+/* $Id$ */
+
+//
+// 
+//
+
+#include "../CreateESDChain.C"
+#include "../PWG0Helper.C"
+
+void runROCESDAnalysis(Char_t* data, Int_t nRuns=20, Int_t offset=0, Bool_t aDebug = kFALSE, Bool_t aProof = kFALSE, 
+    const char* option = "", const char* proofServer = "jgrosseo@lxb6046")
+{
+  if (aProof)
+    connectProof(proofServer);
+
+  TString libraries("libEG;libGeom;libESD");
+  TString packages;
+
+  if (!prepareQuery(libraries, packages, 2))
+    return;
+
+  TChain* chain = CreateESDChain(data, nRuns, offset, kTRUE, kTRUE);
+
+  TList inputList;
+
+  TString selectorName = "AliROCESDAnalysisSelector";
+  AliLog::SetClassDebugLevel(selectorName, AliLog::kInfo);
+
+  if (aDebug != kFALSE)
+  {
+    AliLog::SetClassDebugLevel(selectorName, AliLog::kDebug);
+    selectorName += ".cxx+g";
+  }
+  else
+    selectorName += ".cxx+";
+
+  Int_t result = executeQuery(chain, &inputList, selectorName, option);
+}

diff --git a/PWG0/libPWG0selectors.pkg b/PWG0/libPWG0selectors.pkg
index cf20add7431259d3b82f86c562be44865779bb89..c498d0c3f3170e890bb9d33c5f53c0c300a6e1ad 100644 (file)
--- a/PWG0/libPWG0selectors.pkg
+++ b/PWG0/libPWG0selectors.pkg
@@ -9,7 +9,8 @@ HDRS = dNdEta/AlidNdEtaCorrectionSelector.h \
        dNdEta/AliMultiplicityESDSelector.h \
        dNdEta/AliMultiplicityMCSelector.h \
        dNdEta/AlidNdEtaSystematicsSelector.h \
-       esdTrackCuts/AliTestESDtrackCutsSelector.h
+       esdTrackCuts/AliTestESDtrackCutsSelector.h \
+       TPC/AliROCESDAnalysisSelector.h
 
 SRCS = $(HDRS:.h=.cxx)