Private inizialization of magnetic field elimitated. AliGRPManager used properly
authorgvolpe <gvolpe@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Oct 2009 17:43:08 +0000 (17:43 +0000)
committergvolpe <gvolpe@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 19 Oct 2009 17:43:08 +0000 (17:43 +0000)
HMPID/HESDfromKin.C

index e034633..f0694ad 100644 (file)
@@ -55,10 +55,11 @@ void SimEsd(AliLoader *pHL,AliESDEvent *pEsd,Bool_t htaCheck)
     Printf("-----------------------------------------------");
 }
 
-  InitGRP();
-//  AliMagF *magFieldMap = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
+  AliGRPManager *grpMan = new AliGRPManager();
+  
+  grpMan->ReadGRPEntry();
+  grpMan->SetMagField();
 
-//  AliHMPIDTracker::SetFieldMap(gAL->GetAliRun()->Field(),kTRUE);
   AliHMPIDTracker pTracker;
   AliHMPID *pH=(AliHMPID*)gAL->GetAliRun()->GetDetector("HMPID");
   Int_t iNevt=gAL->GetNumberOfEvents();
@@ -132,203 +133,3 @@ Bool_t OpenCalib()
   if(!pQthre || !pNmean) return kFALSE;  
   return kTRUE;
 }
-//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-Bool_t InitGRP() {
-  //------------------------------------
-  // Initialization of the GRP entry 
-  //------------------------------------
-  
-  AliGRPObject *fGRPData;
-  
-  AliCDBEntry* entry = AliCDBManager::Instance()->Get("GRP/GRP/Data");
-
-  if (entry) {
-
-    TMap* m = dynamic_cast<TMap*>(entry->GetObject());  // old GRP entry
-
-    if (!m) {
-       Printf("Found a TMap in GRP/GRP/Data, converting it into an AliGRPObject");
-       m->Print();
-       fGRPData = new AliGRPObject();
-       fGRPData->ReadValuesFromMap(m);
-    }
-
-    else {
-       Printf("Found an AliGRPObject in GRP/GRP/Data, reading it");
-       fGRPData = dynamic_cast<AliGRPObject*>(entry->GetObject());  // new GRP entry
-       entry->SetOwner(0);
-    }
-
-    AliCDBManager::Instance()->UnloadFromCache("GRP/GRP/Data");
-  }
-
-  if (!fGRPData) {
-     Printf("No GRP entry found in OCDB!");
-     return kFALSE;
-  }
-
-  TString lhcState = fGRPData->GetLHCState();
-  if (lhcState==AliGRPObject::GetInvalidString()) {
-    Printf("GRP/GRP/Data entry:  missing value for the LHC state ! Using UNKNOWN");
-    lhcState = "UNKNOWN";
-  }
-
-  TString beamType = fGRPData->GetBeamType();
-  if (beamType==AliGRPObject::GetInvalidString()) {
-    Printf("GRP/GRP/Data entry:  missing value for the beam type ! Using UNKNOWN");
-    beamType = "UNKNOWN";
-  }
-
-  Float_t beamEnergy = fGRPData->GetBeamEnergy();
-  if (beamEnergy==AliGRPObject::GetInvalidFloat()) {
-    Printf("GRP/GRP/Data entry:  missing value for the beam energy ! Using 0");
-    beamEnergy = 0;
-  }
-  // energy is provided in MeV*120
-  beamEnergy /= 120E3;
-
-  TString runType = fGRPData->GetRunType();
-  if (runType==AliGRPObject::GetInvalidString()) {
-    Printf("GRP/GRP/Data entry:  missing value for the run type ! Using UNKNOWN");
-    runType = "UNKNOWN";
-  }
-
-  Int_t activeDetectors = fGRPData->GetDetectorMask();
-  if (activeDetectors==AliGRPObject::GetInvalidUInt()) {
-    Printf("GRP/GRP/Data entry:  missing value for the detector mask ! Using 1074790399");
-    activeDetectors = 1074790399;
-  }
-
-  fRunInfo = new AliRunInfo(lhcState, beamType, beamEnergy, runType, activeDetectors);
-  printf("qqqqqqqqqqqqqqqqqqqqqqq %s %s %f %s %d\n", lhcState.Data(), beamType.Data(), beamEnergy, runType.Data(), activeDetectors);
-  fRunInfo->Dump();
-
-  //*** Dealing with the magnetic field map
-  if ( TGeoGlobalMagField::Instance()->IsLocked() ) {Printf("Running with the externally locked B field !");}
-  else {
-    // Construct the field map out of the information retrieved from GRP.
-    Bool_t ok = kTRUE;
-    // L3
-    Float_t l3Current = fGRPData->GetL3Current((AliGRPObject::Stats)0);
-    if (l3Current == AliGRPObject::GetInvalidFloat()) {
-      Prtinf("GRP/GRP/Data entry:  missing value for the L3 current !");
-      ok = kFALSE;
-    }
-    
-    Char_t l3Polarity = fGRPData->GetL3Polarity();
-    if (l3Polarity == AliGRPObject::GetInvalidChar()) {
-      Printf("GRP/GRP/Data entry:  missing value for the L3 polarity !");
-      ok = kFALSE;
-    }
-
-    // Dipole
-    Float_t diCurrent = fGRPData->GetDipoleCurrent((AliGRPObject::Stats)0);
-    if (diCurrent == AliGRPObject::GetInvalidFloat()) {
-      Printf("GRP/GRP/Data entry:  missing value for the dipole current !");
-      ok = kFALSE;
-    }
-
-    Char_t diPolarity = fGRPData->GetDipolePolarity();
-    if (diPolarity == AliGRPObject::GetInvalidChar()) {
-      Printf("GRP/GRP/Data entry:  missing value for the dipole polarity !");
-      ok = kFALSE;
-    }
-
-    if (ok) { 
-      if ( !SetFieldMap(l3Current, diCurrent, l3Polarity ? -1:1, diPolarity ? -1:1) )
-       AliFatal("Failed to creat a B field map ! Exiting...");
-      Printf("Running with the B field constructed out of GRP !");
-    }
-    else AliFatal("B field is neither set nor constructed from GRP ! Exitig...");
-    
-  }
-  
-  return kTRUE;
-} 
-
-//_____________________________________________________________________________
-//_____________________________________________________________________________
-  Bool_t SetFieldMap(Float_t l3Cur=30000., Float_t diCur=6000., 
-                    Float_t l3Pol=1., Float_t diPol=1., Float_t beamenergy=7000., 
-                    const Char_t* beamtype="pp",  
-                    const Char_t* path="$(ALICE_ROOT)/data/maps/mfchebKGI_sym.root")
-{
-  //------------------------------------------------
-  // The magnetic field map, defined externally...
-  // L3 current 30000 A  -> 0.5 T
-  // L3 current 12000 A  -> 0.2 T
-  // dipole current 6000 A
-  // The polarities must be the same
-  //------------------------------------------------
-  const Float_t l3NominalCurrent1=30000.; // (A)
-  const Float_t l3NominalCurrent2=12000.; // (A)
-  const Float_t diNominalCurrent =6000. ; // (A)
-
-  const Float_t tolerance=0.03; // relative current tolerance
-  const Float_t zero=77.;       // "zero" current (A)
-  //
-  TString s=(l3Pol < 0) ? "L3: -" : "L3: +";
-  //
-  AliMagF::BMap_t map = AliMagF::k5kG;
-  //
-  double fcL3,fcDip;
-  //
-  l3Cur = TMath::Abs(l3Cur);
-  if (TMath::Abs(l3Cur-l3NominalCurrent1)/l3NominalCurrent1 < tolerance) {
-    fcL3 = l3Cur/l3NominalCurrent1;
-    map  = AliMagF::k5kG;
-    s   += "0.5 T;  ";
-  } else if (TMath::Abs(l3Cur-l3NominalCurrent2)/l3NominalCurrent2 < tolerance) {
-    fcL3 = l3Cur/l3NominalCurrent2;
-    map  = AliMagF::k2kG;
-    s   += "0.2 T;  ";
-  } else if (l3Cur <= zero) {
-    fcL3 = 0;
-    map  = AliMagF::k5kGUniform;
-    s   += "0.0 T;  ";
-    fUniformField=kTRUE;        // track with the uniform (zero) B field
-  } else {
-    AliError(Form("Wrong L3 current (%f A)!",l3Cur));
-    return kFALSE;
-  }
-  //
-  diCur = TMath::Abs(diCur);
-  if (TMath::Abs(diCur-diNominalCurrent)/diNominalCurrent < tolerance) {
-    // 3% current tolerance...
-    fcDip = diCur/diNominalCurrent;
-    s    += "Dipole ON";
-  } else if (diCur <= zero) { // some small current..
-    fcDip = 0.;
-    s    += "Dipole OFF";
-  } else {
-    AliError(Form("Wrong dipole current (%f A)!",diCur));
-    return kFALSE;
-  }
-  //
-  if (l3Pol!=diPol && (map==AliMagF::k5kG || map==AliMagF::k2kG) && fcDip!=0) {
-    AliError("L3 and Dipole polarities must be the same");
-    return kFALSE;
-  }
-  //
-  if (l3Pol<0) fcL3  = -fcL3;
-  if (diPol<0) fcDip = -fcDip;
-  //
-  AliMagF::BeamType_t btype = AliMagF::kNoBeamField;
-  TString btypestr = beamtype;
-  btypestr.ToLower();
-  TPRegexp protonBeam("(proton|p)\\s*-?\\s*\\1");
-  TPRegexp ionBeam("(lead|pb|ion|a)\\s*-?\\s*\\1");
-  if (btypestr.Contains(ionBeam)) btype = AliMagF::kBeamTypeAA;
-  else if (btypestr.Contains(protonBeam)) btype = AliMagF::kBeamTypepp;
-  else {
-    Printf(Form("Cannot determine the beam type from %s, assume no LHC magnet field",beamtype));
-  }
-  Printf("------------------------------");
-  Printf(" Summary for B: %s",s.Data());
-  Printf("------------------------------");
-  AliMagF* fld = new AliMagF("MagneticFieldMap", s.Data(), fcL3, fcDip, map, btype,beamenergy,path);
-  TGeoGlobalMagField::Instance()->SetField( fld );
-  TGeoGlobalMagField::Instance()->Lock();
-  //
-  return kTRUE;
-}