Changeset 62 in svn

Timestamp:

Dec 1, 2008, 11:59:19 AM (16 years ago)

Author:

severine ovyn

Message:

Read BField info from datacard

Location:

trunk

Files:

• data/DataCardDet.dat (modified) (1 diff)
• interface/BFieldProp.h (modified) (1 diff)
• interface/SmearUtil.h (modified) (1 diff)
• src/BFieldProp.cc (modified) (2 diffs)
• src/SmearUtil.cc (modified) (3 diffs)

trunk/data/DataCardDet.dat

@@ -12 +12 @@
 FP420_S          420                    // distance of the RP to the IP, in meters
 FP420_X          0.004                  // distance of the RP to the beam, in meters
+
+#Magnetic Field 
+TRACKING_RADIUS   129                   //radius of the BField coverage
+TRACKING_LENGTH   300                   //length of the BField coverage
+BFIELD_X          0
+BFIELD_Y          0
+BFIELD_Z         3.8
+
 
 #energy resolution for electron/photon

trunk/interface/BFieldProp.h

@@ -29 +29 @@
     void Propagation(const TRootGenParticle *Part,TLorentzVector &genMomentum);
 
-    int TRACKING_RADIUS;
-    int TRACKING_LENGTH;
     int MAXITERATION;
     int MINSEGLENGTH;

trunk/interface/SmearUtil.h

@@ -45 +45 @@
   float FP420_S; // distance of the RP to the IP, in meters
   float FP420_X; // distance of the RP to the beam, in meters
+
+  //Magnetic Field information
+  int TRACKING_RADIUS;                      //radius of the BField coverage
+  int TRACKING_LENGTH;                      //length of the BField coverage
+  float BFIELD_X;          
+  float BFIELD_Y;          
+  float BFIELD_Z;         
 
   

trunk/src/BFieldProp.cc

@@ -27 +27 @@
 TrackPropagation::TrackPropagation() {
 
- TRACKING_RADIUS = 129;
- TRACKING_LENGTH = 300;
  MAXITERATION = 20000;
  MINSEGLENGTH = 70;
@@ -71 +69 @@
      double vz = pz/M;
 
-     double Bx = 0;
-     double By = 0;
-     double Bz = 3.8;
+     double Bx = BFIELD_X;
+     double By = BFIELD_Y;
+     double Bz = BFIELD_Z;
 
      double ax =  (q/M)*(Bz*vy - By*vz);

trunk/src/SmearUtil.cc

@@ -43 +43 @@
 FP420_X      =  0.004;// distance of the RP to the beam, in meters
 
+TRACKING_RADIUS  = 129;                   //radius of the BField coverage
+TRACKING_LENGTH  = 300;                   //length of the BField coverage
+BFIELD_X         = 0.0;
+BFIELD_Y         = 0.0;
+BFIELD_Z         = 3.8;
 
 ELG_Scen =        0.05;                  // S term for central ECAL
@@ -124 +129 @@
     else if(strstr(temp_string.c_str(),"MAX_CALO_FWD")){curstring >> varname >> value; MAX_CALO_FWD = value;}
     else if(strstr(temp_string.c_str(),"MAX_MU")){curstring >> varname >> value; MAX_MU = value;}
+    else if(strstr(temp_string.c_str(),"TRACKING_RADIUS")){curstring >> varname >> value; TRACKING_RADIUS = (int)value;}
+    else if(strstr(temp_string.c_str(),"TRACKING_LENGTH")){curstring >> varname >> value; TRACKING_LENGTH = (int)value;}
+    else if(strstr(temp_string.c_str(),"BFIELD_X")){curstring >> varname >> value; BFIELD_X = value;}
+    else if(strstr(temp_string.c_str(),"BFIELD_Y")){curstring >> varname >> value; BFIELD_Y = value;}
+    else if(strstr(temp_string.c_str(),"BFIELD_Z")){curstring >> varname >> value; BFIELD_Z = value;}
     else if(strstr(temp_string.c_str(),"ELG_Scen")){curstring >> varname >> value; ELG_Scen = value;}
     else if(strstr(temp_string.c_str(),"ELG_Ncen")){curstring >> varname >> value; ELG_Ncen = value;}
@@ -241 +251 @@
         << left << setw(5) <<FP420_X       <<""<< right << setw(10)<<"*"<<"\n";
   f_out<<"*                                                                    *"<<"\n";
+  f_out<<"#***********************************                                 *"<<"\n";
+  f_out<<"# Magnetic field needed informations                                 *"<<"\n";
+  f_out<<"#***********************************                                 *"<<"\n";
+  f_out<<"*                                                                    *"<<"\n";
+  f_out << left << setw(55) <<"* Radius of the BField coverage:              "<<""
+        << left << setw(5) <<TRACKING_RADIUS <<""<< right << setw(10)<<"*"<<"\n";
+  f_out << left << setw(55) <<"* Length of the BField coverage:              "<<""
+        << left << setw(5) <<TRACKING_LENGTH <<""<< right << setw(10)<<"*"<<"\n";
+  f_out << left << setw(55) <<"* BField X component:                         "<<""
+        << left << setw(5) <<BFIELD_X <<""<< right << setw(10)<<"*"<<"\n";
+  f_out << left << setw(55) <<"* BField Y component:                         "<<""
+        << left << setw(5) <<BFIELD_Y <<""<< right << setw(10)<<"*"<<"\n";
+  f_out << left << setw(55) <<"* BField Z component:                         "<<""
+        << left << setw(5) <<BFIELD_Z <<""<< right << setw(10)<<"*"<<"\n";
+  f_out<<"*                                                                    *"<<"\n";
   f_out<<"#************************************                                *"<<"\n";
   f_out<<"# Electromagnetic smearing parameters                                *"<<"\n";