Bug fix - integer rounding (Marian)
authormarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 26 May 2008 19:25:09 +0000 (19:25 +0000)
committermarian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 26 May 2008 19:25:09 +0000 (19:25 +0000)
TPC/AliTPCCalibKr.cxx
TPC/AliTPCCalibKr.h

index 07ab8473fa6fc3016953df7b46f053411ed247a0..3f25c81474151ac306eb8cb36c3d09092c0b92e2 100644 (file)
 // The AliTPCCalibKr class description (TPC Kr calibration).
 //
 //
 // The AliTPCCalibKr class description (TPC Kr calibration).
 //
 //
-// The AliTPCCalibKr fills the array of TH3F histograms (TPC_max_padraw,TPC_max_pad,TPC_ADC_cluster),
-// its data memebers.   
+// The AliTPCCalibKr keeps the array of TH3F histograms (TPC_max_padraw,TPC_max_pad,TPC_ADC_cluster),
+// its data memebers and is filled by AliTPCCalibKrTask under conditions (Accept()).   
 // 
 // 
-// As the input it requires the tree with reconstructed Kr clusters (AliTPCclusterKr objects). 
-// The AliTPCCalibKr objects containing an array of TH3F histograms are stored (by default) in the 
-// ouptut (outHistFile.root) file.
-//
 // The ouput TH3F histograms are later used to determine the calibration parameters of TPC chambers.
 // These calculations are done by using AliTPCCalibKr::Analyse() function. The ouput calibration 
 // parameters (details in AliTPCCalibKr::Analyse()) are stored in the calibKr.root file for each TPC pad.
 // In addition the debugCalibKr.root file with debug information is created.
 //
 // The ouput TH3F histograms are later used to determine the calibration parameters of TPC chambers.
 // These calculations are done by using AliTPCCalibKr::Analyse() function. The ouput calibration 
 // parameters (details in AliTPCCalibKr::Analyse()) are stored in the calibKr.root file for each TPC pad.
 // In addition the debugCalibKr.root file with debug information is created.
 //
+
+/*
 // Usage example:
 //
 // Usage example:
 //
-// 1. Create outHistFile.root histogram file:
-//
-// -- Load libXrdClient.so if data on Xrd cluster e.g. (GSI)
-// gSystem->Load("/usr/local/grid/XRootd/GSI/lib64/libXrdClient.so");
-//
-// -- Load toolkit
-// gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros");
-// gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+");
-// AliXRDPROOFtoolkit tool;
-//
-// -- Make chain of files
-// TChain * chain = tool.MakeChain("KrClusters.txt","Kr","",1000,0);
-//
-// -- Run AliTPCCalibKr task (Only TPC C side)
-// AliTPCCalibKr *task = new AliTPCCalibKr;
-// task->SetInputChain(chain);
-// task->SetASide(kFALSE);
-//
-// task->Process();
-// 
-// 2. Analyse output histograms:
-//
-// TFile f("outHistFile.root");
-// AliTPCCalibKr.Analyse();
-//
-// 3. See calibration parameters e.g.:
-//
-// TFile f("calibKr.root");
-// spectrMean->GetCalROC(70)->GetValue(40,40);
-// fitMean->GetCalROC(70)->GetValue(40,40);
-//
-// 4. See debug information e.g.:
-//
-// TFile f("debugCalibKr.root");
-// .ls;
-//
+
+// 1. Analyse output histograms:
+TFile f("outHistFile.root");
+AliTPCCalibKr *obj = (AliTPCCalibKr*) cOutput.FindObject("AliTPCCalibKr")
+obj->Analyse();
+
+// 2. See calibration parameters e.g.:
+TFile f("calibKr.root");
+spectrMean->GetCalROC(70)->GetValue(40,40);
+fitMean->GetCalROC(70)->GetValue(40,40);
+
+// 3. See debug information e.g.:
+TFile f("debugCalibKr.root");
+.ls;
+
 // -- Print calibKr TTree content 
 // -- Print calibKr TTree content 
-// calibKr->Print();
-//
+calibKr->Print();
+
 // -- Draw calibKr TTree variables
 // -- Draw calibKr TTree variables
-// calibKr.Draw("fitMean");
-//
+calibKr.Draw("fitMean");
+
+*/
+
+
 //
 // Author: Jacek Otwinowski (J.Otwinowski@gsi.de) and Stafan Geartner (S.Gaertner@gsi.de)
 //-----------------------------------------------------------------------------
 //
 // Author: Jacek Otwinowski (J.Otwinowski@gsi.de) and Stafan Geartner (S.Gaertner@gsi.de)
 //-----------------------------------------------------------------------------
@@ -113,29 +94,24 @@ ClassImp(AliTPCCalibKr)
 AliTPCCalibKr::AliTPCCalibKr() : 
   TObject(),
   
 AliTPCCalibKr::AliTPCCalibKr() : 
   TObject(),
   
-  bOutputHisto(kTRUE),
   bASide(kTRUE),
   bCSide(kTRUE),
   bASide(kTRUE),
   bCSide(kTRUE),
-  fClusters(0),
-  fClustKr(0),
-  fTree(0),
   fHistoKrArray(72)
 {
   //
   // default constructor
   //
   fHistoKrArray(72)
 {
   //
   // default constructor
   //
+
+  // init histograms
+  Init();
 }
 
 //_____________________________________________________________________
 AliTPCCalibKr::AliTPCCalibKr(const AliTPCCalibKr& pad) : 
   TObject(pad),
   
 }
 
 //_____________________________________________________________________
 AliTPCCalibKr::AliTPCCalibKr(const AliTPCCalibKr& pad) : 
   TObject(pad),
   
-  bOutputHisto(pad.bOutputHisto),
   bASide(pad.bASide),
   bCSide(pad.bCSide),
   bASide(pad.bASide),
   bCSide(pad.bCSide),
-  fClusters(pad.fClusters),
-  fClustKr(pad.fClustKr),
-  fTree(pad.fTree),
   fHistoKrArray(72)
 {
   // copy constructor
   fHistoKrArray(72)
 {
   // copy constructor
@@ -156,9 +132,12 @@ AliTPCCalibKr::~AliTPCCalibKr()
   //
   // destructor
   //
   //
   // destructor
   //
-  if(fClustKr)  delete fClustKr; fClustKr = 0;
-  if(fClusters) delete fClusters; fClusters = 0;
-  if(fTree)     delete fTree; fTree = 0;
+
+  for (Int_t iSec = 0; iSec < 72; ++iSec) 
+  {
+      TH3F *hNew = (TH3F*)fHistoKrArray.At(iSec);
+      if(hNew) delete hNew; hNew = 0; 
+  }
   fHistoKrArray.Delete();
 }
 
   fHistoKrArray.Delete();
 }
 
@@ -179,24 +158,7 @@ void AliTPCCalibKr::Init()
   // 
   // init input tree and output histograms 
   //
   // 
   // init input tree and output histograms 
   //
-
-  // set input tree
-  if(!fTree) { 
-   Printf("ERROR: Could not read chain from input");
-  }
-  else {
-   fTree->SetBranchStatus("*",1); 
-  }
-
-  // set branch address
-  fClusters = new TClonesArray("AliTPCclusterKr");
-
-  if(!fTree->GetBranch("fClusters")) {
-    Printf("ERROR: Could not get fClusters branch from input");
-  } else {
-   fTree->GetBranch("fClusters")->SetAddress(&fClusters);
-  }
-  
   // create output TObjArray
   fHistoKrArray.Clear();
 
   // create output TObjArray
   fHistoKrArray.Clear();
 
@@ -214,68 +176,19 @@ void AliTPCCalibKr::Init()
       TH3F *hist = CreateHisto(i);
       if(hist) fHistoKrArray.AddAt(hist,i);
        }
       TH3F *hist = CreateHisto(i);
       if(hist) fHistoKrArray.AddAt(hist,i);
        }
-
   }
 }
   }
 }
-
-//_____________________________________________________________________
-Bool_t AliTPCCalibKr::ReadEntry(Int_t evt)
-{
-  // 
-  // read entry from the tree
-  //
-  Long64_t centry = fTree->LoadTree(evt);
-  if(centry < 0) return kFALSE;
-
-  if(!fTree->GetBranch("fClusters")) 
-  {
-    Printf("ERROR: Could not get fClusters branch from input");
-       return kFALSE;
-  } else {
-   fTree->GetBranch("fClusters")->SetAddress(&fClusters);
-  }
-
-  fTree->GetEntry(evt);
-
-return kTRUE;
-}
  
 //_____________________________________________________________________
  
 //_____________________________________________________________________
-Bool_t AliTPCCalibKr::Process()
+Bool_t AliTPCCalibKr::Process(AliTPCclusterKr *cluster)
 {
   //
   // process events 
   // call event by event
   //
 
 {
   //
   // process events 
   // call event by event
   //
 
-  // init tree
-  Init();
-
-  // get events
-  if(!fTree) return kFALSE;
-  Int_t nEvents = fTree->GetEntries();
-
-  // fill histograms 
-  for(Int_t i=0; i<nEvents; ++i)
-  {
-    if(ReadEntry(i) == kFALSE) return kFALSE;
-
-    if(!(i%10000)) cout << "evt: " << i << endl; 
-
-    // get TClonesArray entries
-    fClustKr = 0;
-    Int_t entries = fClusters->GetEntries();
-    for(Int_t j=0; j < entries; ++j)
-       {
-         fClustKr = (AliTPCclusterKr*)fClusters->At(j);
-
-      if(fClustKr) Update(fClustKr);
-         else return kFALSE;
-       }
-  }
-
-  // write output 
-  return Terminate();
+  if(cluster) Update(cluster);
+  else return kFALSE;
 }
 
 //_____________________________________________________________________
 }
 
 //_____________________________________________________________________
