X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=HMPID%2FHconfig.C;h=510b059f3de9a6139d47adc6d1a18d2c6de92056;hb=f88320157f76762bbd0e0c5f25f60b799dafbea4;hp=f29efa535d062aef1d8589bedd713d4a69c7f5c2;hpb=f7a1cc68313147ec921d4c82df1890abe00e4032;p=u%2Fmrichter%2FAliRoot.git diff --git a/HMPID/Hconfig.C b/HMPID/Hconfig.C index f29efa535d0..510b059f3de 100644 --- a/HMPID/Hconfig.C +++ b/HMPID/Hconfig.C @@ -1,10 +1,8 @@ //+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++// -//Change log: 21st June 2007 by Levente Molnar -// Detector description upgrade: ITS,ToF,Absorber, Dipole, V0, Emcal -//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++// #include +#include "AliHMPIDParam.h" Bool_t gGun1 =kFALSE; Bool_t gGun7 =kFALSE; @@ -202,13 +200,13 @@ void HmpConfig::GuiGen(TGCompositeFrame *pMainF) fGenF->AddFrame(fGenBG=new TGButtonGroup(fGenF,"Gene"));//add type button group to vertical generator frame fGenBG->Connect("Pressed(Int_t)","HmpConfig",this,"GenAddSlot(Int_t)"); fGenBG->Connect("Released(Int_t)","HmpConfig",this,"GenRemSlot(Int_t)"); new TGCheckButton(fGenBG,"gun along z" ,kGunZ); - new TGCheckButton(fGenBG,"gun to 1 chamber" ,kGun1); - new TGCheckButton(fGenBG,"gun to 7 chambers" ,kGun7); - new TGCheckButton(fGenBG,"box HMPID phase space" ,kBox ); + new TGCheckButton(fGenBG,"box in 1 chamber" ,kGun1); + new TGCheckButton(fGenBG,"box in 7 chambers" ,kGun7); + new TGCheckButton(fGenBG,"box HMPID phase space" ,kBox ); new TGCheckButton(fGenBG,"HIJING" ,kHijing); new TGCheckButton(fGenBG,"HIJING para" ,kHijingPara); new TGCheckButton(fGenBG,"Pythia" ,kPythia); - new TGCheckButton(fGenBG,"HMPID lib" ,kHmpLib); + new TGCheckButton(fGenBG,"HMPID lib" ,kHmpLib); //N prims for box and Hijing para fGenF->AddFrame(fGenNprimCO=new TGComboBox(fGenF,100)); //add N prims combo to generator vertical frame fGenNprimCO->AddEntry("not used" ,kNotUsed); @@ -248,7 +246,7 @@ void HmpConfig::GuiGen(TGCompositeFrame *pMainF) //Chamber number fGenF->AddFrame(fGenChamCO=new TGComboBox(fGenF,100)); //add chamber number combo to generator vertical frame fGenChamCO->AddEntry("not used",kNotUsed); - for(int i=1;i<=7;i++) fGenChamCO->AddEntry(Form("Chamber %i",i),i); fGenChamCO->Resize(160,20); fGenChamCO->Select(kNotUsed); + for(int i=1;i<=7;i++) fGenChamCO->AddEntry(Form("Chamber %i",i-1),i); fGenChamCO->Resize(160,20); fGenChamCO->Select(kNotUsed); }//GuiGen() //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ void HmpConfig::GenAddSlot(Int_t id) @@ -340,6 +338,11 @@ void HmpConfig::GenRemSlot(Int_t id) //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ void HmpConfig::WriteGen(FILE *pF) { + + CheckGeometry(); + + AliHMPIDParam *pP = AliHMPIDParam::Instance(); + Int_t pid=fGenPidCO->GetSelected(); Float_t pmin=0.1*fGenPminCO->GetSelected(); //particle momentum, GeV Float_t pmax=0.1*fGenPmaxCO->GetSelected(); //particle momentum, GeV @@ -353,39 +356,18 @@ void HmpConfig::WriteGen(FILE *pF) fprintf(pF," AliGenFixed *pGz=new AliGenFixed(1); pGz->SetPart(%i); pGz->SetMomentum(%.1f); pGz->SetOrigin(0,0,-200); pG->AddGenerator(pGz,\"Gz\",1);\n",pid,pmin); if(fGenBG->GetButton(kGun1)->GetState()==kButtonDown){//1 gun towards 1 HMPID chamber - switch(fGenChamCO->GetSelected()){ - case 1: fprintf(pF," AliGenFixed *pG1=new AliGenFixed(1); pG1->SetPart(%i); pG1->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG1->SetTheta(109.5); pG1->SetPhi(10); pG->AddGenerator(pG1,\"g1\",1);\n"); break; - case 2: fprintf(pF," AliGenFixed *pG2=new AliGenFixed(1); pG2->SetPart(%i); pG2->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG2->SetTheta( 90.0); pG2->SetPhi(10); pG->AddGenerator(pG2,\"g2\",1);\n"); break; - case 3: fprintf(pF," AliGenFixed *pG3=new AliGenFixed(1); pG3->SetPart(%i); pG3->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG3->SetTheta(109.5); pG3->SetPhi(30); pG->AddGenerator(pG3,\"g3\",1);\n"); break; - case 4: fprintf(pF," AliGenFixed *pG4=new AliGenFixed(1); pG4->SetPart(%i); pG4->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG4->SetTheta( 87.0); pG4->SetPhi(30); pG->AddGenerator(pG4,\"g4\",1);\n"); break; - case 5: fprintf(pF," AliGenFixed *pG5=new AliGenFixed(1); pG5->SetPart(%i); pG5->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG5->SetTheta( 70.5); pG5->SetPhi(30); pG->AddGenerator(pG5,\"g5\",1);\n"); break; - case 6: fprintf(pF," AliGenFixed *pG6=new AliGenFixed(1); pG6->SetPart(%i); pG6->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG6->SetTheta( 90.0); pG6->SetPhi(50); pG->AddGenerator(pG6,\"g6\",1);\n"); break; - case 7: fprintf(pF," AliGenFixed *pG7=new AliGenFixed(1); pG7->SetPart(%i); pG7->SetMomentum(%.1f);\n",pid,pmin); - fprintf(pF," pG7->SetTheta( 70.5); pG7->SetPhi(50); pG->AddGenerator(pG7,\"g7\",1);\n"); break; - } + Int_t ch = fGenChamCO->GetSelected()-1; + fprintf(pF," AliGenBox *pGB=new AliGenBox(1); pGB->SetPart(%i); pGB->SetMomentumRange(%.1f,%.1f);\n",pid,pmin,pmin+0.5); + fprintf(pF," pGB->SetThetaRange(%.1f,%.1f); pGB->SetPhiRange(%.1f,%.1f); pG->AddGenerator(pGB,\"g\",1);//Chamber %i\n", + pP->ChThMin(ch),pP->ChThMax(ch),pP->ChPhiMin(ch),pP->ChPhiMax(ch),ch); } if(fGenBG->GetButton(kGun7)->GetState()==kButtonDown){//7 guns towards 7 HMPID chambers - fprintf(pF," AliGenFixed *pG1=new AliGenFixed(1); pG1->SetPart(%i); pG1->SetMomentum(%.1f);pG1->SetTheta(109.5-3); pG1->SetPhi(10);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG1,\"g1\",1);\n"); - fprintf(pF," AliGenFixed *pG2=new AliGenFixed(1); pG2->SetPart(%i); pG2->SetMomentum(%.1f);pG2->SetTheta( 90.0-3); pG2->SetPhi(10);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG2,\"g2\",1);\n"); - fprintf(pF," AliGenFixed *pG3=new AliGenFixed(1); pG3->SetPart(%i); pG3->SetMomentum(%.1f);pG3->SetTheta(109.5-3); pG3->SetPhi(30);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG3,\"g3\",1);\n"); - fprintf(pF," AliGenFixed *pG4=new AliGenFixed(1); pG4->SetPart(%i); pG4->SetMomentum(%.1f);pG4->SetTheta( 90.0-3); pG4->SetPhi(30);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG4,\"g4\",1);\n"); - fprintf(pF," AliGenFixed *pG5=new AliGenFixed(1); pG5->SetPart(%i); pG5->SetMomentum(%.1f);pG5->SetTheta( 70.0-3); pG5->SetPhi(30);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG5,\"g5\",1);\n"); - fprintf(pF," AliGenFixed *pG6=new AliGenFixed(1); pG6->SetPart(%i); pG6->SetMomentum(%.1f);pG6->SetTheta( 90.0-3); pG6->SetPhi(50);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG6,\"g6\",1);\n"); - fprintf(pF," AliGenFixed *pG7=new AliGenFixed(1); pG7->SetPart(%i); pG7->SetMomentum(%.1f);pG7->SetTheta( 70.0-3); pG7->SetPhi(50);\n",pid,pmin); - fprintf(pF," pG->AddGenerator(pG7,\"g7\",1);\n"); + for(Int_t ich=0;ich<7;ich++) { + fprintf(pF," AliGenBox *pGB%i=new AliGenBox(1); pGB%i->SetPart(%i); pGB%i->SetMomentumRange(%.1f,%.1f);\n",ich,ich,pid,ich,pmin,pmin+0.5); + fprintf(pF," pGB%i->SetThetaRange(%.1f,%.1f); pGB%i->SetPhiRange(%.1f,%.1f); pG->AddGenerator(pGB%i,\"g%i\",1);//Chamber %i\n", + ich,pP->ChThMin(ich),pP->ChThMax(ich),ich,pP->ChPhiMin(ich),pP->ChPhiMax(ich),ich,ich,ich); + } } if(fGenBG->GetButton(kBox)->GetState()==kButtonDown){// box towards HMPID phase space @@ -471,18 +453,7 @@ void HmpConfig::WriteDet(FILE *pF) //CORE detectors if(fDetBG->GetButton(kPIPE )->GetState()) fprintf(pF,"\n new AliPIPEv3(\"PIPE\",\"Beam Pipe\");\n"); if(fDetBG->GetButton(kSHILD)->GetState()) fprintf(pF,"\n new AliSHILv3(\"SHIL\",\"Shielding Version 2\");\n"); - if(fDetBG->GetButton(kITS )->GetState()){ - fprintf(pF,"\n AliITSvPPRasymmFMD *ITS = new AliITSvPPRasymmFMD(\"ITS\",\"ITS PPR detailed version with asymmetric services\");\n"); - fprintf(pF,"\n ITS->SetMinorVersion(2);\n"); - fprintf(pF,"\n ITS->SetReadDet(kFALSE);\n"); - fprintf(pF,"\n ITS->SetThicknessDet1(200.);\n"); - fprintf(pF,"\n ITS->SetThicknessDet2(200.); \n"); - fprintf(pF,"\n ITS->SetThicknessChip1(150.);\n"); - fprintf(pF,"\n ITS->SetThicknessChip2(150.);\n"); - fprintf(pF,"\n ITS->SetRails(0);\n"); - fprintf(pF,"\n ITS->SetCoolingFluid(1);\n"); - fprintf(pF,"\n ITS->SetEUCLID(0);\n"); - } + if(fDetBG->GetButton(kITS )->GetState()) fprintf(pF,"\n new AliITSv11(\"ITS\",\"ITS v11\");\n"); if(fDetBG->GetButton(kTPC )->GetState()) fprintf(pF,"\n new AliTPCv2(\"TPC\",\"Default\");\n"); if(fDetBG->GetButton(kFRAME)->GetState()) fprintf(pF,"\n AliFRAMEv2 *pFrame=new AliFRAMEv2(\"FRAME\",\"Space Frame\"); pFrame->SetHoles(1);\n"); @@ -543,7 +514,7 @@ void HmpConfig::GuiBatch(TGHorizontalFrame *pMainF) new TGRadioButton(fCluBG, "Clusters CORE" ,kAll ); new TGRadioButton(fCluBG, "Clusters HMPID" ,kHmp ); - new TGCheckButton(fTrkBG, "Apply RecoParam" ,kRecoPar); + new TGCheckButton(fTrkBG, "Apply RecoParamV1" ,kRecoPar); new TGCheckButton(fTrkBG, "Load Align data" ,kAln ); new TGCheckButton(fTrkBG, "Prim vertex" ,kVtx ); @@ -689,36 +660,31 @@ void HmpConfig::WriteBatch() fprintf(fp,"void %s()\n{\n",sBatchName); fprintf(fp," gSystem->Exec(\"rm -rf *RRR \"); //remove garbage\n"); - if(fRecB->GetState()){ - - if(fMagBG->GetButton(kFld0)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,0,1,10.,2);//no field\n"); - else if(fMagBG->GetButton(kFld2)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,1,1,10,0);//0.2 Tesla field\n"); - else if(fMagBG->GetButton(kFld4)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,1,1,10,1);//0.4 Tesla field\n"); - else if(fMagBG->GetButton(kFld5)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,1,1,10,2);//0.5 Tesla field\n"); - else if(fMagBG->GetButton(kFld_2)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,0);//-0.2 Tesla field\n"); - else if(fMagBG->GetButton(kFld_4)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,1);//-0.4 Tesla field\n"); - else if(fMagBG->GetButton(kFld_5)->GetState()) fprintf(fp," AliMagF *field=new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,2);//-0.5 Tesla field\n"); - - fprintf(fp," Bool_t uniform=kTRUE;\n"); - fprintf(fp," AliTracker::SetFieldMap(field,uniform);\n\n"); - +/* + if(fMagBG->GetButton(kFld0)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 0,1,10.,2));//no field\n"); + else if(fMagBG->GetButton(kFld2)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 1,1,10,AliMagF::k2kG));//0.2 Tesla field\n"); + else if(fMagBG->GetButton(kFld4)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 1,1,10,AliMagF::k4kG));//0.4 Tesla field\n"); + else if(fMagBG->GetButton(kFld5)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2, 1,1,10,AliMagF::k5kG));//0.5 Tesla field\n"); + else if(fMagBG->GetButton(kFld_2)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,AliMagF::k2kG));//-0.2 Tesla field\n"); + else if(fMagBG->GetButton(kFld_4)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,AliMagF::k4kG));//-0.4 Tesla field\n"); + else if(fMagBG->GetButton(kFld_5)->GetState()) fprintf(fp," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,AliMagF::k5kG));//-0.5 Tesla field\n"); +*/ fprintf(fp," AliReconstruction *pRec=new AliReconstruction;\n"); fprintf(fp," gBenchmark->Start(\"ALICE\"); TDatime time; //start benchmarking\n\n"); //--------------------------------------------- if (fTrkBG->GetButton(kRecoPar)->GetState()) { - fprintf(fp," AliHMPIDRecoParam * hmpidRecoParam = AliHMPIDRecoParam::GetUserModeParam(); //Get the HMPID reco param\n"); - fprintf(fp," hmpidRecoParam->SetUserCutMode(kTRUE); //Switch to RecoParam\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(0,4); //eg cut for UserCutSigma (Values: ch0)\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(1,4); //eg cut for UserCutSigma (Values: ch1)\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(2,4); //eg cut for UserCutSigma (Values: ch2)\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(3,4); //eg cut for UserCutSigma (Values: ch3)\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(4,4); //eg cut for UserCutSigma (Values: ch4)\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(5,4); //eg cut for UserCutSigma (Values: ch5)\n"); - fprintf(fp," hmpidRecoParam->SetUserCut(6,4); //eg cut for UserCutSigma (Values: ch6)\n"); - fprintf(fp," AliHMPIDReconstructor::SetRecoParam(hmpidRecoParam); //Pass the RecoPar to the Reconstructor\n"); + fprintf(fp," AliHMPIDRecoParamV1 *hmpidRecoParam=AliHMPIDReconstructor::GetRecoParam(); //Get the HMPID reco param\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(0,4); //eg cut for UserCutSigma (Values: ch0)\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(1,4); //eg cut for UserCutSigma (Values: ch1)\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(2,4); //eg cut for UserCutSigma (Values: ch2)\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(3,4); //eg cut for UserCutSigma (Values: ch3)\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(4,4); //eg cut for UserCutSigma (Values: ch4)\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(5,4); //eg cut for UserCutSigma (Values: ch5)\n"); + fprintf(fp," hmpidRecoParam->SetHmpUserCut(6,4); //eg cut for UserCutSigma (Values: ch6)\n"); + fprintf(fp," AliHMPIDReconstructor::SetRecoParam(hmpidRecoParam); //Pass the RecoPar to the Reconstructor\n"); } //--------------------------------------------- @@ -797,13 +763,13 @@ void HmpConfig::WriteConfig() WritePhys(pF); //physics processes //Field - if(fMagBG->GetButton(kFld0)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,0,0,10,2)); //no field\n\n"); - else if(fMagBG->GetButton(kFld2)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,1,1,10,0));//0.2 Tesla field\n\n"); - else if(fMagBG->GetButton(kFld4)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,1,1,10,1));//0.4 Tesla field\n\n"); - else if(fMagBG->GetButton(kFld5)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,1,1,10,2));//0.5 Tesla field\n\n"); - else if(fMagBG->GetButton(kFld_2)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,0));//-0.2 Tesla field\n\n"); - else if(fMagBG->GetButton(kFld_4)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,1));//-0.4 Tesla field\n\n"); - else if(fMagBG->GetButton(kFld_5)->GetState()) fprintf(pF," gAlice->SetField(new AliMagF(\"Maps\",\"Maps\",2,-1,1,10,2));//-0.5 Tesla field\n\n"); + if(fMagBG->GetButton(kFld0)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",0,0,AliMagF::k5kGUniform));// NO field\n"); + else if(fMagBG->GetButton(kFld2)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k2kG));//0.2 Tesla field\n"); + // else if(fMagBG->GetButton(kFld4)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k4kG));//0.4 Tesla field\n"); + else if(fMagBG->GetButton(kFld5)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",1,1,AliMagF::k5kG));//0.5 Tesla field\n"); + else if(fMagBG->GetButton(kFld_2)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,AliMagF::k2kG));//-0.2 Tesla field\n"); + // else if(fMagBG->GetButton(kFld_4)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,AliMagF::k4kG));//-0.4 Tesla field\n"); + else if(fMagBG->GetButton(kFld_5)->GetState()) fprintf(pF," TGeoGlobalMagField::Instance()->SetField(new AliMagF(\"Maps\",\"Maps\",-1,1,AliMagF::k5kG));//-0.5 Tesla field\n"); fprintf(pF," pAL->CdGAFile();\n\n"); //???? //Generator @@ -824,3 +790,13 @@ void Hconfig() { new HmpConfig("Config.C"); } +//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ +void CheckGeometry() +{ + if(!(gSystem->IsFileInIncludePath("geometry.root"))) { + Printf(""); + Printf("HConfig needs for the geometry.root file. Please put it in the current dir!. Bye."); + Printf(""); + gApplication->Terminate(0); + } +}