c->SetX(0, fSeg2[0]->GetAnod(fInput->DetElemId(), c->GetX(0)));
c->SetX(1, fSeg2[1]->GetAnod(fInput->DetElemId(), c->GetX(1)));
}
- // c->SetDetElementID(fInput->DetElemId());
+ // c->SetDetElemId(fInput->DetElemId());
// If reasonable chi^2 add result to the list of rawclusters
if (chi2 < 0.3) {
AddRawCluster(*c);
FillCluster(&cnew,cath);
}
cnew.SetClusterType(cnew.PhysicsContribution());
- // cnew.SetDetElementID(fInput->DetElemId());
+ // cnew.SetDetElemId(fInput->DetElemId());
AddRawCluster(cnew);
fNPeaks++;
}
// Setting detection element in raw cluster for alignment
// BB 19/05/05
- c.SetDetElementID(fInput->DetElemId());
+ c.SetDetElemId(fInput->DetElemId());
TClonesArray &lrawcl = *fRawClusters;
new(lrawcl[fNRawClusters++]) AliMUONRawCluster(c);
}
fNcluster[0]=fNcluster[1]=-1;
fGhost=0;
- fDetElementID = 0;
+ fDetElemId = 0;
}
//____________________________________________________
Int_t AliMUONRawCluster::Compare(const TObject *obj) const
Int_t GetOffset(Int_t i, Int_t j) const;
Float_t GetContrib(Int_t i, Int_t j) const;
Int_t GetPhysics(Int_t i) const;
- Int_t GetDetElementID() const ; // BB added
+ Int_t GetDetElemId() const ;
Int_t SetCharge(Int_t i,Int_t Q);
Int_t SetX(Int_t i, Float_t X);
Int_t SetY(Int_t i, Float_t Y);
Int_t SetZ(Int_t i, Float_t Z);
- void SetDetElementID(Int_t ID); // BB Added
+ void SetDetElemId(Int_t Id);
Int_t SetTrack(Int_t i, Int_t track);
Int_t SetPeakSignal(Int_t i, Int_t peaksignal);
Int_t SetMultiplicity(Int_t i, Int_t mul);
// 2 none give satisfactory chi2
Int_t fNcluster[2]; // Number of clusters
Float_t fChi2[2]; // Chi**2 of fit
- Int_t fDetElementID; // ID number of the detection element (slat) on which the cluster is found.
+ Int_t fDetElemId; // ID number of the detection element (slat) on which the cluster is found.
ClassDef(AliMUONRawCluster,1) //Cluster class for MUON
};
inline void AliMUONRawCluster::SetPhysics(Int_t i, Int_t physics)
{ fPhysicsMap[i] = physics; }
-inline void AliMUONRawCluster::SetDetElementID(Int_t ID)
-{ fDetElementID = ID; }
+inline void AliMUONRawCluster::SetDetElemId(Int_t Id)
+{ fDetElemId = Id; }
-inline Int_t AliMUONRawCluster::GetDetElementID() const
-{ return fDetElementID;}
+inline Int_t AliMUONRawCluster::GetDetElemId() const
+{ return fDetElemId;}
#endif
{
//
// Config file for MUON test
- // Gines MARTINEZ, Subatech, mai 2003, august 2003
- //
+ //
//=====================================================================
// Libraries required by geant321
gSystem->Load("libgeant321.so");
- new TGeant3("C++ Interface to Geant3");
+ new TGeant3TGeo("C++ Interface to Geant3");
//=======================================================================
// Create the output file
Text_t filename[100];
// New MUONv1 version (geometry defined via builders)
AliMUON *MUON = new AliMUONv1("MUON", "default");
+ // If align = true, the detection elements transformations
+ // are taken from the input files and not from the code
//MUON->SetAlign(true);
- // If align = true, the detection elements transformations
- // are taken from the input files and not from the code
-
-// ((AliMUONv1*)MUON)->SetStepManagerVersionDE(true);
- // To save detection element Id in the hits
- MUON->SetSegmentationType(1);// default wise to old (1), new (2)
-
- //MUON->AddGeometryBuilder(new AliMUONSt1GeometryBuilder(MUON));
+
+ // To save detection element Id in the hits
+ ((AliMUONv1*)MUON)->SetStepManagerVersionDE(true);
+
+ MUON->SetSegmentationType(2);// default wise to old (1), new (2)
MUON->AddGeometryBuilder(new AliMUONSt1GeometryBuilderV2(MUON));
- MUON->AddGeometryBuilder(new AliMUONSt2GeometryBuilder(MUON));
+ MUON->AddGeometryBuilder(new AliMUONSt2GeometryBuilderV2(MUON));
MUON->AddGeometryBuilder(new AliMUONSlatGeometryBuilder(MUON));
MUON->AddGeometryBuilder(new AliMUONTriggerGeometryBuilder(MUON));
}
Float_t z1 = mRecPoint->GetZ(1);
Float_t chi2_0 = mRecPoint->GetChi2(0);
Float_t chi2_1 = mRecPoint->GetChi2(1);
- Int_t de = mRecPoint->GetDetElementID();
- // printf(">>> RecPoint %4d x %6.3f %6.3f y %6.3f %6.3f z %6.3f %6.3f Q0 %4d Q1 %4d Hit %4d Track1 %4d Track2 %4d Chi2 %6.3f %6.3f \n",
- //irecpoint, x0, x1, y0, y1, z0, z1, Q0, Q1, Track0, Track1, Track2, chi2_0, chi2_1);
- cout << mRecPoint->GetDetElemId() << endl;
- printf(">>> RecPoint %4d x %6.3f y %6.3f z %6.3f DetElem %2d \n",irecpoint,x0,y0,z0,de);
+ Int_t de = mRecPoint->GetDetElemId();
+ printf(">>> RecPoint %4d DetElem %2d x %6.3f %6.3f y %6.3f %6.3f z %6.3f %6.3f Q0 %4d Q1 %4d Hit %4d Track1 %4d Track2 %4d Chi2 %6.3f %6.3f \n",
+ irecpoint,de,x0,x1,y0,y1,z0,z1,Q0,Q1,Track0, Track1, Track2, chi2_0, chi2_1);
} // end recpoint loop
} // end chamber loop
muondata.ResetRawClusters();
for hits, you have to do the following before:
root [0] gAlice->SetConfigFunction("Config( \"/home/martinez/aliroot/work_NewIO/test/\" , \"box\" );");
-
============================================================
How to run MUONCheck macro
============================================================
============================================================
- How to check the Geometry
+ How to check the Geometry
============================================================
gAlice->Init("$ALICE_ROOT/MUON/Config.C");
.L $ALICE/geant3/TGeant3/G3GUI.C
G3GUI()
+============================================================
+ How to check the Geometry with the new Geometrical modeler
+ ftp://root.cern.ch/root/doc/chapter16.pdf
+ http://agenda.cern.ch/fullAgenda.php?ida=a05212
+============================================================
+gAlice->Init("$ALICE_ROOT/MUON/Config.C");
+gGeoManager->GetMasterVolume()->Draw();
+
============================================================
How to run MUONdisplay
===========================================================
csh Script for the full reconstruction with raw data generator
===========================================================
-
-
-aliroot -b << EOF
-AliSimulation MuonSim("YourConfig.C")
-MuonSim.SetWriteRawData("MUON");
-MuonSim.Run(XXX)
-.q
-EOF
-
-aliroot -b << EOF
-TPluginManager* pluginManager = gROOT->GetPluginManager();
-pluginManager->AddHandler("AliReconstructor", "MUON","AliMUONReconstructor", "MUON","AliMUONReconstructor()")
-AliReconstruction MuonRec("galice.root")
-MuonRec.SetRunTracking("")
-MuonRec.SetRunVertexFinder(kFALSE)
-MuonRec.SetRunLocalReconstruction("MUON")
-MuonRec.SetFillESD("MUON")
-MuonRec.Run()
-.q
-EOF
-
-aliroot -b << EOF
-.includepath $ALICE_ROOT/STEER
-.includepath $ALICE_ROOT/MUON
-.L $ALICE_ROOT/MUON/MUONmassPlot_ESD.C++
-MUONmassPlot("galice.root",0,99999);
-.q
-EOF
+See file AlirootRun_MUONtest.script
+in $ALICE_ROOT/MUON
============================================================