Moving test task to TPC folder
authorzampolli <chiara.zampolli@cern.ch>
Thu, 31 Jul 2014 13:05:03 +0000 (15:05 +0200)
committerzampolli <chiara.zampolli@cern.ch>
Thu, 31 Jul 2014 13:05:22 +0000 (15:05 +0200)
ANALYSIS/ANALYSISaliceLinkDef.h
ANALYSIS/AliAnalysisTaskPt.cxx [deleted file]
ANALYSIS/AliAnalysisTaskPt.h [deleted file]
ANALYSIS/CMakelibANALYSISalice.pkg
TPC/CMakelibTPCcalib.pkg
TPC/TPCcalibLinkDef.h

index 9a35105..3e2e6c6 100644 (file)
@@ -50,7 +50,6 @@
 
 #pragma link C++ class AliAnalysisUtils+;
 #pragma link C++ class AliAnalysisTaskBadChunkID+;
-#pragma link C++ class  AliAnalysisTaskPt+;
 
 #ifdef WITHXML
 #pragma link C++ class AliTagAnalysis+;
diff --git a/ANALYSIS/AliAnalysisTaskPt.cxx b/ANALYSIS/AliAnalysisTaskPt.cxx
deleted file mode 100644 (file)
index 3f67553..0000000
+++ /dev/null
@@ -1,195 +0,0 @@
-#include "TChain.h"
-#include "TTree.h"
-#include "TH1F.h"
-#include "TCanvas.h"
-#include "TObjArray.h"
-
-#include "AliAnalysisTask.h"
-#include "AliAnalysisManager.h"
-
-#include "AliESDEvent.h"
-#include "AliVVevent.h"
-#include "AliVVtrack.h"
-#include "AliESDtrackCuts.h"
-#include "AliVEventHandler.h"
-#include "../TPC/Rec/AliTPCseed.h"
-#include "../TPC/Rec/AliTPCclusterMI.h"
-
-#include "AliAnalysisTaskPt.h"
-
-// example of an analysis task creating a p_t spectrum
-// Authors: Panos Cristakoglou, Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing
-
-ClassImp(AliAnalysisTaskPt)
-
-//________________________________________________________________________
-AliAnalysisTaskPt::AliAnalysisTaskPt(const char *name) 
-: AliAnalysisTask(name, ""), fESD(0), fESDfriend(0), fHistPt(0), fCuts(0), fEv(0), fHistQ(0), fListOut(0)
-{
-  // Constructor
-
-  // Define input and output slots here
-  // Input slot #0 works with a TChain
-  DefineInput(0, TChain::Class());
-  // Output slot #0 writes into a TH1 container
-  DefineOutput(0, TList::Class());
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskPt::ConnectInputData(Option_t *) 
-{
-  // Connect ESD or AOD here
-  // Called once
-
-  printf("----> AliAnalysisTaskPt::ConnectInputData\n");
-  TTree* tree = dynamic_cast<TTree*> (GetInputData(0));
-  if (!tree) {
-    Printf("ERROR: Could not read chain from input slot 0");
-  } else {
-    // Disable all branches and enable only the needed ones
-    // The next two lines are different when data produced as AliESDEvent is read
-    /*
-    tree->SetBranchStatus("*", kFALSE);
-    tree->SetBranchStatus("fTracks.*", kTRUE);
-    */
-
-    AliVEventHandler *esdH = AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler();
-
-    if (!esdH) {
-      Printf("ERROR: Could not get ESDInputHandler");
-    } else {
-      Printf("----> AliAnalysisTaskPt::ConnectInputData Getting the event from handler %p", esdH);
-      //fESD = dynamic_cast<AliESDEvent*>(esdH->GetEvent());
-      fESD = esdH->GetEvent();
-      fESDfriend = esdH->GetFriendEvent();
-    }
-    if (!fESD) {
-      Printf("ERROR, no ESD event");
-    }
-    if (!fESDfriend) {
-      Printf("ERROR, no ESD friend");
-    }
-  }
-
-  Printf("fESD = %p, fESDfriend = %p", fESD, fESDfriend);
-  printf("<---- AliAnalysisTaskPt::ConnectInputData\n");
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskPt::CreateOutputObjects()
-{
-  // Create histograms
-  // Called once
-
-  fListOut = new TList();
-  fListOut->SetOwner();
-  fListOut->SetName("listHistos");
-
-  fHistPt = new TH1F("fHistPt", "P_{T} distribution", 15, 0.1, 3.1);
-  fHistPt->GetXaxis()->SetTitle("P_{T} (GeV/c)");
-  fHistPt->GetYaxis()->SetTitle("dN/dP_{T} (c/GeV)");
-  fHistPt->SetMarkerStyle(kFullCircle);
-
-  fHistQ = new TH1F("fHistQ", "TPC clusters: Q distribution", 1000, 0, 10000);
-  fHistQ->GetXaxis()->SetTitle("Q");
-  fHistQ->GetYaxis()->SetTitle("dN/dQ");
-  fHistQ->SetMarkerStyle(kFullCircle);
-
-  fListOut->Add(fHistPt);
-  fListOut->Add(fHistQ);
-
-  PostData(0, fListOut);
-
-  fCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(1);
-}
-
-//________________________________________________________________________
-void AliAnalysisTaskPt::Exec(Option_t *) 
-{
-  // Main loop
-  // Called for each event
-
-  if (!fESD) {
-    Printf("ERROR: fESD not available");
-    return;
-  }
-  if (!fESDfriend) {
-    Printf("ERROR: fESDfriend not available");
-    return;
-  }
-
-  Printf("There are %d tracks in this event", fESD->GetNumberOfTracks());
-  Printf("... and there are %d friends in this event", fESDfriend->GetNumberOfTracks());
-
-  // Track loop to fill a pT spectrum
-  for (Int_t iTracks = 0; iTracks < fESD->GetNumberOfTracks(); iTracks++) {
-    AliVVtrack* track = fESD->GetTrack(iTracks);
-    if (!track) {
-      Printf("ERROR: Could not receive track %d", iTracks);
-      continue;
-    }
-
-    fHistPt->Fill(track->Pt());
-  } //track loop 
-
-
-  // Friend Track loop
-  for (Int_t iFriend = 0; iFriend < fESDfriend->GetNumberOfTracks(); iFriend++) {
-    AliESDfriendTrack* friendTrack = fESDfriend->GetTrack(iFriend);
-    if (!friendTrack) {
-      Printf("ERROR: Could not receive track %d", iFriend);
-      continue;
-    } 
-    TObject* calibObject;
-    AliTPCseed* seed = NULL;
-    for (Int_t idx = 0; (calibObject = friendTrack->GetCalibObject(idx)); ++idx) {
-      Printf(" |Cal %d = %p", idx, calibObject); 
-      if ((seed = dynamic_cast<AliTPCseed*>(calibObject))) {
-       Printf("Found TPC seed");
-       for (Int_t irow = 0; irow < 160; irow++){
-         AliTPCclusterMI* cluMI = seed->GetClusterPointer(irow);
-         if (cluMI){
-           Printf("Found cluster at row %d", irow);
-           Float_t q = cluMI->GetQ();
-           Printf("Q = %f", q);
-           fHistQ->Fill(q);
-         }
-         else {
-           Printf("Row %d does not contain clusters", irow);
-         }
-       }        
-      }
-    }    
-  }
-
-  // Post output data.
-  PostData(0, fListOut);
-  fEv++;
-}      
-
-//________________________________________________________________________
-void AliAnalysisTaskPt::Terminate(Option_t *) 
-{
-  // Draw result to the screen
-  // Called once at the end of the query
-
-  Printf("Terminate called: fESD = %p", fESD);
-
-  fListOut = dynamic_cast<TList*> (GetOutputData(0)); 
-
-  if (fListOut) {
-    fHistPt = dynamic_cast<TH1F*>(fListOut->FindObject("fHistPt")); 
-    if (!fHistPt) {
-      Printf("ERROR: fHistPt not available");
-      return;
-    }
-   
-    TCanvas *c1 = new TCanvas("AliAnalysisTaskPt","Pt",10,10,510,510);
-    c1->cd(1)->SetLogy();
-    fHistPt->DrawCopy("E");
-  }
-  else {
-    Printf("In Terminate: no TList found");
-  }
-
-}
diff --git a/ANALYSIS/AliAnalysisTaskPt.h b/ANALYSIS/AliAnalysisTaskPt.h
deleted file mode 100644 (file)
index 153b1a2..0000000
+++ /dev/null
@@ -1,42 +0,0 @@
-#ifndef AliAnalysisTaskPt_cxx
-#define AliAnalysisTaskPt_cxx
-
-// example of an analysis task creating a p_t spectrum
-// Authors: Panos Cristakoglou, Jan Fiete Grosse-Oetringhaus, Christian Klein-Boesing
-
-class TH1F;
-class AliESDEvent;
-class AliESDfriend;
-class AliVVevent;
-class AliESDtrackCuts;
-class TList;
-
-#include "AliAnalysisTask.h"
-
-class AliAnalysisTaskPt : public AliAnalysisTask {
- public:
- AliAnalysisTaskPt() : AliAnalysisTask(), fESD(0), fESDfriend(0), fHistPt(0), fCuts(0), fEv(0), fHistQ(0), fListOut(0) {}
-  AliAnalysisTaskPt(const char *name);
-  virtual ~AliAnalysisTaskPt() {}
-  
-  virtual void   ConnectInputData(Option_t *);
-  virtual void   CreateOutputObjects();
-  virtual void   Exec(Option_t *option);
-  virtual void   Terminate(Option_t *);
-
- private:
-  AliVVevent*       fESD;          // ESD object
-  AliESDfriend*     fESDfriend;    // ESD friend object
-  TH1F*             fHistPt;       // Pt spectrum
-  AliESDtrackCuts*  fCuts;         // cuts
-  Int_t fEv;
-  TH1F*             fHistQ;        // TPC clusters Q spectrum
-  TList*            fListOut;      // output list
-
-  AliAnalysisTaskPt(const AliAnalysisTaskPt&); // not implemented
-  AliAnalysisTaskPt& operator=(const AliAnalysisTaskPt&); // not implemented
-  
-  ClassDef(AliAnalysisTaskPt, 1); // example of analysis
-};
-
-#endif
index 59b8170..9d0adbe 100644 (file)
@@ -64,7 +64,6 @@ set ( SRCS
     AliAnalysisUtils.cxx
     BadChunkFilter/AliAnalysisTaskBadChunkID.cxx
     AliAODv0KineCuts.cxx
-    AliAnalysisTaskPt.cxx
     )
 
 if( ROOTHASALIEN STREQUAL "yes")
index b5ece9b..e508501 100644 (file)
@@ -49,6 +49,7 @@ set ( SRCS
     Calib/AliTPCkalmanAlign.cxx 
     Calib/AliTPCcalibSummary.cxx 
     Calib/AliTPCcalibGainMult.cxx 
+    Calib/AliAnalysisTaskPt.cxx
     )
 
 string ( REPLACE ".cxx" ".h" HDRS "${SRCS}" )
index 584e2ab..19bf0f2 100644 (file)
@@ -49,6 +49,7 @@
                                                    // --- move functionality to AliTPCCorrectionFit
 
 #pragma link C++ class AliTPCcalibSummary;         // Tree creation of calibration parameters
+#pragma link C++ class  AliAnalysisTaskPt+;
 
 #endif