1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 #include "AliMUONDataInterface.h"
19 #include "AliMUONGeometryTransformer.h"
20 #include "AliMUONVDigit.h"
21 #include "AliMUONVCluster.h"
22 #include "AliMUONLocalTrigger.h"
23 #include "AliMUONRegionalTrigger.h"
24 #include "AliMUONGlobalTrigger.h"
25 #include "AliMUONTriggerCircuit.h"
26 #include "AliMUONVClusterStore.h"
27 #include "AliMUONVDigitStore.h"
28 #include "AliMUONVTriggerStore.h"
31 #include "AliMpEncodePair.h"
32 #include "AliMpDEManager.h"
33 #include "AliMpConstants.h"
36 #include "AliLoader.h"
37 #include "AliRunLoader.h"
38 #include "AliHeader.h"
39 #include "AliCDBManager.h"
41 #include <Riostream.h>
46 #include <TIterator.h>
50 //-----------------------------------------------------------------------------
51 /// \class AliMUONDataInterface
53 /// An easy to use interface to the MUON data data stored in
54 /// TreeS, TreeD and TreeR.
56 /// For MC related information (i.e. TreeH, TreeK, TreeTR), see
57 /// AliMUONMCDataInterface.
60 /// This interface in not necessarily the fastest way to fetch the data but
61 /// it is the easiest.
63 /// \author Laurent Aphecetche, Subatech & Artur Szostak <artursz@iafrica.com> (University of Cape Town)
64 //-----------------------------------------------------------------------------
67 ClassImp(AliMUONDataInterface)
71 Int_t AliMUONDataInterface::fgInstanceCounter(0);
73 //______________________________________________________________________________
74 AliMUONDataInterface::AliMUONDataInterface(const char* filename)
83 fCurrentIteratorType(kNoIterator),
90 /// @param filename should be the full path to a valid galice.root file
94 if ( AliCDBManager::Instance() != NULL &&
95 AliCDBManager::Instance()->GetDefaultStorage() == NULL ) {
96 AliFatal("CDB default storage not defined.");
102 if ( ! AliMpCDB::LoadDDLStore() ) {
103 AliFatal("Could not access mapping from OCDB !");
107 //______________________________________________________________________________
108 AliMUONDataInterface::~AliMUONDataInterface()
112 if ( fLoader != 0x0 )
114 delete fLoader->GetRunLoader();
119 //______________________________________________________________________________
121 AliMUONDataInterface::DigitStore(Int_t event)
123 /// Return digitStore for a given event.
124 /// Return 0x0 if event not found.
125 /// Returned pointer should not be deleted
127 /// \note If a previous store has been retrieved by one of the methods of
128 /// this class, but for a different event number, then those stores will
129 /// be deleted and no longer valid.
130 /// If you require access to the data for the earlier retrieved store,
131 /// but for different events, then you should deep copy / clone the object.
133 if (not IsValid()) return 0x0;
135 if (event == fCurrentEvent)
137 if (fDigitStore != 0x0)
143 if ( not LoadEvent(event) ) return 0x0;
146 fLoader->LoadDigits();
148 TTree* treeD = fLoader->TreeD();
151 AliError("Could not get treeD");
155 fDigitStore = AliMUONVDigitStore::Create(*treeD);
156 if ( fDigitStore != 0x0 )
158 fDigitStore->Clear();
159 fDigitStore->Connect(*treeD);
163 fLoader->UnloadDigits();
168 //______________________________________________________________________________
169 AliMUONVClusterStore*
170 AliMUONDataInterface::ClusterStore(Int_t event)
172 /// Return clusterStore for a given event.
173 /// Return 0x0 if event not found.
174 /// Returned pointer should not be deleted
176 /// \note If a previous store has been retrieved by one of the methods of
177 /// this class, but for a different event number, then those stores will
178 /// be deleted and no longer valid.
179 /// If you require access to the data for the earlier retrieved store,
180 /// but for different events, then you should deep copy / clone the object.
182 if (not IsValid()) return 0x0;
184 if (event == fCurrentEvent)
186 if (fClusterStore != 0x0)
187 return fClusterStore;
192 if ( not LoadEvent(event) ) return 0x0;
195 fLoader->LoadRecPoints();
197 TTree* treeR = fLoader->TreeR();
200 AliError("Could not get treeR");
204 fClusterStore = AliMUONVClusterStore::Create(*treeR);
205 if ( fClusterStore != 0x0 )
207 fClusterStore->Clear();
208 fClusterStore->Connect(*treeR);
212 fLoader->UnloadRecPoints();
214 return fClusterStore;
217 //_____________________________________________________________________________
218 AliMUONVTriggerStore*
219 AliMUONDataInterface::TriggerStore(Int_t event, const char* treeLetter)
221 /// Return the triggerStore for a given event.
222 /// Return 0x0 if event not found.
223 /// Returned pointer should not be deleted
224 /// treeLetter can be R or D to tell from which tree to read the information
226 /// \note If a previous store has been retrieved by one of the methods of
227 /// this class, but for a different event number, then those stores will
228 /// be deleted and no longer valid.
229 /// If you require access to the data for the earlier retrieved store,
230 /// but for different events, then you should deep copy / clone the object.
232 if (not IsValid()) return 0x0;
234 if (event == fCurrentEvent)
236 if (fTreeLetter == treeLetter)
238 if (fTriggerStore != 0x0)
239 return fTriggerStore;
243 // Reset only the fTriggerStore since the others might still be valid
244 // for the same event.
245 if (fTriggerStore != 0x0)
247 delete fTriggerStore;
254 // Event has changed so reset all the stores.
256 if ( not LoadEvent(event) ) return 0x0;
261 TString stree(treeLetter);
266 fLoader->LoadDigits();
267 tree = fLoader->TreeD();
269 else if ( stree == "R" )
271 fLoader->LoadRecPoints();
272 tree = fLoader->TreeR();
277 AliError(Form("Could not get tree%s",treeLetter));
281 fTriggerStore = AliMUONVTriggerStore::Create(*tree);
282 if ( fTriggerStore != 0x0 )
284 fTriggerStore->Clear();
285 fTriggerStore->Connect(*tree);
291 fLoader->UnloadDigits();
293 else if ( stree == "R" )
295 fLoader->UnloadRecPoints();
299 return fTriggerStore;
302 //______________________________________________________________________________
304 AliMUONDataInterface::DumpDigits(Int_t event, Bool_t sorted)
306 /// Dump the digits for a given event, sorted if so required
308 if ( fDigitStore != 0x0 )
312 DumpSorted(*fDigitStore);
316 fDigitStore->Print();
321 //______________________________________________________________________________
323 AliMUONDataInterface::DumpRecPoints(Int_t event, Bool_t sorted)
325 /// Dump the recpoints for a given event, sorted if so required
327 if ( fClusterStore != 0x0 )
331 DumpSorted(*fClusterStore);
335 fClusterStore->Print();
340 //_____________________________________________________________________________
342 AliMUONDataInterface::DumpSorted(const AliMUONVStore& store) const
344 /// Dump the given store, in sorted order
346 TIter next(store.CreateIterator());
349 list.SetOwner(kFALSE);
351 while ( ( object = next() ) )
361 //_____________________________________________________________________________
363 AliMUONDataInterface::DumpTrigger(Int_t event, const char* treeLetter)
365 /// Dump trigger for a given event from a given tree (if event>=0)
366 /// or loop over all events and build a trigger ntuple if event<0
367 /// treeLetter can be R or D to tell from which tree to read the information
371 NtupleTrigger(treeLetter);
375 TriggerStore(event,treeLetter);
377 if ( fTriggerStore != 0x0 )
379 fTriggerStore->Print();
384 //_____________________________________________________________________________
386 AliMUONDataInterface::NtupleTrigger(const char* treeLetter)
388 //// Loop over events to build trigger ntuples
391 TString sTreeLetter(treeLetter);
392 sTreeLetter.ToUpper();
394 if ( sTreeLetter != "R" && sTreeLetter != "D" )
396 AliError(Form("Cannot handle tree%s. Use D or R",treeLetter));
401 TNtuple tupleGlo("TgtupleGlo","Global Trigger Ntuple",
402 "ev:slpt:shpt:uplpt:uphpt:lplpt:lplpt");
403 TNtuple tupleLoc("TgtupleLoc","Local Trigger Ntuple",
404 "ev:LoCircuit:LoStripX:LoDev:StripY:LoLpt:LoHpt:y11:y21:x11");
406 // initialize counters
414 AliMUONGeometryTransformer transformer;
415 transformer.LoadGeometryData(Form("%s/geometry.root",
416 gSystem->DirName(fLoader->GetRunLoader()->GetFileName())));
417 AliMUONTriggerCircuit triggerCircuit(&transformer);
419 // select output file name from selected Tree
420 Char_t fileNameOut[30];
421 if (sTreeLetter == "D")
423 AliInfo(Form("reading from Digits\n"));
424 sprintf(fileNameOut,"TriggerCheckFromDigits.root");
426 else if (sTreeLetter == "R")
428 AliInfo(Form("reading from RecPoints\n"));
429 sprintf(fileNameOut,"TriggerCheckFromRP.root");
433 Int_t nevents = NumberOfEvents();
435 for (Int_t ievent=0; ievent<nevents; ++ievent)
437 if (ievent%100==0) AliInfo(Form("Processing event %d\n",ievent));
439 AliMUONVTriggerStore* triggerStore = TriggerStore(ievent,treeLetter);
443 AliError(Form("Could not read %s from tree%s","Trigger",treeLetter));
447 // get global trigger info
448 AliMUONGlobalTrigger* gloTrg = triggerStore->Global();
449 sLowpt+=gloTrg->SingleLpt();
450 sHighpt+=gloTrg->SingleHpt();
451 uSLowpt+=gloTrg->PairUnlikeLpt();
452 uSHighpt+=gloTrg->PairUnlikeHpt();
453 lSLowpt+=gloTrg->PairLikeLpt();
454 lSHighpt+=gloTrg->PairLikeHpt();
456 // loop on local triggers
457 TIter next(triggerStore->CreateIterator());
458 AliMUONLocalTrigger* locTrg(0x0);
459 while ( ( locTrg = static_cast<AliMUONLocalTrigger*>(next()) ) )
461 Bool_t xTrig=locTrg->IsTrigX();
462 Bool_t yTrig=locTrg->IsTrigY();
465 { // fill ntuple if trigger in X and Y
466 tupleLoc.Fill(ievent,locTrg->LoCircuit(),
472 triggerCircuit.GetY11Pos(locTrg->LoCircuit(),locTrg->LoStripX()),
473 triggerCircuit.GetY21Pos(locTrg->LoCircuit(),locTrg->LoStripX()+locTrg->LoDev()+1),
474 triggerCircuit.GetX11Pos(locTrg->LoCircuit(),locTrg->LoStripY()));
476 tupleGlo.Fill(ievent,gloTrg->SingleLpt(),gloTrg->SingleHpt(),
477 gloTrg->PairUnlikeLpt(),gloTrg->PairUnlikeHpt(),
478 gloTrg->PairLikeLpt(),gloTrg->PairLikeHpt());
479 } // end of loop on local triggers
480 } // end of loop on events
482 // print info and store ntuples
484 printf("=============================================\n");
485 printf("================ SUMMARY ==================\n");
487 printf("Total number of events processed %d \n",nevents);
489 printf(" Global Trigger output Low pt High pt\n");
490 printf(" number of Single :\t");
491 printf("%i\t%i\t",sLowpt,sHighpt);
493 printf(" number of UnlikeSign pair :\t");
494 printf("%i\t%i\t",uSLowpt,uSHighpt);
496 printf(" number of LikeSign pair :\t");
497 printf("%i\t%i\t",lSLowpt,lSHighpt);
499 printf("=============================================\n");
502 TFile myFile(fileNameOut, "RECREATE");
508 //_____________________________________________________________________________
510 AliMUONDataInterface::LoadEvent(Int_t event)
512 /// Load event if different from the current one.
513 /// Returns kFALSE on error and kTRUE if the event was loaded.
517 AliDebug(1,Form("Loading event %d using runLoader %p",event,fLoader->GetRunLoader()));
518 if (fLoader->GetRunLoader()->GetEvent(event) == 0)
520 fCurrentEvent = event;
527 //______________________________________________________________________________
529 AliMUONDataInterface::NumberOfEvents() const
531 /// Number of events in the current galice.root file we're attached to
532 if (not IsValid()) return -1;
533 return fLoader->GetRunLoader()->GetNumberOfEvents();
536 //_____________________________________________________________________________
538 AliMUONDataInterface::Open(const char* filename)
540 /// Connect to a given galice.root file
546 if ( fLoader != 0x0 )
548 delete fLoader->GetRunLoader();
555 TString foldername(Form("%s-%d",ClassName(),fgInstanceCounter));
557 while (AliRunLoader::GetRunLoader(foldername) != 0x0)
559 delete AliRunLoader::GetRunLoader(foldername);
562 AliRunLoader* runLoader = AliRunLoader::Open(filename,foldername);
563 if (runLoader == 0x0)
565 AliError(Form("Cannot open file %s",filename));
569 runLoader->LoadHeader();
570 if ( ! runLoader->GetHeader() ) {
571 AliError("Cannot load header.");
575 Int_t runNumber = runLoader->GetHeader()->GetRun();
576 AliCDBManager::Instance()->SetRun(runNumber);
578 runLoader->UnloadHeader();
580 fLoader = runLoader->GetDetectorLoader("MUON");
583 AliError("Cannot get AliMUONLoader");
589 AliError(Form("Could not access %s filename. Object is unuseable",filename));
593 //_____________________________________________________________________________
594 Bool_t AliMUONDataInterface::GetEvent(Int_t event)
596 /// Loads all reconstructed data for the given event.
598 if (DigitStore(event) == 0x0) return kFALSE;
599 if (ClusterStore(event) == 0x0) return kFALSE;
600 if (TriggerStore(event) == 0x0) return kFALSE;
604 //_____________________________________________________________________________
605 Int_t AliMUONDataInterface::NumberOfDigits(Int_t detElemId)
607 /// Returns the number of digits to be found on a given detector element.
608 /// @param detElemId The detector element ID number to search on.
610 TIterator* iter = GetIterator(kDigitIteratorByDetectorElement, detElemId);
611 return CountObjects(iter);
614 //_____________________________________________________________________________
615 AliMUONVDigit* AliMUONDataInterface::Digit(Int_t detElemId, Int_t index)
617 /// Returns the a pointer to the index'th digit on the specified detector element.
618 /// @param detElemId The detector element ID number to search on.
619 /// @param index The index number of the digit to fetch in the range [0 .. N-1],
620 /// where N = NumberOfDigits(detElemId)
622 TIterator* iter = GetIterator(kDigitIteratorByDetectorElement, detElemId);
623 return static_cast<AliMUONVDigit*>( FetchObject(iter, index) );
626 //_____________________________________________________________________________
627 Int_t AliMUONDataInterface::NumberOfDigits(Int_t chamber, Int_t cathode)
629 /// Returns the number of digits to be found on a specific chamber and cathode.
630 /// @param chamber The chamber number in the range [0 .. 13].
631 /// @param cathode The cathode in the range [0 .. 1], where 0 is the bending and
632 /// 1 is the non-bending plane.
634 TIterator* iter = GetIterator(kDigitIteratorByChamberAndCathode, chamber, cathode);
635 return CountObjects(iter);
638 //_____________________________________________________________________________
639 AliMUONVDigit* AliMUONDataInterface::Digit(Int_t chamber, Int_t cathode, Int_t index)
641 /// Returns the a pointer to the index'th digit on the specified chamber and cathode.
642 /// @param chamber The chamber number in the range [0 .. 13].
643 /// @param cathode The cathode in the range [0 .. 1], where 0 is the bending and
644 /// 1 is the non-bending plane.
645 /// @param index The index number of the digit to fetch in the range [0 .. N-1],
646 /// where N = NumberOfDigits(chamber, cathode)
648 TIterator* iter = GetIterator(kDigitIteratorByChamberAndCathode, chamber, cathode);
649 return static_cast<AliMUONVDigit*>( FetchObject(iter, index) );
652 //_____________________________________________________________________________
653 Int_t AliMUONDataInterface::NumberOfRawClusters(Int_t chamber)
655 /// Returns the number of reconstructed raw clusters on the specified chamber.
656 /// @param chamber The chamber number in the range [0 .. 13].
658 TIterator* iter = GetIterator(kRawClusterIterator, chamber);
659 return CountObjects(iter);
662 //_____________________________________________________________________________
663 AliMUONVCluster* AliMUONDataInterface::RawCluster(Int_t chamber, Int_t index)
665 /// Returns a pointer to the index'th raw cluster on the specified chamber.
666 /// @param chamber The chamber number in the range [0 .. 13].
667 /// @param index The index number of the raw cluster to fetch in the range [0 .. N-1],
668 /// where N = NumberOfRawClusters(chamber)
670 TIterator* iter = GetIterator(kRawClusterIterator, chamber);
671 return static_cast<AliMUONVCluster*>( FetchObject(iter, index) );
674 //_____________________________________________________________________________
675 Int_t AliMUONDataInterface::NumberOfLocalTriggers()
677 /// Returns the number of reconstructed local trigger objects.
679 TIterator* iter = GetIterator(kLocalTriggerIterator);
680 return CountObjects(iter);
683 //_____________________________________________________________________________
684 AliMUONLocalTrigger* AliMUONDataInterface::LocalTrigger(Int_t index)
686 /// Returns a pointer to the index'th local trigger object.
687 /// @param index The index number of the local trigger object to fetch in the range [0 .. N-1],
688 /// where N = NumberOfLocalTriggers()
690 TIterator* iter = GetIterator(kLocalTriggerIterator);
691 return static_cast<AliMUONLocalTrigger*>( FetchObject(iter, index) );
694 //_____________________________________________________________________________
695 Int_t AliMUONDataInterface::NumberOfRegionalTriggers()
697 /// Returns the number of regional trigger objects reconstructed.
699 TIterator* iter = GetIterator(kRegionalTriggerIterator);
700 return CountObjects(iter);
703 //_____________________________________________________________________________
704 AliMUONRegionalTrigger* AliMUONDataInterface::RegionalTrigger(Int_t index)
706 /// Returns a pointer to the index'th regional trigger object.
707 /// @param index The index number of the regional trigger object to fetch in the range [0 .. N-1],
708 /// where N = NumberOfRegionalTriggers()
710 TIterator* iter = GetIterator(kRegionalTriggerIterator);
711 return static_cast<AliMUONRegionalTrigger*>( FetchObject(iter, index) );
714 //_____________________________________________________________________________
715 AliMUONGlobalTrigger* AliMUONDataInterface::GlobalTrigger()
717 /// Returns a pointer to the reconstructed global trigger object for the event.
719 AliMUONVTriggerStore* store = TriggerStore(fCurrentEvent);
720 if (store == 0x0) return 0x0;
721 return store->Global();
724 //_____________________________________________________________________________
725 void AliMUONDataInterface::ResetStores()
727 /// Deletes all the store objects that have been created and resets the pointers to 0x0.
728 /// The temporary iterator object is automatically reset. See ResetIterator for more details.
731 if (fDigitStore != 0x0)
736 if (fTriggerStore != 0x0)
738 delete fTriggerStore;
741 if (fClusterStore != 0x0)
743 delete fClusterStore;
748 //_____________________________________________________________________________
749 TIterator* AliMUONDataInterface::GetIterator(IteratorType type, Int_t x, Int_t y)
751 /// Creates an appropriate iterator object and returns it.
752 /// If the iterator has already been created then that one is returned otherwise
753 /// a new object is created.
754 /// Depending on the value of 'type' the semantics of parameters x and y can change.
755 /// @param type The type of iterator to create.
756 /// @param x This is the detector element ID if type == kDigitIteratorByDetectorElement
757 /// If type equals kDigitIteratorByChamberAndCathode or kRawClusterIterator
758 /// then this is the chamber number. In all other cases this parameter is
760 /// @param y If type == kDigitIteratorByChamberAndCathode then this parameter is the
761 /// cathode number. In all other cases this parameter is
764 if (type == fCurrentIteratorType and fDataX == x and fDataY == y)
767 if (fCurrentEvent == -1)
769 AliError("No event was selected. Try first using GetEvent().");
777 case kDigitIteratorByDetectorElement:
780 AliMUONVDigitStore* store = DigitStore(fCurrentEvent);
781 if (store == 0x0) return 0x0;
782 fIterator = store->CreateIterator(detElem, detElem, 2);
783 if (fIterator == 0x0) return 0x0;
784 fCurrentIteratorType = kDigitIteratorByDetectorElement;
789 case kDigitIteratorByChamberAndCathode:
793 if (chamber < 0 or AliMpConstants::NofChambers() <= chamber)
796 "Must have give a chamber value in the range [0..%d], but got a value of: %d",
797 AliMpConstants::NofChambers() - 1,
802 if (cathode < 0 or 1 < cathode)
804 AliError(Form("Must have give a cathode value in the range [0..1], but got a value of: %d", cathode));
808 AliMUONVDigitStore* store = DigitStore(fCurrentEvent);
809 if (store == 0x0) return 0x0;
810 MpPair_t pair = AliMpDEManager::GetDetElemIdRange(chamber);
811 fIterator = store->CreateIterator(AliMp::PairFirst(pair), AliMp::PairSecond(pair), cathode);
812 if (fIterator == 0x0) return 0x0;
813 fCurrentIteratorType = kDigitIteratorByChamberAndCathode;
819 case kRawClusterIterator:
822 AliMUONVClusterStore* store = ClusterStore(fCurrentEvent);
823 if (store == 0x0) return 0x0;
824 fIterator = store->CreateChamberIterator(chamber, chamber);
825 if (fIterator == 0x0) return 0x0;
826 fCurrentIteratorType = kRawClusterIterator;
831 case kLocalTriggerIterator:
833 AliMUONVTriggerStore* store = TriggerStore(fCurrentEvent);
834 if (store == 0x0) return 0x0;
835 fIterator = store->CreateLocalIterator();
836 if (fIterator == 0x0) return 0x0;
837 fCurrentIteratorType = kLocalTriggerIterator;
841 case kRegionalTriggerIterator:
843 AliMUONVTriggerStore* store = TriggerStore(fCurrentEvent);
844 if (store == 0x0) return 0x0;
845 fIterator = store->CreateRegionalIterator();
846 if (fIterator == 0x0) return 0x0;
847 fCurrentIteratorType = kRegionalTriggerIterator;
856 //_____________________________________________________________________________
857 void AliMUONDataInterface::ResetIterator()
859 /// The temporary iterator object is deleted if it exists and the pointer reset to 0x0.
860 /// The iterator type and temporary data indicating the state of the iterator are
863 if (fIterator != 0x0) delete fIterator;
864 fCurrentIteratorType = kNoIterator;
865 fCurrentIndex = fDataX = fDataY = -1;
869 //_____________________________________________________________________________
870 Int_t AliMUONDataInterface::CountObjects(TIterator* iter)
872 /// Counts the number of objects in the iterator and resets it.
873 /// @return The number of objects in 'iter'.
875 if (iter == 0x0) return -1;
878 while ( iter->Next() != 0x0 ) count++;
884 //_____________________________________________________________________________
885 TObject* AliMUONDataInterface::FetchObject(TIterator* iter, Int_t index)
887 /// Fetches the index'th object from the iterator counting the first object
888 /// returned by iterator after it is reset as index == 0. The next object
889 /// has index == 1 and so on where the last object returned by the iterator
890 /// has index == N-1 where N = CountObjects(iter)
891 /// This method will only reset the iterator if index is smaller than
892 /// fCurrentIndex, which is used to track the iteration progress and is
893 /// updated when a new object if returned by this method.
894 /// @param iter The iterator to fetch an object from.
895 /// @param index The index number of the object to fetch in the range [0 .. N-1]
896 /// where N = CountObjects(iter)
900 AliError(Form("Index is out of bounds. Got a value of %d.", index));
904 if (iter == 0x0) return 0x0;
905 if (index <= fCurrentIndex)
911 TObject* object = 0x0;
912 while (fCurrentIndex < index)
914 object = iter->Next();
917 AliError(Form("Index is out of bounds. Got a value of %d.", index));