@@ -340,6 +253,7 @@ Bool_t AliTPCCalibKr::Update(AliTPCclusterKr  *cl)
   return kTRUE;
 }
 
   return kTRUE;
 }
 
+//_____________________________________________________________________
 Bool_t AliTPCCalibKr::Accept(AliTPCclusterKr  *cl){
   //
   // cuts
 Bool_t AliTPCCalibKr::Accept(AliTPCclusterKr  *cl){
   //
   // cuts
@@ -353,13 +267,13 @@ Bool_t AliTPCCalibKr::Accept(AliTPCclusterKr  *cl){
     TCut cutAll = cutR0+cutR1+cutR2+cutR3+cutS1;
   */
   //R0
     TCut cutAll = cutR0+cutR1+cutR2+cutR3+cutS1;
   */
   //R0
-  if (cl->GetADCcluster()/ cl->GetSize() >200)        return kFALSE;
+  if (cl->GetADCcluster()/ Float_t(cl->GetSize()) >200)        return kFALSE;
   // R1
   // R1
-  if (cl->GetADCcluster()/ cl->GetSize() <7)          return kFALSE;
+  if (cl->GetADCcluster()/ Float_t(cl->GetSize()) <7)          return kFALSE;
   //R2
   //R2
-  if (cl->GetMax().GetAdc()/ cl->GetADCcluster() >0.4)  return kFALSE;
+  if (cl->GetMax().GetAdc()/ Float_t(cl->GetADCcluster()) >0.4)  return kFALSE;
   //R3
   //R3
-  if (cl->GetMax().GetAdc()/ cl->GetADCcluster() <0.01) return kFALSE;
+  if (cl->GetMax().GetAdc()/ Float_t(cl->GetADCcluster()) <0.01) return kFALSE;
   //S1
   if (cl->GetSize()>200) return kFALSE;
   if (cl->GetSize()<6)  return kFALSE;
   //S1
   if (cl->GetSize()>200) return kFALSE;
   if (cl->GetSize()<6)  return kFALSE;
@@ -367,46 +281,12 @@ Bool_t AliTPCCalibKr::Accept(AliTPCclusterKr  *cl){
 
 }
 
 
 }
 
-
-
 //_____________________________________________________________________
 TH3F* AliTPCCalibKr::GetHistoKr(Int_t chamber) const
 {
   // get histograms from fHistoKrArray
   return (TH3F*) fHistoKrArray.At(chamber);
 }
 //_____________________________________________________________________
 TH3F* AliTPCCalibKr::GetHistoKr(Int_t chamber) const
 {
   // get histograms from fHistoKrArray
   return (TH3F*) fHistoKrArray.At(chamber);
 }
-
-//_____________________________________________________________________
-Bool_t AliTPCCalibKr::Terminate() 
-{
-  //
-  // store AliTPCCalibKr in the output file 
-  //
-  if(bOutputHisto) {
-    TFile *outFile = new TFile("outHistFile.root","RECREATE"); 
-   
-    if(outFile) 
-       {
-         outFile->cd();
-
-         for(int i=0; i<72; ++i) {
-            if( IsCSide(i) == kTRUE && bCSide == kTRUE)
-              printf("C side chamber: %d, 3D histo entries: %10.f \n",i,((TH3F*)fHistoKrArray.At(i))->GetEntries());
-
-            if( IsCSide(i) == kFALSE && bASide == kTRUE)
-              printf("A side chamber: %d, 3D histo entries: %10.f \n",i,((TH3F*)fHistoKrArray.At(i))->GetEntries());
-         }
-         this->Write();
-         outFile->Close();
-
-         return kTRUE;
-       }
-       else 
-         return kFALSE;
-  }
-
-return kFALSE;
-}
  
 //_____________________________________________________________________
 void AliTPCCalibKr::Analyse() 
  
 //_____________________________________________________________________
 void AliTPCCalibKr::Analyse() 
@@ -498,7 +378,7 @@ void AliTPCCalibKr::Analyse()
         delete gausFit;
         delete projH;
         if (fitMean <= 0) continue;
         delete gausFit;
         delete projH;
         if (fitMean <= 0) continue;
-        printf("[ch%i r%i, p%i] entries = %f, maxEntries = %f, fitMean = %f, fitRMS = %f\n", chamber, iRow, iPad, entries, maxEntries, fitMean, fitRMS);
+        //printf("[ch%i r%i, p%i] entries = %f, maxEntries = %f, fitMean = %f, fitRMS = %f\n", chamber, iRow, iPad, entries, maxEntries, fitMean, fitRMS);
     
         // write the calibration parameters for each pad that the 3d histogram was projected onto
         // (with considering the step size) to the CalPads
     
         // write the calibration parameters for each pad that the 3d histogram was projected onto
         // (with considering the step size) to the CalPads
@@ -582,3 +462,41 @@ TH1D* AliTPCCalibKr::ProjectHisto(TH3F* histo3D, const char* name, Int_t xMin, I
   projH->SetEntries((Long64_t)entries);
   return projH;
 }
   projH->SetEntries((Long64_t)entries);
   return projH;
 }
+
+//_____________________________________________________________________
+Long64_t AliTPCCalibKr::Merge(TCollection* list) {
+// merge function 
+//
+cout << "Merge " << endl;
+
+if (!list)
+return 0;
+
+if (list->IsEmpty())
+return 1;
+
+TIterator* iter = list->MakeIterator();
+TObject* obj = 0;
+
+    iter->Reset();
+    Int_t count=0;
+       while((obj = iter->Next()) != 0)
+       {
+         AliTPCCalibKr* entry = dynamic_cast<AliTPCCalibKr*>(obj);
+         if (entry == 0) continue; 
+
+               for(int i=0; i<72; ++i) { 
+                 if(IsCSide(i) == kTRUE && bCSide == kTRUE) { 
+                 ((TH3F*)fHistoKrArray.At(i))->Add( ((TH3F*)entry->fHistoKrArray.At(i)) );  
+                 } 
+
+                 if(IsCSide(i) == kFALSE && bASide == kTRUE) { 
+                   ((TH3F*)fHistoKrArray.At(i))->Add( ((TH3F*)entry->fHistoKrArray.At(i)) );  
+                 } 
+               } 
+
+         count++;
+       }
+
+return count;
+}
index 43d65e5cfdb525333660e552083e571798692335..7425eafa7bdfbc6fab2b2d6df2144f1ac6207899 100644 (file)
@@ -25,8 +25,7 @@ public:
 
   //
   void Init();
 
   //
   void Init();
-  Bool_t Process();
-  Bool_t ReadEntry(Int_t evt);
+  Bool_t Process(AliTPCclusterKr *cluster);
   Bool_t Accept(AliTPCclusterKr *cluster);
   Bool_t Update(AliTPCclusterKr *cluster);
   TH3F*  CreateHisto(Int_t chamber);
   Bool_t Accept(AliTPCclusterKr *cluster);
   Bool_t Update(AliTPCclusterKr *cluster);
   TH3F*  CreateHisto(Int_t chamber);
@@ -37,25 +36,21 @@ public:
   Bool_t IsCSide(Int_t chamber);
   Bool_t IsIROC(Int_t chamber);
 
   Bool_t IsCSide(Int_t chamber);
   Bool_t IsIROC(Int_t chamber);
 
-  Bool_t Terminate();
   void Analyse();
   static TH1D* ProjectHisto(TH3F* histo3D, const char* name = "_pz", Int_t firstxbin = 0, Int_t lastxbin = 0, Int_t firstybin = 0, Int_t lastybin = 0);
 
   void Analyse();
   static TH1D* ProjectHisto(TH3F* histo3D, const char* name = "_pz", Int_t firstxbin = 0, Int_t lastxbin = 0, Int_t firstybin = 0, Int_t lastybin = 0);
 
-  void SetInputChain(TChain *inChain) {fTree = (TTree*) inChain;}
-  void SetOutputHisto(Bool_t bHisto = kTRUE) {bOutputHisto = bHisto;}
-
   void SetASide(Bool_t bA = kTRUE) {bASide = bA;} // fill histo only A TPC side
   void SetBSide(Bool_t bC = kTRUE) {bCSide = bC;} // fill histo only C TPC side
 
   void SetASide(Bool_t bA = kTRUE) {bASide = bA;} // fill histo only A TPC side
   void SetBSide(Bool_t bC = kTRUE) {bCSide = bC;} // fill histo only C TPC side
 
+  //Merge output objects (needed by PROOF)
+  virtual Long64_t Merge(TCollection* list);
+  
 private:
 
   Bool_t bOutputHisto;        //!
   Bool_t bASide;              //!
   Bool_t bCSide;              //!
 
 private:
 
   Bool_t bOutputHisto;        //!
   Bool_t bASide;              //!
   Bool_t bCSide;              //!
 
-  TClonesArray *fClusters;    //!
-  AliTPCclusterKr *fClustKr;  //!
-  TTree *fTree;               //!
   TObjArray fHistoKrArray;    //  Calibration histograms for Kr distribution
 
 public:
   TObjArray fHistoKrArray;    //  Calibration histograms for Kr distribution
 
 public: