May 2012 Data, Zero Bias (221K events):

May 2012 Data, Minimum Bias (838K events):

High-beta pp Run at 8 TeV - July 2012

pA HI Run 2013.

cmsrel CMSSW_5_X_Y
cd CMSSW_5_X_Y/src
cmsenv

cd $CMSSW_BASE/src
cvs co -r [cmstotem-V00-02-02|cmstotem-pp2760-V00-01-04] -d UATree UserCode/UATree

addpkg CommonTools/RecoAlgos
mv UATree/UABaseTree/macros/RecoTrack* CommonTools/RecoAlgos/interface/
mv UATree/UABaseTree/macros/PFClusterToRefCandidate.h CommonTools/RecoAlgos/src/

cvs co UserCode/FerencSiklerVertexing
cd UserCode/FerencSiklerVertexing
scram b -c
cd $CMSSW_BASE/src

cd UATree/UADataFormat/src
bash dictionary.sh
mv LinkDef.cc LinkDefcc ; mv LinkDef.h LinkDefh

cd $CMSSW_BASE/src
scram b -r -j8

Example configuration files

Data

from UATree.UABaseTree.UABaseTree_forward_cfg import process
process.source.fileNames = ['file1','file2',...]
process.maxEvents.input = 1000
process.GlobalTag.globaltag = '[GR_R_53_V7|GR_R_53_V21|...]::All'

process.uabasetree.outputfilename = "UABaseTree_CMS-TOTEM.root"

process.uabasetree.hlt_paths = (
    'HLT_<...>',
    'HLT_<...>',
    ...
)

Monte Carlo:

from UATree.UABaseTree.UABaseTree_forward_MC_cfg import process
process.source.fileNames = ['file1','file2',...]
process.maxEvents.input = 1000
process.GlobalTag.globaltag = 'START53_V7E::All'

process.uabasetree.outputfilename = "UABaseTree_CMS-TOTEM-MC.root"

