How to do a susy D3PD based Analysis

Creating a class from a D3PD.

irst note that the all the susy d3pds have a TTree called 'susy'. The first thing is open a root sesion and do:

TChain c("susy");
c.Add("/path/to/data/file_d3pd_00001.root")
c.Add("/path/to/data/file_d3pd_00002.root")
c.MakeClass("analysis");
.q

Now you can see two generated files: analysis.h and analysis.C. The method MakeClass can receives any string. Note that treename is the name of the tree stored in the D3PD, and it should be the same for all the D3PSs intended to be added to the chain.

Now a little modifications to the analysis.h file. The first thing is to add a header include:

#include <vector>

and before the definition of the class add:

using namespace std;

This is because the generated code was designed to run on lxplus, not in an stand alone application.

My good runs lists

    bool grl[70];
    // period A
    grl[0]  = ( RunNumber == 152166 ) and ( (lbn >= 206) and (lbn <= 300) );
    grl[1]  = ( RunNumber == 152214 ) and ( (lbn >= 159) and (lbn <= 201) );
    grl[2]  = ( RunNumber == 152220 ) and ( (lbn >=   1) and (lbn <=   3) );
    grl[3]  = ( RunNumber == 152221 ) and ( (lbn >=   5) and (lbn <= 167) );
    grl[4]  = ( RunNumber == 152345 ) and ( (lbn >= 128) and (lbn <= 207) );
    grl[5]  = ( RunNumber == 152409 ) and ( ( (lbn >= 126) and (lbn <= 209) ) or ( (lbn >= 286) and (lbn <= 716) ) );
    grl[6]  = ( RunNumber == 152441 ) and ( ( (lbn >= 317) and (lbn <= 370) ) or ( (lbn >= 372) and (lbn <= 405) ) or ( (lbn >= 411) and (lbn <= 671) ) );
    grl[7]  = ( RunNumber == 152508 ) and ( ( (lbn >= 196) and (lbn <= 204) ) or ( (lbn >= 241) and (lbn <= 261) ) );
    grl[8]  = ( RunNumber == 152777 ) and ( ( (lbn >=  63) and (lbn <=  80) ) or ( (lbn >=  89) and (lbn <= 299) ) or ( (lbn >= 321) and (lbn <= 329) ) );
    grl[9]  = ( RunNumber == 152844 ) and ( (lbn >= 195) and (lbn <= 234) );
    grl[10] = ( RunNumber == 152845 ) and ( ( (lbn >= 111) and (lbn <= 134) ) or ( (lbn >= 145) and (lbn <= 169) ) or ( (lbn >= 186) and (lbn <= 349) ) );
    grl[11] = ( RunNumber == 152878 ) and ( ( (lbn >=  99) and (lbn <= 159) ) or ( (lbn >= 191) and (lbn <= 214) ) );
    grl[12] = ( RunNumber == 152933 ) and ( (lbn >=  46) and (lbn <= 173) );
    grl[13] = ( RunNumber == 152994 ) and ( (lbn >= 293) and (lbn <= 353) );
    grl[14] = ( RunNumber == 153030 ) and ( ( (lbn >= 120) and (lbn <= 164) ) or ( (lbn >= 186) and (lbn <= 203) ) );
    ( RunNumber == 153134 )? grl[15] = false : grl[15] = false;
    ( RunNumber == 153136 )? grl[16] = false : grl[16] = false ;
    grl[17] = ( RunNumber == 153159 ) and ( ( (lbn >=  87) and (lbn <= 104) ) or ( (lbn >= 126) and (lbn <= 177) ) );
    grl[18] = ( RunNumber == 153200 ) and ( (lbn >= 156) and (lbn <= 178) );
    // period B
    grl[19] = ( RunNumber == 153565 ) and ( ( (lbn >=  258) and (lbn <= 269) ) or ( (lbn >= 281) and (lbn <= 953) ) );
    ( RunNumber == 153599 )? grl[20] = false : grl[20] = false;
    grl[21] = ( RunNumber == 154810 ) and (lbn == 158);
    ( RunNumber == 154813 )? grl[22] = false : grl[22] = false;
    grl[23] = ( RunNumber == 154815 ) and ( (lbn >=   7) and (lbn <=   8) );
    ( RunNumber == 154817 )? grl[24] = false : grl[24] = false;
    ( RunNumber == 154822 )? grl[25] = false : grl[25] = false;
    grl[26] = ( RunNumber == 155073 ) and ( (lbn >=  82) and (lbn <= 407) );
    grl[27] = ( RunNumber == 155112 ) and ( ( (lbn >= 127) and (lbn <= 166) ) or ( (lbn >= 186) and (lbn <= 625) ) );
    grl[28] = ( RunNumber == 155116 ) and ( ( (lbn >=   8) and (lbn <=  54) ) or ( (lbn >=  61) and (lbn <=  93) ) );
    ( RunNumber == 155118 )? grl[29] = false : grl[29] = false;
    grl[30] = ( RunNumber == 155160 ) and ( (lbn >= 240) and (lbn <= 503) );
    // period C
    grl[31] = ( RunNumber == 155228 ) and ( (lbn >= 131) and (lbn <= 148) );
    grl[32] = ( RunNumber == 155280 ) and ( (lbn >= 360) and (lbn <= 375) );
    grl[33] = ( RunNumber == 155569 ) and ( (lbn >= 228) and (lbn <= 467) );
    grl[34] = ( RunNumber == 155634 ) and ( ( (lbn >= 144) and (lbn <= 147) ) or ( (lbn >= 149) and (lbn <= 333) ) );
    grl[35] = ( RunNumber == 155669 ) and ( (lbn >= 264) and (lbn <= 311) );
    grl[36] = ( RunNumber == 155678 ) and ( (lbn >= 241) and (lbn <= 308) );
    grl[37] = ( RunNumber == 155697 ) and ( (lbn >= 264) and (lbn <= 506) );
    grl[38] = ( RunNumber == 156682 ) and ( (lbn >= 405) and (lbn <= 512) );
    // period D
    ( RunNumber == 158045 )? grl[39] = false : grl[39] = false;
    ( RunNumber == 158116 )? grl[40] = false : grl[40] = false;
    grl[41] = ( RunNumber == 158269 ) and ( ( (lbn >=   6) and (lbn <=  59) ) or ( (lbn >=  68) and (lbn <=  81) ) );
    grl[42] = ( RunNumber == 158299 ) and ( (lbn >= 378) and (lbn <= 434) );
    ( RunNumber == 158392 )? grl[43] = false : grl[43] = false;
    grl[44] = ( RunNumber == 158443 ) and ( ( (lbn >= 210) and (lbn <= 213) ) or ( (lbn >= 216) and (lbn <= 230) ) );
    grl[45] = ( RunNumber == 158466 ) and ( ( (lbn >= 265) and (lbn <= 268) ) or ( (lbn >= 273) and (lbn <= 277) ) );
    grl[46] = ( RunNumber == 158545 ) and ( (lbn >= 236) and (lbn <= 261) );
    grl[47] = ( RunNumber == 158548 ) and ( (lbn >= 111) and (lbn <= 304) );
    ( RunNumber == 158549 )? grl[48] = false : grl[48] = false;
    grl[49] = ( RunNumber == 158582 ) and ( (lbn >=  56) and (lbn <= 372) );
    grl[50] = ( RunNumber == 158632 ) and ( ( (lbn >= 163) and (lbn <= 178) ) or ( (lbn >= 182) and (lbn <= 272) ) );
    grl[51] = ( RunNumber == 158801 ) and ( ( (lbn >= 135) and (lbn <= 150) ) or ( (lbn >= 160) and (lbn <= 311) ) );
    grl[52] = ( RunNumber == 158975 ) and ( ( (lbn >= 105) and (lbn <= 244) ) or ( (lbn >= 253) and (lbn <= 282) ) or ( (lbn >= 315) and (lbn <= 352) ) );
    grl[53] = ( RunNumber == 159041 ) and ( ( (lbn >= 105) and (lbn <= 144) ) or ( (lbn >= 164) and (lbn <= 392) ) );
    grl[54] = ( RunNumber == 159086 ) and ( ( (lbn >= 116) and (lbn <= 168) ) or ( (lbn >= 179) and (lbn <= 227) ) or ( (lbn >= 235) and (lbn <= 515) ) or ( (lbn >= 522) and (lbn <= 556) ) );
    grl[55] = ( RunNumber == 159113 ) and ( (lbn >= 259) and (lbn <= 519) );
    grl[56] = ( RunNumber == 159179 ) and ( (lbn >= 221) and (lbn <= 334) );
    grl[57] = ( RunNumber == 159202 ) and ( (lbn >= 186) and (lbn <= 269) );
    grl[58] = ( RunNumber == 159203 ) and ( (lbn >=   3) and (lbn <=  82) );
    grl[59] = ( RunNumber == 159224 ) and ( ( (lbn >= 249) and (lbn <= 423) ) or ( (lbn >= 425) and (lbn <= 430) ) or ( (lbn >= 448) and (lbn <= 612) ) or ( (lbn >= 614) and (lbn <= 813) ) );

    if( grl[0]  or grl[1]  or grl[2]  or grl[3]  or grl[4]  or grl[5]  or grl[6]  or grl[7]  or grl[8]  or grl[9]  or
        grl[10] or grl[11] or grl[12] or grl[13] or grl[14] or grl[15] or grl[16] or grl[17] or grl[18] or grl[19] or
        grl[20] or grl[21] or grl[22] or grl[23] or grl[24] or grl[25] or grl[26] or grl[27] or grl[28] or grl[29] or
        grl[30] or grl[31] or grl[32] or grl[33] or grl[34] or grl[35] or grl[36] or grl[37] or grl[38] or grl[39] or
        grl[40] or grl[41] or grl[42] or grl[43] or grl[44] or grl[45] or grl[46] or grl[47] or grl[48] or grl[49] or
        grl[50] or grl[51] or grl[52] or grl[53] or grl[54] or grl[55] or grl[56] or grl[57] or grl[58] or grl[59]
      )
    {

-- FernandoAndresQuinonezGranados - 18-Jul-2010