// providing simulation data management
//-----------------------------------------------------------------------------
+#include <TTree.h>
#include "AliMUON.h"
#include "AliMUONSDigitizerV2.h"
#include "AliMUONDigitizerV3.h"
#include "AliMUONDigitMaker.h"
+#include "AliMUONDigit.h"
#include "AliMUONCalibrationData.h"
#include "AliMUONDigitStoreV1.h"
-#include "AliMUONVTriggerStore.h"
+#include "AliMUONTriggerStoreV1.h"
#include "AliMUONHitStoreV1.h"
#include "AliMUONChamberTrigger.h"
#include "AliMUONSt2GeometryBuilderV2.h"
#include "AliMUONSlatGeometryBuilder.h"
#include "AliMUONTriggerGeometryBuilder.h"
+#include "AliMUONDigitCalibrator.h"
+#include "AliMUONRecoParam.h"
+#include "AliCDBManager.h"
+#include "AliCDBEntry.h"
#include "AliMUONRawWriter.h"
fDigitizerWithNoise(1),
fDigitizerNSigmas(4.0),
fIsTailEffect(kTRUE),
+ fConvertTrigger(kFALSE),
fRawWriter(0x0),
fDigitMaker(0x0),
fHitStore(0x0),
fDigitStoreConcreteClassName(),
- fCalibrationData(0x0)
+ fCalibrationData(0x0),
+ fDigitCalibrator(0x0)
{
/// Default Constructor
fDigitizerWithNoise(1),
fDigitizerNSigmas(4.0),
fIsTailEffect(kTRUE),
+ fConvertTrigger(kFALSE),
fRawWriter(0x0),
fDigitMaker(new AliMUONDigitMaker),
fHitStore(0x0),
fDigitStoreConcreteClassName("AliMUONDigitStoreV2S"),
- fCalibrationData()
+ fCalibrationData(),
+ fDigitCalibrator(0x0)
{
/// Standard constructor
delete fDigitMaker;
delete fHitStore;
delete fCalibrationData;
+ delete fDigitCalibrator;
}
//_____________________________________________________________________________
//_____________________________________________________________________
Bool_t AliMUON::Raw2SDigits(AliRawReader* rawReader)
{
-/// Convert raw data to SDigit
-
- fLoader->LoadDigits("READ");
+ /// Convert raw data to SDigit
+
if (!fLoader->TreeS()) fLoader->MakeSDigitsContainer();
-
+
TTree* treeS = fLoader->TreeS();
- AliMUONVDigitStore* sDigitStore = new AliMUONDigitStoreV1;
+ AliMUONVDigitStore* sDigitStore = AliMUONVDigitStore::Create(DigitStoreClassName());
+ AliMUONVTriggerStore* triggerStore = 0x0;
+
sDigitStore->Connect(*treeS);
+
+ if (!fDigitMaker) fDigitMaker = new AliMUONDigitMaker;
+
+ if (fConvertTrigger) {
+ triggerStore = new AliMUONTriggerStoreV1;
+ triggerStore->Connect(*treeS,true);
+ fDigitMaker->SetMakeTriggerDigits(true);
+ }
+
+ if (!fDigitCalibrator)
+ {
+ AliMUONRecoParam* recoParam = 0x0;
+
+ AliCDBEntry* entry = AliCDBManager::Instance()->Get("MUON/Calib/RecoParam");
+
+ if (entry)
+ {
+ // load recoParam according OCDB content (single or array)
+ if (!(recoParam = dynamic_cast<AliMUONRecoParam*>(entry->GetObject())))
+ {
+ TObjArray* recoParamArray = static_cast<TObjArray*>(entry->GetObject());
+
+ for(Int_t i = 0; i < recoParamArray->GetEntriesFast(); ++i)
+ {
+ recoParam = static_cast<AliMUONRecoParam*>(recoParamArray->UncheckedAt(i));
+ if (recoParam && recoParam->IsDefault()) break;
+ recoParam = 0x0;
+ }
+ }
+ }
+
+ if (!recoParam)
+ {
+ AliFatal("Cannot work without recoparams !");
+ }
+
+ TString calibMode = recoParam->GetCalibrationMode();
+
+ fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,recoParam,calibMode.Data());
+ }
+
+ fDigitMaker->Raw2Digits(rawReader,sDigitStore,triggerStore);
- fDigitMaker->Raw2Digits(rawReader,sDigitStore,0x0);
+ fDigitCalibrator->Calibrate(*sDigitStore);
+ TIter next(sDigitStore->CreateIterator());
+ AliMUONDigit* sdigit;
+
+ // now tweak the digits to make them "as fresh as possible", i.e.
+ // reset their calibrated status, as they'll be calibrated again
+ // once embedded.
+ while ( ( sdigit = static_cast<AliMUONDigit*>(next()) ) )
+ {
+ sdigit->Calibrated(kFALSE);
+ }
+
+ treeS->Fill();
+
fLoader->WriteSDigits("OVERWRITE");
-
+
fLoader->UnloadSDigits();
-
+
delete sDigitStore;
+
+ delete triggerStore;
return kTRUE;
}