RAW Data (Saving ZDC and FSC Digi's - pA/pp Run 2013):

from UATree.UABaseTree.UABaseTree_forward_RAW_cfg import process
process.source.fileNames = ['file1','file2',...]
process.maxEvents.input = 1000
process.GlobalTag.globaltag = 'GR_R_53_V21::All'

process.uabasetree.outputfilename = "UABaseTree_CMS-TOTEM.root"

pA HI Run 2013

pA HI Run 2013.

Low energy pp Run at 2.76 TeV - January 2013

High-beta pp Run at 8 TeV - July 2012

Monte Carlo

Setting TTree branches

  MyEvtId* evtId = NULL;
  MyL1TrigOld* l1Trig = NULL;
  MyHLTrig* hltTrig = NULL;

  vector<MyTracks>* track_coll = NULL;
  vector<MyVertex>* vertex_coll = NULL;
  vector<MyPFJet>* pfJet_coll = NULL;
  vector<MyPFCand>* pFlow_coll = NULL;  
 
  tree->SetBranchAddress("cmsEvtUA",&evtId);
  tree->SetBranchAddress("cmsTrigUA",&l1Trig);
  tree->SetBranchAddress("cmsHLTTrigUA",&hltTrig);
  tree->SetBranchAddress("cmsTracksUA",&track_coll);
  tree->SetBranchAddress("cmsVerticesUA",&vertex_coll);
  tree->SetBranchAddress("cmsPFJetsUA",&pfJet_coll);
  tree->SetBranchAddress("cmsParticleFlowUA",&pFlow_coll);

Accessing L1 physics and technical trigger decision

      for (int itrig = 0 ; itrig < 128 ; ++itrig){
         if( l1Trig->PhysTrigWord[itrig] == 1)
            histosTH1F["decisionPhysTrig"]->Fill( itrig, event_weight );
      }

      for (int itrig = 0 ; itrig < 64 ; ++itrig){
         if( l1Trig->TechTrigWord[itrig] == 1 )
            histosTH1F["decisionTechTrig"]->Fill( itrig, event_weight );
      }

Accessing HLT result

      bool accept_HLT = (*hltTrig).HLTmap[ "HLT_L1DoubleJet20part1_v1" ];

Accessing vertex information

      for(vector<MyVertex>::iterator it_vtx = vertex_coll->begin() ; it_vtx != vertex_coll->end() ; ++it_vtx){
         histosTH1F["vtx_zpos"]->Fill( it_vtx->z, event_weight );
         histosTH1F["vtx_xpos"]->Fill( it_vtx->x, event_weight );
         histosTH1F["vtx_ypos"]->Fill( it_vtx->y, event_weight );
      }

Accessing tracking information

      for(vector<MyTracks>::iterator it_trk = track_coll->begin() ; it_trk != track_coll->end() ; ++it_trk){
         histosTH1F["track_pt"]->Fill( it_trk->Pt(), event_weight );
         histosTH1F["track_eta"]->Fill( it_trk->Eta(), event_weight );
         histosTH1F["track_phi"]->Fill( it_trk->Phi(), event_weight );
      }

Accessing jet information

      for(vector<MyPFJet>::iterator it_jet = pfJet_coll->begin() ; it_jet != pfJet_coll->end() ; ++it_jet){
         map<string,MyBaseJet>::iterator it_map = it_jet->mapjet.begin();
         for(; it_map != it_jet->mapjet.end(); ++it_map)
            if(verbose) cout << it_map->first << endl;

         MyBaseJet const& basejet = it_jet->mapjet[ "ak5PFL2L3Residual" ];
         histosTH1F["jet_pt"]->Fill( basejet.Pt(), event_weight );
         if(basejet.Pt() > 0.) histosTH1F["jet_eta"]->Fill( basejet.Eta(), event_weight );
         histosTH1F["jet_phi"]->Fill( basejet.Phi(), event_weight );
      }

Accessing particle-flow information

      for(vector<MyPFCand>::iterator it_pfcand = pFlow_coll->begin();
                                     it_pfcand != pFlow_coll->end(); ++it_pfcand){
         int partType = it_pfcand->particleId;
         double eta = it_pfcand->Eta();
         double energy = it_pfcand->Energy();
         histosTH2F["energyVsEtaAllTypes"]->Fill( eta, energy, event_weight );
      }

double GetMass(double P, double I, double C, double K){
  return sqrt((I-C)/K)*P;
}

TF1* GetMassLine(double M, double K, double C, double MinI)
{
   double BetaMax = 0.9;
   double PMax = sqrt((BetaMax*BetaMax*M*M)/(1-BetaMax*BetaMax));

   double BetaMin = 0.1;
   double PMin = sqrt((BetaMin*BetaMin*M*M)/(1-BetaMin*BetaMin));

   if(MinI>0)PMax = std::min(PMax,sqrt((M*M*K)/(MinI-C)));

   TF1* MassLine = new TF1("MassLine","[2] + ([0]*[0]*[1])/(x*x)", PMin, PMax);
   MassLine->SetParName  (0,"M");
   MassLine->SetParName  (1,"K");
   MassLine->SetParName  (2,"C");
   MassLine->SetParameter(0, M);
   MassLine->SetParameter(1, K);
   MassLine->SetParameter(2, C);
   MassLine->SetLineWidth(2);
   return MassLine;
}

double dEdxEstimOnTheFly(const fwlite::ChainEvent& ev, const reco::TrackRef&   track, double scaleFactor=1.0, bool usePixel=false, bool useClusterCleaning=true){
     fwlite::Handle<HSCPDeDxInfoValueMap> dEdxHitsH;
     dEdxHitsH.getByLabel(ev, "dedxHitInfo");
     if(!dEdxHitsH.isValid()){printf("Invalid dEdxHitInfo\n");return NULL;}
     const ValueMap<HSCPDeDxInfo>& dEdxHitMap = *dEdxHitsH.product();

     const HSCPDeDxInfo& hscpHitsInfo = dEdxHitMap.get((size_t)track.key());

     std::vector<double> vect_charge;
     for(unsigned int h=0;h<hscpHitsInfo.charge.size();h++){
        DetId detid(hscpHitsInfo.detIds[h]);  
        if(!usePixel && detid.subdetId()<3)continue; // skip pixels
        if(useClusterCleaning && !hscpHitsInfo.shapetest[h])continue; //skip strip hits failling the shape test.

        double Norm = (detid.subdetId()<3)?3.61e-06:3.61e-06*265;  //compute the normalization factor to get the energy in MeV/cm
        Norm*=10.0; //mm --> cm

        //save the dE/dx in MeV/cm to a vector.  The scalefactor can be used to correct for data/MC agreement.
        vect_charge.push_back(scaleFactor*Norm*hscpHitsInfo.charge[h]/hscpHitsInfo.pathlength[h]);
        //printf("%f ", Norm*hscpHitsInfo.charge[h]/hscpHitsInfo.pathlength[h]);
     }
     int size = vect_charge.size();

     double result=0;

     //build the harmonic 2 dE/dx estimator
     double expo = -2;
     for(int i = 0; i< size; i ++){
        result+=pow(vect_charge[i],expo); 
     }
     result = (size>0)?pow(result/size,1./expo):0.;
     //printf(" : %f\n",result);

     return result;
}