Reverting commit 44732
authorhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 29 Oct 2010 15:52:51 +0000 (15:52 +0000)
committerhristov <hristov@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 29 Oct 2010 15:52:51 +0000 (15:52 +0000)
ZDC/AliZDCTowerCalibTask.cxx [deleted file]
ZDC/AliZDCTowerCalibTask.h [deleted file]
ZDC/CMakelibZDCrec.pkg
ZDC/RunZDCTowerCalib.C [deleted file]
ZDC/ZDCrecLinkDef.h
ZDC/libZDCrec.pkg

diff --git a/ZDC/AliZDCTowerCalibTask.cxx b/ZDC/AliZDCTowerCalibTask.cxx
deleted file mode 100644 (file)
index bcb2a7e..0000000
+++ /dev/null
@@ -1,208 +0,0 @@
-/**************************************************************************
- * 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.                  *
- **************************************************************************/
-
-#include "TChain.h"
-#include "TTree.h"
-#include "TH1F.h"
-#include "TCanvas.h"
-#include "TDecompLU.h" 
-
-#include "AliAnalysisTask.h"
-#include "AliAnalysisManager.h"
-
-#include "AliESDEvent.h"
-#include "AliESDInputHandler.h"
-#include "AliESDZDC.h"
-#include "AliZDCTowerCalibTask.h"
-#include "AliZDCTowerCalib.h"
-#include "AliCDBManager.h"
-#include "AliCDBId.h"
-#include "AliCDBRunRange.h"
-#include "AliCDBStorage.h"
-
-ClassImp(AliZDCTowerCalibTask)
-
-//________________________________________________________________________
-AliZDCTowerCalibTask::AliZDCTowerCalibTask() 
-  : AliAnalysisTask(), 
-    fESD(0),
-    fAZNA(),
-    fAZNC(),
-    fBZNA(),
-    fBZNC(),
-    fADCMin(0)
-{
-}
-
-//________________________________________________________________________
-AliZDCTowerCalibTask::AliZDCTowerCalibTask(const char *name) 
-  : AliAnalysisTask(name, ""), 
-    fESD(0),
-    fAZNA(),
-    fAZNC(),
-    fBZNA(),
-    fBZNC(),
-    fADCMin(0)
-{
-  // Default 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, TH1F::Class());
-  
-  fAZNA.ResizeTo(4,4); 
-  fAZNC.ResizeTo(4,4); 
-  fBZNA.ResizeTo(4); 
-  fBZNC.ResizeTo(4); 
-  
-  fAZNA.Zero(); 
-  fAZNC.Zero(); 
-  fBZNA.Zero(); 
-  fBZNC.Zero(); 
-  fADCMin = 0; 
-}
-
-//________________________________________________________________________
-void AliZDCTowerCalibTask::ConnectInputData(Option_t *) 
-{
-  // Connect ESD or AOD here
-  // Called once
-  
-  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);
-    */
-    
-    AliESDInputHandler *esdH = dynamic_cast<AliESDInputHandler*> (AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
-    if (!esdH) Printf("ERROR: Could not get ESDInputHandler");
-    else fESD = esdH->GetEvent();
-  }
-}
-
-//________________________________________________________________________
-void AliZDCTowerCalibTask::CreateOutputObjects(){
-  // Create histograms
-  // Called once
-}
-
-//________________________________________________________________________
-void AliZDCTowerCalibTask::Exec(Option_t *) 
-{
-  // Main loop
-  // Called for each event
-
-  if (!fESD) {
-    Printf("ERROR: fESD not available");
-    return;
-  }
-
-  AliESDZDC *zdcEvent = fESD->GetESDZDC();
-  Double_t *ezn1 = 0, *ezn2 = 0; 
-  if(fESD->GetEventType()==7){
-    ezn1 = (Double_t *) zdcEvent->GetZN1TowerEnergy(); // znc
-    ezn2 = (Double_t *) zdcEvent->GetZN2TowerEnergy(); // zna
-  }
-  Double_t ezn1sum = ezn1[1] + ezn1[2] + ezn1[3] + ezn1[4]; 
-  Double_t ezn2sum = ezn2[1] + ezn2[2] + ezn2[3] + ezn2[4]; 
-  // select event
-  Bool_t isZNAAccepted = kTRUE; 
-  if (ezn2[0] < fADCMin || ezn2sum < fADCMin) isZNAAccepted = kFALSE; 
-  Bool_t isZNCAccepted = kTRUE; 
-  if (ezn1[0] < fADCMin || ezn1sum < fADCMin) isZNCAccepted = kFALSE; 
-  // calculate coefficient matrix and known terms 
-  for (Int_t i=0; i<4; i++) { 
-    for (Int_t j=0; j<4; j++) { 
-      if (isZNAAccepted) fAZNA[i][j] += ezn2[ i + 1 ] * ezn2[ j + 1 ] / ezn2[0]; 
-      if (isZNCAccepted) fAZNC[i][j] += ezn1[ i + 1 ] * ezn1[ j + 1 ] / ezn1[0]; 
-    }
-    if (isZNAAccepted) fBZNA[i] += ezn2[ i + 1 ]; 
-    if (isZNCAccepted) fBZNC[i] += ezn1[ i + 1 ];     
-  }
-  // Post output data.
-  //  PostData(0, fHistPt);
-}      
-
-//________________________________________________________________________
-void AliZDCTowerCalibTask::Terminate(Option_t *) 
-{
-  // solve the system of linear equations giving the calibration coefficients 
-  // for zna and znc
-
-  TDecompLU luZNA(fAZNA);
-  Bool_t isZNAOk;
-  TVectorD coeffZNA = luZNA.Solve(fBZNA,isZNAOk);
-  if (isZNAOk) { 
-    printf ("ZNA: coefficient matrix:\n"); 
-    fAZNA.Print(); 
-    printf ("ZNA: known terms:\n"); 
-    fBZNA.Print(); 
-    printf ("ZNA: Fitted calibration coefficients:\n"); 
-    coeffZNA.Print(); 
-  }
-  else { 
-    printf ("Singular coefficient matrix for ZNA! Setting all coefficients to 1\n"); 
-    for (Int_t i=0; i<4; i++) coeffZNA[i] = 1; 
-  } 
-
-  TDecompLU luZNC(fAZNC);
-  Bool_t isZNCOk;
-  TVectorD coeffZNC = luZNC.Solve(fBZNC,isZNCOk);
-  if (isZNCOk) { 
-    printf ("ZNC: coefficient matrix:\n"); 
-    fAZNC.Print(); 
-    printf ("ZNC: known terms:\n"); 
-    fBZNC.Print(); 
-    printf ("ZNC: Fitted calibration coefficients:\n"); 
-    coeffZNC.Print(); 
-  }
-  else {
-    printf ("Singular coefficient matrix for ZNC! Setting all coefficients to 1\n"); 
-    for (Int_t i=0; i<4; i++) coeffZNC[i] = 1; 
-  } 
-  // write to OCDB
-
-
-
-  AliZDCTowerCalib *towerCalib = new AliZDCTowerCalib();
-  
-  towerCalib->SetZN1EqualCoeff(0, 1.);
-  towerCalib->SetZN2EqualCoeff(0, 1.);
-  
-  for(Int_t j=1; j<5; j++){  
-    towerCalib->SetZN1EqualCoeff(j, coeffZNC[j+1]);
-    towerCalib->SetZN2EqualCoeff(j, coeffZNA[j+1]);
-  }
-  towerCalib->Print("");
-  
-  AliCDBManager *manager = AliCDBManager::Instance();
-  manager->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
-  
-  AliCDBMetaData *md = new AliCDBMetaData();
-  md->SetResponsible("Chiara Oppedisano");
-  md->SetComment("Calibration object for ZDC written by a macro");
-  md->SetObjectClassName("AliZDCTowerCalib");
-  md->SetAliRootVersion(gSystem->Getenv("ARVERSION"));
-  md->SetBeamPeriod(0);
-  
-  AliCDBId id("ZDC/Calib/TowerCalib",fESD->GetRunNumber(),AliCDBRunRange::Infinity());
-  manager->GetDefaultStorage()->Put(towerCalib,id, md);
-}
diff --git a/ZDC/AliZDCTowerCalibTask.h b/ZDC/AliZDCTowerCalibTask.h
deleted file mode 100644 (file)
index cc90b7c..0000000
+++ /dev/null
@@ -1,41 +0,0 @@
-#ifndef AliZDCTowerCalibTask_h
-#define AliZDCTowerCalibTask_h
-
-// analysis task for ZN towers intercalibration
-// Author: Alessandro De Falco, INFN Cagliari
-
-#include "TMatrixD.h"
-#include "TVectorD.h"
-
-class TH1F;
-class AliESDEvent;
-
-#include "AliAnalysisTask.h"
-
-class AliZDCTowerCalibTask : public AliAnalysisTask {
- public:
-  AliZDCTowerCalibTask();
-  AliZDCTowerCalibTask(const char *name);
-  virtual ~AliZDCTowerCalibTask() {}
-  
-  virtual void   ConnectInputData(Option_t *);
-  virtual void   CreateOutputObjects();
-  virtual void   Exec(Option_t *option);
-  virtual void   Terminate(Option_t *);
-  virtual void   SetADCMin(Double_t adcmin) { fADCMin = adcmin; } 
-
-  
- private:
-  AliESDEvent *fESD;   //ESD object
-  TMatrixD fAZNA;      // coefficient matrix for ZNA calorimeter
-  TMatrixD fAZNC;      // coefficient matrix for ZNC calorimeter
-  TVectorD fBZNA;      // vector of known terms for ZNA calorimeter
-  TVectorD fBZNC;      // vector of known terms for ZNC calorimeter
-  Double_t fADCMin; 
-  AliZDCTowerCalibTask(const AliZDCTowerCalibTask&); // not implemented
-  AliZDCTowerCalibTask& operator=(const AliZDCTowerCalibTask&); // not implemented
-  
-  ClassDef(AliZDCTowerCalibTask, 1); 
-};
-
-#endif
index d03d755..a0e7cc4 100644 (file)
@@ -25,7 +25,7 @@
 # SHLIBS - Shared Libraries and objects for linking (Executables only)           #
 #--------------------------------------------------------------------------------#
 
-set ( SRCS  AliZDCReco.cxx AliZDCRecoParam.cxx AliZDCRecoParampp.cxx AliZDCRecoParamPbPb.cxx AliZDCReconstructor.cxx AliZDCQADataMakerRec.cxx AliZDCTowerCalibTask.cxx)
+set ( SRCS  AliZDCReco.cxx AliZDCRecoParam.cxx AliZDCRecoParampp.cxx AliZDCRecoParamPbPb.cxx AliZDCReconstructor.cxx AliZDCQADataMakerRec.cxx)
 
 string ( REPLACE ".cxx" ".h" HDRS "${SRCS}" )
 
diff --git a/ZDC/RunZDCTowerCalib.C b/ZDC/RunZDCTowerCalib.C
deleted file mode 100644 (file)
index 087674b..0000000
+++ /dev/null
@@ -1,64 +0,0 @@
-void RunZDCTowerCalib(char *path="/alice/cern.ch/user/d/defalco/data", 
-                   Bool_t isLocal=kFALSE, 
-                   Double_t adcmin=7){
-  // load analysis framework
-  gSystem->Load("libANALYSIS");
-  gSystem->Load("libANALYSISalice") ;
-
-  gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
-  AliAnalysisGrid *alienHandler = 0;    // only used in GRID analysis   
-  TChain* chain = 0;
-  if (!isLocal) {
-    gROOT->LoadMacro("CreateAlienHandler.C"); // Create and configure the alien handler plugin
-    alienHandler = CreateAlienHandler(path);  
-    if (!alienHandler) return;
-  }
-  else { 
-    char filename[500]; 
-    sprintf (filename,"%s/AliESDs.root",path); 
-    chain = new TChain("esdTree"); 
-    chain->Add(filename); 
-    chain->Add(filename); 
-  }
-
-  gSystem->AddIncludePath("-I$ALICE_ROOT/include");
-  gSystem->AddIncludePath("-I$ALICE_ROOT/ZDC");
-
-  // Create the analysis manager
-
-  AliAnalysisManager *mgr = new AliAnalysisManager("ZNCalibManager");
-  if (!isLocal) mgr->SetGridHandler(alienHandler);
-
-  AliVEventHandler* esdH = new AliESDInputHandler();
-  mgr->SetInputEventHandler(esdH);
-
-  // Create task
-
-  gROOT->LoadMacro("AliZDCTowerCalibTask.cxx+g");
-  AliAnalysisTask *task = new AliZDCTowerCalibTask("TaskZDCCalib");
-  ((AliZDCTowerCalibTask*) task)->SetADCMin(adcmin); 
-
-  // Add task
-  mgr->AddTask(task);
-
-  // Create containers for input/output
-
-  AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
-//   AliAnalysisDataContainer *coutput = mgr->CreateContainer("chist", TH1::Class(),    AliAnalysisManager::kOutputContainer, "Pt.ESD.1.root");
-
-  // Connect input/output
-  mgr->ConnectInput(task, 0, cinput);
-//   mgr->ConnectOutput(task, 0, coutput);
-
-  // Enable debug printouts
-  mgr->SetDebugLevel(2);
-
-  if (!mgr->InitAnalysis())
-    return;
-
-  mgr->PrintStatus();
-
-  if (isLocal) mgr->StartAnalysis("local", chain);
-  else mgr->StartAnalysis("grid");
-
-}
index 4775adf..920e308 100644 (file)
@@ -14,6 +14,5 @@
 #pragma link C++ class  AliZDCRecoParamPbPb+;
 #pragma link C++ class  AliZDCReconstructor+;
 #pragma link C++ class  AliZDCQADataMakerRec+;
-#pragma link C++ class  AliZDCTowerCalibTask+;
 
 #endif
index 319243b..a914a8f 100644 (file)
@@ -3,7 +3,7 @@
 
 SRCS= AliZDCReco.cxx AliZDCRecoParam.cxx AliZDCRecoParampp.cxx \
 AliZDCRecoParamPbPb.cxx AliZDCReconstructor.cxx \
-AliZDCQADataMakerRec.cxx AliZDCTowerCalibTask.cxx
+AliZDCQADataMakerRec.cxx
 
 HDRS= $(SRCS:.cxx=.h)