From: schutz Date: Tue, 14 Mar 2000 18:27:54 +0000 (+0000) Subject: The web becomes richer X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=c2583e52368d428e3188ca9019b5d77c673e272a;p=u%2Fmrichter%2FAliRoot.git The web becomes richer --- diff --git a/PHOS/PHOS-HTML/CVS_name.gif b/PHOS/PHOS-HTML/CVS_name.gif new file mode 100644 index 00000000000..0abd4fb884e Binary files /dev/null and b/PHOS/PHOS-HTML/CVS_name.gif differ diff --git a/PHOS/PHOS-HTML/ClassIndex.html b/PHOS/PHOS-HTML/ClassIndex.html index 1ed00d0d1d9..d3d23a5b491 100644 --- a/PHOS/PHOS-HTML/ClassIndex.html +++ b/PHOS/PHOS-HTML/ClassIndex.html @@ -1,4 +1,4 @@ - + @@ -19,101 +19,6 @@
diff --git a/PHOS/PHOS-HTML/PHOSConfig_Central.C b/PHOS/PHOS-HTML/PHOSConfig_Central.C new file mode 100644 index 00000000000..b8797ed3e14 --- /dev/null +++ b/PHOS/PHOS-HTML/PHOSConfig_Central.C @@ -0,0 +1,660 @@ +void Config(Int_t thickness) +{ + Text_t filename[50]; + + new AliGeant3("C++ Interface to Geant3"); + + //======================================================================= + // Create the output file + sprintf(filename,"junk.root"); + printf(">>> Output file is %s \n",filename); + + TFile *rootfile = new TFile(filename,"recreate"); + rootfile->SetCompressionLevel(2); + TGeant3 *geant3 = (TGeant3*)gMC; + + //======================================================================= + // ******* GEANT STEERING parameters FOR ALICE SIMULATION ******* + geant3->SetTRIG(1); //Number of events to be processed + geant3->SetSWIT(4,10); + geant3->SetDEBU(0,0,1); + //geant3->SetSWIT(2,2); + geant3->SetERAN(1.e-6,1.e3,90); + geant3->SetDCAY(1); + geant3->SetPAIR(1); + geant3->SetCOMP(1); + geant3->SetPHOT(1); + geant3->SetPFIS(0); + geant3->SetDRAY(0); + geant3->SetANNI(1); + geant3->SetBREM(1); + geant3->SetMUNU(1); + geant3->SetCKOV(1); + geant3->SetHADR(1); //Select pure GEANH (HADR 1) or GEANH/NUCRIN (HADR 3) + geant3->SetLOSS(2); + geant3->SetMULS(1); + geant3->SetRAYL(1); + geant3->SetAUTO(1); //Select automatic STMIN etc... calc. (AUTO 1) or manual (AUTO 0) + geant3->SetABAN(0); //Restore 3.16 behaviour for abandoned tracks + geant3->SetOPTI(2); //Select optimisation level for GEANT geometry searches (0,1,2) + Float_t cut = 1.e-3; // 1MeV cut by default + Float_t tofmax = 1.e10; + // GAM ELEC NHAD CHAD MUON EBREM MUHAB EDEL MUDEL MUPA TOFMAX + geant3->SetCUTS(cut,cut, cut, cut, cut, cut, cut, cut, cut, cut, tofmax); + // + //======================================================================= + // --- Specify event type to be tracked through the ALICE setup + // --- All positions are in cm, angles in degrees, and P and E in GeV + // + // The following Cocktail generator is defined to simulate the neutral and + // charged background in the ALICE detector. This background is important + // in the case of photon detector as PHOS. We simulated a cocktail of + // pions (pi+, pi- and pi0) , kaons (K+, K-, Kshort and Klong), eta mesons, + // omega mesons and main baryons (protons, antiprotons, neutrons and + // antineutrons) + // + // 1-Nov-1999 Gines MARTINEZ, GPS @ SUBATECH, Nantes, France + // + Int_t ParticleDensity = 100; // Number of particles created in the selected rapidity range + Int_t NumberOfPions = 0.77 * ParticleDensity ; // ~77% of pions: pi^+, pi^- and pi^0 with equal probability + Int_t NumberOfKaons = 0.11 * ParticleDensity ; // ~11% of kaons: K^+, K^-, K_short and K_long with equal probability + Int_t NumberOfEtas = 0.04 * ParticleDensity ; // ~4.0% of eta mesons + Int_t NumberOfOmegas = 0.01 * ParticleDensity ; // ~1.0% of omega mesons + Int_t NumberOfBaryons = 0.07 * ParticleDensity ; // ~7.0% of baryons: p, pbar, n and nbar with equal probability + + + // P I O N S + AliGenParam *generpion = new AliGenParam(NumberOfPions,Pion); + // AliGenPHOSlib::GetPt(Pion), AliGenPHOSlib::GetY(Pion), AliGenPHOSlib::GetIp(Pion) ; + generpion->SetWeighting(non_analog); // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt + generpion->SetForceDecay(nodecay); // Requiring that the generated particle is directly tracked by GEANT + // + // K A O N S + AliGenParam *generkaon = new AliGenParam(NumberOfKaons,Kaon); + // AliGenPHOSlib::GetPt(Kaon), AliGenPHOSlib::GetY(Kaon), AliGenPHOSlib::GetIp(Kaon) ); + generkaon->SetWeighting(non_analog); // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt + generkaon->SetForceDecay(nodecay); // Requiring that the generated particle is directly tracked by GEANT + // + // E T A S + AliGenParam *genereta = new AliGenParam(NumberOfEtas,Eta); + // AliGenPHOSlib::GetPt(Eta), AliGenPHOSlib::GetY(Eta), AliGenPHOSlib::GetIp(Eta) ); + genereta->SetWeighting(non_analog); // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt + genereta->SetForceDecay(nodecay); // Requiring that the generated particle is directly tracked by GEANT + // + // O M E G A S + AliGenParam *generomega = new AliGenParam(NumberOfOmegas,Omega); + // AliGenPHOSlib::GetPt(Omega), AliGenPHOSlib::GetY(Omega), AliGenPHOSlib::GetIp(Omega) ); + generomega->SetWeighting(non_analog); // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt + generomega->SetForceDecay(nodecay); // Requiring that the generated particle is directly tracked by GEANT + // + // B A R Y O N S + AliGenParam *generbaryon = new AliGenParam(NumberOfBaryons,Baryon); + // AliGenPHOSlib::GetPt(Baryon), AliGenPHOSlib::GetY(Baryon), AliGenPHOSlib::GetIp(Baryon) ); + generbaryon->SetWeighting(non_analog); // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt + generbaryon->SetForceDecay(nodecay); // Requiring that the generated particle is directly tracked by GEANT + + AliGenCocktail *gener = new AliGenCocktail(); // Cocktail class + gener->SetPtRange(.02,10.00); // Transverse momentum range + gener->SetPhiRange(180.,360.); // Azimuthal angle range + gener->SetYRange(-0.25,0.25); // Pseudorapidity range + gener->SetOrigin(0,0,0); // Vertex position + gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position + // + gener->AddGenerator(generpion,"pion",1.); + gener->AddGenerator(generkaon,"kaon",1.); + gener->AddGenerator(genereta,"eta",1.); + gener->AddGenerator(generomega,"omega",1.); + gener->AddGenerator(generbaryon,"baryon",1.); + gener->Init(); + + gAlice->SetField(0,0); //Specify maximum magnetic field in Tesla (neg. ==> default field) + +Int_t iMAG=1; +Int_t iITS=0; +Int_t iTPC=0; +Int_t iTOF=0; +Int_t iRICH=0; +Int_t iZDC=0; +Int_t iCASTOR=0; +Int_t iTRD=0; +Int_t iABSO=0; +Int_t iDIPO=0; +Int_t iHALL=1; +Int_t iFRAME=0; +Int_t iSHIL=0; +Int_t iPIPE=0; +Int_t iFMD=0; +Int_t iMUON=0; +Int_t iPHOS=1; +Int_t iPMD=0; +Int_t iSTART=0; + + + //=================== Alice BODY parameters ============================= + AliBODY *BODY = new AliBODY("BODY","Alice envelop"); + + +if(iMAG) { +//=================== MAG parameters ============================ +// --- Start with Magnet since detector layouts may be depending --- +// --- on the selected Magnet dimensions --- +AliMAG *MAG = new AliMAG("MAG","Magnet"); +} + + +if(iABSO) { +//=================== ABSO parameters ============================ +AliABSO *ABSO = new AliABSOv0("ABSO","Muon Absorber"); +} + +if(iDIPO) { +//=================== DIPO parameters ============================ + +AliDIPO *DIPO = new AliDIPOv2("DIPO","Dipole version 2"); +} + +if(iHALL) { +//=================== HALL parameters ============================ + +AliHALL *HALL = new AliHALL("HALL","Alice Hall"); +} + + +if(iFRAME) { +//=================== FRAME parameters ============================ + +AliFRAME *FRAME = new AliFRAMEv0("FRAME","Space Frame"); + +} + +if(iSHIL) { +//=================== SHIL parameters ============================ + +AliSHIL *SHIL = new AliSHILv0("SHIL","Shielding"); +} + + +if(iPIPE) { +//=================== PIPE parameters ============================ + +AliPIPE *PIPE = new AliPIPEv0("PIPE","Beam Pipe"); +} + + +if(iITS) { +//=================== ITS parameters ============================ +// +// EUCLID is a flag to output (=1) both geometry and media to two ASCII files +// (called by default ITSgeometry.euc and ITSgeometry.tme) in a format +// understandable to the CAD system EUCLID. The default (=0) means that you +// dont want to use this facility. +// +AliITS *ITS = new AliITSv5("ITS","normal ITS"); +ITS->SetEUCLID(0); +} + +if(iTPC) { +//============================ TPC parameters ================================ +// --- This allows the user to specify sectors for the SLOW (TPC geometry 2) +// --- Simulator. SecAL (SecAU) <0 means that ALL lower (upper) +// --- sectors are specified, any value other than that requires at least one +// --- sector (lower or upper)to be specified! +// --- Reminder: sectors 1-24 are lower sectors (1-12 -> z>0, 13-24 -> z<0) +// --- sectors 25-72 are the upper ones (25-48 -> z>0, 49-72 -> z<0) +// --- SecLows - number of lower sectors specified (up to 6) +// --- SecUps - number of upper sectors specified (up to 12) +// --- Sens - sensitive strips for the Slow Simulator !!! +// --- This does NOT work if all S or L-sectors are specified, i.e. +// --- if SecAL or SecAU < 0 +// +// +//----------------------------------------------------------------------------- + +AliTPC *TPC = new AliTPCv1("TPC","Normal TPC"); +AliTPCD *paramd = TPC->GetDigParam(); +AliTPCParam *param = &(paramd->GetParam()); + +// Set geometrical parameters + +param->SetSectorAngles(20.,10.,20.,10.); +param->SetInnerRadiusLow(83.9); +param->SetInnerRadiusUp(141.3); +param->SetOuterRadiusLow(146.9); +param->SetOuterRadiusUp(249.4); +param->SetInSecLowEdge(81.6); +param->SetInSecUpEdge(143.6); +param->SetOuSecLowEdge(144.2); +param->SetOuSecUpEdge(252.1); +param->SetEdge(1.5); +param->SetDeadZone(1.15); +param->SetPadLength(2.0); +param->SetPadWidth(0.3); +param->SetPadPitchLength(2.05); +param->SetPadPitchWidth(0.35); +param->Update(); + +if (TPC->IsVersion() != 2) paramd->Write("Param1"); + +// set gas mixture + +TPC->SetGasMixt(2,20,10,-1,0.9,0.1,0.); +TPC->SetSecAL(1); +TPC->SetSecAU(1); +// Meaningless with versions other than 2 +TPC->SetSecLows(1, 2, 3, 1+18, 2+18, 3+18); +TPC->SetSecUps(1+36, 2+36, 3+36, 1+38+18, 2+38+18, 3+38+18, -1,-1,-1,-1,-1,-1); +TPC->SetSens(1); +} + +if(iTOF) { +//=================== TOF parameters ============================ +AliTOF *TOF = new AliTOFv1("TOF","normal TOF"); +} + +if(iRICH) { +//=================== RICH parameters =========================== + + AliRICH *RICH = new AliRICHv0("RICH","normal RICH"); + + RICH->SetSMAXAR(0.03); + RICH->SetSMAXAL(-1); +// +// Version 0 +// Default Segmentation + AliRICHsegmentationV0* RsegV0 = new AliRICHsegmentationV0; + RsegV0->SetPADSIZ(.8, .8); + RsegV0->SetDAnod(0.8/3); +// Default response + AliRICHresponseV0* Rresponse0 = new AliRICHresponseV0; + AliRICHresponseCkv* RresponseCkv = new AliRICHresponseCkv; + +//------------------------Chambers 0-6 ---------------------------- + for (Int_t i=0; i<7; i++) { + RICH->SetSegmentationModel(i, 1, RsegV0); + RICH->SetResponseModel(i, mip , Rresponse0); + RICH->SetResponseModel(i, cerenkov, RresponseCkv); + RICH->Chamber(i).SetRSIGM(5.); + RICH->Chamber(i).SetMUCHSP(43.); + RICH->Chamber(i).SetMUSIGM(0.18, 0.18); + RICH->Chamber(i).SetMAXADC( 1024); + RICH->Chamber(i).SetSqrtKx3(0.77459667); + RICH->Chamber(i).SetKx2(0.962); + RICH->Chamber(i).SetKx4(0.379); + RICH->Chamber(i).SetSqrtKy3(0.77459667); + RICH->Chamber(i).SetKy2(0.962); + RICH->Chamber(i).SetKy4(0.379); + RICH->Chamber(i).SetPitch(0.25); + RICH->SetNsec(i,1); + } +} + +if(iZDC) { +//=================== ZDC parameters ============================ + +AliZDC *ZDC = new AliZDCv1("ZDC","normal ZDC"); +} + +if(iCASTOR) { +//=================== CASTOR parameters ============================ + +AliCASTOR *CASTOR = new AliCASTORv1("CASTOR","normal CASTOR"); +} + +if(iTRD) { +//=================== TRD parameters ============================ + +AliTRD *TRD = new AliTRDv1("TRD","TRD version 0"); +// Select the gas mixture (0: 97% Xe + 3% isobutane, 1: 90% Xe + 10% CO2) +TRD->SetGasMix(0); +} + +if(iFMD) { +//=================== FMD parameters ============================ + +AliFMD *FMD = new AliFMDv1("FMD","normal FMD"); +} + +if(iMUON) { +//=================== MUON parameters =========================== + +AliMUON *MUON = new AliMUONv0("MUON","normal MUON"); + + MUON->SetMaxStepGas(0.1); + MUON->SetMaxStepAlu(0.1); +// +// Version 0 +// +// First define the number of planes that are segmented (1 or 2) by a call +// to SetNsec. +// Then chose for each chamber (chamber plane) the segmentation +// and response model. +// They should be equal for the two chambers of each station. In a future +// version this will be enforced. +// +// + Int_t chamber; + Int_t station; +// Default response + AliMUONresponseV0* response0 = new AliMUONresponseV0; + response0->SetSqrtKx3(0.7131); + response0->SetKx2(1.0107); + response0->SetKx4(0.4036); + response0->SetSqrtKy3(0.7642); + response0->SetKy2(0.9706); + response0->SetKy4(0.3831); + response0->SetPitch(0.25); + response0->SetSigmaIntegration(10.); + response0->SetChargeSlope(50); + response0->SetChargeSpread(0.18, 0.18); + response0->SetMaxAdc(4096); +//-------------------------------------------------------- +// Configuration for Chamber TC1/2 (Station 1) ---------- +//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ +// Float_t rseg1[4]={17.5, 55.2, 71.3, 95.5}; + Float_t rseg1[4]={15.5, 55.2, 71.3, 95.5}; + Int_t nseg1[4]={4, 4, 2, 1}; +// + chamber=1; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV01 *seg11=new AliMUONsegmentationV01; + + seg11->SetSegRadii(rseg1); + seg11->SetPADSIZ(3, 0.5); + seg11->SetDAnod(3.0/3./4); + seg11->SetPadDivision(nseg1); + + MUON->SetSegmentationModel(chamber-1, 1, seg11); +// + AliMUONsegmentationV02 *seg12=new AliMUONsegmentationV02; + seg12->SetSegRadii(rseg1); + seg12->SetPADSIZ(0.75, 2.0); + seg12->SetDAnod(3.0/3./4); + seg12->SetPadDivision(nseg1); + + MUON->SetSegmentationModel(chamber-1, 2, seg12); + + MUON->SetResponseModel(chamber-1, response0); + + chamber=2; +//^^^^^^^^^ +// + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV01 *seg21=new AliMUONsegmentationV01; + seg21->SetSegRadii(rseg1); + seg21->SetPADSIZ(3, 0.5); + seg21->SetDAnod(3.0/3./4); + seg21->SetPadDivision(nseg1); + MUON->SetSegmentationModel(chamber-1, 1, seg21); +// + AliMUONsegmentationV02 *seg22=new AliMUONsegmentationV02; + seg22->SetSegRadii(rseg1); + seg22->SetPADSIZ(0.75, 2.); + seg22->SetDAnod(3.0/3./4); + seg22->SetPadDivision(nseg1); + MUON->SetSegmentationModel(chamber-1, 2, seg22); + + MUON->SetResponseModel(chamber-1, response0); +// +//-------------------------------------------------------- +// Configuration for Chamber TC3/4 ----------------------- +//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ +// Float_t rseg2[4]={23.5, 47.1, 87.7, 122.5}; + Float_t rseg2[4]={21.5, 47.1, 87.7, 122.5}; + Int_t nseg2[4]={4, 4, 2, 1}; +// + chamber=3; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV01 *seg31=new AliMUONsegmentationV01; + seg31->SetSegRadii(rseg2); + seg31->SetPADSIZ(3, 0.5); + seg31->SetDAnod(3.0/3./4); + seg31->SetPadDivision(nseg2); + MUON->SetSegmentationModel(chamber-1, 1, seg31); +// + AliMUONsegmentationV02 *seg32=new AliMUONsegmentationV02; + seg32->SetSegRadii(rseg2); + seg32->SetPADSIZ(0.75, 2.); + seg32->SetPadDivision(nseg2); + seg32->SetDAnod(3.0/3./4); + + MUON->SetSegmentationModel(chamber-1, 2, seg32); + + MUON->SetResponseModel(chamber-1, response0); + + chamber=4; +//^^^^^^^^^ +// + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV01 *seg41=new AliMUONsegmentationV01; + seg41->SetSegRadii(rseg2); + seg41->SetPADSIZ(3, 0.5); + seg41->SetDAnod(3.0/3./4); + seg41->SetPadDivision(nseg2); + MUON->SetSegmentationModel(chamber-1, 1, seg41); +// + AliMUONsegmentationV02 *seg42=new AliMUONsegmentationV02; + seg42->SetSegRadii(rseg2); + seg42->SetPADSIZ(0.75, 2.); + seg42->SetPadDivision(nseg2); + seg42->SetDAnod(3.0/3./4); + + MUON->SetSegmentationModel(chamber-1, 2, seg42); + + MUON->SetResponseModel(chamber-1, response0); + + +//-------------------------------------------------------- +// Configuration for Chamber TC5/6 ----------------------- +//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ +/* + seg5 = new AliMUONsegmentationV1; + AliMUONresponseV0* response5 = new AliMUONresponseV0; + // K3 = 0.62 + response5->SetSqrtKx3(0.78740079); + response5->SetKx2(0.95237319); // 0.5 * kPI * (1- 0.5*sqrtky3 ) + response5->SetKx4(0.37480633); // 0.25/TMath::ATan(sqrtkx3) + // K3 = 0.55 + response5->SetSqrtKy3(0.74161985); + response5->SetKy2(0.98832946); + response5->SetKy4(0.39177817); + response5->SetPitch(0.325); + response5->SetSigmaIntegration(10.); + response5->SetChargeSlope(50); + response5->SetChargeSpread(0.4, 0.4); + response5->SetMaxAdc(4096); + + chamber=5; + MUON->SetNsec(chamber-1,1); + MUON->SetSegmentationModel(chamber-1, 1, seg5); + MUON->SetResponseModel(chamber-1, response5); + + chamber=6; + MUON->SetNsec(chamber-1,1); + MUON->SetSegmentationModel(chamber-1, 1, seg5); + MUON->SetResponseModel(chamber-1, response5); +// +// Station 3 + station=3; + MUON->SetPADSIZ(station, 1, 0.975, 0.55); +*/ + + chamber=5; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV0 *seg51=new AliMUONsegmentationV0; + seg51->SetPADSIZ(0.75, 0.5); + seg51->SetDAnod(3.0/3./4); + MUON->SetSegmentationModel(chamber-1, 1, seg51); +// + AliMUONsegmentationV0 *seg52=new AliMUONsegmentationV0; + seg52->SetPADSIZ(0.5,0.75); + seg52->SetDAnod(3.0/3./4); + MUON->SetSegmentationModel(chamber-1, 2, seg52); + + MUON->SetResponseModel(chamber-1, response0); + + chamber=6; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV0 *seg61=new AliMUONsegmentationV0; + seg61->SetPADSIZ(0.75, 0.5); + seg61->SetDAnod(3.0/3./4); + MUON->SetSegmentationModel(chamber-1, 1, seg61); +// + AliMUONsegmentationV0 *seg62=new AliMUONsegmentationV0; + seg62->SetPADSIZ(0.5,0.75); + seg62->SetDAnod(3.0/3./4); + MUON->SetSegmentationModel(chamber-1, 2, seg62); + + MUON->SetResponseModel(chamber-1, response0); + +//-------------------------------------------------------- +// Configuration for Chamber TC7/8 (Station 4) ---------- +//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ + + Int_t nseg4[4]={4, 4, 2, 1}; + + chamber=7; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV04 *seg71=new AliMUONsegmentationV04; + seg71->SetPADSIZ(10.,0.5); + seg71->SetDAnod(0.25); + seg71->SetPadDivision(nseg4); + MUON->SetSegmentationModel(chamber-1, 1, seg71); + + AliMUONsegmentationV05 *seg72=new AliMUONsegmentationV05; + seg72->SetPADSIZ(1,10); + seg72->SetDAnod(0.25); + seg72->SetPadDivision(nseg4); + MUON->SetSegmentationModel(chamber-1, 2, seg72); + + MUON->SetResponseModel(chamber-1, response0); + + chamber=8; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); + AliMUONsegmentationV04 *seg81=new AliMUONsegmentationV04; + seg81->SetPADSIZ(10., 0.5); + seg81->SetPadDivision(nseg4); + seg81->SetDAnod(0.25); + MUON->SetSegmentationModel(chamber-1, 1, seg81); + + AliMUONsegmentationV05 *seg82=new AliMUONsegmentationV05; + seg82->SetPADSIZ(1, 10); + seg82->SetPadDivision(nseg4); + seg82->SetDAnod(0.25); + MUON->SetSegmentationModel(chamber-1, 2, seg82); + + MUON->SetResponseModel(chamber-1, response0); +//-------------------------------------------------------- +// Configuration for Chamber TC9/10 (Station 5) --------- +//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ + chamber=9; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); +// + AliMUONsegmentationV04 *seg91=new AliMUONsegmentationV04; + seg91->SetPADSIZ(10.,0.5); + seg91->SetDAnod(0.25); + seg91->SetPadDivision(nseg4); + MUON->SetSegmentationModel(chamber-1, 1, seg91); + + AliMUONsegmentationV05 *seg92=new AliMUONsegmentationV05; + seg92->SetPADSIZ(1,10); + seg92->SetDAnod(0.25); + seg92->SetPadDivision(nseg4); + + MUON->SetSegmentationModel(chamber-1, 2, seg92); + + MUON->SetResponseModel(chamber-1, response0); + + chamber=10; +//^^^^^^^^^ + MUON->SetNsec(chamber-1,2); + AliMUONsegmentationV04 *seg101=new AliMUONsegmentationV04; + seg101->SetPADSIZ(10., 0.5); + seg101->SetPadDivision(nseg4); + seg101->SetDAnod(0.25); + MUON->SetSegmentationModel(chamber-1, 1, seg101); + + AliMUONsegmentationV05 *seg102=new AliMUONsegmentationV05; + seg102->SetPADSIZ(1,10); + seg102->SetPadDivision(nseg4); + seg102->SetDAnod(0.25); + MUON->SetSegmentationModel(chamber-1, 2, seg102); + + MUON->SetResponseModel(chamber-1, response0); +//-------------------------------------------------------- +// Configuration for Trigger staions --------------------- +// (not yet used/implemented) ---------------------------- +//^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ + + chamber=11; + MUON->SetNsec(chamber-1,1); + AliMUONsegmentationV0 *seg1112=new AliMUONsegmentationV0; + seg1112->SetDAnod(0.51/3.); + + MUON->SetSegmentationModel(chamber-1, 1, seg1112); + MUON->SetResponseModel(chamber-1, response0); + + chamber=12; + MUON->SetNsec(chamber-1,1); + MUON->SetSegmentationModel(chamber-1, 1, seg1112); + MUON->SetResponseModel(chamber-1, response0); +// +// Trigger Station 1 + station=6; + MUON->SetPADSIZ(station, 1, 0.75, 0.5); + + chamber=13; + MUON->SetNsec(chamber-1,1); + AliMUONsegmentationV0 *seg1314=new AliMUONsegmentationV0; + seg1314->SetDAnod(0.51/3.); + + MUON->SetSegmentationModel(chamber-1, 1, seg1314); + MUON->SetResponseModel(chamber-1, response0); + + chamber=14; + MUON->SetNsec(chamber-1,1); + MUON->SetSegmentationModel(chamber-1, 1, seg1314); + MUON->SetResponseModel(chamber-1, response0); +// +// Trigger Station 2 + station=7; + MUON->SetPADSIZ(station, 1, 0.75, 0.5); +} + if(iPHOS) { + //=================== PHOS parameters =========================== + + AliPHOS *PHOS = new AliPHOSv0("PHOS","GPS2"); + Float_t x = (Float_t) thickness / 100. ; + if (thickness == 0) x = 0.00001 ; + PHOS->GetGeometry()-> SetLeadConverterThickness(x); + } +if(iPMD) { +//=================== PMD parameters ============================ + +AliPMD *PMD = new AliPMDv0("PMD","normal PMD"); +PMD->SetPAR(1., 1., 0.8, 0.02); +PMD->SetIN(6., 18., -580., 27., 27.); +PMD->SetGEO(0.0, 0.2, 4.); +PMD->SetPadSize(0.8, 1.0, 1.0, 1.5); + +} + +if(iSTART) { +//=================== START parameters ============================ +AliSTART *START = new AliSTARTv0("START","START Detector"); +} + +} diff --git a/PHOS/PHOS-HTML/geometry.html b/PHOS/PHOS-HTML/geometry.html index 02df3d934db..6db433ecf17 100644 --- a/PHOS/PHOS-HTML/geometry.html +++ b/PHOS/PHOS-HTML/geometry.html @@ -1,91 +1,99 @@ - - - - PHOS Geometry in AliRoot - + + + + PHOS Geometry in AliRoot + -

-Geometry & Materials -

- - -This document result from the effort lead by the Photon Group at SUBATECH -to organize within the PHOS collaboration the software needed for the performance -and physics simulation of PHOS and for the reconstruction of the raw data -(presently simulated one's and real data in the future). The PHOS software -so far distributed with the AliRoot -V3.02 package is still rather shaky, lacks documentation and organisation. -It was therefore necessary to establish a few rules to achieve a coherent -software package usable by any interested user. -

-Rules

- -
    -
  1. -SUBATECH -represented by Yves Schutz is the -librarian of the PHOS contribution to AliRoot
    2. - -
  2. -Any modification to the existing software or any new piece of software -must be approved by the librarian who is in charge of the updates
    3. - -
  3. -Developpers must used CVS to download the official software from the CERN -repository
    4. - -
  4. -Developpers must comply with the AliRoot coding convention
    5. - -
  5. -FORTRAN is a strictly prohibited language
    6. - -
  6. -Any new piece of code must be documented along the line: what is the purpose, -how is it done, how to use it
    7. -
-

-The PHOS package is organized along two directions which we will try to -keep as independent as possible. The simulation is the first one, the reconstruction -the second one. -

-Simulation

-The simulation is done in two steps, one which describes the detector geometry -and material and one which does the particle tracking and stores the hits -and digits in a TTree itself written on a disk file. These two steps are -steered through a macro file launched at the AliRoot prompt. -

-Geometry

- -

-Definition

-PHOS consists of two distinct parts. The calorimeter, named EMCA, consists -of the PbW04 crystals within their housing. The crystals are assembled -in fNModules modules of fNPhi rows along the x-axis direction -and fNZ modules along the z-axis direction. Since the total number -of crystals and their layout within ALICE is presently (11/11/1999) not -final it was decided to fully parametrize the geometry. A change in any -of the three previously defined parameters will provide automatically a -new layout of the modules. -
The second part of PHOS is the charged particle identifier. Since until -end of 2000 no decision will be taken on the final design of this element -two or more versions can be considered. Today (11/11/1999) only one version -is implemented and is named PPSD for PHOS Pre-Shower Detector. It consists -of a first layer of fNumberOfModulesPhi x fNumberOfModulesZ gas -detectors per PHOS module, each detectors being subdivided into fNumberOfPadsPhi -x fNumberOfPadsZ gas cells, a Lead converter and a second layer -of gas detectors identical to the first layer. This second part is positionned -on top of the EMCA module. -

-The class AliPHOSv0 ( inheritance tree ) describes this geometry. We have defined also a class - AliPHOSv1 ( inheritance tree ) which in addition contains the photodiode response and the light transport through the crystal. -

-The parameters of the geometry are given in class AliPHOSGeometry ( inheritance tree ) -

+

+ Geometry & Materials +

+

+ Simulation +

+ The simulation is done in two steps, one which describes the detector geometry + and material and one which does the particle tracking and stores the hits + and digits in a TTree itself written on a disk file. These two steps are + steered through the Config.C file(example), + a root macro launched at the AliRoot prompt. To customize this file you must: +

+
    +
  1. + select a file name for the root output: + 
    +	//=======================================================================
+	//  Create the output file
+	sprintf(filename,"myfilename.root");
+
    +
    2. +
  2. + Define your event generator and simulations parameters: + 
    +         AliGenCocktail *gener = new AliGenCocktail(); // Cocktail class
+         gener->SetPtRange(.02,10.00);                 // Transverse momentum range   
+         gener->SetPhiRange(180.,360.);                // Azimuthal angle range  
+         gener->SetYRange(-0.25,0.25);                 // Pseudorapidity range
+         gener->SetOrigin(0,0,0);                      // Vertex position
+         gener->SetSigma(0,0,5.6);                     // Sigma in (X,Y,Z) (cm) on IP position
+         gener->Init();                                // Initialize the generator
+
    +
    3. +
  3. + Select the detector seen by the tracking: + 
    +         Int_t iPHOS=1;    // PHOS is in
+         Int_t iPMD=0;     // PMD  is out
+
    +
    4. +
  4. + Create a PHOS object (GPS2 is so far the unique configuration option, see below): + 
     
+         if(iPHOS) {
+         //=================== PHOS parameters =========================== 
+         AliPHOS *PHOS  = new AliPHOSv0("PHOS","GPS2");
+
    +
    5. +
+

+ You van now run the simulation by entering the following command at the aliroot prompt: + 

+       aliroot> gAlice->Run(1,"PHOSConfig_Central.C");  
+
+

+

+ Geometry +

+

+ Definition +

+ PHOS consists of two distinct parts. The calorimeter, named EMCA, consists + of the PbW04 crystals within their housing. The crystals are assembled + in fNModules modules of fNPhi rows along the x-axis direction + and fNZ modules along the z-axis direction. Since the total number + of crystals and their layout within ALICE is presently (11/11/1999) not + final it was decided to fully parametrize the geometry. A change in any + of the three previously defined parameters will provide automatically a + new layout of the modules. +
The second part of PHOS is the charged particle identifier. Since until + end of 2000 no decision will be taken on the final design of this element + two or more versions can be considered. Today (11/11/1999) only one version + is implemented and is named PPSD for PHOS Pre-Shower Detector. It consists + of a first layer of fNumberOfModulesPhi x fNumberOfModulesZ gas + detectors per PHOS module, each detectors being subdivided into fNumberOfPadsPhi + x fNumberOfPadsZ gas cells, a Lead converter and a second layer + of gas detectors identical to the first layer. This second part is positionned + on top of the EMCA module. +

+ The class AliPHOSv0 ( inheritance tree ) + describes this geometry. We have defined also a class + AliPHOSv1 ( inheritance tree ) + which in addition contains the photodiode response and the light transport through the crystal. +

+ The parameters of the geometry are given in class AliPHOSGeometry ( inheritance tree ) +

Geant Tree for PHOS

@@ -97,110 +105,120 @@ The parameters of the geometry are given in class

Figure 2.: Class diagram for the geometry/simulation package
-

-AliPHOS : This is the base class. It derives from AliDetector. -It's only purpose so far (11/11/1999) is to describe the materials -(à la AliMC) needed for the EMCA and PPSD construction. -

-AliPHOSv4: It derives from AliPHOS. -
-

-

-AliPHOSGeometry: It derives from TObject to make it -persistent. It is a singleton, i.e., a pointer to the unique instance of -this class is obtained by: -

AliPHOSGeometry * Geom = AliPHOSGeometry::GetInstance()
- - -

-AliPHOSHit: It derives from AliHit. It stores the -hits as the pair (absolute Id, deposited energy) in the hit TTree when -ordered by the step manager of AliPHOSv4 -

-AliPHOSDigit : It derives from AliDigit. It stores -the digits as the pair (absolute Id, energy) in the  digit TTree when -ordered by the finish event  of AliPHOSv4 -

-Usage

-

-

-Int_t iPHOS=1;
-if(iPHOS) {
-AliPHOS * phos = new AliPHOSv4("PHOS", char * name) ;
-}
-
-
where name can take the following values: - +
+ EMCA consists of 5 modules of 64x64 modules each, positionned at + azimuthal angles -40, -20, 0, +20, +40 (see figures). PPSD consists per + EMCA module of 2 layers of 4x4 gas detectors each, each detector having + 24x24 gas cells.
© >Groupe Photons Subatech (Yves SCHUTZ) @@ -209,7 +227,7 @@ EMCA module of 2 layers of 4x4 gas detectors each, each detector having
-Last modified: Mon Mar 13 10:46:53 CET 2000 +Last modified: Tue Mar 14 18:38:42 CET 2000
Valid HTML 4.0!Valid CSS!
diff --git a/PHOS/PHOS-HTML/index.html b/PHOS/PHOS-HTML/index.html index 4b2de288874..20d972805a3 100644 --- a/PHOS/PHOS-HTML/index.html +++ b/PHOS/PHOS-HTML/index.html @@ -9,6 +9,8 @@ + image9 = new Image(356,240); + image9.src = "rootlogo.gif"; + image10 = new Image(356,240); + image10.src = "CVS_name.gif"; + image11 = new Image(356,240); + image11.src = "usecase.gif"; +