* provided "as is" without express or implied warranty. *
**************************************************************************/
-// Macro MUONTrigger.C (TO BE COMPILED)
-// for testing the C++ trigger code
+// This macro is to be used to check the trigger algorithm w/o having to
+// (re-)perform simulation and digitalization.
+// It loads the digits, erase TreeR and store the current trigger output in
+// TreeR.
+// The different trigger outputs can be compared by looking at the GLT branch
+// of TreeD (filled during simulation) and the TC branch of TreeR (filled from
+// a copy of TreeD during reconstruction or with this macro).
+// Note: rec points from tracking chamber will be lost.
+//
+// usage: (to be compiled)
+// MUONTrigger("galice.root",0) -> Default Trigger Code
+// MUONTrigger("galice.root",1) -> New Trigger Code
#if !defined(__CINT__) || defined(__MAKECINT__)
#include "AliRun.h"
#include "AliMUONLoader.h"
#include "AliMUONData.h"
#include "AliMUONDigit.h"
-#include "AliMUONTriggerDecision.h"
+#include "AliMUONTriggerDecisionV1.h"
+#include "AliMUONTriggerElectronics.h"
+#include "AliMUONCalibrationData.h"
+#include "AliCDBManager.h"
#endif
-void MUONTrigger(char * FileName="galice.root", Int_t nevents=1, Int_t idebug=2)
+void MUONTrigger(char * FileName="galice.root", Int_t NewTriggerCode=0)
{
// Creating Run Loader and openning file containing Digits
AliRunLoader * RunLoader = AliRunLoader::Open(FileName,"MUONLoader","UPDATE");
if (RunLoader ==0x0) {
- printf(">>> Error : Error Opening %s file \n",FileName);
- return;
+ printf(">>> Error : Error Opening %s file \n",FileName);
+ return;
}
// Loading AliRun master
if (RunLoader->GetAliRun() == 0x0) RunLoader->LoadgAlice();
gAlice = RunLoader->GetAliRun();
// Loading MUON subsystem
-// AliMUONLoader *MUONLoader = (AliMUONLoader*) RunLoader->GetLoader("MUONLoader");
AliLoader * MUONLoader = RunLoader->GetLoader("MUONLoader");
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");
-
+
// Creating MUONTriggerDecision
- AliMUONTriggerDecision* decision = new AliMUONTriggerDecision(MUONLoader ,idebug,MUONData);
-
- Int_t nEvents = RunLoader->GetNumberOfEvents();
- AliMUONDigit * mDigit;
-
- Int_t tracks[10];
- Int_t charges[10];
+ TTask *TriggerProcessor;
+ if (NewTriggerCode == 0) {
+ cout << " using default trigger code " << "\n";
+ TriggerProcessor = new AliMUONTriggerDecisionV1(MUONData);
+ } else {
+ cout << " using new trigger code " << "\n";
+ 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");
+ }
+ }
+
+ AliMUONDigit * mDigit;
Int_t digits[7];
for(Int_t ievent = 0; ievent < nevents; ievent++) {
RunLoader->GetEvent(ievent);
MUONData->SetTreeAddress("D");
- for(Int_t icathode=0; icathode<2; icathode++) {
- MUONData->GetCathode(icathode);
- for(Int_t ichamber=10; ichamber<14; ichamber++) {
- Int_t idigit, ndigits;
- ndigits = (Int_t) MUONData->Digits(ichamber)->GetEntriesFast();
-// printf(">>> Chamber Cathode ndigits %d %d %d\n",ichamber,icathode,ndigits);
- for(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->Signal();
- digits[4] = mDigit->Physics();
- digits[5] = mDigit->Hit();
- digits[6] = mDigit->DetElemId();
-/*
- Int_t nptracks = mDigit->GetNTracks();
- for (Int_t i = 0; i < nptracks; i++) {
- tracks[i] = mDigit->GetTrack(i);
- charges[i] = mDigit->GetCharge(i);
- }
-*/
-// printf("ichamber icathode ix iy %d %d %d %d \n",ichamber,icathode,mDigit->PadX(),mDigit->PadY());
-
- decision->AddDigit(ichamber, tracks, charges, digits);
- } // loop on digits
- } // loop on chambers
- } // loop on cathodes
- MUONData->ResetDigits();
- decision->Trigger();
- decision->ClearDigits();
+ MUONData->GetDigits();
+ for(Int_t ichamber=10; ichamber<14; ichamber++) {
+ Int_t idigit, ndigits;
+ ndigits = (Int_t) MUONData->Digits(ichamber)->GetEntriesFast();
+
+ for(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->Signal();
+ 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
+
+ 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";
+ }
+ }
+
+ MUONData->MakeBranch("TC");
+ MUONData->SetTreeAddress("TC");
+ TriggerProcessor->ExecuteTask();
+
+ MUONData->Fill("TC");
+ MUONLoader->WriteRecPoints("OVERWRITE");
+ MUONData->ResetDigits();
+
} // loop on events
MUONLoader->UnloadDigits();
+ MUONLoader->UnloadRecPoints();
}
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff