// 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"
#include "AliLoader.h"
#include "AliCDBManager.h"
-#include "AliRunDigitizer.h"
+#include "AliDigitizationInput.h"
#include "AliMC.h"
#include "AliRun.h"
#include "AliRawDataHeaderSim.h"
fCurIterPad(0),
fIsMaxStep(kTRUE),
fTriggerScalerEvent(kFALSE),
- fTriggerResponseV1(kFALSE),
+ fTriggerResponseV1(0),
fTriggerCoinc44(0),
- fTriggerEffCells(0),
+ fTriggerEffCells(kTRUE),
fDigitizerWithNoise(1),
fDigitizerNSigmas(4.0),
fIsTailEffect(kTRUE),
+ fConvertTrigger(kFALSE),
fRawWriter(0x0),
fDigitMaker(0x0),
fHitStore(0x0),
fDigitStoreConcreteClassName(),
- fCalibrationData(0x0)
+ fCalibrationData(0x0),
+ fDigitCalibrator(0x0)
{
/// Default Constructor
fCurIterPad(0),
fIsMaxStep(kTRUE),
fTriggerScalerEvent(kFALSE),
- fTriggerResponseV1(kFALSE),
+ fTriggerResponseV1(0),
fTriggerCoinc44(0),
- fTriggerEffCells(0),
+ fTriggerEffCells(kTRUE),
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;
}
//_____________________________________________________________________________
}
//____________________________________________________________________
-AliDigitizer* AliMUON::CreateDigitizer(AliRunDigitizer* manager) const
+AliDigitizer* AliMUON::CreateDigitizer(AliDigitizationInput* digInput) const
{
/// Return digitizer
- AliMUONDigitizerV3* digitizer = new AliMUONDigitizerV3(manager, fDigitizerWithNoise);
+ AliMUONDigitizerV3* digitizer = new AliMUONDigitizerV3(digInput, fDigitizerWithNoise);
AliMUONDigitizerV3::SetNSigmas(fDigitizerNSigmas);
- digitizer->SetCalibrationData(fCalibrationData);
+ digitizer->SetCalibrationData(fCalibrationData,GetRecoParam());
return digitizer;
}
/// Perform Hits2Digits using SDigitizerV2
AliMUONSDigitizerV2 sdigitizer;
- sdigitizer.ExecuteTask();
+ sdigitizer.Digitize();
}
//_____________________________________________________________________
//_____________________________________________________________________
Bool_t AliMUON::Raw2SDigits(AliRawReader* rawReader)
{
-/// Convert raw data to SDigit
-
+ /// 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 = GetRecoParam();
+
+ if (!recoParam)
+ {
+ AliFatal("Cannot work without recoparams !");
+ }
+
+ fDigitCalibrator = new AliMUONDigitCalibrator(*fCalibrationData,recoParam);
+ }
+
+ fDigitMaker->Raw2Digits(rawReader,sDigitStore,triggerStore);
+
+ 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);
+ }
- fDigitMaker->Raw2Digits(rawReader,sDigitStore,0x0);
treeS->Fill();
+
fLoader->WriteSDigits("OVERWRITE");
-
+
fLoader->UnloadSDigits();
-
+
delete sDigitStore;
+
+ delete triggerStore;
return kTRUE;
}
}
//____________________________________________________________________
-Bool_t AliMUON::GetTriggerResponseV1() const
+Int_t AliMUON::GetTriggerResponseV1() const
{
///
/// Returns fTriggerResponseV1
}
+//____________________________________________________________________
+AliMUONRecoParam* AliMUON::GetRecoParam() const
+{
+ 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;
+ }
+ }
+ }
+ return recoParam;
+}
+