/* $Id$ */
-#include <TClonesArray.h>
-#include <TMinuit.h>
+// ----------------------------
+// Class AliMUONClusterInput
+// ----------------------------
+// Global data service for hit reconstruction
+// Author: to be added
-#include "AliRun.h"
-#include "AliMUON.h"
-#include "AliMUONChamber.h"
-#include "AliMUONConstants.h"
#include "AliMUONClusterInput.h"
+
+#include "AliMUONGeometryTransformer.h"
+#include "AliMUONGeometrySegmentation.h"
+#include "AliMUONSegFactory.h"
+#include "AliMUONSegmentation.h"
+#include "AliMUONConstants.h"
#include "AliMUONMathieson.h"
#include "AliMUONRawCluster.h"
#include "AliMUONDigit.h"
+
#include "AliLog.h"
+#include <TClonesArray.h>
+#include <TClass.h>
+#include <TMinuit.h>
+#include <TGeoManager.h>
+
+/// \cond CLASSIMP
ClassImp(AliMUONClusterInput)
+/// \endcond
-AliMUONClusterInput* AliMUONClusterInput::fgClusterInput = 0;
-TMinuit* AliMUONClusterInput::fgMinuit = 0;
-AliMUONMathieson* AliMUONClusterInput::fgMathieson = 0;
+AliMUONClusterInput* AliMUONClusterInput::fgClusterInput = 0;
+TMinuit* AliMUONClusterInput::fgMinuit = 0;
+AliMUONMathieson* AliMUONClusterInput::fgMathieson = 0;
+AliMUONGeometryTransformer* AliMUONClusterInput::fgTransformer = 0;
+AliMUONSegmentation* AliMUONClusterInput::fgSegmentation = 0;
+//______________________________________________________________________________
AliMUONClusterInput::AliMUONClusterInput()
: TObject(),
+ fNseg(0),
+ fChamber(0),
fCluster(0),
+ fZ(0.),
fChargeCorrel(1.),
- fSegmentationType(1),
fDetElemId(0)
-
+
{
+/// Default constructor
+
fDigits[0]=0;
fDigits[1]=0;
- fSegmentation[0]=0;
- fSegmentation[1]=0;
fSegmentation2[0]=0;
fSegmentation2[1]=0;
}
+//______________________________________________________________________________
AliMUONClusterInput* AliMUONClusterInput::Instance()
{
-// return pointer to the singleton instance
+/// return pointer to the singleton instance
if (fgClusterInput == 0) {
fgClusterInput = new AliMUONClusterInput();
fgMinuit = new TMinuit(8);
+
+ // Create segmentation with activated Root geometry
+ if ( ! gGeoManager ) {
+ AliFatalClass("Geometry not loaded.");
+ return fgClusterInput;
+ }
+ fgTransformer = new AliMUONGeometryTransformer(true);
+ fgTransformer->ReadGeometryData("volpath.dat", gGeoManager);
+ AliMUONSegFactory factory(fgTransformer);
+ fgSegmentation = factory.CreateSegmentation();
}
return fgClusterInput;
}
+//______________________________________________________________________________
AliMUONClusterInput::~AliMUONClusterInput()
{
-// Destructor
+/// Destructor
delete fgMinuit;
delete fgMathieson;
+ delete fgTransformer;
+ delete fgSegmentation;
+ fgMinuit = 0;
+ fgMathieson = 0;
}
-AliMUONClusterInput::AliMUONClusterInput(const AliMUONClusterInput& clusterInput):TObject(clusterInput)
-{
-// Protected copy constructor
-
- AliFatal("Not implemented.");
-}
-
-void AliMUONClusterInput::SetDigits(Int_t chamber, TClonesArray* dig1, TClonesArray* dig2)
-{
-// Set pointer to digits with corresponding segmentations and responses (two cathode planes)
- fChamber=chamber;
- fDigits[0]=dig1;
- fDigits[1]=dig2;
- fNDigits[0]=dig1->GetEntriesFast();
- fNDigits[1]=dig2->GetEntriesFast();
-
- AliMUON *pMUON;
- AliMUONChamber* iChamber;
-
- pMUON = (AliMUON*) gAlice->GetModule("MUON");
- if ((fSegmentationType = pMUON->WhichSegmentation()) != 1)
- AliFatal("Wrong segmentation type");
-
- iChamber = &(pMUON->Chamber(chamber));
- fgMathieson = new AliMUONMathieson();
-
- fSegmentation[0]=iChamber->SegmentationModel(1);
- fSegmentation[1]=iChamber->SegmentationModel(2);
-
- fNseg = 2;
- if (chamber < AliMUONConstants::NTrackingCh()) {
- if (chamber > 1 ) {
- fgMathieson->SetPitch(AliMUONConstants::Pitch());
- fgMathieson->SetSqrtKx3AndDeriveKx2Kx4(AliMUONConstants::SqrtKx3());
- fgMathieson->SetSqrtKy3AndDeriveKy2Ky4(AliMUONConstants::SqrtKy3());
- fChargeCorrel = AliMUONConstants::ChargeCorrel();
- } else {
- fgMathieson->SetPitch(AliMUONConstants::PitchSt1());
- fgMathieson->SetSqrtKx3AndDeriveKx2Kx4(AliMUONConstants::SqrtKx3St1());
- fgMathieson->SetSqrtKy3AndDeriveKy2Ky4(AliMUONConstants::SqrtKy3St1());
- fChargeCorrel = AliMUONConstants::ChargeCorrelSt1();
- }
- }
-}
-
-void AliMUONClusterInput::SetDigits(Int_t chamber, TClonesArray* dig)
-{
-// Set pointer to digits with corresponding segmentations and responses (one cathode plane)
- fDigits[0]=dig;
- AliMUON *pMUON;
- AliMUONChamber* iChamber;
-
- pMUON = (AliMUON*) gAlice->GetModule("MUON");
- if ((fSegmentationType = pMUON->WhichSegmentation()) != 1)
- AliFatal("Wrong segmentation type");
-
- iChamber = &(pMUON->Chamber(chamber));
- fSegmentation[0]=iChamber->SegmentationModel(1);
-
- fNseg=1;
-}
+//______________________________________________________________________________
void AliMUONClusterInput::SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig1, TClonesArray* dig2)
{
- // Set pointer to digits with corresponding segmentations and responses (two cathode planes)
+ /// Set pointer to digits with corresponding segmentations and responses (two cathode planes)
fChamber = chamber;
fDetElemId = idDE;
fDigits[0] = dig1;
fNDigits[0] = dig1->GetEntriesFast();
fNDigits[1] = dig2->GetEntriesFast();
- AliMUON *pMUON;
- AliMUONChamber* iChamber;
-
- pMUON = (AliMUON*) gAlice->GetModule("MUON");
- iChamber = &(pMUON->Chamber(chamber));
-
+ delete fgMathieson;
fgMathieson = new AliMUONMathieson();
- if ((fSegmentationType = pMUON->WhichSegmentation()) != 2)
- AliFatal("Wrong segmentation type");
- fSegmentation2[0]=iChamber->SegmentationModel2(1);
- fSegmentation2[1]=iChamber->SegmentationModel2(2);
+ fSegmentation2[0]= fgSegmentation->GetModuleSegmentationByDEId(fDetElemId, 0);
+ fSegmentation2[1]= fgSegmentation->GetModuleSegmentationByDEId(fDetElemId, 1);
fNseg = 2;
if (chamber < AliMUONConstants::NTrackingCh()) {
}
}
+//______________________________________________________________________________
void AliMUONClusterInput::SetDigits(Int_t chamber, Int_t idDE, TClonesArray* dig)
{
-// Set pointer to digits with corresponding segmentations and responses (one cathode plane)
+/// Set pointer to digits with corresponding segmentations and responses (one cathode plane)
fChamber = chamber;
fDetElemId = idDE;
fDigits[0] = dig;
- AliMUON *pMUON;
- AliMUONChamber* iChamber;
-
- pMUON = (AliMUON*) gAlice->GetModule("MUON");
- iChamber = &(pMUON->Chamber(chamber));
- if ((fSegmentationType = pMUON->WhichSegmentation()) != 2)
- AliFatal("Wrong segmentation type");
-
- fSegmentation2[0]=iChamber->SegmentationModel2(1);
-
+ fSegmentation2[0]= fgSegmentation->GetModuleSegmentationByDEId(fDetElemId, 0);
fNseg=1;
}
+//______________________________________________________________________________
void AliMUONClusterInput::SetCluster(AliMUONRawCluster* cluster)
{
-// Set the current cluster
- //PH printf("\n %p \n", cluster);
- fCluster=cluster;
- Float_t qtot;
- Int_t i, cath, ix, iy;
- AliMUONDigit* digit;
- fNmul[0]=cluster->GetMultiplicity(0);
- fNmul[1]=cluster->GetMultiplicity(1);
- //PH printf("\n %p %p ", fDigits[0], fDigits[1]);
-
- for (cath=0; cath<2; cath++) {
- qtot=0;
- for (i=0; i<fNmul[cath]; i++) {
- // pointer to digit
- digit =(AliMUONDigit*)
+/// Set the current cluster
+ //PH printf("\n %p \n", cluster);
+ fCluster=cluster;
+ Float_t qtot;
+ Int_t i, cath, ix, iy;
+ AliMUONDigit* digit;
+ fNmul[0]=cluster->GetMultiplicity(0);
+ fNmul[1]=cluster->GetMultiplicity(1);
+ //PH printf("\n %p %p ", fDigits[0], fDigits[1]);
+
+ for (cath=0; cath<2; cath++) {
+ qtot=0;
+ for (i=0; i<fNmul[cath]; i++) {
+ // pointer to digit
+ digit =(AliMUONDigit*)
(fDigits[cath]->UncheckedAt(cluster->GetIndex(i,cath)));
// pad coordinates
ix = digit->PadX();
qtot+=fCharge[i][cath];
// Current z
Float_t xc, yc;
- if (fSegmentationType == 1)
- fSegmentation[cath]->GetPadC(ix,iy,xc,yc,fZ);
- else
- fSegmentation2[cath]->GetPadC(fDetElemId,ix,iy,xc,yc,fZ);
-
+ fSegmentation2[cath]->GetPadC(fDetElemId,ix,iy,xc,yc,fZ);
} // loop over cluster digits
fQtot[cath]=qtot;
fChargeTot[cath]=Int_t(qtot);
} // loop over cathodes
}
-
-
+//______________________________________________________________________________
Float_t AliMUONClusterInput::DiscrChargeS1(Int_t i,Double_t *par)
{
-// Compute the charge on first cathod only.
+/// Compute the charge on first cathod only.
return DiscrChargeCombiS1(i,par,0);
}
+//______________________________________________________________________________
Float_t AliMUONClusterInput::DiscrChargeCombiS1(Int_t i,Double_t *par, Int_t cath)
{
-// par[0] x-position of cluster
-// par[1] y-position of cluster
-
- Float_t q1;
- if (fSegmentationType == 1) {
-
- fSegmentation[cath]->SetPad(fix[i][cath], fiy[i][cath]);
- // First Cluster
- fSegmentation[cath]->SetHit(par[0],par[1],fZ);
- q1 = fgMathieson->IntXY(fSegmentation[cath]);
-
- } else {
-
- fSegmentation2[cath]->SetPad(fDetElemId, fix[i][cath], fiy[i][cath]);
- // First Cluster
- fSegmentation2[cath]->SetHit(fDetElemId, par[0],par[1],fZ);
- q1 = fgMathieson->IntXY(fDetElemId, fSegmentation2[cath]);
- }
-
+/// \todo add comment
+/// - par[0] x-position of cluster
+/// - param par[1] y-position of cluster
+
+ Float_t q1;
+ fSegmentation2[cath]-> SetPad(fDetElemId, fix[i][cath], fiy[i][cath]);
+ // First Cluster
+ fSegmentation2[cath]-> SetHit(fDetElemId, par[0],par[1],fZ);
+ q1 = fgMathieson->IntXY(fDetElemId, fSegmentation2[cath]);
+
Float_t value = fQtot[cath]*q1;
return value;
}
+//______________________________________________________________________________
Float_t AliMUONClusterInput::DiscrChargeS2(Int_t i,Double_t *par)
{
-// par[0] x-position of first cluster
-// par[1] y-position of first cluster
-// par[2] x-position of second cluster
-// par[3] y-position of second cluster
-// par[4] charge fraction of first cluster
-// 1-par[4] charge fraction of second cluster
+/// \todo add comment
+/// - par[0] x-position of first cluster
+/// - par[1] y-position of first cluster
+/// - par[2] x-position of second cluster
+/// - par[3] y-position of second cluster
+/// - par[4] charge fraction of first cluster
+/// - 1-par[4] charge fraction of second cluster
Float_t q1, q2;
- if (fSegmentationType == 1) {
-
- fSegmentation[0]->SetPad(fix[i][0], fiy[i][0]);
- // First Cluster
- fSegmentation[0]->SetHit(par[0],par[1],fZ);
- q1 = fgMathieson->IntXY(fSegmentation[0]);
-
- // Second Cluster
- fSegmentation[0]->SetHit(par[2],par[3],fZ);
- q2 = fgMathieson->IntXY(fSegmentation[0]);
-
- } else {
-
- fSegmentation2[0]->SetPad(fDetElemId, fix[i][0], fiy[i][0]);
- // First Cluster
- fSegmentation2[0]->SetHit(fDetElemId, par[0],par[1],fZ);
- q1 = fgMathieson->IntXY(fSegmentation[0]);
-
- // Second Cluster
- fSegmentation2[0]->SetHit(fDetElemId,par[2],par[3],fZ);
- q2 = fgMathieson->IntXY(fDetElemId, fSegmentation2[0]);
- }
-
+ fSegmentation2[0]->SetPad(fDetElemId, fix[i][0], fiy[i][0]);
+ // First Cluster
+ fSegmentation2[0]->SetHit(fDetElemId, par[0],par[1],fZ);
+ q1 = fgMathieson->IntXY(fDetElemId, fSegmentation2[0]);
+
+ // Second Cluster
+ fSegmentation2[0]->SetHit(fDetElemId,par[2],par[3],fZ);
+ q2 = fgMathieson->IntXY(fDetElemId, fSegmentation2[0]);
+
Float_t value = fQtot[0]*(par[4]*q1+(1.-par[4])*q2);
return value;
}
+//______________________________________________________________________________
Float_t AliMUONClusterInput::DiscrChargeCombiS2(Int_t i,Double_t *par, Int_t cath)
{
-// par[0] x-position of first cluster
-// par[1] y-position of first cluster
-// par[2] x-position of second cluster
-// par[3] y-position of second cluster
-// par[4] charge fraction of first cluster - first cathode
-// 1-par[4] charge fraction of second cluster
-// par[5] charge fraction of first cluster - second cathode
+/// \todo add comment
+/// - par[0] x-position of first cluster
+/// - par[1] y-position of first cluster
+/// - par[2] x-position of second cluster
+/// - par[3] y-position of second cluster
+/// - par[4] charge fraction of first cluster - first cathode
+/// - 1-par[4] charge fraction of second cluster
+/// - par[5] charge fraction of first cluster - second cathode
Float_t q1, q2;
- if (fSegmentationType == 1) {
-
- fSegmentation[cath]->SetPad(fix[i][cath], fiy[i][cath]);
- // First Cluster
- fSegmentation[cath]->SetHit(par[0],par[1],fZ);
- q1 = fgMathieson->IntXY(fSegmentation[cath]);
-
- // Second Cluster
- fSegmentation[cath]->SetHit(par[2],par[3],fZ);
- q2 = fgMathieson->IntXY(fSegmentation[cath]);
-
- } else {
- fSegmentation2[cath]->SetPad(fDetElemId,fix[i][cath], fiy[i][cath]);
- // First Cluster
- fSegmentation2[cath]->SetHit(fDetElemId,par[0],par[1],fZ);
- q1 = fgMathieson->IntXY(fDetElemId, fSegmentation2[cath]);
-
- // Second Cluster
- fSegmentation2[cath]->SetHit(fDetElemId,par[2],par[3],fZ);
- q2 = fgMathieson->IntXY(fDetElemId, fSegmentation2[cath]);
- }
+ fSegmentation2[cath]->SetPad(fDetElemId,fix[i][cath], fiy[i][cath]);
+ // First Cluster
+ fSegmentation2[cath]->SetHit(fDetElemId,par[0],par[1],fZ);
+ q1 = fgMathieson->IntXY(fDetElemId, fSegmentation2[cath]);
+
+ // Second Cluster
+ fSegmentation2[cath]->SetHit(fDetElemId,par[2],par[3],fZ);
+ q2 = fgMathieson->IntXY(fDetElemId, fSegmentation2[cath]);
+
Float_t value;
if (cath==0) {
value = fQtot[0]*(par[4]*q1+(1.-par[4])*q2);
}
return value;
}
-
-AliMUONClusterInput& AliMUONClusterInput
-::operator = (const AliMUONClusterInput& rhs)
-{
-// Protected assignement operator
-
- if (this == &rhs) return *this;
-
- AliFatal("Not implemented.");
-
- return *this;
-}