MUONTrigger.C restored thanks to Laurent (Philippe C.)
authorpcrochet <pcrochet@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 26 Oct 2007 16:16:43 +0000 (16:16 +0000)
committerpcrochet <pcrochet@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 26 Oct 2007 16:16:43 +0000 (16:16 +0000)
MUON/MUONTrigger.C
MUON/READMEtrigger.txt
MUON/loadmacros.C

index 6fb68fd..91a2519 100644 (file)
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include "AliRun.h"
 #include "AliMUON.h"
-#include "AliMUONLoader.h"
-#include "AliMUONData.h"
 #include "AliMUONDigit.h"
 #include "AliMUONTriggerElectronics.h"
 #include "AliMUONCalibrationData.h"
 #include "AliCDBManager.h"
+#include "AliMUONDataInterface.h"
+#include "AliMUONMCDataInterface.h"
+#include "AliMUONVTriggerStore.h"
+#include "AliMUONDigitStoreV1.h"
 #include <TClonesArray.h>
+#include "AliMpCDB.h"
+#include <TFile.h>
 #endif
-void MUONTrigger(char * FileNameSim="galice_sim.root", char * FileName="galice.root")
+
+void MUONTrigger(const char* filename)
 {
     // Creating Run Loader and openning file containing Digits 
-    AliRunLoader * RunLoader = AliRunLoader::Open(FileName,"MUONLoader","UPDATE");
+    AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONLoader","UPDATE");
     if (RunLoader ==0x0) {
-        printf(">>> Error : Error Opening %s file \n",FileName);
+        printf(">>> Error : Error Opening %s file \n",filename);
         return;
     }
     // Loading AliRun master
@@ -44,82 +49,96 @@ void MUONTrigger(char * FileNameSim="galice_sim.root", char * FileName="galice.r
     gAlice = RunLoader->GetAliRun();
     
     // Loading MUON subsystem
-    AliLoader * MUONLoader = RunLoader->GetLoader("MUONLoader");
+    AliLoader* MUONLoader = RunLoader->GetDetectorLoader("MUON");
     MUONLoader->LoadDigits("READ");
-    MUONLoader->LoadRecPoints("UPDATE"); // absolutely essential !!!
-    
-    Int_t nevents;
-    nevents = RunLoader->GetNumberOfEvents();
-    
-    // Creating MUON data container
-    AliMUONData* MUONData = new AliMUONData(MUONLoader,"MUON","MUON");
+    MUONLoader->LoadRecPoints("UPDATE"); // absolutely essential !!!    
     
     // Creating MUONTriggerDecision
-    TTask *TriggerProcessor;
     AliCDBManager* cdbManager = AliCDBManager::Instance();
     cdbManager->SetDefaultStorage("local://$ALICE_ROOT");
-    Int_t runnumber = gAlice->GetRunNumber();
-    AliMUONCalibrationData *CalibrationData = new AliMUONCalibrationData(runnumber);
-    TriggerProcessor = new AliMUONTriggerElectronics(MUONData,CalibrationData);
     
-    // Testing if Trigger has already been done
-    RunLoader->GetEvent(0);    
-    if (MUONLoader->TreeR()) {
-       if (MUONData->IsTriggerBranchesInTree()) {
-           MUONLoader->UnloadRecPoints();
-           MUONLoader->LoadRecPoints("RECREATE");
-           printf("Recreating recpoints files\n");
-       }
-    }
+    Int_t runnumber = cdbManager->GetRun();
+    AliMpCDB::LoadDDLStore();
     
-    AliMUONDigit * mDigit;    
-    Float_t digits[7];
+    AliMUONCalibrationData *CalibrationData = new AliMUONCalibrationData(runnumber);
+    AliMUONTriggerElectronics *TriggerProcessor        = new AliMUONTriggerElectronics(CalibrationData);
+    
+    Int_t nevents = RunLoader->GetNumberOfEvents();
+    AliMUONVDigitStore* digitStore=0x0;
+    AliMUONVTriggerStore* triggerStore=0x0;
     
     for(Int_t ievent = 0; ievent < nevents; ievent++) {
        printf(">>> Event %i out of %i \n",ievent,nevents);
-       RunLoader->GetEvent(ievent);
-       MUONData->SetTreeAddress("D");
+       RunLoader->GetRunLoader()->GetEvent(ievent);
        
-       MUONData->GetDigits();
-       for(Int_t ichamber=10; ichamber<14; ichamber++) {           
-           Int_t ndigits = (Int_t) MUONData->Digits(ichamber)->GetEntriesFast();
-           for(Int_t idigit=0; idigit<ndigits; idigit++) {
-               mDigit = static_cast<AliMUONDigit*>(MUONData->Digits(ichamber)->At(idigit));             
-               
-               digits[0] = mDigit->PadX();
-               digits[1] = mDigit->PadY();
-               digits[2] = mDigit->Cathode();
-               digits[3] = mDigit->Charge();
-               digits[4] = mDigit->Physics();
-               digits[5] = mDigit->Hit();
-               digits[6] = mDigit->DetElemId();
-               
-//             printf("ichamber ix iy %d %d %d \n",ichamber,mDigit->PadX(),mDigit->PadY());
-               
-           } // loop on digits
-       } // loop on chambers
+       MUONLoader->LoadRecPoints("update");
+       MUONLoader->CleanRecPoints();
+       MUONLoader->MakeRecPointsContainer();
+       TTree* clustersTree = MUONLoader->TreeR();
+       TFile* cfile = clustersTree->GetCurrentFile();
+       if ( !cfile ) 
+       {
+           cout << " could not find Cluster file " << endl;
+           return;
+       }
        
+       MUONLoader->LoadDigits("read");
+       TTree* digitsTree = MUONLoader->TreeD();
+       TFile* dfile = digitsTree->GetCurrentFile();
+       if ( !dfile ) 
+       {
+           cout << " could not find Digit file " << endl;
+           return;
+       }
        
-       if (MUONLoader->TreeR() == 0x0) {       
-           MUONLoader->MakeRecPointsContainer();
-       } else {
-           if (MUONData->IsTriggerBranchesInTree()){ 
-               if (ievent==0) MUONLoader->UnloadRecPoints();
-               MUONLoader->MakeRecPointsContainer();
-               cout << "Recreating RecPointsContainer and deleting previous ones" << "\n";
+// here start reconstruction   
+       if (!digitStore) digitStore = AliMUONVDigitStore::Create(*digitsTree);  
+       if (!triggerStore) triggerStore = AliMUONVTriggerStore::Create(*digitsTree);
+       // insure we start with empty stores
+       if ( digitStore ) 
+       {
+           digitStore->Clear(); 
+           Bool_t alone = ( triggerStore ? kFALSE : kTRUE );
+           Bool_t ok = digitStore->Connect(*digitsTree,alone);
+           if (!ok)
+           {
+               cerr << "Could not connect digitStore to digitsTree \n";
+               return;
            }
-       }       
+       } else {
+           cerr << "digitStore does not exist " << "\n";
+           return;
+       }
+       
+       digitsTree->GetEvent(0);
        
-       MUONData->MakeBranch("TC");     
-       MUONData->SetTreeAddress("TC");
-       TriggerProcessor->ExecuteTask();
+// process trigger response
+       TriggerProcessor->Digits2Trigger(*digitStore,*triggerStore);
+       
+       //triggerStore->Print();
+       
+       Bool_t ok(kFALSE);
+       if ( triggerStore ) {
+           ok = triggerStore->Connect(*clustersTree,kTRUE);
+           if (!ok)
+           {
+               cerr << "Could not create triggerStore branches in TreeR " << "\n";
+               return;
+           }
+       } else {
+           cerr << "triggerStore does not exist " << "\n";
+           return;
+       }
+
+// fill TreeR
+       clustersTree->Fill();
+       MUONLoader->UnloadDigits();
+       MUONLoader->WriteRecPoints("OVERWRITE");
+       MUONLoader->UnloadRecPoints();
 
-       MUONData->Fill("TC");
-       MUONLoader->WriteRecPoints("OVERWRITE");  
-        MUONData->ResetDigits();
+    }  // loop on events
 
-    } // loop on events
-    MUONLoader->UnloadDigits();
-    MUONLoader->UnloadRecPoints();
 }
 
+
+
index 65c69eb..0dcc6f0 100644 (file)
@@ -20,7 +20,7 @@ Note: rec points from tracking chamber will be lost.
 Usage:
 <pre>
 root [0] .L $ALICE_ROOT/MUON/MUONTrigger.C+
-root [1] MUONTrigger()
+root [1] MUONTrigger("galice.root")
 </pre>
 
 
index 110efe5..8fe3d0b 100644 (file)
@@ -123,7 +123,7 @@ void loadmacros ()
   gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONRecoCheck.C++");           
   // gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONResoEffChamber.C++"); 
   gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONStatusMap.C++");        
-  // gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONTrigger.C++");        
+  gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONTrigger.C++");        
   gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONTriggerEfficiency.C++");   
   gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONTriggerEfficiencyPt.C++"); 
   gROOT->LoadMacro("$ALICE_ROOT/MUON/MUONTriggerChamberEfficiency.C++");