Changeset 512 in svn for trunk

Timestamp:

Jul 23, 2009, 3:42:42 PM (15 years ago)

Author:

Xavier Rouby

Message:

new commented lines for the misalignment studies

File:

• trunk/src/VeryForward.cc (modified) (4 diffs)

trunk/src/VeryForward.cc

@@ -124 +124 @@
   beamline1->add(rp220_1); 
   beamline1->add(rp420_1);
+  //beamline1->alignElement("\"MQXA.1R5\"",+0.0005,0);
+  //beamline1->alignElement("\"MQM.9R5.B1\"",-0.0001,0);
+  //beamline1->alignElement("\"MQML.5R5.B1\"",+0.0001,0);
   
   beamline2->fill(DET->RP_beam2Card,-1,DET->RP_IP_name);                              
@@ -135 +138 @@
   beamline2->add(rp220_2);
   beamline2->add(rp420_2); 
+  //beamline2->alignElement("\"MQXA.1L5\"",+0.0005,0);
   // rp220_1, rp220_2, rp420_1 and rp420_2 will be deallocated in ~H_AbstractBeamLine
   // do not put explicit delete
@@ -239 +243 @@
   // so if RP/FP too far away, the particle must be a proton or a mu+
   if( std::min(DET->RP_420_s,DET->RP_220_s) > 17 && (particle->PID != pP && particle->PID != 13)) return;
-
   if( fabs(genMomentum.Eta()) > DET->CEN_max_calo_fwd )
     {
@@ -255 +258 @@
       // here below, p1.getX(), p1.getY(), p1.getTX() and p1.getTY() =0 unless some smearing is done
       // all in micrometers or microradians
+
+      // for checking purposes
+      /*p1.smearAng();   p1.smearPos(); tx=0; ty=0;      // to be removed !!!! test only !!!!
+      ofstream xdist("xdist",ios::ios_base::app);
+      xdist << endl;
+      xdist << (1E3)*particle->X + p1.getX() + DET->RP_cross_x << " = " << (1E3)*particle->X << " + " << p1.getX() << " + " << DET->RP_cross_x << endl;
+      xdist << (1E3)*particle->Y + p1.getY() + DET->RP_cross_y << " = " << (1E3)*particle->Y << " + " << p1.getY() << " + " << DET->RP_cross_y << endl;
+      xdist << tx + p1.getTX()- kickers_on*DET->RP_cross_ang_x << " = " << tx << " + " << p1.getTX() << " - " << kickers_on << "*" << DET->RP_cross_ang_x<< endl;
+      xdist << ty + p1.getTY()- kickers_on*DET->RP_cross_ang_y << " = " << ty << " + " << p1.getTY() << " - " << kickers_on << "*" << DET->RP_cross_ang_y<< endl;
+      xdist.close();
+      */
+
+      
       p1.setPosition((1E3)*particle->X + p1.getX() + DET->RP_cross_x, 
                      (1E3)*particle->Y + p1.getY() + DET->RP_cross_y,